Research Papers: Ocean Engineering

Amplitude Induced Nonlinearity in Piston Mode Resonant Flow: A Fully Viscous Numerical Analysis

[+] Author and Article Information
Luca Bonfiglio

Department of Mechanical Engineering,
MIT-Sea Grant College,
Massachusetts Institute of Technology,
Cambridge, MA 02139
e-mail: bonfi@mit.edu

Stefano Brizzolara

Kevin T. Crofton Department of Aerospace and
Ocean Engineering,
Virginia Tech,
Blacksburg, VA 24061
e-mail: stebriz@vt.edu

1Corresponding author.

Contributed by the Ocean, Offshore, and Arctic Engineering Division of ASME for publication in the JOURNAL OF OFFSHORE MECHANICS AND ARCTIC ENGINEERING. Manuscript received April 8, 2016; final manuscript received July 13, 2017; published online September 29, 2017. Assoc. Editor: Muk Chen Ong.

J. Offshore Mech. Arct. Eng 140(1), 011101 (Sep 29, 2017) (11 pages) Paper No: OMAE-16-1038; doi: 10.1115/1.4037487 History: Received April 08, 2016; Revised July 13, 2017

Near field flow characteristics around catamarans close to resonant conditions involve violent viscous flow such as energetic vortex shedding and steep wave making. This paper presents a systematic and comprehensive numerical investigation of these phenomena at various oscillating frequencies and separation distances of twin sections. The numerical model is based on the solution of Navier–Stokes equations assuming laminar-flow conditions with a volume of fluid (VOF) approach which has proven to be particularly effective in predicting strongly nonlinear radiated waves which directly affect the magnitude of the hydrodynamic forces around resonant frequencies. Considered nonlinear effects include wave breaking, vortex shedding and wave-body wave-wave interactions. The method is first validated using available experiments on twin circular sections: the agreement in a very wide frequency range is improved over traditional linear potential flow based solutions. Particular attention is given to the prediction of added mass and damping coefficients at resonant conditions where linear potential flow methods fail, if empirical viscous corrections are not included. The results of the systematic investigation show for the first time how the so-called piston-mode motion characteristics are nonlinearly dependent on the gap width and motion amplitude. At low oscillation amplitudes, flow velocity reduces and so does the energy lost for viscous effects. On the other hand for higher oscillation amplitude, the internal free surface breaks dissipating energy hence reducing the piston mode amplitude. These effects cannot be numerically demonstrated without a computational technique able to capture free surface nonlinearity and viscous effects.

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


