Numerical Analysis of Unsteady Flow around a Transversely Oscillating Circular Cylinder

International Journal of Aeronautical and Space Sciences.
2012.
Mar,
13(1):
27-33

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/bync/ 3.0/) which permits unrestricted non-commercial use, distribution and reproduction in any medium, provided the original work is properly cited.

- Received : February 01, 2012
- Accepted : March 16, 2012
- Published : March 30, 2012

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

The relationship between the excitation frequency and the vortex shedding frequency is analyzed during the oscillation of the circular cylinder. Two-dimension unsteady Navier-Stoke’s equation is calculated by using the Optimized High Order Compact (OHOC) scheme. The flow condition is Mach number 0.3 and Reynold’s number 1000. From the results acquired by calculation, it can be inferred that, when the excitation frequency is near the vortex shedding frequency at the fixed cylinder wake, the oscillation frequency of lift and drag coefficients appears to lock-on. The lock-on refers to a phenomenon in which the aerodynamic coefficient appears as one primary oscillation frequency through excitation and its amplitude is amplified. In the non-lock-on zone, the excitation frequency is not in the lock-on mode anymore and beat is formed in which two or more primary oscillation frequencies of the aerodynamic coefficient are mixed together.
Non-dimension parameters are cylinder diameter(D), freestream velocity(U) and freestream density(
ρ
). The characters of t, ξ and η represent time, and a generalized coordinate system. The
Q
is a non- dimensional flux vector that has transformed with the generalized coordinates,
and
represents the nonviscous flux vector and the viscous flux vector in each direction . In order to deal with the oscillating airfoil, a moving mesh had to be used and the mesh movements were taken into consideration for the conversion matrices
[8]
. The term
s
denotes the geometric conservation law and it considers the mesh movement as shown in Equation (2).
In order to maintain high order and high resolution for the time differential, Runge-Kutta’s 4th order scheme was used
[8]
. For the space differential, OHOC (Optimized High- Order Compact), the higher-order and high resolution scheme of 4th-order accuracy, was used
[9]
. The implicit method that uses the following 7 points was used to obtain the numerical differential of each physical property. For the case of α ≠ 0 and β = 0 or α ≠ 0 and β = 0, the tri-diagonal or penta-diagonal matrix was used to obtain the differential. Due to the characteristics of the central difference scheme, there occurred truncation errors and phase errors. In order to minimize such errors, the α, β and an coefficient proposed by Kim and Lee
[9]
was used.
With the OHOC scheme, it is possible to achieve high resolution for space, but it is impossible to exactly reproduce the characteristics of the wave propagation as with the upwind scheme. Due to these characteristics of the central difference scheme, the dissipation and diffusion errors are generated and they significantly affect the numerical stability. In this study, the Adaptive Nonlinear Artificial Dissipation model proposed by kim and Lee
[10]
was used to enhance numerical stability.
For mesh, we used O-type meshes according to the characteristics of the circular cylinder shape and concentrated them around the cylinder where the change of flow was expected to be large. As for the number of meshes, the meshes shown in
Fig. 1
where 250x200 (50,000) meshes were employed and they were used, based on the computational results of the several mesh systems. The moving mesh was used based on the algebraic computation of the meshes that move along with moving cylinder surface within the zone 20 times the diameter. For other zones, fixed meshes were used.
When acoustic waves reach the boundary surface at subsonic flow, physical reflected waves and non-physical reflected waves that occur numerically are generated. While passing through the boundary, the physical reflected waves have long wavelength but they do not affect the numerical stability. However, the non-physical reflected waves have short wavelength and they generate unnecessary numerical waves due to the numerical differentiation at the boundary. This as a result causes numerical instability in the entire computational domain. Particularly, when the turbulent flow or vortexes pass through the boundary, highly nonphysical reflected waves occur. Generally, in order to inhibit non-physical reflection that occurs at the boundary, inflow and outflow conditions according to the eigenvalues are used. However, the characteristic boundary conditions are not sufficient enough to inhibit non-physical reflection. In this study, the zonal characteristic boundary conditions proposed by Sandberg and Sandham
[11]
were applied to minimize the effect of non-physical reflected waves on the
Computaional grids around a circular cylinder
computational domain. The zonal characteristic boundary condition was developed on the basis of the general local characteristic boundary conditions, and it is easily applicable to the general characteristic boundary conditions.
S_{f}
) and amplitude for the flow conditions of
R_{e∞}
= 1000 and
M
= 0.3. The amplitudes of oscillation were set to 0.1D, 0.2D and 0.3D respectively. The frequency in the range of 0.16 ≤
S_{f}
≤ 0.25 was applied. Before analyzing the oscillating cylinder, a fixed cylinder was analyzed to compare the results with other results.
By analyzing the flow around a stationary cylinder, the natural Karman vortex shedding frequency, that is a dimensionless Strouhal number (
S_{S}
), was computed by using the FFT analysis of the lift coefficient.
With respect to the changes in the lift and drag coefficients over time, it is seen that a periodic steady state was maintained and it is as shown in
Fig. 2
(a).
Fig. 2
(b) shows the
(a) Lift and drag coefficients
(b) PSD of lift coefficients
FFT analysis of the lift coefficient and it can be seen that the PSD was strong only at one point of the frequency. For the log (PSD), the
S_{S}
number is clearly seen without any other frequency disturbance. The result of
S_{S}
= 0.22 is similar to
S_{S}
= 0.21 by Blake
[12]
,
S_{S}
= 0.206 by Williamson
[13]
and
S_{S}
= 0.200 by Jordan and Fromm
[14]
.
The excitation frequency of a circular cylinder was applied in the range of 0.16 ≤
S_{f}
≤ 0.25 for the amplitude of 0.2D in order to check the lock-on zone on the basis of the Strouhal
(a) Lift and drag coefficients for S_{f} = 0.16
(b) PSD of the lift coefficients for S_{f} = 0.16
(c) Lift and drag coefficients for S_{f} = 0.21
number of the flow that is around the stationary cylinder.
Fig. 3
(a),
Fig. 3
(b),
(e)
and
(f)
show the non-lock-on phenomenon at an amplitude of 0.2D and the excitation frequencies were 0.16 and 0.24. It shows that two or more frequencies are mixed when the lift and the drag coefficients change and that there is a beat between the frequencies.
Figs. 3
(c) and
(d)
show the lock-on at the amplitude of 0.2 and the excitation frequency of 0.22. This is when Strouhal number of the cylinder agrees along with the excitation frequency.
(d) PSD of lift coefficients for S_{f} = 0.21
(e) Lift and drag coefficients for S_{f} = 0.24
(f) PSD of lift coefficients for S_{f} = 0.24
(a) Lift and drag coefficients for S_{f} = 0.20
(b) PSD of the lift coefficients for S_{f} = 0.20
(c) Lift and drag coefficients for S_{f} = 0.22
(d) PSD of lift coefficients for S_{f} = 0.22
When the lift and drag coefficients change, one primary oscillation frequency appears without the interference of other frequencies. It appears that the extreme value of PSD becomes very strong at one point frequency and log (PSD) also becomes strong at one point frequency without any other frequency disturbance. As shown in the PSD value, one strong extreme value happens to be larger compared to that at the fixed cylinder. In other words, the maximum lift coefficient also increases and this is due to the characteristic Strouhal frequency of the cylinder. The excitation frequency creates resonance and thus it increases the force applied to the cylinder.
Fig. 4
(a) ~
(d)
show the non-lock-on case of frequency of
S_{f}
=0.20 and the lock-on case of
S_{f}
=0.22 with an amplitude of 0.1D, respectively.
Figs. 5
(a) and (b) show that the lock-on phenomenon. In this the vortex shedding is shown during a cycle with the amplitude of 0.2D and the excitation frequency of 0.21. In
Fig.5
of the case of lock-on, it is shown that the distance between vortices is regular respect to resonance. However, for the case of non-lock-on, the distance slightly changes and this is due to the influence of beat frequency. In order to show the irregularity between vortices that is caused by
Pressure contour when the amplitude was 0.2D and S_{f} = 0.21
the beat frequency, the contour lines of pressure distribution during a cycle with the amplitude of 0.3 and the excitation frequencies of 0.23 are presented in the
Figs. 6
(a) and (b). The figures make the irregularity of the vortex distribution clearer.
Fig. 7
shows that the lock-on zone and the non-lockon zone according to the amplitude and the excitation frequency.
Fig.7
t the lock-on zone becomes wider as the amplitude becomes larger. This phenomenon appears as the cylinder oscillation affects the fluid flow more, as the amplitude becomes larger.
Pressure contour when the amplitude was 0.3D and S_{f} = 0.23
Lock-on zone dependent on the amplitude and the excitation frequency

