Research Papers: Multiphase Flows

Finite Volume, Computational Fluid Dynamics-Based Investigation of Supercavity Pulsations

[+] Author and Article Information
Grant M. Skidmore

Penn State Applied Research Laboratory,
P.O. Box 30,
State College, PA 16804-0030;
Melbourne School of Engineering,
University of Melbourne,
Parkville, VIC 3052, Australia
e-mails: skidmore.grant@gmail.com;

Jules W. Lindau

Penn State Applied Research Laboratory,
P.O. Box 30,
State College, PA 16804-0030
e-mail: jwl10@arl.psu.edu

Timothy A. Brungart

Penn State Applied Research Laboratory,
P.O. Box 30,
State College, PA 16804-0030
e-mail: tab7@arl.psu.edu

Michael J. Moeny

Penn State Applied Research Laboratory,
P.O. Box 30,
State College, PA 16804-0030
e-mail: mjm369@arl.psu.edu

Michael P. Kinzel

Penn State Applied Research Laboratory,
P.O. Box 30,
State College, PA 16804-0030
e-mail: mpk176@arl.psu.edu

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received May 11, 2016; final manuscript received March 29, 2017; published online June 20, 2017. Assoc. Editor: Riccardo Mereu.The United States Government retains, and by accepting the article for publication, the publisher acknowledges that the United States Government retains, a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this work, or allow others to do so, for United States Government purposes.

J. Fluids Eng 139(9), 091301 (Jun 20, 2017) (10 pages) Paper No: FE-16-1300; doi: 10.1115/1.4036596 History: Received May 11, 2016; Revised March 29, 2017

Computations of pulsating supercavity flows behind axisymmetric disk cavitators are presented. The method of computation is a finite volume discretization of the equations of mixture fluid motion. The gas phase is treated as compressible, the liquid phase as incompressible, and the interface accuracy enhanced using a volume of fluid (VOF) approach. The re-entrant, pulsating, and twin vortex modes of cavity closure are delineated and computationally resolved, including the expected hysteresis. A phase diagram of cavitation number versus ventilation rate at three Froude conditions is computationally constructed. Sample re-entrant, pulsation, and twin vortex snapshots are presented. Pulsation results are compared with stability criterion from the literature as well as examined for their expected character. Computations appear to capture the complete spectrum of cavity closure conditions. A detailed comparison of computational simulation and physical experiment at similar conditions is also included as a means to validate the computational results.

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


Song, C. , 1962, “ Pulsation of Ventilated Cavities,” J. Ship Res., 5(4), pp. 8–20. http://www.sname.org/HigherLogic/System/DownloadDocumentFile.ashx?DocumentFileKey=c7be08ee-f3aa-4474-bdfa-9a55f8dd73c6
Paryshev, E. , 2006, “ Approximate Mathematical Models in High-Speed Hydrodynamics,” J. Eng. Math., 55(1–4), pp. 41–64. [CrossRef]
Savchenko, Y. , 2001, “ Supercavitation-Problems and Perspectives,” 4th International Symposium on Cavitation (CAV), Pasadena, CA, June 20–23, Paper No. 66. http://caltechconf.library.caltech.edu/66/1/Savchenko.pdf
Campbell, I. , and Hilborne, D. , 1958, “ Air Entrainment Behind Artificially Ventilated Cavities,” 2nd Symposium on Naval Hydrodynamics, Washington, DC, Aug. 25–29, pp. 467–481.
Michel, J. M. , 1984, “ Some Features of Water Flows With Ventilated Cavities,” ASME J. Fluids Eng., 106(3), pp. 319–326. [CrossRef]
Skidmore, G. M. , Brungart, T. A. , Lindau, J. W. , and Moeny, M. J. , 2015, “ Noise Generated by Ventilated Supercavities,” Noise Control Eng. J., 63(1), pp. 94–101. [CrossRef]
Skidmore, G. , 2013, “ The Pulsation of Ventilated Supercavities,” M.S. thesis, The Pennsylvania State University, University Park, PA. https://etda.libraries.psu.edu/catalog/18680
Karn, A. , Arndt, R. A. , and Hong, J. , 2016, “ An Experimental Investigation Into Supercavity Closure Mechanisms,” J. Fluid Mech., 789, pp. 259–284. [CrossRef]
Kunz, R. F. , Boger, D. A. , Stinebring, D. R. , Chyczewski, T. S. , Lindau, J. W. , Gibeling, H. J. , Venkateswaran, S. , and Govindan, T. R. , 2000, “ A Preconditioned Navier–Stokes Method for Two-Phase Flows With Application to Cavitation Prediction,” Comput. Fluids, 29(8), pp. 849–875. [CrossRef]
Venkateswaran, S. , Lindau, J. W. , Kunz, R. F. , and Merkle, C. L. , 2002, “ Computation of Multiphase Mixture Flows With Compressibility Effects,” J. Comput. Phys., 180(1), pp. 54–77. [CrossRef]
Senocak, I. , and Shyy, W. , 2002, “ A Pressure-Based Method for Turbulent Cavitating Flow Computations,” J. Comput. Phys., 176(2), pp. 363–383. [CrossRef]
Venkateswaran, S. , and Merkle, C. L. , 1998, “ Evaluation of Artificial Dissipation Models and Their Relationship to the Accuracy of Euler and Navier–Stokes Computations,” 16th International Conference on Numerical Methods in Fluid Dynamics, Arcachon, France, July 6–10, pp. 427–432.
CD-Adapco, 2014, “ STAR-CCM+ 9.02.005 User Guide,” CD-Adapco, Melville, NY.
Kinzel, M. P. , 2008, “ Computational Techniques and Analysis of Cavitating-Fluid Flows,” Ph.D. thesis, The Pennsylvania State University, State College, PA. https://etda.libraries.psu.edu/catalog/8854
Wilcox, D. C. , 1998, Turbulence Modeling for CFD, 2nd ed., DCW Industries, La Canada, CA.
Demirdzic, I. , Muzaferija, S., Peric, M., Schreck, H., and Seidl, V., 1999, “ Computation of Flows with Free Surfaces,” Scientific Computing in Chemical Engineering II, F. Keil, W. Mackens, H. Vob, and J. Werther, eds., Springer, Berlin, pp. 360–367.
Phillips, T. S. , and Roy, C. J. , 2014, “ Richardson Extrapolation-Based Discretization Uncertainty Estimation for Computational Fluid Dynamics,” ASME J. Fluids Eng., 136(12), p. 121401. [CrossRef]
Logvinovich, G. V. , 1972, “ Hydrodynamics of Free-Boundary Flows,” Israel Program for Scientific Translations, Jerusalem, Israel, Report No. NASA-TT-F-658. https://archive.org/details/nasa_techdoc_19720006655
May, A. , 1975, “ Water Entry and the Cavity-Running Behavior of Missiles,” NAVSEA Hydroballistics Advisory Committee, Silver Spring, MD, Technical Report No. SEAHAC/TR 75-2. http://www.dtic.mil/dtic/tr/fulltext/u2/a020429.pdf
Spurk, J. , 2002, “ On the Gas Loss From Ventilated Supercavities,” Acta Mech., 155(3), pp. 125–135. [CrossRef]
Kinzel, M. , Lindau, J. , and Kunz, R. , 2009, “ Air Entrainment Mechanisms From Artificial Supercavities: Insight Based on Numerical Simulations,” 7th International Symposium on Cavitation (CAV), Ann Arbor, MI, Aug. 16–20, pp. 562–575. https://deepblue.lib.umich.edu/bitstream/handle/2027.42/84310/CAV2009-final136.pdf%3Bsequence=1
Lindau, J. W. , Kunz, R. F. , Sankaran, V. , Stinebring, D. R. , and Peltier, L. J. , 2005, “ Development and Application of Turbulent, Multiphase CFD to Supercavitation,” Second International Symposium on Seawater Drag Reduction (ISSDR), Busan, South Korea, May 23–26, pp. 221–238.


