Upon completing this tutorial, the user will be familiar with performing a simulation of external, laminar flow around a 2-D geometry. The specific geometry chosen for the tutorial is a cylinder. Consequently, the following capabilities of SU2 will be showcased:
- Steady, 2-D Laminar Navier-Stokes equations
- Roe (Second-Order) numerical scheme in space
- Euler implicit time integration
- Navier-Stokes wall and far-field boundary conditions
In this tutorial, we discuss some numerical method options, including how to activate a slope limiter for upwind methods.
The resources for this tutorial can be found in the TestCases/navierstokes/cylinder directory. You will need the configuration file (lam_cylinder.cfg) and the mesh file (mesh_cylinder_lam.su2).
Experimental results for drag over a cylinder at low Reynolds numbers are reported in the following article:
D. J. Tritton, "Experiments on the flow past a circular cylinder at low Reynolds numbers," Journal of Fluid Mechanics, Vol. 6, No. 4, pp. 547-567, 1959.
The following tutorial will walk you through the steps required when solving for the flow around a cylinder using SU2. It is assumed you have already obtained and compiled the SU2_CFD code for a serial computation. If you have yet to complete these requirements, please see the Download and Installation pages.
The flow around a (geometrically) two-dimensional circular cylinder is case that has been used both as a validation case and as a legitimate research case. At very low Reynolds numbers of less than 46, the flow is steady and symmetrical. As the Reynolds number is increased, asymmetries and time-dependence develop, eventually resulting in the Von Karmann vortex street, and then on to turbulence.
This problem will solve the for the external flow over the cylinder with these conditions:
- Freestream temperature = 273.15 K
Freestream Mach number = 0.1
- Angle of attack (AoA) = 0.0 degrees
- Reynolds number = 40 for a cylinder radius of 0.3048 m (1 ft)
The problem geometry is two-dimensional.The mesh has 26,192 triangular elements. It is fine near the surface of the cylinder to resolve the boundary layer. The exterior boundary is approximately 15 diameters away from the cylinder surface to avoid interaction between the boundary conditions. Far-field boundary conditions are used at the outer boundary. No-slip boundary conditions are placed on the surface of the cylinder.
Figure (1): The computational mesh for the 2-D cylinder test case. The outer boundary in red is the far-field, and the small circle in the center is the cylinder which uses the Navier-Stokes Wall boundary condition.
Configuration File Option
Several of the key configuration file options for this simulation are highlighted here. Please see the configuration file page for a general description of all options.
Some options concerning numerics:
For laminar flow around the cylinder, the 2nd-order Roe upwinding method showed good performance when the Venkatakrishnan limiter was used. Without the limiter, the computation is much less stable and may not converge. The viscous terms are computed with the corrected average of gradients method.
The cylinder simulation for the 13,336 node mesh is small and will execute relatively quickly on a single workstation or laptop in serial. To run this test case, follow these steps at a terminal command line:
- Move to the directory containing the config file (lam_cylinder.cfg) and the mesh file (mesh_cylinder_lam.su2). Make sure that the SU2 tools were compiled, installed, and that their install location was added to your path.
- Run the executable by entering "SU2_CFD lam_cylinder.cfg" at the command line.
- SU2 will print residual updates with each iteration of the flow solver, and the simulation will terminate after meeting the specified convergence criteria.
- Files containing the results will be written upon exiting SU2. The flow solution can be visualized in ParaView (.vtk) or Tecplot (.dat for ASCII).
The following results show the flow around the cylinder as calculated by SU2 (note that these were for a slightly higher Mach number of 0.3).
Figure (2): Pressure contours around the cylinder.
Figure (3): Laminar viscosity contours for this steady, low Reynolds number flow.
Figure (4): Mach number contours around the cylinder with streamlines. Note the large laminar separation region behind the cylinder at Re = 40.
Please take a minute and complete a brief survey about your experience with this tutorial. Your feedback is anonymous and goes toward improving your experience with the software.