Research Papers: Flows in Complex Systems

Effects of Submergence on Low and Moderate Reynolds Number Free-Surface Flow Around a Matrix of Cubes

[+] Author and Article Information
Z. Ikram

School of Engineering and Materials,
Queen Mary University of London,
327 Mile End Road,
London E1 4NS, UK
e-mail: zaheerikram@gmail.com

E. J. Avital

School of Engineering and Materials,
Queen Mary University of London,
327 Mile End Road,
London E1 4NS, UK
e-mail: e.avital@qmul.ac.uk

J. J. R. Williams

School of Engineering and Materials,
Queen Mary University of London,
327 Mile End Road,
London E1 4NS, UK
e-mail: j.j.r.williams@qmul.ac.uk

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received July 4, 2014; final manuscript received October 12, 2015; published online January 4, 2016. Assoc. Editor: Zhongquan Charlie Zheng.

J. Fluids Eng 138(5), 051102 (Jan 04, 2016) (11 pages) Paper No: FE-14-1354; doi: 10.1115/1.4031852 History: Received July 04, 2014; Revised October 12, 2015

The effect of reducing submergence depth at a low and moderate Reynolds number flow is investigated using large eddy simulation (LES) around a matrix of cubes. The submerged body is modeled using an immersed boundary method, while the free-surface is accounted for using a moving mesh. Results show that for reducing the submergence depth, the forces acting on the cube reduce as the force variation increased. Variation in depth is also found to influence the level of damping and redistribution of turbulence near the free-surface boundary. Both submergence depth and Reynolds number are also found to influence the dominant free-surface signature and shedding frequencies from the cube. In the interobstacle region (IOR), the variation of Reynolds number and submergence depth is found to have little effect.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


Shen, Z. , Wang, B. , Cui, G. , and Zhang, Z. , 2015, “ Flow Pattern and Pollutant Dispersion Over Three Dimensional Building Arrays,” Atmos. Environ., 116, pp. 202–215. [CrossRef]
Huang, C. J. , and Dong, C. M. , 2001, “ On the Interaction of a Solitary Wave and a Submerged Dike,” Coastal Eng., 43(3–4), pp. 265–286. [CrossRef]
Hussein, H. J. , and Martinuzzi, R. J. , 1996, “ Energy Balance for Turbulent Flow Around a Surface Mounted Cube Placed in a Channel,” Phys. Fluids, 8(3), pp. 764–780. [CrossRef]
Shah, K. B. , and Ferziger, J. H. , 1997, “ A Fluid Mechanicians View of Wind Engineering: Large Eddy Simulation of Flow Past a Cubic Obstacle,” J. Wind Eng. Ind. Aerodyn., 67–68, pp. 211–224. [CrossRef]
Martinuzzi, R. J. , and Havel, B. , 2004, “ Vortex Shedding From Two Surface-Mounted Cubes in Tandem,” Int. J. Heat Fluid Flow, 25(3), pp. 364–372. [CrossRef]
Pinarbasi, A. , Pinar, E. , Akilli, H. , and Ince, E. , 2015, “ Shallow Water Experiments of Flow Past Two Identical Square Cylinders in Tandem,” Eur. J. Mech.-B/Fluids, 49, pp. 100–107. [CrossRef]
Meinders, E. R. , and Hanjalic, K. , 1999, “ Vortex Structure and Heat Transfer in Turbulent Flow Over a Wall-Mounted Matrix of Cubes,” Int. J. Heat Fluid Flow, 20(3), pp. 255–267. [CrossRef]
Niceno, B. , Dronkers, A. D. T. , and Hanjalic, K. , 2002, “ Turbulent Heat Transfer From a Multi-Layered Wall-Mounted Cube Matrix: A Large Eddy Simulation,” Int. J. Heat Fluid Flow, 23(2), pp. 173–185. [CrossRef]
Martinuzzi, R. , and Tropea, C. , 1993, “ Flow Around Surface-Mounted, Prismatic Obstacles Placed in a Fully Developed Channel Flow,” ASME J. Fluids Eng., 115(1), pp. 85–92. [CrossRef]
Thomas, T. G. , and Williams, J. J. R. , 1999, “ Simulation of Skewed Turbulent Flow Past a Surface Mounted Cube,” J. Wind Eng. Ind. Aerodyn., 81(1–4), pp. 347–360. [CrossRef]
Xie, Z. , and Castro, I. P. , 2006, “ LES and RANS for Turbulent Flow Over Arrays of Wall-Mounted Obstacles,” Flow, Turbul. Combust., 76(3), pp. 291–312. [CrossRef]
Okamoto, S. , Tsunoda, K. , Katsumata, T. , Abe, N. , and Kijima, M. , 1996, “ Turbulent Near-Wakes of Periodic Array of Square Blocks on a Plate,” Int. J. Heat Fluid Flow, 17(3), pp. 211–218. [CrossRef]
Williams, J. J. R. , 2005, “ Curvilinear Turbulence Modeling of Open Channel Flow,” J. Hydraul. Res., 43(2), pp. 158–164. [CrossRef]
Ikram, Z. , Avital, E. J. , and Williams, J. J. R. , 2012, “ Detached Eddy Simulation of Free-Surface Flow Around a Submerged Submarine Fairwater,” ASME J. Fluids Eng., 134(6), p. 061103. [CrossRef]
Jasak, H. , Weller, H. G. , and Gosman, A. D. , 1999, “ High Resolution NVD Differencing Scheme for Arbitrarily Unstructured Meshes,” Int. J. Numer. Methods Fluids, 31(2), pp. 431–449. [CrossRef]
Bensow, R. E. , Persson, T. , Fureby, C. , Svennberg, U. , and Alin, N. , 2004, “ Large Eddy Simulation of the Viscous Flow Around Submarine Hulls,” 25th Symposium on Naval Hydrodynamics, St. John's, NL, Canada, Aug. 8–13, pp. 1–16.
Le, H. , and Moin, P. , 1991, “ An Improvement of Fractional Step Methods for the Incompressible Navier–Stokes Equations,” J. Comput. Phys., 92(2), pp. 369–379. [CrossRef]
Thomas, T. G. , Leslie, D. C. , and Williams, J. J. R. , 1995, “ Free-Surface Simulations Using a Conservative 3D Code,” J. Comput. Phys., 116(1), pp. 52–68. [CrossRef]
Chan, R. K.-C. , and Street, R. L. , 1970, “ A Computer Study of Finite-Amplitude Water Waves,” J. Comput. Phys., 6(1), pp. 68–94. [CrossRef]
Apsley, D. , and Hu, W. , 2003, “ CFD Simulation of Two- and Three-Dimensional Free-Surface Flow,” Int. J. Numer. Methods Fluids, 42(5), pp. 465–491. [CrossRef]
Yang, C. , and Lohner, R. , 2003, “ Prediction of Flows Over an Axisymmetric Body With Appendages,” 8th International Conference on Numerical Ship Hydrodynamics, Busan, Korea, Sept. 22–25, pp. 1–15.
Shi, J. , Thomas, T. G. , and Williams, J. J. R. , 2000, “ Free-Surface Effects in Open Channel Flow at Moderate Froude and Reynold's Numbers,” J. Hydraul. Res., 38(6), pp. 465–474. [CrossRef]
Werner, H. , and Wengle, H. , 1991, “ Large-Eddy Simulation of Turbulent Flow Over and Around a Cube in a Plate Channel,” 8th International Symposium on Turbulent Shear Flows, Munich, Germany, Sept. 9–11, pp. 155–168.
Tseng, Y. H. , and Ferziger, J. H. , 2003, “ A Ghost-Cell Immersed Boundary Method for Flow in Complex Geometry,” J. Comput. Phys., 192(2), pp. 593–623. [CrossRef]
Majumdar, S. , Iaccarino, G. , and Durbin, P. , 2001, “ RANS Solvers With Adaptive Structured Boundary Non-Conforming Grids,” Annual Research Briefs 2001, Stanford, CA, pp. 353–366.
Schmidt, S. , and Thiele, F. , 2002, “ Comparison of Numerical Methods Applied to the Flow Over Wall-Mounted Cubes,” Int. J. Heat Fluid Flow, 23(3), pp. 330–339. [CrossRef]
Hunt, J. C. R. , Wray, A. A. , and Moin, P. , 1988, “ Eddies, Streams, and Convergence Zones in Turbulent Flows,” Summer Program 1988, Center for Turbulence Research, Stanford, CA, pp. 193–208.


Grahic Jump Location
Fig. 1

(a) Simplified representation of the uniformly spaced matrix of cubes, as studied by Meinders and Hanjalic [7] and (b) XY plane view of the grid mesh

Grahic Jump Location
Fig. 2

Comparison of the computed streamwise velocity to the experimental measurements of Meinders and Hanjalic [7]

Grahic Jump Location
Fig. 3

Comparison of computed Reynolds stresses against experimental measurements of Meinders and Hanjalic [7], for Reh = 3584

Grahic Jump Location
Fig. 4

Time-averaged streamline trace around the cube take about the XZ plane for y/h = 2.0 for Reh = 3584 (left) and d/h = 2

Grahic Jump Location
Fig. 5

Time-averaged streamwise and vertical velocity components and mean pressure distributions around the cube along the XZ plane for y/h = 2.0. Results for all cases are plotted.

Grahic Jump Location
Fig. 6

Time-averaged normal Reynolds stress and < u′w′> cross Reynolds stress component along the XZ plane for y/h = 2.0 for Reh = 3584

Grahic Jump Location
Fig. 7

Time-averaged distribution of the fluctuation kinetic energy around the cube about the symmetric plane (left), and XY plane z/h = 0.4 (right) for Reh = 3584

Grahic Jump Location
Fig. 8

Turbulent kinetic production and dissipation along the center XZ plane for Reh = 3584

Grahic Jump Location
Fig. 9

Energy spectrums for Reh = 40,000 (left) and Reh = 3584 (right) plotted at position 2.54x/h in the streamwise direction for d/h = 2.0

Grahic Jump Location
Fig. 10

Instantaneous streamline traces along the frontal face (60% of the face) of a cube in a matrix configuration at various time (dimensionless) intervals for Reh = 3584 and Reh = 40,000

Grahic Jump Location
Fig. 11

Instantaneous streamlines in XZ midplane and fixed time intervals for both Reh = 3584 and Reh = 40,000 at d/h = 2.0

Grahic Jump Location
Fig. 12

Instantaneous vortical structures in the near-wall region (z/h < 0.3) visualized for Q = 130 and for d/h = 2.0

Grahic Jump Location
Fig. 13

Time variation and free-surface interaction of vortical structures viewed about the XZ axes for d/h = 2.0, visualized isosurface value of Q = 130

Grahic Jump Location
Fig. 14

Vortex data acquired over 20 shedding cycles for time that Q > 0, the data are taken across the IOR for Reh = 3584

Grahic Jump Location
Fig. 15

Instantaneous surface plots shown at fixed intervals for d/h = 2.0 at Reh = 3584 and Reh = 40,000




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In