Grahic Jump Location
Fig. 1

Snapshots from water tunnel testing and corresponding CFD-based representations of supercavity closure modes. For CFD results, volume fraction isosurfaces of 50% gas, opaque, and 10% gas, translucent, are shown.

Grahic Jump Location
Fig. 2

Series of computational meshes applied to convergence study

Grahic Jump Location
Fig. 3

Mesh convergence study. Time histories of pressure at probe located 0.67 cavitator diameters aft of forward face.

Grahic Jump Location
Fig. 4

Pressure histories from probes distributed along length of pulsing, Fr = 24.5, CFD cavity. A 50% volume fraction gas isosurface is shown.

Grahic Jump Location
Fig. 5

Computational and experimental internal cavity pressure for twin vortex, re-entrant jet, and pulsating supercavities

Grahic Jump Location
Fig. 10

CFD representation of supercavities near the limit of linear Paryshev stability. Cavities illustrated with isosurface of gas volume fraction equal to 0.5.

Grahic Jump Location
Fig. 11

Drawings and photographs of water tunnel test configuration (dimensions in mm)

Grahic Jump Location
Fig. 12

Illustration of computational mesh employed to generate the Froude number matched pulsating supercavity simulation

Grahic Jump Location
Fig. 6

Ventilation rate versus cavitation number for CFD results at three Froude numbers

Grahic Jump Location
Fig. 7

Highly resolved CFD representation of ventilated twin vortex cavity (Fr = 12.3, CQ = 1.07, σ = 0.055). Cavity illustrated with gas volume fraction equal to 0.5 translucent isosurface. Mesh density illustrated on isosurface. Disk cavitator colored by pressure.

Grahic Jump Location
Fig. 8

Resolution of supercavities with incompressible and compressible finite volume computational physics. Cavities illustrated with 0.5 gas volume fraction isosurfaces.

Grahic Jump Location
Fig. 9

CFD representation of mode two pulsing cavity (Fr = 24.5, CQ = 0.268, σ = 0.045). Cavity illustrated with isosurface of gas volume fraction equal to 0.5. Progress of wave corresponding to pulsation illustrated with arrows.

Grahic Jump Location
Fig. 13

Froude 4.5, mode 2, pulsating supercavities: (a) image from ARL-PSU 0.305 m water tunnel, (b) CFD cavity illustrated with 0.5 gas volume fraction isosurface, and (c) CFD cavity illustrated with 0.25 gas volume fraction isosurface

Grahic Jump Location
Fig. 14

Pressure histories, Froude 4.5, mode 2, pulsating supercavity, water tunnel, and CFD

Grahic Jump Location
Fig. 15

Pressure spectra, Froude 4.5, mode 2, pulsating supercavity, water tunnel, and CFD



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