Molin, B. , 1999, “ On the Piston Mode in Moonpools,” 14th International Workshop on Water Waves and Floating Bodies (IWWWFB), Port Huron, MI, Apr. 11–14, pp. 103–106.
McIver, P. , McIver, M. , and Zhang, J. , 2003, “ Excitation of Trapped Water Waves by the Forced Motion of Structures,” J. Fluid Mech., 494, pp. 141–162. [CrossRef]
McIver, M. , 1996, “ An Example of Non-Uniqueness in the Two-Dimensional Linear Water Wave Problem,” J. Fluid Mech., 315, pp. 257–266. [CrossRef]
McIver, P. , 2005, “ Complex Resonances in the Water-Wave Problem for a Floating Structure,” J. Fluid Mech., 536, pp. 423–443. [CrossRef]
McIver, P. , and Newman, J. , 2003, “ Trapping Structures in the Three Dimensional Water-Wave Problem,” J. Fluid Mech., 484, pp. 283–301. [CrossRef]
Kristiansen, T. , and Faltinsen, O. , 2010, “ A Two-Dimensional Numerical and Experimental Study of Resonant Coupled Ship and Piston Mode Motion,” Appl. Ocean Res., 32(2), pp. 158–176. [CrossRef]
Molin, B. , 2001, “ On the Piston and Sloshing Modes in Moonpools,” J. Fluid Mech., 430, pp. 27–50. [CrossRef]
Faltinsen, O. , Rognebakke, O. , and Timokha, A. , 2007, “ Two Dimensional Resonant Piston-Like Sloshing in a Moonpool,” J. Fluid Mech., 575, pp. 359–397. [CrossRef]
Yeung, R. , and Seah, R. , 2007, “ On Helmholtz and Higher-Order Resonance of Twin Floating Bodies,” J. Eng. Math., 58(1–4), pp. 251–265. [CrossRef]
Porter, R. , and Evans, D. , 2008, “ Examples of Trapped Modes in the Presence of Freely Floating Structures,” J. Fluid Mech., 606, pp. 189–207. [CrossRef]
Newman, J. , 1997, “ Resonant Diffraction Problems,” 12th International Workshop on Water Waves and Floating Bodies (IWWWFB), Carry Le Rouet, France, Mar. 16–19, pp. 307–308.
Faltinsen, O. , and Timokha, A. , 2015, “ On Damping of Two-Dimensional Piston-Mode Sloshing in a Rectangular Moonpool Under Forced Heave Motions,” J. Fluid Mech., 772, p. R1. [CrossRef]
Buchner, B. , Dijk, A. , and de Wilde, J. , 2001, “ Numerical Multiple-Body Simulations of Side-by-Side Mooring to an FPSO,” 11th International Offshore and Polar Engineering Conference, Stavanger, Norway, June 17–22, SPE Paper No. ISOPE-I-01-053.
Newman, J. , 2004, “ Progress in Wave Load Computations on Offshore Structure,” 23rd International Conference on Offshore Mechanics and Arctic Engineering, Vancouver, BC, Canada, June 20–24.
Chen, X. , Duan, W. , and Dai, Y. , 2005, “ Accurate Computation of Second-Order Low-Frequency Loads,” 19th National Conference on Hydrodynamics/Seventh National Congress on Hydrodynamic, Harbin, China, Aug., pp. 1–10.
Pauw, W. , Huijsmans, R. , and Voogt, A. , 2007, “ Advances in the Hydrodynamic of Side-by-Side Moored Vessels,” ASME Paper No. OMAE2007-29374.
Centeno, R. , Fonseca, N. , and Soares, C. G. , 2000, “ Prediction of Motion of Catamarans Accounting for Viscous Effects,” Int. Shipbuild. Prog., 47(451), pp. 303–323.
Thwaites, B. , 1960, Incompressible Aerodynamics, Oxford University Press, New York, Chap. 5.
Lu, L. , Liang, C. , Bin, T. , and Liang, S. , 2010, “ Comparison of Potential Flow and Viscous Fluid Models in Gap Resonance,” 25th International Workshop on Water Waves and Floating Bodies (IWWWFB), Harbin, China, May 9–12, pp. 1–4.
Kistiansen, T. , and Faltinsen, O. , 2012, “ Gap Resonances Analyzed by a New Domain-Decomposition Method Combining Potential and Viscous Flow Draft,” Appl. Ocean Res., 34, pp. 198–208. [CrossRef]
Kristiansen, T. , Sauder, T. , and Firoozkoohi, R. , 2013, “ Validation of a Hybrid Code Combining Potential and Viscous Flow With Application to 3D Moonpool,” ASME Paper No. OMAE2013-10748.
Fredriksen, A. G. , Kristiansen, T. , and Faltinsen, O. M. , 2014, “ Experimental and Numerical Investigation of Wave Resonance in Moonpools at Low Forward Speed,” Appl. Ocean Res., 47, pp. 28–46. [CrossRef]
Fredriksen, A. G. , Kristiansen, T. , and Faltinsen, O. M. , 2015, “ Wave-Induced Response of a Floating Two-Dimensional Body With a Moonpool,” Philos. Trans. R. Soc. London A, 373(2033), pp. 1–17.
Sun, L. , Eatock Taylor, R. , and Taylor, P. H. , 2010, “ First-and Second-Order Analysis of Resonant Waves Between Adjacent Barges,” J. Fluids Struct., 26(6), pp. 954–978. [CrossRef]
Ohkusu, M. , 1969, “ On the Heaving Motion of a Circular Cylinder on the Surface of a Fluid,” Rep. Res. Inst. Appl. Mech., 17(58), pp. 167–185.
Wang, S. , and Wahab, R. , 1971, “ Heaving Oscillations of Twin Cylinders in a Free Surface,” J. Ship Res., 15(1), pp. 33–48.
Lee, C. , Jones, H. , and Bedel, J. , 1971, “ Added Mass and Damping Coefficients of Heaving Twin Cylinders in a Free Surface,” Naval Ship Research and Development Center, Bethesda, MD, Technical Report No. 3695.
Bonfiglio, L. , Brizzolara, S. , and Chryssostomidis, C. , 2012, “ Added Mass and Damping of Oscillating Bodies: A Fully Viscous Numerical Approach,” Nineth WSEAS International Conference on Fluid Mechanics (FLUIDS 1́2), Cambridge, MA, Jan. 25–27, pp. 210–215.
Bonfiglio, L. , and Brizzolara, S. , 2014, “ Unsteady Viscous Flow With Non Linear Free Surface Around Oscillating SWATH Ship Sections,” WSEAS Trans. Fluid Mech., 9, pp. 49–57.
Bonfiglio, L. , and Brizzolara, S. , 2013, “ Influence of Viscosity on Radiation Forces: A Comparison Between Monohull, Catamaran and SWATH,” 23rd International Offshore and Polar Engineering, Anchorage, AK, June 30–July 5, SPE Paper No. ISOPE-I-13-389.
Fletcher, R. , 1975, “ Conjugate Gradient Methods for Indefinite Systems,” Lecture Notes in Mathematics, Vol. 506, Springer, Berlin, pp. 773–789. [CrossRef]
Jasak, H. , and Tukovic, Z. , 2007, “ Automatic Mesh Motion for the Unstructured Finite Volume Method,” Trans. FAMENA, 30(2), pp. 1–20.
Jasak, H. , 2009, “ Dynamic Mesh Handling in OpenFOAM,” AIAA Paper No. 2009-341.
Issa, R. , 1985, “ Solution of the Implicitly Discretized Fluid Flow Equations by Operator Splitting,” J. Comput. Phys., 62(1), pp. 40–65. [CrossRef]
Ohkusu, M. , 1970, “ On the Motions of Multihull Ships in Waves,” Rep. Res. Inst. Appl. Mech., 18(60), pp. 33–60.
Frank, W. , 1967, “ Oscillation of Cylinders in or Below the Free Surface of Deep Fluids,” Naval Ship Research and Development Center, Washington, DC, Technical Report No. 2375.


