Challenge

Flow over a Cylinder(Ansys Fluent)

Objective: Simulate the flow over a cylinder and explain the phenomenon of Karman vortex street.

PART-I

  1. Simulate the flow with the steady and unsteady case and calculate the Strouhal Number for Re= 100.

PART-II

  1. Calculate the coefficient of drag and lift over a cylinder by setting the Reynolds number to 10,100,1000,10000 & 100000. (Run with steady solver)
  2. Discuss the effect of Reynolds number on the coefficient of drag.

Introduction:

Karman Vortex Street

In Fluid dynamic, a Karman vortex street (or a von Karman vortex street) is a repeating pattern of swirling vortices, caused by a process known as Vortex shedding, which is responsible for the unsteady separation of flow of a fluid around blunt bodies.

A vortex street will form only at a certain range of flow velocities, specified by a range of Reynolds numbers (Re), typically above a limiting Re value of about 90. The (global) Reynolds number for a flow is a measure of the ratio of inertia to viscous force in the flow of a fluid around a body or in a channel, and may be defined as a non-dimensional parameter of the global speed of the whole fluid flow:

where:

= the free stream flow speed(i.e. the flow speed far from the fluid boundaries like the body speed relative to the fluid at rest, or an inviscid flow speed, computed through the Bernoulli equation), which is the original global flow parameter, i.e. the target to be non-dimensionalized.

= a characteristic length parameter of the body or channel

= the free stream kinematic viscosity parameter of the fluid, which in turn is the ratio:

Strouhal number:
In dimensional analysis, a strouhal number is a dimensionless number describing oscillating flow mechanics.
For a large Strouhal number (order of 1), viscosity dominates fluid flow, resulting in a collective oscillating movement of the fluid. For the low Strouhal number (order of `10^-4`and below), the high speed, quasi-steady state portion of the movement dominates the oscillation. Oscillation for intermediate Strouhal number is characterized by the buildup and rapidly subsequent shedding of vortices.
The Strouhal number is given by;
`S_t= (f D)/U^'`
where,
f= frequency of vortex shedding
D= characteristic diameter of the cylinder
`U^'`= fluid flow velocity
Drag Co-efficient:
In fluid dynamics, the drag coefficient is a dimensionless quantity that is used to quantify the drag or resistance of an object in a fluid environment, such as air or water. It is used in the drag equation in which a lower drag coefficient indicates the object will have less aerodynamics or hydrodynamic drag. The drag coefficient is always associated with a particular surface area.
The drag coefficient is given by;
`C_d= (2*F_d)/(rho u^2 A`
where,
`F_d` is the drag force, which is by definition the force is acting opposite to the relative motion of any object moving with respect to a surrounding fluid.
`rho` is the mass density of the fluid
u is the flow speed of the object relative to the fluid.
A is the area
Lift Co-efficient:
The lift coefficient is a dimensionless quantity that relates the lift generated by the lifting body to the fluid density around the body, the fluid velocity, and associated reference area.
It is given by;
`C_L=(L)/(qS)=(L)/(1/2rho u^2S)`
where,
L is the lift force
S is the relevant surface area
q is the fluid dynamic pressure. in turn, linked to fluid density to the flow speed.
Procedure:
1. Geometry: Created the geometry by using the line and circle command. Circle diameter of 2m and dimension having 20*60m and the circle is placed 20m from the inlet.
2. Mesh: After geometry is loaded we will mesh it.
  • But in the meshing nodes, the shape is more quadrilateral than the triangle so we need to change it by changing the method to the triangle in the insert option.
  • After changing the node's shape into triangles we should improve the mesh around the cylinder by sizing in the mesh which is used to control precisely the mesh.
  • Then we should add inflation which is used to create a body-fitted mesh.

3. Setting up Ansys Fluent:

Fluent details/Procedure

  1. Laminar model
  2. Fluid: Air
  3. Method of discretization: SIMPLE
  4. Density/kinematic viscosity: 1/0.05
  5. Creating a monitor point behind the cylinder by selecting points under surfaces.
  6. After creating the monitor point, create reports/plots for drag, lift velocity w.r.t cylinder.
  7. Create plots for the velocity at the monitor point to analyze the vortex shedding.
  8. Use hybrid initialization for steady-state and standard for transient state simulations.

PART-I

Simulate the flow with the steady and unsteady case and calculate the Strouhal Number for Re= 100 and inlet velocity 2.5 m/s.

Steady-state:

  • Velocity and Pressure Contour

  • Residual plot

  • Drag and Lift coefficient

  • Velocity at the monitor point

  • Strouhal Number plot

Unsteady-state:

  • Velocity and Pressure Contour

  • Residual plot

  • Drag and Lift coefficient

  • Velocity at the monitor point

  • Strouhal Number plot

PART-II

a. Simulating the flow with a steady-state solver for Re=10 and inlet velocity=0.25 m/s.

  • Velocity and Pressure Contour

  • Residual plot

  • Drag and Lift coefficient

  • Velocity at the monitor point

b. Simulating the flow with a steady-state solver for Re=100 and inlet velocity=2.5 m/s.

  • Velocity and Pressure Contour

  • Residual plot

  • Drag and Lift coefficient

  • Velocity at the monitor point

c. Simulating the flow with a steady-state solver for Re=1000 and inlet velocity=25 m/s.

  • Velocity and Pressure Contour

  • Residual plot

  • Drag and Lift coefficient

  • Velocity at the monitor point

d. Simulating the flow with a steady-state solver for Re=10000 and inlet velocity=250 m/s.

  • Velocity and Pressure Contour

  • Residual plot

  • Drag and Lift coefficient

  • Velocity at the monitor point

e. Simulating the flow with a steady-state solver for Re=100000 and inlet velocity=2500 m/s.

  • Velocity and Pressure Contour

  • Residual plot

  • Velocity at the monitor point

Conclusion:

  • In PART-I, the flow over a cylinder is simulated for both Steady-state and unsteady state with Re=100, and with results generated the phenomenon of Von Karman vortex street was observed.
  • Necessary plots were plotted in order to observe that flow is converged and the vortex shedding is happening.
  • In PART-II, five different simulations (i.e Re=10,100,1000,10000 & 100000) were simulated to observe the behaviour of drag and lift coefficients by keeping other variable constant.
  • As we refine the mesh, we need to run more iteration in order to get results.

References:

  1. https://skill-lync.com/knowledgebase/how-to-create-an-animation-in-ansys-fluent
  2. https://skill-lync.com/knowledgebase/how-to-perform-a-fft-using-ansys-fluent
  3. https://skill-lync.com/knowledgebase/what-are-the-types-of-dimensionless-numbers
Your portfolio, next

Ship work like this
in your own course.

Talk to a career counsellor

Free 20-minute call. We'll map the right course to where you want to work.

By submitting you agree to our Terms and Privacy Policy.

Talk to a career counsellor

Free 20-minute call. We'll map the right course to where you want to work.

By submitting you agree to our Terms and Privacy Policy.