DFG flow around cylinder, laminar case (Re=20)
General information and references
The DFG-benchmark in general simulates flow in a pipe around a fixed circular object. Different inflow speeds decide about the behaviour of the flow. The benchmark is set up in 2D.
Please consider also the following literature concerning this benchmark:
|||On High Order Methods for the Stationary Incompressible Navier-Stokes Equations; University of Heidelberg; Preprint 42/98, 1998; Download|
|||Higher-order finite element discretizations in a benchmark problem for incompressible flows; John, V.; Matthies, G.; Int. J. Num. Meth. Fluids 2001; 37: 885-903|
Geometry and flow configuration:
The geometry of the benchmark consists of a simple channel of length 2.2 and height 0.41. At (x,y)=(0.2,0.2) a cylinder with diameter L=0.1 is placed.
- Figure: The CFD benchmark geometry
The fluid is characterized by a viscosity-parameter of ν=0.001 and a density of ρ=1. The problem to solve is the stationary Navier-Stokes equation
On the upper wall, the lower wall and on the cylinger, Dirichlet-0-boundary with velocity
is imposed. The left wall is set to a parabolic inflow profile with maximum inflow velocity U=0.3:
The boundary conditions on the right wall are Neumann boundary conditions (i.e. do-nothing boundary conditions).
The defined benchmark setup results in a Reynolds number of
To verify the results visually, check the following pictures, which show the velocity, the pressure and the streamfunction in the domain:
The DFG benchmark is knwon for many years and had been extensively analyzed. In the following, we describe how to check the computed results by numerical numbers, independent of the above pictures.
In the stationary case, the correctness of the results are measured by three numbers:
- The drag coefficient
- The lift coefficient
- The pressure drop
with A=(0.15,0.2) and B=(0.25,0.2)being the leftmost and rightmost point of the cylinder, respectively.
- Figure: Measuring the pressure
The surface of the cylinder is denoted by S. Let , denote the tangential and normal vector of S, respectively. Then, the forces are calculated using the formula
The definition of the stress tensor σ in this situation is
In this situation, the stress tensor can also be defined as the following surface integral, which gives usually better values in line integration along the surface S:
This leads to the following formula for calculating the forces:
The reference values for this situation were calculated in  as follows: