0
Research Papers: CFD and VIV

RANS Simulation Versus Experiments of Flow Induced Motion of Circular Cylinder With Passive Turbulence Control at 35,000 < RE < 130,000

[+] Author and Article Information
Wei Wu

Department of Naval Architecture and
Marine Engineering,
University of Michigan,
Ann Arbor, MI
e-mail: wuwei@umich.edu

Michael M. Bernitsas

Professor
CTO of Vortex Hydro Energy,
Department of Naval Architecture and
Marine Engineering,
University of Michigan,
Ann Arbor, MI
e-mail: michaelb@umich.edu

Kevin Maki

Assistant Professor
Department of Naval Architecture and
Marine Engineering,
University of Michigan,
Ann Arbor, MI
e-mail: kjmaki@umich.edu

Contributed by the Ocean, Offshore, and Arctic Engineering Division of ASME for publication in the JOURNAL OF OFFSHORE MECHANICS AND ARCTIC ENGINEERING. Manuscript received August 17, 2011; final manuscript received June 19, 2014; published online July 29, 2014. Assoc. Editor: Solomon Yim.

J. Offshore Mech. Arct. Eng 136(4), 041802 (Jul 29, 2014) (10 pages) Paper No: OMAE-11-1076; doi: 10.1115/1.4027895 History: Received August 17, 2011; Revised June 19, 2014

Two-dimensional (2D) Unsteady Reynolds-Averaged Navier–Stokes equations (URANS) equations with the Spalart–Allmaras turbulence model are used to simulate the flow and body kinematics of the transverse motion of spring-mounted circular cylinder. The flow is in the high-lift TrSL3 regime of a Reynolds number in the range 35,000 < Re < 130,000. Passive turbulence control (PTC) in the form of selectively distributed surface roughness is used to alter the cylinder flow induced motion (FIM). Simulation is performed using a solver based on the open source Computational Fluid Dynamics (CFD) tool OpenFOAM, which solves continuum mechanics problems with a finite-volume discretization method. Roughness parameters of PTC are chosen based on tests conducted in the Marine Renewable Energy Lab (MRELab) of the University of Michigan. The numerical tool is first tested on smooth cylinder in vortex-induced vibration (VIV) and results are compared with available experimental measurements and URANS simulations. For the cylinder with PTC cases, the sandpaper grit on the cylinder wall is modeled as a rough-wall boundary condition. Two sets of cases with different system parameters (spring, damping) are simulated and the results are compared with experimental data measured in the MRELab. The amplitude ratio curve shows clearly three different branches, including the VIV initial and upper branches, and a galloping branch. The numerical branches are similar to those observed experimentally. Frequency ratio, vortex patterns, transitional behavior, and lift are also predicted well for PTC cylinders at such high Reynolds numbers.

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

References

Bernitsas, M. M., Raghavan, K., Ben-Simon, Y., and Garcia, E. M. H., 2008, “VIVACE (Vortex Induced Vibration Aquatic Clean Energy): A New Concept in Generation of Clean and Renewable Energy From Fluid Flow,” ASME J. Offshore Mech. Arct. Eng., 130(4), p. 041101. [CrossRef]
Bernitsas, M. M., and Raghavan, K., 2009, “Converter of Current, Tide, or Wave Energy,” U.S. Patent and Trademark Office Patent No. 7,493,759 B2.
Bernitsas, M. M., Ben-Simon, Y., Raghavan, K., and Garcia, E. M. H., 2009, “The VIVACE Converter: Model Tests at High Damping and Reynolds Numbers Around 105,” ASME J. Offshore Mech. Arct. Eng., 131(1), p. 011102. [CrossRef]
Bernitsas, M. M., Raghavan, K., and Duchene, G., 2008, “Induced Separation and Vorticity Using Roughness in VIV of Circular Cylinders at 8 × 103 < Re < 2.0 × 105,” ASME 27th International OMAE Conference, Estoril, Portugal, June 15–20. [CrossRef]
Raghavan, K., and Bernitsas, M. M., 2011, “Experimental Investigation of Reynolds Number Effect on Vortex Induced Vibration of Rigid Cylinder on Elastic Supports,” Ocean Eng., 38(5–6), pp. 719–731. [CrossRef]
Chang, C. C., Kumar, R. A., and Bernitsas, M. M., 2010, “VIV and Galloping of Single Circular Cylinder With Surface Roughness at 3.0 × 104 ≤ Re ≤ 1.2 × 105,” Ocean Eng., 38(16), pp. 1713–1732. [CrossRef]
Chang, C. C., 2010, “Hydrokinetic Energy Harnessing by Enhancement of Flow Induced Motion Using Passive Turbulence Control,” Ph.D. thesis, University of Michigan, Ann Arbor, MI.
Zdravkovich, M. M., 1997, Flow Around Circular Cylinders Vol. 1: Fundamentals, Oxford University, Oxford, UK.
Khalak, A., and Williamson, C. H. K., 1996, “Dynamics of a Hydroelastic Cylinder With Very Low Mass and Damping,” J. Fluids Struct., 10(5), pp. 455–472. [CrossRef]
Khalak, A., and Wliamson, C. H. K., 1999, “Motions Forces and Mode Transitions in Vortex-Induced Vibrations at Low Mass-Damping,” J. Fluids Struct., 13(7–8), pp. 813–851. [CrossRef]
Guilmineau, E., and Queutey, P., 2004, “Numerical Simulation of Vortex-Induced Vibration of a Circular Cylinder With Low Mass-Damping in a Turbulent Flow,” J. Fluids Struct., 19(4), pp. 449–466. [CrossRef]
Spalart, P. R., and Allmaras, S. R., 1994, “A One-Equation Turbulence Model for Aerodynamic Flows,” La Recherche Aerospatiale, pp. 5–21.
Travin, A., Shur, M., Strelets, M., and Spalart, P., 2000, “Detached-Eddy Simulations Past a Circular Cylinder,” Flow, Turbul. Combust., 63(1–4), pp. 293–310. [CrossRef]
Shur, M. L., Spalart, P. R., Strelets, M. Kh., and Travin, A. K., 1996, “Navier–Stokes Simulation of Shedding Turbulent Flow Past a Circular Cylinder and a Cylinder With a Backward Splitter Plate,” Third ECCOMAS CFD Conference, G. A.Desideri, C.Hirsch, P.Le Tallec, M.Pandolfi, and J.Periaux, eds., Wiley, Chichester, Paris, Sept. 9–13, pp. 676–682.
Morgan, P. E., and Visbal, M. R., 2006, “Application of Hybrid RANS/ILES to Geometries With Separated Flows,” 3rd AIAA Flow Control Conference, San Francisco, CA, June 5–7.
Breuer, M., 2000, “A Challenging Test Case for Large Eddy Simulation: High Reynolds Number Circular Cylinder Flow,” Int. J. Heat Fluid Flow, 21(5), pp. 648–654. [CrossRef]
Cantwell, B., and Coles, D., 1983, “An Experimental Study of Entrainment and Transport in the Turbulent Near Wake of a Circular Cylinder,” J. Fluid Mech., 136, pp. 321–374. [CrossRef]
Wieselsberger, C., 1921, “New Data on the Law of Hydro and Aerodynamic Resistance,” Physikalische Zeitschrift, (in German).
Roshko, A., 1961, “Experiments on the Flow Past a Circular Cylinder at Very High Reynolds Number,” J. Fluid Mech., 10(3), pp. 345–356. [CrossRef]
Wanderley, J. B. V., Sphaier, S. H., and Levi, C., 2008, “A Numerical Investigation of Vortex Induced Vibration on an Elastically Mounted Rigid Cylinder,” ASME 27th International Conference, (OMAE2008), Estoril, Portugal, June 15–20, pp. 703–711. [CrossRef]
Wanderley, J. B. V., Sphaier, S. H., and Levi, C., 2009, “A Numerical Investigation of the Hysteresis Effect on Vortex Induced Vibration on an Elastically Mounted Rigid Cylinder,” ASME 28th International Conference, (OMAE2009), Honolulu, HI, May 31–June 5, pp. 653–660. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

Close-up of the structured-medium grid S2 for smooth cylinder

Grahic Jump Location
Fig. 2

Configuration of sand strips (PTC) along the cylinder [7]

Grahic Jump Location
Fig. 3

Grid R1 for a cylinder with PTC

Grahic Jump Location
Fig. 4

Amplitude ratio A* (AMax/D) (m*ζ = 0.013, m* = 2.4) compared with Ref. [9]

Grahic Jump Location
Fig. 5

Frequency ratio fo/fn,water (m*ζ=0.013, m* = 2.4) compared with Ref. [10]

Grahic Jump Location
Fig. 6

(a) Double-peaked spectrum of amplitude Re = 3000 (Ur,air = 3, Ur,water = 3.57) in Quasi-Periodic sub-branch and (b) single-peaked spectrum of amplitude Re = 3800 (Ur,air = 3.8, Ur,water = 4.52) in periodic sub-branch

Grahic Jump Location
Fig. 7

Different vortex patterns for smooth cylinder in VIV (m*ζ = 0.013, m* = 2.4)

Grahic Jump Location
Fig. 8

Time traces of displacement ratio and lift coefficient for smooth cylinder: (a) Re = 5000 (Ur,air = 5, Ur,water = 5.95) and (b) Re = 9000 (Ur,air = 9, Ur,water = 10.71)

Grahic Jump Location
Fig. 9

Force comparison between current study and Ref. [11] (2D-URANS k-ω): (a) Maximum lift coefficient and (b) mean drag coefficient

Grahic Jump Location
Fig. 10

Amplitude ratio comparison between present study and Ref. [7], (K = 1600 N/m, ζharn = 0.08)

Grahic Jump Location
Fig. 11

Displacement ratio (y(t)/D) and lift coefficient for different Re (K = 1600 N/m, ζharn = 0.08)

Grahic Jump Location
Fig. 12

Frequency ratio comparison between present study and Ref. [7], University of Michigan (K = 1600 N/m, ζharn = 0.08)

Grahic Jump Location
Fig. 13

Amplitude ratio (ARMS*) comparison between present study and Ref. [7], University of Michigan (K = 2000 N/m, ζharn = 0.08)

Grahic Jump Location
Fig. 14

Frequency ratio comparison between present study and Ref. [7], University of Michigan (K = 2000 N/m, ζharn = 0.08)

Grahic Jump Location
Fig. 15

Different vortex patterns for PTC cylinder in VIV and galloping (m*ζ = 0.013, m* = 2.4)

Grahic Jump Location
Fig. 16

FFT analysis of displacement ratio (y(t)/D) and lift coefficient for different Re (K = 1600 N/m, ζharn = 0.08)

Tables

Errata

Discussions

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