Optimized High Order Compact Scheme
;
Transversely Oscillating Circular Cylinder
;
Lock-on and non-Lock-on
;
Vortex Shedding

1. Introduction

The analysis of the flow around the moving object accounts for a significant part in the recent computational fluid dynamics. Even though, numerous studies have been carried out until now, there is still much to be identified on the flow around a circular cylinder. Despite its simple shape, a circular cylinder is known to cause complex flow phenomena such as separation flow and vortex shedding. This is due to its large pressure gradient which makes numerical computation and interpretation difficult. As fluid flow around a blunt body is similar to that around a circular cylinder, analysis of flow around a circular cylinder is the basis for analyzing flow around blunt bodies. Structural oscillation is one of the important factors associated with the service life of a structure. Moreover, aero-acoustic noise has a great impact on the noise environment around a structure. Therefore, many people have attempted various methods to analyze and control the cause of vortex shedding. A method for controlling vortex shedding is to attach a plate at the back of a circular cylinder or to control wakes by using another cylinder, to disturb the inflow or absorb and spout fluid on the cylinder surface, to periodically oscillate the cylinder, etc. Accordingly, flow around an oscillating circular cylinder has been studied and it is yet an object of important numerical computation
[1
,
2
,
3]
.
The most important phenomenon for analysis of the flow around an oscillating object is the “lock-on.” The lock-on is a phenomenon in which the oscillation of an object coincides with the vortex shedding that is generated from the object. The circular cylinder has its characteristic Karman vortex shedding frequency which is dependent on the Reynold’s number. The lock-on happens when the oscillation frequency of the cylinder is approximately same or similar to the Karman vortex shedding frequency
[1
,
2
,
3]
.
Many researchers have carried out numerous numerical or experimental studies for the fluid flow around the oscillating cylinder. Williamson and Roshko
[1]
carried out excitation oscillation test in order to analyze the wakes that are generated depending on the changes in oscillation amplitudes and oscillating frequencies. Kim and Yoo
[2]
, measured the changes in dynamic properties of wakes behind a circular cylinder by using a high-speed PIV system to compare and examine it with natural shedding. Atluri, et al
[3]
. Examined wake flow of a circular cylinder oscillation which is dependent on the Reynold’s Numbers and frequencies. Koopman
[4]
examined the lock-on phenomenon for a circular cylinder oscillating with low Reynold’s numbers. Gerrand
[5]
examined the vortex shedding mechanism that occurs in a fixed blunt body and its various patterns. Griffin and Ramberg
[6]
examined the transverse oscillation effect in a lock-on zone and observed horizontal vortex intervals by using the visualization technique. Kang, et al
[7]
. numerically computed the flow of a circular cylinder oscillation by using the DNS of the central difference scheme of 4th order accuracy and the upwind scheme of 7th order accuracy.
For this study, we identified the lock-on zone according to the changes of the transverse amplitude and the excitation frequency of the cylinder, and carried out a comparative analysis for flow phenomena in the lock-on and the non lock-on conditions. In other words, we analyzed the correlation between the excitation frequency of the circular cylinder which oscillates in the transverse direction and the vortex shedding frequency that is generated around the cylinder. As for the governing equation, we used two-dimensional unsteady compressible Navier-Stoke’s equation but we did not use the turbulence model for the computation of low Reynold’s number. We used the fourthorder accuracy Runge-Kutta’s scheme for time differential and the Optimized High-Order Compact (OHOC) scheme with the fourth-order accuracy for the space differential. For the zone, we used the O-type moving mesh up to 40 times the diameter of the cylinder. The flow condition was set to and the transverse amplitudes were set to 0.1, 0.2 and 0.3 times the diameter respectively. The excitation frequency was applied at the zone of . By performing PSD analysis for the lift coefficient for each flow condition, we analyzed the characteristics of the primary oscillation frequency and resonance frequency. When the excitation frequency is near the natural vortex shedding frequency of the fixed cylinder, the lift coefficient and the vortex shedding were amplified by the excitation oscillation, and the the one primary oscillation frequency appeared in the lock-on zone. In the non-lock-on zone beyond the zone, the beat was formed in which two or more primary oscillation frequencies were mixed together. Furthermore, as the amplitude increased, the frequency zone where lock-on appeared got wider.
2. Governing Equation and Numerical Method

