Week 2 - Flow over a Cylinder.

AIM:

Simulate the flow over a cylinder and explain the phenomenon of Karman vortex street. Understanding the vortex shedding for different Reynolds numbers by changing the inlet velocities accordingly. Create a monitor point behind the cylinder at a distance 4 times the diameter which can be used to calculate and analyze the vortex shedding.

OBJECTIVE:

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

2. 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)

3. Discuss the effect of Reynolds number on the coefficient of drag. [ Results should be validated with any standard literature and error should be within 5 %]

4. In both approaches, show that the flow has converged. That is, what quantities or plots need to be looked at to determine that the flow has converged?

5. Plots of the coefficient of drag and lift.

6. Plots to show vortex shedding behind the cylinder.

7. Mention the preferred material at the end of the challenge submission page.

Introduction:

Vortex shadding:

vortex shedding is an oscillating flow that takes place when a fluid such as air or water flows past a bluff (as opposed to streamlined) body at certain velocities, depending on the size and shape of the body. In this flow, vortices are created at the back of the body and detach periodically from either side of the body forming a Von Karman vortex street. The fluid flow past the object creates alternating low-pressure vortices on the downstream side of the object. The object will tend to move toward the low-pressure zone.

If the bluff structure is not mounted rigidly and the frequency of vortex shedding matches the resonance frequency of the structure, then the structure can begin to resonate, vibrating with harmonic oscillations driven by the energy of the flow. This vibration is the cause for overhead power line wires humming in the wind, and for the fluttering of automobile whip radio antennas at some speeds. Tall chimneys constructed of thin-walled steel tubes can be sufficiently flexible that, in air flow with a speed in the critical range, vortex shedding can drive the chimney into violent oscillations that can damage or destroy the chimney.

Reynolds Number (Re):

Reynolds number for a flow is a measure of the ratio of inertial to viscous forces in the flow of a fluid around a body.

Re = rho*V*d/Dynamic Viscosity

The range of Re values will vary with the size and shape of the body from which the eddies are being shed, as well as with the kinematic viscosity of the fluid. Over a large Re_d range (47<Re_d<10⁵ for circular cylinders; reference length is d: diameter of the circular cylinder) eddies are shed continuously from each side of the circle boundary, forming rows of vortices in its wake. The alternation leads to the core of a vortex in one row being opposite the point midway between two vortex cores in the other row, giving rise to the distinctive pattern shown in the picture. Ultimately, the energy of the vortices is consumed by viscosity as they move further downstream, and the regular pattern disappears.

Strouhal number (St):

Strouhal number is a dimensionless number describing oscillating flow mechanisms.

The Strouhal number is often given as

St = f*d/u

Where,

f is the frequency of vortex shading,

L is the characteristic length and

U is the flow velocity.

3. Co-efficient of drag:

In fluid dynamics, the drag coefficient`(C_d)` 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 a less aerodynamic or hydrodynamic drag. The drag coefficient is always associated with a particular surface area.

The drag coefficient of any object comprises the effects of the two basic contributors to fluid dynamic drag: skin friction and form drag. The drag coefficient of a lifting airfoil or hydrofoil also includes the effects of lift-induced drag. The drag coefficient of a complete structure such as an aircraft also includes the effects of interference drag.

The drag coefficient `(C_d)` is defined as

`C_d=(2*F_d)/(rhou^2A)`

where:

`F_d`is the drag force, which is by definition the force component in the direction of the flow velocity,

`rho` is the mass density of the fluid,

`u`is the flow speed of the object relative to the fluid,

`A` is the reference area.

4. Co-efficient of lift:

The lift coefficient `(C_L)` is a dimensionless coefficient that relates the lift generated by a lifting body to the fluid density around the body, the fluid velocity, and an associated reference area. A lifting body is a foil or a complete foil-bearing body such as a fixed-wing aircraft. C_L is a function of the angle of the body to the flow, its Reynolds number, and its Mach number. The section lift coefficient c_l refers to the dynamic lift characteristics of a two-dimensional foil section, with the reference area replaced by the foil chord.

where

is the lift force,

is the relevant surface area and

q is the fluid dynamic pressure, in turn, linked to the fluid density, and to the flow speed.

A brief explanation to solve the flow over a cylinder

1. Geometry

Creating the geometry by sketching using the line and circle command, circle diameter of 2m and the dimensions having 20X60m and the circle is placed 20m from the inlet.

2.Mesh

After the geometry is loaded we should mesh it and here I used 0.19m as elemental size

But in the meshing nodes, the shape is more quadrilateral than the triangle so we should this 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 body-fitted mesh

And we should name the boundaries as given below

where cylinder edge is named as wall

3.Setting up Ansys fluent

Before opening the set up we should update the mesh. Then we should change the viscous flow to laminar

Open Materials and create a fluid that is to be used for the simulation for this case, it is air and also inputs the density and kinematic viscosity values.

Open Zones and input the inlet velocity as per the reynolds number and check if the fluid is updated.

Update the reference values which is located inside set up at the right-hand side.

Fluent details/procedure:

a. laminar model

b. fluid: Air

c. Method of discretization: Simple

d. Density/kinematic viscosity: 1/0.05

e. Creating a monitor point behind the cylinder by selecting point under surfaces in the results bar.

f. After creating the monitor point, create reports/plots for drag, lift, velocity wrt cylinder.

g. Create plots for the velocity at the monitor point to analyze the vortex shedding.

h. Use hybrid initialization for steady-state and standard for transient state simulations.

Different Cases

(i) Simulate the flow with the steady-state solver for Re=10.

For Re=10 velocity v=0.25m/s

• Velocity and pressure contours

• Residuals plot

• Drag and lift coefficient

• velocity at the monitor point

(ii) Simulate the flow with the steady-state solver for Re=100.

For Re=100 velocity v=2.5m/s

• Velocity and pressure contours

• Residuals plot

• Drag and lift coefficient

• velocity at the monitor point

(iii) Simulate the flow with the steady-state solver for Re= 1000.

For Re=1000 velocity v=25m/s

• Velocity and pressure contours

• Residuals plot

• Drag and lift coefficient

• velocity at the monitor point

(iv) Simulate the flow with the steady-state solver for Re= 10000.

For Re=10000 velocity v=250m/s

• Velocity and pressure contours

• residual plot

• Drag and lift Coefficient

• velocity at the monitor point

(v) Simulate the flow with the steady-state solver for Re= 100000.

For Re=100000 velocity v=2500m/s

• Velocity and pressure contours

• residual plot

• Drag and lift Coefficient

• velocity at the monitor point