Reynolds and froude number effect on the flow past an interface-piercing circular cylinder

International Journal of Naval Architecture and Ocean Engineering.
2014.
Sep,
6(3):
529-561

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

- Published : September 30, 2014

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

Large-eddy simulation
;
Interface-piercing circular cylinder
;
Run-up
;
Wave breaking
;
Vortical structures

INTRODUCTION

The turbulent flow past a circular cylinder has been investigated extensively for a long time due to its importance in many engineering applications. It shows different features at different Reynolds numbers (
COMPUTATIONAL METHODS

CFDShip-Iowa version 6.2, a sharp interface orthogonal curvilinear grid solver for two-phase incompressible flows recently developed at IIHR (
Yang and Stern, 2009
;
Suh et al., 2011
) is used for all computations. The interface is tracked by a Coupled Level Set and Volume-of-Fluid (CLSVOF) method (
Wang et al., 2009
;
Wang et al., 2012
) instead of the level set method used in
Suh et al. (2011)
to capture detailed breaking waves.
- Mathematical model

An orthogonal curvilinear coordinate system is adopted for the governing equations. The Navier-Stokes equations for twophase, immiscible, incompressible flows have been written in the orthogonal curvilinear coordinates. The continuity equation is given as follows:
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- Numerical method

The finite-difference method is used to discretize the Navier-Stokes equations on a non-uniform orthogonal curvilinear grid. A staggered variable arrangement is adopted, i.e., the contravariant velocity components
SIMULATION CONDITIONS, VERIFICATION AND VALIDATION

- Grid, computational domain, and geometry

Body-fitted cylindrical grids were used for all cases and
Table 1
shows the number of grid points and distance of the first grid point from the cylinder wall (
Simulation conditions.

PPT Slide

Lager Image

PPT Slide

Lager Image

- Flow conditions, boundary conditions, and initial conditions

The diameter of cylinder,
PPT Slide

Lager Image

PPT Slide

Lager Image

- Verification and validation of integral variables

The time histories of the drag coefficient
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

Verification study with caseRe= 458,000 /Fr= 1.64

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

Iterative uncertainty using drag coefficient (CD).

PPT Slide

Lager Image

Validation study with caseRe= 458,000 /Fr= 1.64.

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- Flow field validation forRe=234,000/Fr=0.84

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

REYNOLDS NUMBER EFFECT WITH FIXED FROUDE NUMBER = 0.84

Three different Reynolds numbers from sub-critical to post-critical (27,000, 234,000, and 458,000) with
- Mean separation pattern

Fig. 13
shows the separation pattern of the mean flow with vortex core lines, which is obtained using the approach discussed in
Kandasamy et al. (2009)
and
Sujudi and Haimes (1995)
. The separated shear layer was visualized approximately using the iso-surfaces of the stagnation
PPT Slide

Lager Image

- Vorticity

For a single-phase flow past a circular cylinder, the regular Karman vortex street and organized wake patterns are prominent for low
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- Reynolds stresses

The streamwise Reynolds stress (
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- Vorticity transport

Suh et al. (2011)
studied the mechanism of the vorticity generation at sub-critical
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- Overall Froude number effect with differentRe

There is little influence of Reynolds number at the air-water interface based on variables related to the air-water interface such as mean interface elevation, run-up height, and Kelvin wave patterns, etc. The instantaneous air-water interface elevations for all different Fr with different Re are shown in
Fig. 25
. For the lowest
PPT Slide

Lager Image

PPT Slide

Lager Image

FROUDE NUMBER EFFECT ON FLOW STRUCTURES WITH RE = 458,000

Three different Froude numbers (0.84, 1.24, and 1.64) with
- Air-water interface structures at differentFr

The mean air-water interface elevation contours for all three
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- Vorticity and velocity at the interface

Three-dimensional instantaneous coherent vortical structures identified by the second invariant of the velocity gradient tensor (Hunt et al., 1998) for the three Fr are shown in
Fig. 31
. At
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- Wave breaking at differentFr

Fig. 34
shows the iso-surface of the air-water interface with slices of the wave profile orthogonal to the cylinder wall for the three
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

CONCLUSIONS AND FUTURE WORK

The two-phase turbulent flow past an interface-piercing circular cylinder at different Reynolds and Froude numbers is studied numerically using large-eddy simulation with a Lagrangian dynamic subgrid-scale model. Verification and validation studies are performed to show the accuracy of the simulations with comparisons against available experimental data.
It is found that the variation in sub-critical
Acknowledgements

This work was sponsored by the Office of Naval Research (ONR) under Grant No. N000141-01-00-1-7, with Dr. Ki-Han Kim as the program manager. The simulations were performed at the Department of Defense (DoD) Supercomputing Resource Centers through the High Performance Computing Modernization Program (HPCMP).

Achenbach E.
1968
Distribution of local pressure and skin friction around a circular cylinder in cross-flow up to Re = 5×106
Journal of Fluid Mechanics
34
(4)
625 -
639

Bonmarin P.
1989
Geometric properties of deep-water breaking waves
Journal of Fluid Mechanics
209
405 -
433

Bhushan S.
,
Stern F.
,
Doctors L.J.
2010
Verification and validation of URANS wave resistance for air cushion vehicles, and comparison with linear theory
Journal of Ship Research
55
249 -
267

Cantwell B.J.
,
Coles D.
1983
An experimental study of entrainment and transport in turbulent near wake of a circular cylinder
Journal of Fluid Mechanics
136
321 -
374

Chaplin J.R.
,
Teigen P.
2003
Steady flow past a vertical surface-piercing circular cylinder
Journal of Fluids and Structures
18
(3-4)
271 -
285

Falgout R.D.
,
Jones J.E.
,
Yang U.M.
2006
Numerical solution of partial differential equations on parallel computers
Springer
Berlin

Galvin C.J.
1968
Breaker type classification on three laboratory beaches
Journal of Geophysical Research
73
(12)
3651 -
3659

Grue J.
,
Jensen A.
2006
Experimental velocities and accelerations in very steep wave events in deep water
European Journal of Mechanics B/Fluids
25
(5)
554 -
564

Hay A.D.
Flow about semi-submerged cylinders of finite length, Princeton University Report
Princeton University
Princeton, NJ

Hunt J.
,
Wray A.
,
Moin P.
1988
Eddies, stream, and convergence zones in turbulent flows
Center for Turbulence Research
In: Proceedings CTR Summer Program
Stanford, CA
193 -
208

Inoue M.
,
Bara N.
,
Himeno Y.
1993
Experimental and numerical study of viscous flow field around an advancing vertical circular cylinder piercing a free-surface
Journal of Kansai Society Naval Architecture
220
57 -
64

Jiang G.S.
,
Shu C.W.
1996
Efficient implementation of weighted ENO schemes
Journal of Computational Physics
126
202 -
228

Kandasamy M.
,
Xing T.
,
Stern F.
2009
Unsteady free surface wave-induced separation: Vortical structures and instabilities
Journal of Fluids and Structures
25
(2)
343 -
363

Kawamura T.
,
Mayer S.
,
Garapon A.
,
Sørensen L.
2002
Large eddy simulation of a flow past free surface piercing circular cylinder
Journal of Fluids Engineering
124
91 -
101

Kimmoun O.
,
Branger H.
2007
A particle image velocimetry investigation on laboratory surf-zone breaking waves over a sloping beach
Journal of Fluid Mechanics
588
353 -
397

Koo B.
,
Wang Z.
,
Yang J.
,
Stern F.
2012
Impulsive plunging wave breaking downstream of a bump in a shallow water flume-Part II: Numerical simulations
Journal of Fluids and Structures
32
121 -
134

Kravchenko A.G.
,
Moin P.
2000
Numerical studies of flow over a circular cylinder at ReD=3900
Physics of Fluids
12
403 -
417

Norberg C.
1992
Pressure forces on a circular cylinder in cross flow
Dynamicsand Instabilities
Göttingen, Germany
In: IUTAM Symposium Bluff-Body Wake
7-11 September 1992
275 -
278

Peregrine D.H.
1983
Breaking waves on beaches
Annual Review of Fluid Mechanics
15
147 -
178

Sarghini F.
,
Piomelli U.
,
Balaras E.
1999
Scale-similar models for large-eddy simulations
Physics of Fluids
11
(6)
1596 -
1607

Shakeri M.
,
Tavakolinejad M.
,
Duncan J.H.
2009
An experimental investigation of divergent bow waves simulated by a two-dimensional plus temporal wave marker technique
Journal of Fluid Mechanics
634
217 -
243

Shu C.W.
,
Osher S.
1988
Efficient implementation of essentially non oscillatory shock-capturing schemes
Journal of Computational Physics
77
(2)
439 -
471

Singh S.P.
,
Mittal S.
2005
Flow past a cylinder: shear layer instability and drag crisis
International Journal for Numerical Methods in Fluids
47
(1)
75 -
98

Stern F.
,
Wilson R.
,
Shao J.
2006
Quantitative V&V of CFD simulations and certification of CFD codes
International Journal for Numerical Methods in Fluids
50
(11)
1335 -
1355

Suh J.
,
Yang J.
,
Stern F.
2011
The effect of air-water interface on the vortex shedding from a vertical circular cylinder
Journal of Fluids and Structure
27
(1)
1 -
22

Sujudi D.
,
Haimess R.
1995
Identification of swirling flow in 3D vector fields, Technical report
Department of Aeronautics and Astronautics

Tallent J.R.
,
Yamashita T.
,
Tsuchiya Y.
1990
Transformation characteristics of breaking water waves
Water Wave Kinematics
178
509 -
523

Wang Z.
,
Yang J.
,
Koo B.G.
,
Stern F.
2009
A coupled level set and volume-of-fluid method for sharp interface simulation of plunging breaking waves
International Journal of Multiphase Flow
35
(3)
227 -
246

Wang Z.
,
Yang J.
,
Stern F.
2010
Numerical simulations of wave breakings around a wedge-shaped bow
Proceedings of 28th Symposium on Naval Hydrodynamics
Pasadena, California
12-17 September 2010

Wang Z.
,
Yang J.
,
Stern F.
2012
A new volume-of-fluid method with a constructed distance function on general structured grids
Journal of Computational Physics
231
(9)
3703 -
3722

Wissink J.G.
,
Rodi W.
2008
Numerical study of the near wake of a circular cylinder
International Journal of Heat and Fluid Flow
29
(4)
1060 -
1070

Xing T.
,
Stern F.
2010
Factors of safety for Richardson extrapolation
ASME Journal of Fluids Engineering
132
1 -
13

Yang J.
,
Stern F.
2009
Sharp interface immersed-boundary/level-set method for wave-body interactions
Journal of Computational Physics
228
(17)
6590 -
6616

Yeon S.
,
Yang J.
,
Stern F.
2013
Large eddy simulation of drag crisis in turbulent flow past a circular cylinder
ITTC workshop on wave run-up and vortex sheddin
Nantes, France
17-18 October 2013

Yoon S.H.
,
Kim D.H.
,
Sadat-Hosseini H.
,
Yang J.
,
Stern F.
2013
High-fidelity CFD simulation of wave run-up around vertical cylinders in monochromatic waves
ITTC workshop on wave run-up and vortex shedding
Nantes, France
17-18 October 2013

Yu G.
,
Avital E.J.
,
Williams J.J.R.
2008
Large eddy simulation of flow past free surface piercing circular cylinders
Journal of Fluids Engineering
130
10134.1 -
9

Citing 'Reynolds and froude number effect on the flow past an interface-piercing circular cylinder
'

@article{ E1JSE6_2014_v6n3_529}
,title={Reynolds and froude number effect on the flow past an interface-piercing circular cylinder}
,volume={3}
, number= {3}
, journal={International Journal of Naval Architecture and Ocean Engineering}
, publisher={The Society of Naval Architects of Korea}
, author={Koo, Bonguk
and
Yang, Jianming
and
Yeon, Seong Mo
and
Stern, Frederick}
, year={2014}
, month={Sep}