The 2-dimensional unsteady compressible Navier- Stokes equation non-dimensionalized in the generalized coordinates was used, as in Equation (1).
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

3. Results and discussion

In this study, the phenomenon of flow was analyzed depending on the magnitude of excitation frequency (
Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Lager Image

Conclusion

This study numerically analyzed the flow around the two-dimensional circular cylinder which oscillated in the transverse direction. As for the governing equation, we used the two-dimensional unsteady compressible Navier-Stokes equation. For the space difference, we used OHOC scheme to ensure high-order and high-resolution accuracy. For the time difference, we used the 4th order Runge-Kutta’s scheme. Moreover, in order to suppress the nonlinear discontinuous wave that causes numerical instability, ANAD scheme was used for artificial dissipation. As for the boundary condition, the non-reflecting zonal characteristic boundary condition was used. This suppressed the unphysical reflected wave and the reflected wave that is generated by the vortex with high energy density when it passed through the boundary. The cylinder had the characteristic vortex shedding frequency and this was dependant on the Reynold’s number. When the oscillation frequency of the cylinder was excited at the characteristic vortex shedding frequency of the cylinder or at a similar range, the cylinder went into the lock-on condition. In this, one primary oscillation frequency was generated and the resonance was created to increase the amplitude of the lift. In the non-lock-on zone, two or more primary oscillation frequencies were generated by the difference between the excitation frequency and the natural vortex shedding frequency. There was a beat among these frequencies. It was confirmed that the lock-on zone increased as the excitation amplitude increased.
Acknowledgements

This study was supported by research fund from Chosun University, 2011 .

Williamson C.H.K
,
Roshko A.
"Vortex Formation in the wake of an Oscillating Cylinder"
Journal of Fluids and Structures
2
(355)
355 -
381
** DOI : 10.1016/S0889-9746(88)90058-8**

Kim W
,
Yoo J. Y.
2003
"Change of Vortex Dynamics in the Cylinder Wake by the Lock on to Oscillatory Incident Flow"
Korean Society of Mechanical Engineers
B, 27
(11)
1645 -
1654
** DOI : 10.3795/KSME-B.2003.27.11.1645**

Atluri S.
,
Rao V.K.
,
Dalton C.
2008
"A numerical investigation of the near-wake structure in the variable frequency excitation oscillation of a circular cylinder"
Journal of Fluids and Structures Fluids and Structures
25
229 -
244

Koopmann G.H.
1967
"The vortex wakes of vibrating cylinders at low Reynolds numbers"
J Fluid Mech
28
501 -
512
** DOI : 10.1017/S0022112067002253**

Gerrand J.H.
1966
"The mechanics of the formation region of the vortices behind bluff bodies"
J. Fluid Mech
25
401 -
413
** DOI : 10.1017/S0022112066001721**

Griffin O.M.
,
Ramberg S.E.
1976
"Vortex shedding from a circular cylinder vibrating inline with an incident uniform flow"
J Fluid Mech
75
257 -
276
** DOI : 10.1017/S0022112076000207**

Kang S.J.
,
Tanahashi M.
,
Miyuchi T.
,
Lee Y.H.
2001
"Direct Numerical Simulation of the Flow past an Oscillating Circular Cylinder"
Korean Society of Computational Fluids Engineering
6
(4)
26 -
34

Hoffmann K.C.
,
Chiang S.T.
1993
"Computational Fluid Dynamics for Engineers"
Engineering Education System USA
1
26 -
34

Kim J.W.
,
Lee D.J.
1996
"Optimized compact finite difference schemes with maximum resolution"
AIAA
34
(5)
887 -
893
** DOI : 10.2514/3.13164**

Kim J.W.
,
Lee D.J.
1999
"Adaptive nonlinear artificial dissipation model for computational aeroacoustics"
3rd CAA Workshop on Benchmark Problems
USA

Sandberg R.D.
,
Sandham N.D.
2006
"Nonreflecting zonal characteristic boundary condition for direct numerical simulation of aerodynamic sound"
AIAA J.
44
(2)
402 -
405
** DOI : 10.2514/1.19169**

Blake W.K.
1986
"Dipole sound from Cylinders", Mechanics of flow induced sound and vibration, vol.1
1st ed.
Academic Press
New York
219 -
287

Williamson C.H.K.
1996
"Three dimensional wake transition behind a cylinder"
Journal of Fluid Mechanics
328
345 -
407
** DOI : 10.1017/S0022112096008750**

Jordan S.K.
,
Fromm J.E.
1972
"Oscillatory Drag, Lift and Torque on a circular cylinder in a uniform flow"
Phisics of Fluids
15
(3)
371 -
376
** DOI : 10.1063/1.1693918**

Citing 'Numerical Analysis of Unsteady Flow around a Transversely Oscillating Circular Cylinder
'

@article{ HGJHC0_2012_v13n1_27}
,title={Numerical Analysis of Unsteady Flow around a Transversely Oscillating Circular Cylinder}
,volume={1}
, url={http://dx.doi.org/10.5139/IJASS.2012.13.1.27}, DOI={10.5139/IJASS.2012.13.1.27}
, number= {1}
, journal={International Journal of Aeronautical and Space Sciences}
, publisher={The Korean Society for Aeronautical & Space Sciences}
, author={Moon, Ji-soo
and
Kim, Jae-soo}
, year={2012}
, month={Mar}