Research Papers: Fundamental Issues and Canonical Flows

Time Domain Simulation of Lifting Bodies Acting at or Near the Free Surface With Vortex Particle Wakes

[+] Author and Article Information
Rachel Gouveia

Navatek, Ltd.,
1080 Kingstown Road,
Wakefield, RI 02879
e-mail: rgouveia@navatekltd.com

Stephanie Fitzpatrick

Navatek, Ltd.,
1080 Kingstown Road,
Wakefield, RI 02879
e-mail: sfitzpatrick@navatekltd.com

Amanda Costa

Navatek, Ltd.,
1080 Kingstown Road,
Wakefield, RI 02879
e-mail: acosta@navatekltd.com

David Kring

Navatek, Ltd.,
1080 Kingstown Road,
Wakefield, RI 02879
e-mail: dkring@navatekltd.com

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received February 14, 2018; final manuscript received November 27, 2018; published online January 14, 2019. Assoc. Editor: Olivier Coutier-Delgosha.

J. Fluids Eng 141(4), 041201 (Jan 14, 2019) (8 pages) Paper No: FE-18-1101; doi: 10.1115/1.4042147 History: Received February 14, 2018; Revised November 27, 2018

Boundary element method (BEM) potential-flow solvers are regularly used in industrial applications due to their quick setup and computational time. In aerodynamics, vortex particle methods (VPM) are widely used with BEM potential-flow solvers for modeling lift. However, they are seldom applied to the ocean environment. This paper discusses the implementation of a VPM into Aegir, an existing time-domain, seakeeping, medium-fidelity, BEM potential-flow solver. The wake in the VPM is modeled using both a small dipole buffer wake sheet and vortex particles. It has been observed that this method captures both the details of complex wake patterns behind lift-producing surfaces and the expected lift force, thus improving the accuracy of the solution. Two new contributions presented in this paper include the extension of the VPM from previous source-based methods to a potential formulation and full interaction with free surface waves.

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


Willis, D. J. , Peraire, J. , and White, J. K. , 2007, “ A Combined pFFT-Multipole Tree Code, Unsteady Panel Method With Vortex Particle Wakes,” Int. J. Numer. Methods Fluids, 53(8), pp. 1399–1422. [CrossRef]
Ploumhans, P. , Daeninck, G. , and Winckelmans, G. , 2004, “ Simulation of Three-Dimensional Bluff-Body Flows Using the Vortex Particle and Boundary Element Methods,” Flow, Turbul. Combust., 73(2), pp. 117–131. [CrossRef]
Daeninck, G. , Chatelain, P. , Rubel, M. , Winckelmans, G. , and Leonard, A. , 2004, “ Simulation of Vehicle Aerodynamics Using a Vortex Element Method,” The Aerodynamics of Heavy Vehicles: Trucks, Buses, and Trains (Lecture Notes in Applied and Computational Mechanics), Springer, Berlin, pp. 353–366.
Rood, E. P. , 1994, “ Interpreting Vortex Interactions With a Free Surface,” ASME J. Fluids Eng., 116(1), pp. 91–94. [CrossRef]
Nakos, D. E. , and Sclavounos, P. D. , 1990, “ Steady and Unsteady Ship Wave Patterns,” J. Fluid Mech., 215, pp. 263–288. [CrossRef]
Costa, A. J. , Kowalyshyn, D. , Tuil, K. , Young, Y. L. , Milewski, W. M. , and Kring, D. C. , 2015, “ Influence of Free Surface on Hydrodynamic Lift Using a High-Order Boundary Element Method,” 13th International Conference on Fast Sea Transportation (Fast 2015), Washington, DC, Sept. 1–4. http://docplayer.net/47339106-Influence-of-free-surface-on-hydrodynamic-lift-using-a-high-order-boundary-element-method.html
Winckelmans, G. S. , and Leonard, A. , 1993, “ Contributions to Vortex Particle Methods for the Computation of Three-Dimensional Incompressible Unsteady Flows,” J. Comput. Phys., 109(2), pp. 247–273. [CrossRef]
Kring, D. C. , 1994, “ Time Domain Ship Motions by a Three-Dimensional Rankine Panel Method,” Ph.D. dissertation, Massachusetts Institute of Technology, Cambridge, MA. https://dspace.mit.edu/handle/1721.1/11939
Kring, D. C. , Milewski, W. M. , and Fine, N. E. , 2004, “ Validation of a NURBS-Based BEM for Multihull Ship Seakeeping,” 25th Symposium on Naval Hydrodynamics, St. John's, NL, Canada.
Abbott, I. H. , and von Doenhoff, A. E. , 1959, Theory of Wing Sections: Including Summary of Airfoil Data, Dover, Mineola, NY.
Theodorsen, T. , 1949, “ General Theory of Aerodynamic Instability and the Mechanism of Flutter,” National Advisory Committee for Aeronautics, Washington, DC, Report No. 496.
Kennedy, R. , Helfers, D. , and Young, Y. L. , 2015, “ A Reduced-Order Model for an Oscillating Hydrofoil Near the Free Surface,” 13th International Conference on Fast Sea Transportation, Washington, DC, Sept. 1–4.
Muscari, R. , Felli, M. , and Di Mascio, A. , 2011, “ Analysis of the Flow Past a Fully Appended Hull With Propellers by Computational and Experimental Fluid Dynamics,” ASME J. Fluids Eng., 133(6), p. 061104. [CrossRef]
Mook, D. T. , and Dong, B. , 1994, “ Perspective: Numerical Simulations of Wakes and Blade-Vortex Interaction,” ASME J. Fluids Eng., 116(1), pp. 5–21. [CrossRef]
Young, Y. L. , and Kinnas, S. A. , 2001, “ A BEM for the Prediction of Unsteady Midchord Face and/or Back Propeller Cavitation,” ASME J. Fluids Eng., 123(2), pp. 311–319. [CrossRef]


Grahic Jump Location
Fig. 2

Converting dipole panels to particles

Grahic Jump Location
Fig. 3

Positions of the vortex particles canceling the end vortex of the dipole wake

Grahic Jump Location
Fig. 8

Pressure distributions at the root of the hydrofoil of instantaneous 1 deg, 3 deg, and 5 deg pitch runs from presented VPM compared to the pre-existing dipole wake method

Grahic Jump Location
Fig. 5

Forced heaving NACA0012 with VPM compared to Theodorsen theoretical 2D lift force

Grahic Jump Location
Fig. 6

Wake pattern of an instantaneous pitch run from presented VPM using the high-order approach

Grahic Jump Location
Fig. 7

Wake patterns of an instantaneous pitch run from the presented VPM for 1 deg (left) and 5 deg (right)

Grahic Jump Location
Fig. 4

Steady forcing convergence study of the VPM compared to Prandtl's lifting line theory

Grahic Jump Location
Fig. 12

NACA0015 forced to pitch at a depth to chord ratio 0.5. Notice that as the particles in the wake of the foil begin to oscillate dramatically near the free surface, they pierce the boundary and cross above the surface.

Grahic Jump Location
Fig. 9

Resulting particle advection from forced heave simulation of aspect ratio 1 hydrofoil

Grahic Jump Location
Fig. 10

Lift coefficient time traces for a 3D NACA0015 foil pitching between ±5 deg with reduced pitching frequencies (a) k = 0.5, (b) k = 1.0, and (c) k = 2.0

Grahic Jump Location
Fig. 11

Forced pitching NACA0015 with VPM at depth to chord ratio: (a) 1.0, (b) 0.75, and (c) 0.5

Grahic Jump Location
Fig. 13

NACA0015 forced to pitch at depth to chord ratio 0.5. The particles are observed to cross above the free surface.

Grahic Jump Location
Fig. 14

Steady forcing VPM time domain simulation with downstream bodies

Grahic Jump Location
Fig. 15

Simulation with free surface and multiple hydrofoils



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