Grahic Jump Location
Fig. 1

Catamaran cross section and probes for wave elevation. The circular cross section center of each hull is placed at a distance b from the catamaran symmetry plane and is characterized by a radius a corresponding to the depth of the catamaran at rest. Dimensions indicated are expressed in meters.

Grahic Jump Location
Fig. 2

Grid sensitivity study: coarsest-medium-finest meshes designed to be simulated at 6 rad/s. Three different mesh regions: structured around the free surface, body-fitted structured around the body and unstructured in the far field. Mesh size at cylinder wall d1=3e−4 m for M1 (coarse), d2=1e−4 m for M2 (medium) and d3=2e−5 m for M3 (fine).

Grahic Jump Location
Fig. 3

Validation study: CFD results are presented through solid and dashed line for ξ=0.635 cm and ξ=1.27 cm, respectively. Dashed-dot line represents potential flow theory. Experiments have been represented through filled dots when stiffeners have been used during the experiments. (a) Nondimensional added mass coefficient: b=1.5a, (b) nondimensional damping coefficient: b=1.5a, (c) nondimensional added mass coefficient: b=2.0a, (d) nondimensional damping coefficient: b=2.0a, (e) nondimensional added mass coefficient: b=4.0a, and (f) nondimensional damping coefficient: b=4.0a.

Grahic Jump Location
Fig. 4

Nondimensional vorticity contour for hull configuration b=1.5a at the piston mode frequency: (a) ξ=0.635 cm, (b) ξ=1.27 cm, (c) ξ=2.54 cm, and (d) ξ=3.81 cm

Grahic Jump Location
Fig. 5

Nondimensional wave elevation at probes located at centerline y = 0 (A), y=b+2a (B), y=b+4a (C), and y=b+8a (D). Hull distance b=1.5a, heave amplitude ξ=0.635 cm, oscillation frequency ω=6.50 rad/s.

Grahic Jump Location
Fig. 6

Nondimensional added mass, damping coefficients and inner wave amplitude (probeA). From left to right b=1.5a, b=2a and b=4a. From top to bottom ξ=0.635 cm, ξ=1.270 cm, ξ=2.540 cm and ξ=3.81 cm. Vertical lines correspond to undamped piston mode frequency resulting from Eq. (8) with λ = 1. Nondimensional added mass values are indicated on the left axis while nondimensional damping and inner wave amplitude values are indicated on the right axis.

Grahic Jump Location
Fig. 7

Hydrodynamic coefficients resulting from a forced heave motion of amplitude ξ=1.27 cm. Symbols on the curves represent computed points: (a) nondimensional added mass coefficient and (b) nondimensional damping coefficient.

Grahic Jump Location
Fig. 8

Nondimensional added mass (left axis), damping, inner (probe A) and outer wave amplitude (probe B, C, D) (right axis) for hull distance b=1.5a and motion amplitude ξ=1.27 cm. Damping experiences a relative minimum at nondimensional frequency of 0.903 corresponding to the trapping mode frequency presented in Fig. 9. (a) ω=6.25 rad/s : piston mode frequency, (b) ω=7.25 rad/s : zero damping frequency, and (c) wave elevation plot comparison (circles: piston mode frequency, crosses: zero damping frequency).

Grahic Jump Location
Fig. 9

Volume of fluid contour and wave elevation plots for ω=6.25 rad/s and ω=7.25 rad/s. Hull distance b=1.5a and oscillation amplitude ξ=1.27 cm.

Grahic Jump Location
Fig. 10

Free decay test of a column of water located at z = a (deck height) for an hull configuration b=1.5a: (a) inner wave elevation and (b) fast Fourier transform (FFT) of inner wave elevation

Grahic Jump Location
Fig. 11

Nondimensional piston-mode frequency trend with respect to heave amplitude for three different gap widths

Grahic Jump Location
Fig. 12

Nondimensional inner wave amplitude trend with respect to heave amplitude for three gap widths

Grahic Jump Location
Fig. 13

Nondimensional inner wave amplitude trend with respect to gap width for four heave amplitudes



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