Large Eddy Simulation of Flow Past Free Surface Piercing Circular Cylinders

Flows past a free surface piercing cylinder are studied numerically by large eddy simulation at Froude numbers up to FrD=3.0 and Reynolds numbers up to ReD=1×105. A two-phase volume of fluid technique is employed to simulate the air-water flow and a flux corrected transport algorithm for transport of the interface. The effect of the free surface on the vortex structure in the near wake is investigated in detail together with the loadings on the cylinder at various Reynolds and Froude numbers. The computational results show that the free surface inhibits the vortex generation in the near wake, and as a result, reduces the vorticity and vortex shedding. At higher Froude numbers, this effect is stronger and vortex structures exhibit a 3D feature. However, the free surface effect is attenuated as Reynolds number increases. The time-averaged drag force on the unit height of a cylinder is shown to vary along the cylinder and the variation depends largely on Froude number. For flows at ReD=2.7×104, a negative pressure zone is developed in both the air and water regions near the free surface leading to a significant increase of drag force on the cylinder in the vicinity of the free surface at about FrD=2.0. The mean value of the overall drag force on the cylinder increases with Reynolds number and decreases with Froude number but the reduction is very small for FrD=1.6–2.0. The dominant Strouhal number of the lift oscillation decreases with Reynolds number but increases with Froude number.

[1]  A. Roshko Experiments on the flow past a circular cylinder at very high Reynolds number , 1961, Journal of Fluid Mechanics.

[2]  M. Rudman INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 24, 671–691 (1997) VOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS , 2022 .

[3]  P. Moin,et al.  Model consistency in large eddy simulation of turbulent channel flows , 1988 .

[4]  Jinhee Jeong,et al.  On the identification of a vortex , 1995, Journal of Fluid Mechanics.

[5]  I. Orlanski A Simple Boundary Condition for Unbounded Hyperbolic Flows , 1976 .

[6]  E. R. V. Driest On Turbulent Flow Near a Wall , 1956 .

[7]  Eugenio Oñate,et al.  BOUNDARY CONDITIONS FOR FINITE ELEMENT SIMULATIONS OF CONVECTIVE FLOWS WITH ARTIFICIAL BOUNDARIES , 1996 .

[8]  C. Peskin Numerical analysis of blood flow in the heart , 1977 .

[9]  S. Balachandar,et al.  Effect of three‐dimensionality on the lift and drag of nominally two‐dimensional cylinders , 1995 .

[10]  J. Gerrard The mechanics of the formation region of vortices behind bluff bodies , 1966, Journal of Fluid Mechanics.

[11]  Exact solutions of the generalized Navier– Stokes equations for benchmarking , 2002 .

[12]  M. Breuer LARGE EDDY SIMULATION OF THE SUBCRITICAL FLOW PAST A CIRCULAR CYLINDER: NUMERICAL AND MODELING ASPECTS , 1998 .

[13]  Ming-Hseng Tseng,et al.  Computation of three‐dimensional flow around square and circular piers , 2000 .

[14]  Stefan Mayer,et al.  Large Eddy Simulation of a Flow Past a Free Surface Piercing Circular Cylinder , 2002 .

[15]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[16]  Y. Himeno,et al.  EXPERIMENTAL AND NUMERICAL STUDY OF VISCOUS FLOW FIELD AROUND AN ADVANCING VERTICAL CIRCULAR CYLINDER PIERCING A FREE-SURFACE , 1993 .

[17]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[18]  Peter W. Bearman,et al.  Aspect ratio and end plate effects on vortex shedding from a circular cylinder , 1992, Journal of Fluid Mechanics.