Models - Mph.cylinder Flow
Models - Mph.cylinder Flow
3a
This model is licensed under the COMSOL Software License Agreement 5.3a.
All trademarks are the property of their respective owners. See www.comsol.com/trademarks.
Introduction
The flow of fluid behind a blunt body such as an automobile is difficult to compute due
to the unsteady flows. The wake behind such a body consists of unordered eddies of all
sizes that create large drag on the body. In contrast, the turbulence in the thin boundary
layers next to the streamlined bodies of aircraft and fish create only weak disturbances of
flow.
An exception to this occurs when you place a slender body at right angles to a slow flow
because the eddies organize. A von Kármán vortex street appears with a predictable
frequency and involves the shedding of eddies from alternating sides. Everyday examples
of this phenomenon include singing telephone wires and an automobile radio antenna
vibrating in an air stream.
Model Definition
To illustrate how you can study such effects, the following model examines unsteady,
incompressible flow past a long cylinder placed in a channel at right angle to the oncoming
fluid. With a symmetric inlet velocity profile, the flow needs some kind of asymmetry to
trigger the vortex production. This can be achieved by placing the cylinder with a small
offset from the center of the flow. In this case, an unstructured mesh is used, and the small
asymmetry in the mesh proves to be enough to trigger the vortex production.
The simulation time necessary for a periodic flow pattern to appear is difficult to predict.
A key predictor is the Reynolds number, which is based on cylinder diameter. For low
values (below 100) the flow is steady. In this simulation, the Reynolds number equals 100,
which gives a developed von Kármán vortex street, but the flow still is not fully turbulent.
The frequency and amplitude of oscillations are stable features, but flow details are
extremely sensitive to perturbations. To gain an appreciation for this sensitivity, you can
compare flow images taken at the same time but with such minor differences as are created
by different tolerances for the time stepping. It is important to note that this sensitivity is
a physical reality and not simply a numerical artifact.
The viscous forces on the cylinder are proportional to the gradient of the velocity field at
the cylinder surface. Evaluating the velocity gradient on the boundary by directly
differentiating the FEM solution is possible but not very accurate. The differentiation
produces 1st-order polynomials when second-order elements are used for the velocity
field. A far better approach is to use a pair of reaction force operators to compute the
integrals of the viscous forces, comparable to second-order accurate integrals of the viscous
forces. An alternative approach would be to use a pair of weak-constraint variables to
enforce the no slip condition. Preferably use the reaction force operator instead of weak
constraints when computing integrals of reaction forces or fluxes in postprocessing.
The drag and lift forces themselves are not as interesting as the dimensionless drag and lift
coefficients. These depend only on the Reynolds number and an object’s shape, not its
size. The coefficients are defined as
2F D
C D = -------------------------
2
-
ρU mean A
2F L
C L = -------------------------
2
-
ρU mean A
Figure 1: A plot of the last time step clearly shows the von Kármán path.
The flow around a cylinder is a common benchmark test for CFD algorithms. Various
research teams have tried their strengths on this problem using different techniques.
Results from some of these experiments have been collected by Schäfer and Turek (Ref. 1),
who also used them to compute a probable value for the “real” answer.
Figure 2 shows how the lift coefficient develops a periodic variation as the von Kármán
vortex structure is formed.
Modeling Instructions
From the File menu, choose New.
NEW
In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D.
2 In the Select Physics tree, select Fluid Flow>Single-Phase Flow>Laminar Flow (spf).
3 Click Add.
4 Click Study.
5 In the Select Study tree, select Preset Studies>Time Dependent.
6 Click Done.
GLOBAL DEFINITIONS
Parameters
1 In the Model Builder window, under Global Definitions click Parameters.
2 In the Settings window for Parameters, locate the Parameters section.
3 In the table, enter the following settings:
Next, create a smoothed step function feature that you will use for ramping up the inflow
velocity.
GEOMETRY 1
Rectangle 1 (r1)
1 On the Geometry toolbar, click Primitives and choose Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type 2.2.
4 In the Height text field, type 0.4.
5 Right-click Rectangle 1 (r1) and choose Build Selected.
Circle 1 (c1)
1 On the Geometry toolbar, click Primitives and choose Circle.
2 In the Settings window for Circle, locate the Position section.
3 In the x text field, type 0.2.
4 In the y text field, type 0.2.
5 Locate the Size and Shape section. In the Radius text field, type 0.05.
6 Right-click Circle 1 (c1) and choose Build Selected.
Difference 1 (dif1)
1 On the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the object r1 only.
3 In the Settings window for Difference, locate the Difference section.
4 Find the Objects to subtract subsection. Select the Active toggle button.
5 Select the object c1 only.
6 On the Geometry toolbar, click Build All.
MATERIALS
Material 1 (mat1)
1 In the Model Builder window, under Component 1 (comp1) right-click Materials and
choose Blank Material.
2 In the Settings window for Material, locate the Material Contents section.
3 In the table, enter the following settings:
Inlet 1
1 On the Physics toolbar, click Boundaries and choose Inlet.
2 Select Boundary 1 only.
To define a parabolic velocity profile use the predefined local curve parameter s. Ramp
up the velocity using the previously defined step function. Append the inverse unit
bracket [1/s] to the time variable t because the step function expects a dimensionless
argument.
Outlet 1
1 On the Physics toolbar, click Boundaries and choose Outlet.
2 Select Boundary 4 only.
MESH 1
1 In the Model Builder window, under Component 1 (comp1) click Mesh 1.
2 In the Settings window for Mesh, locate the Mesh Settings section.
3 From the Element size list, choose Finer.
4 Click Build All.
If you zoom in on the inlet and the cylinder you can see the boundary layers that the
physics-controlled mesh gives.
STUDY 1
RESULTS
Velocity (spf)
Add a Particle Tracing with Mass node to the first default plot group to reproduce the plot
in Figure 1.
To reproduce Figure 2 and Figure 3 of the lift and drag coefficients, first add an Integral
data set for computing the total reaction force on the cylinder.
Selection
1 On the Results toolbar, click More Data Sets and choose Evaluation>Integral.
2 On the Results toolbar, click Selection.
3 In the Settings window for Selection, locate the Geometric Entity Selection section.
4 From the Geometric entity level list, choose Boundary.
5 Select Boundaries 5–8 only.
1D Plot Group 3
1 On the Results toolbar, click 1D Plot Group.
2 In the Settings window for 1D Plot Group, locate the Data section.
3 From the Data set list, choose Integral 1.
Point Graph 1
1 Right-click 1D Plot Group 3 and choose Point Graph.
2 In the Settings window for Point Graph, locate the y-Axis Data section.
3 In the Expression text field, type (-reacf(v)[N]*2/(spf.rho*U_mean^2*
0.1[m^2]))[1/m].
1D Plot Group 4
Finally, visualize the drag coefficient using the following steps:
1 On the Home toolbar, click Add Plot Group and choose 1D Plot Group.
2 In the Settings window for 1D Plot Group, locate the Data section.
3 From the Data set list, choose Integral 1.
Point Graph 1
1 Right-click 1D Plot Group 4 and choose Point Graph.