CAE Software【Femtet】Murata Software Co., Ltd.
Example3 Flow around Cylinder
General
-
The steady-state analysis is applied to the flow around the cylinder.
-
The fluid velocity distribution, the fluid velocity vectors, the streamlines, and the force on the wall are solved.
-
Unless specified in the list below, the default conditions will be applied.
-
Obtain this session's project file. (Save the project file before open)
Analysis Space
|
Item |
Setting |
|
Analysis Space |
2D |
|
Model Unit |
mm |
Show Results
|
Item |
Tab |
Setting |
|
Solver |
Solver selection |
Fluid analysis [Bernoulli] |
|
Analysis Type |
Fluid Analysis |
Steady-state analysis |
|
Laminar Flow/Turbulent Flow |
Fluid Analysis |
Select Laminar Flow |
|
Layer Mesh Setting for Wall Surface (General Settings) |
Fluid Analysis |
Specify mesh height of 1st layer Height of 1st layer mesh: 0.5 [mm] |
|
Meshing Setup |
Mesh |
General mesh size: 10[mm] |
Model
The material of Air (000_Air) is set to a rectangular sheet body. The boundary conditions of inlet and outlet are set on the left edge and the right edge respectively.
The slip wall outer boundary condition is applied to the top and bottom edges where the boundary condition is not set.
The model is a circle sheet body and the material is iron (007_Fe).
The edges surrounding a circle is a boundary of solid and fluid. Solid wall is automatically set to them.
Body Attributes and Materials
|
Body Number/Type |
Body Attribute Name |
Material Name |
|
0/Solid |
air |
000_Air(*) |
|
1/Solid |
column |
007_Fe * |
* Available from the material DB
Boundary Condition
|
Boundary Condition Name/Topology |
Tab |
Boundary Condition Type |
Setting |
|
Inlet/Edge |
Fluid |
Inlet |
Forced Inflow |
|
Outlet/Face |
Fluid |
Outlet |
Natural Outflow |
|
Outer Boundary Condition |
Fluid |
Slip wall |
– |
If Reynolds number exceeds 100, the time dependency of the turbulence occurs. It makes calculation difficult in the steady-state analysis.
The Reynolds number calculated from this model form, material property, and fluid velocity is about 26.5. The steady-state analysis of the laminar flow can be executed.
Viscosity μ=1.82e-5[Pa s]
Density ρ: 1.205[kg/m3]
Kinematic viscosity v=μ/ρ=1.82e-5/1.205=1.510e-5[m2/s]
Fluid velocity V=0.01[m/s]
Diameter of cylinder D=0.04[m]
Reynolds number Re = V*L/ν=0.01*0.04/1.510e-5 = 26.5
Results
The vectors of the fluid velocity distribution around the cylinder are shown below.
The vectors are adjusted to the same length in the graphics setup.
It is known that the so-called twin whirlpools occur behind the cylinder when the Reynolds number is around 30.
Two whirlpools can be observed as below.
The streamlines are shown below.
The pitch of the startpoints is adjusted in the graphics setup.
The streamlines are generated.
The force on the wall is shown by the table.
[RESERVED_inner_fluid_007_Fe] is the force applied to the cylinder by the fluid.
Drag F:The force in the flow direction (X-direction) is 6.698 x 10 ^ -9 [N].
The drag coefficient CD can be calculated from the drag F.
Thickness in depth direction t = 1[mm]
Cross sectional area S = D * t = 0.02 * 0.001 = 4e-5 [m2]
Dynamic pressure Pk = 0.5 * ρ*V^2 = 0.5 * 1.205 * 0.01 * 0.01 = 6.025e-5 [Pa]
Drag coefficient CD = F / Pk / S = 6.698e-9 / 6.025e-5 /2e-5 = 2.78