Numerical simulation of vortex induced vibrations on a circular cylinder at different Reynold’s number
Shivam Yadav1 , Akshaj Kulshreshtha2 , Baij Nath Singh3
1, 2Bennett University, Greater Noida, India
3Indian Institute of Technology Dhanbad, Dhanbad- 826004, Jharkhand, India
Vibroengineering PROCEDIA, Vol. 29, 2019, p. 171-175.
Received 24 October 2019; accepted 1 November 2019; published 28 November 2019
43rd International Conference on Vibroengineering in Greater Noida (Delhi), India, November 28-30, 2019
This research investigates the effect of vortex induced vibrations on flow past a circular cylinder for two-dimensional unsteady incompressible flow at different Reynold’s number. The pressure bases steady solver is used for computation along with standard k-ε turbulence model. The change in the lift and drag coefficient with respect to increase in Reynolds number is studied and contours of vorticity are plotted. The pressure distribution on the fixed cylinder for different Reynolds number is also presented. It is found that drag coefficient reduces with the increasing Re and lift coefficient increases up to Reynold’s number 104. Moreover, the pressure difference on the fixed cylinder increases with the increasing Reynold’s number.
Keywords: vortex induced vibration (VIV), wake flow, standard k-ε turbulence model.
The flow past circular cylinder is gaining attention of researchers due to its primitive significance in the engineering applications. Steady flow past a circular cylinder gives the various flow regimes at different Reynolds numbers and coefficient of drag rely upon the Reynolds number .
In his research, Xiong elaborates that Vortex induced vibrations generates when a structure is placed normally to the direction of flow in a Newtonian fluid so that it can fluctuate due to separation in vortices shedding and both natural and shedding frequency becomes equal .
Vortex shedding can successfully lead to separation of delay, drag reduction, and reduction of vibration and noise. There are many features which results from the interactions between the external disturbance, and the ﬂow around the circular cylinder.
Vortex shedding is an oscillating flow of a fluid like air, and water around the bluff bodies at Reynold’s number based on the characteristic length of the body. On the bluff bodies, the flow separation takes place to causes the pressure difference on different surface of the body. The vortices detach from the body periodically producing the Von-Karman street. Results computed from the experimental studies are most reliable, but due to the cost limitations, nowadays numerical methods are more preferred as they also provide reliable, and acceptable results .
Fluid passing through the model separates the flow at one or more than one sharp corner and makes the contours and as the Reynolds number increases the substantial difficulties starts occurring .
By using laminar separation bubble, and turbulent model, an increase in a Reynolds number we found that there is a rapid reduction in drag force and increases in lift force .
By doing experiment on a circular cylinder the investigation of vortex induced vibrations with Low damping ratio gives the idea that there is a change in velocity .
To do the computational simulation first order SAS-SST turbulent model has been used on a 2-D cylinder near different wake regions and found that shape of mean velocity is related to fluctuation in velocities .
In this paper, we consider the case of flow past a circular cylinder at different Reynolds number. Our research work includes the computational work. The simulations carried out for 2-D Circular Cylinder only. Thus, the present work deals with the lift and drag coefficient in comparison to increase in Reynolds number.
2. Governing equations
Two major governing equations for 2-D incompressible turbulent laminar flow we have used.
Continuity equation is given in Eq. (1):
Navier stokes equations are given in Eq. (2) and Eq. (3).
For -direction component:
For -direction component:
where is the density, is the free stream velocity, is the pressure, is the dynamic viscosity of fluid, and are the velocities in and directions respectively.
2.1. Quantities of interest
is the Reynold’s number, and are the lift and drag coefficient. is the pressure coefficient and is the skin friction coefficient. is the diameter of cylinder, and are the lift and drag force respectively.
3. 3.Computational domain and meshing
The numerical computation of drag and lift coefficient for different Reynold’s number is studied using Ansys Fluent 19.1. The 2-D cylinder of diameter 5 mm is formed and computational domain is created around it having dimensions of 150 mm×70 mm. The fluid domain with boundary conditions is given in Fig. 1. The inlet is assigned as velocity inlet and outlet is having atmospheric pressure i.e. 101325 Pa and all other walls are having no slip condition. The structured 2-D quadrilateral mesh is used for the domain and inflation layers are formed around the cylinder. The mesh with inflation layers is shown in Fig. 2.
The steady pressure-based solver is used to study the aerodynamic effects around the cylinder. The turbulence model used is standard - and SIMPLE scheme is used for pressure and velocity coupling with second order upwind for momentum.
Fig. 1. Domain around the cylinder
Fig. 2. Meshing around the cylinder with inflation layers
4. Results and discussion
The numerical simulation around cylinder is carried out by varying the velocity at inlet to get different Reynold’s number and to analyze the change in lift, drag and pressure coefficient with respect to . Table 1 shows the drag coefficient and lift coefficient values for different Reynold’s number. The plot for the averaged and along the increasing is shown in Fig. 3 and Fig. 4 respectively. It is observed that the drag coefficient decreases on increasing the while lift coefficient increases up to 104 and drops for higher .
Table 1. Coefficient of lift and drag coefficient for different Reynold’s number
Coefficient of drag ()
Coefficient of lift ()
Fig. 3. Coeffcient of drag vs Reynold’s number
Fig. 4. Coeffcient of lift vs Reynold’s number
Fig. 5 represents the coefficient of lift and pressure distribution on fixed cylinder for at 102 ,103 ,104 and 105. The is plotted for the iteration ranging from 800 to 1800. The value increasing as increases till 104 while decreases for further higher values which can be seen from the graph. The value of pressure coefficient is having overlapped curve for at 102 and the difference between pressure distribution increases along the increasing .
Fig. 5. The lift coefficient and pressure coefficient distribution along the fixed cylinder wall at different Reynold’s number
The vorticity contour is plotted in Fig. 6 to check the wake pattern on the cylinder for increasing . The range of vortex shedding is kept from –104 to 104. The oscillations around the cylinder increases with Reynold’s number forming von Karman street around the cylinder. It can be seen that separate shear layers are rolling on the upper and lower wall of fixed cylinder. The red color in the vorticity contour represents the clockwise rotation while blue color represents the counterclockwise rotation of vortices.
To understand the oscillations around the cylinder due to formation of vortex shedding the power spectral density of the lift coefficient is determined using Fast Fourier Transformation. The maximum frequency for the 105 using FFT is 0.018 Hz. Fig. 7 represents the power spectral distribution for the same.
Fig. 6. Velocity contours at different Reynolds numbers
Fig. 7. Power spectral distribution vs frequency
In this paper the numerical computation for the lift and drag coefficient with increasing number of Reynold’s number is represented. The drag coefficient reduces with the increasing and lift coefficient increases up to Reynold’s number 104. Similarly, the pressure difference on the fixed cylinder increases with the increasing Reynold’s number. The vorticity contour and power spectral distribution of lift coefficient is also presented in the results and discussion section.
- Martinez R. C., Sweeney L. G., Finlay W. H. Aerodynamic forces and moment on a sphere or cylinder attached to a wall in a Blasius boundary layer. Engineering Applications of Computational Fluid Mechanics, Vol. 3, Issue 3, 2009, p. 289-295. [Publisher]
- Xiong Y., Peng S., Zhang M., Yang D. Numerical study on the vortex-induced vibration of a circular cylinder in viscoelastic fluids. Journal of Non-Newtonian Fluid Mechanics, Vol. 272, 2019, p. 104170. [Publisher]
- Singha S., Sinhamahapatra K. Flow past a circular cylinder between parallel walls at low Reynolds numbers. Ocean Engineering, Vol. 37, Issues 8-9, 2010, p. 757-769. [Publisher]
- Rajani B. N., Kandasamy A., Majumdar S. Numerical simulation of laminar flow past a circular cylinder. Applied Mathematical Modelling, Vol. 33, Issue 3, 2009, p. 1228-1247. [Publisher]
- Chopra G., Mittal S. Numerical simulations of flow past a circular cylinder. Journal of Physics: Conference Series, Vol. 822, 2017, p. 012019. [Publisher]
- Brika D., Laneville A. Vortex-induced vibrations of a long flexible circular cylinder. Journal of Fluid Mechanics, Vol. 250, 1993, p. 1-481. [Publisher]
- Shim Y. M., Sharma R., Richards P. Numerical study of the flow over a circular cylinder in the near wake at Reynolds number 3900. 39th AIAA Fluid Dynamics Conference, 2009. [Publisher]
Lecture Notes in Mechanical EngineeringAnuj Yadav, N. K. Singh