# Difference between revisions of "Inline oscillating cylinder by Michael Alletto"

Jozsef Nagy (talk | contribs) |
Jozsef Nagy (talk | contribs) |
||

Line 32: | Line 32: | ||

==Set up== | ==Set up== | ||

− | The | + | The figure below shows the setup of the simulation: At the inflow a uniform velocity of U<sub>inf</sub> = 1 m/s in x-direction is applied. The pressure is set to zero gradient at this boundary. At the top and bottom walls a slip boundary condition is applied. At the outlet the pressure is to be set equal to zero and for the velocity a zero gradient boundary condition is applied. The kinematic viscosity &rho = 0.01 m^2/s is set in order that the Reynolds number Re = U<sub>inf</sub> D / &nu based on the cylinder diameter D = 1 m is equal to Re = 100. For this Reynolds number the flow over a stationary cylinder is two dimensional. For this reason the flow is assumed to be two dimensional. This is imposed by using only one cell in span wise direction and setting both boundary conditions in this direction to the type empty. At the cylinder wall a no slip boundary condition is applied for the velocity and zero gradient for the pressure. The cylinder oscillates in x-direction with a prescribed sinosoidal motion: x<sub>c</sub> = A sin(2\pi f<sub>c</sub> t) = A sin (&omega t). Four different frequencies are computed: f<sub>c</sub> = 0, 0.16, 0.33, 0.66 1/s. The amplitude is set to A = 0.14m. The purpose of this tutorial was to check whether the increase of the lift at f<sub>c</sub> = 2f<sub>st</sub> can be captured and if both techniques to simulate moving bodies deliver comparable results. |

− | |||

− | == | ||

− | |||

− | |||

[[File:MovingMesh_alettoflap.png|900px|center|mesh of the flapping wing]] | [[File:MovingMesh_alettoflap.png|900px|center|mesh of the flapping wing]] | ||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

− | |||

==References== | ==References== |

## Revision as of 05:06, 10 April 2020

**contributor**: Michael Alletto**affiliation**: WRD MS**contact**: click here for email address**OpenFOAM versions**: v1912**Published under**: CC BY-NC-ND license (creative commons licenses)

Go back to Collection by authors.

Go back to Meshing.

Go back to Turbulence.

# 2D flapping wing

You can download the case file OscillatingCylinderSend.tar.gz here.

## Introduction

In this tutorial we will examine the inline oscillatory motion of a cylinder in a uniform ow. Inline means that the cylinder is oscillating parallel to the incoming flow. This kind of configuration is often used to validate solvers aiming to simulate large body motions (see e.g. [1, 4, 5]). The reason is that that despite its simplicity this configuration is associated with a variety of complex physical phenomena.

One of this phenomena are frequency lock in: When the frequency of the oscillating cylinder is close to twice the shedding frequency of the stationary cylinder f_{st}, the frequency of the vortices shed became equal to the one of the cylinder. The amplitude of the cylinder motion which is required to maintain the lock in regime is increasing with departure from twice the shedding frequency of a fixed cylinder [2].

Another interesting phenomena is the increase of the lift coefficient C_{l} when the cylinder is oscillating at 2f_{st} [6]. For frequencies which are either lower or higher than 2f_{st} C_{l} is much lower [6]. [3] made simulations with an in line oscillating cylinder at a Reynolds number of Re = 175. For this Re number the ow around a fixed cylinder is two dimensional. For an oscillation frequency f_{c} of 2f_{st} and an amplitude of 65% of the cylinder diameter the frequency spectrum of the lift force increases considerably. For lower amplitudes of the oscillation only a few peaks are present in the frequency spectrum.
In order to simulate this kind of ow congurations, OpenFOAM offers two different approaches: the overset methodology and the morphing mesh method. The latter can be used if the cylinder motion is not to large to deteriorate the mesh quality. Regarding the overset method we have one or more meshes which displace relative to a background mesh. The coupling between the moving meshes and the background mesh is achieved via interpolation of the
ow quantities. Regarding the morphing mesh, the connectivity between the cells remain unchanged.However the position of the points building a cell can move relative to each other.

Since both methods described above can be used to simulate the flow of an inline oscillating cylinder, both methods can be compared against each other besides the validation against experiments and other simulations. The comparison of this two methods and the validation against experiments is the purpose of this tutorial.

## Set up

The figure below shows the setup of the simulation: At the inflow a uniform velocity of U_{inf} = 1 m/s in x-direction is applied. The pressure is set to zero gradient at this boundary. At the top and bottom walls a slip boundary condition is applied. At the outlet the pressure is to be set equal to zero and for the velocity a zero gradient boundary condition is applied. The kinematic viscosity &rho = 0.01 m^2/s is set in order that the Reynolds number Re = U_{inf} D / &nu based on the cylinder diameter D = 1 m is equal to Re = 100. For this Reynolds number the flow over a stationary cylinder is two dimensional. For this reason the flow is assumed to be two dimensional. This is imposed by using only one cell in span wise direction and setting both boundary conditions in this direction to the type empty. At the cylinder wall a no slip boundary condition is applied for the velocity and zero gradient for the pressure. The cylinder oscillates in x-direction with a prescribed sinosoidal motion: x_{c} = A sin(2\pi f_{c} t) = A sin (&omega t). Four different frequencies are computed: f_{c} = 0, 0.16, 0.33, 0.66 1/s. The amplitude is set to A = 0.14m. The purpose of this tutorial was to check whether the increase of the lift at f_{c} = 2f_{st} can be captured and if both techniques to simulate moving bodies deliver comparable results.

## References

[1] Dominic D Chandar, Bharathi Boppana, and Vasanth Kumar. A comparative study of different overset grid solvers between openfoam, star-ccm+ and ansys-fluent. In 2018 AIAA Aerospace Sciences Meeting, page 1564, 2018.

[2] Owen M Griffn and MS Hall. Vortex shedding lock-on and flow control in bluff body wakes. 1991.

[3] Justin S Leontini, David Lo Jacono, and Mark C Thompson. A numerical study of an inline oscillating cylinder in a free stream. Journal of Fluid Mechanics, 688:551-568, 2011.

[4] Chuan-Chieh Liao, Yu-Wei Chang, Chao-An Lin, and JM McDonough. Simulating flows with moving rigid boundary using immersed-boundary method. Computers & Fluids, 39(1):152-167, 2010. 11

[5] Shen-Wei Su, Ming-Chih Lai, and Chao-An Lin. An immersed boundary technique for simulating complex flows with rigid boundary. Computers & Fluids, 36(2):313-324, 2007.

[6] Y Tanida, A Okajima, and Y Watanabe. Stability of a circular cylinder oscillating in uniform ow or in a wake. Journal of Fluid Mechanics, 61(4):769-784, 1973.

[7] Charles HK Williamson. Vortex dynamics in the cylinder wake. Annual review of Fluid mechanics, 28(1):477-539, 1996.