Document Type : Original Article


Department of Aerospace Engineering, K. N. Tusi University of Technology


 One of the most effective ways of high-speed motion in water is the motion in the supercavitation regime. This way provides the possibility to avoid considerable viscose resistance of boundary layer and consequently reach to very small drag coefficient which can be several times smaller than, that of the continuous flow. In this study the numerical simulation of developed and supercavitating flow is performed. The CFX code which served as a platform for the present work is a three-dimensional code that solves the Reynolds-Averaged Navier-Stokes equations with a finite volume method. The cavitation model is implemented based on the use of Rayleigh-Plesset equation to estimate the rate of vapor production. A high Reynolds number form ĸ-ε model is implemented to provide turbulence closure. For steady state flows and poor mesh resolution near the wall (using log-law wall functions), there is a priori no difference between the two equations formulations. For the different case studies, multi-block structured meshes were generated and the numerical simulation is performed in a wide range of cavitation numbers. Results are presented for steady state flows with natural cavitation about various bodies. Comparisons are made with available measurement of surface pressure distribution, cavitation bubble geometry (cavity length and cavity width) and drag coefficient. The simulated results are in a good agreement with the experimental data. Finally, the three-dimensional results are presented for a submerged body running at several angles of attack.


Main Subjects

[1] N.H.Singhal, A.K.Athavale, M.Li and Y.Jiang, Mathematical basis and validation of the full cavitation model, J.Fluids Eng., Vol. 124, 1-8, (2002).
[2] C.L.Merkle, J.Feng, and P.E.O.Buelow, Computational modeling of the dynamics of sheet cavitation, Proceeding of the 3rd International Symposium on Cavitation, (CAV98), Gernoble, France, (1998).
[3] R.F.Kunse, D.A.Boger, D.A.Stinberg, T.S.Chyczewski, J.W.Lindau, and H.J.Gibeling, A preconditioned Navier-Stokes method for two-phase flows with application to cavitation, Computers & Fluids, Vol. 29, 849-875, (2000).
[4] R.F.Kunse, D.A.Boger, D.A.Stinberg, T.S.Chyczewski, H.J.Gibeling, S.Vencateswaran, and T.A.Govindan, A preconditionated Navier-Stokes method for two-phase flows with application to cavitation prediction, Computers & Fluids, 29, 849-875 (2000).
[5] J.W.Lindau, R.F.Kunse, Advancement and application of multiphase CFD modeling to high speed supercavitating flows, Final report SF298 final report, (2004).
[6] J.Sauer and G.H.Schnerr, Unsteady cavitating flow – a new cavitation model based on modified front capturing method and bubble dynamics, In Proc. Of  FEDSM00 4th Fluids Engineering Summer Conference, (2000).
[7] I.Senocak, and W.Shyy, A pressure-based method for turbulent cavitating flow computations, J, Comput. Phys, 176, 363–383, (2002).
[8] I.Senocak, Computational Methodology for the Simulation Turbulent Cavitating Flows. PhD thesis, University of Florida, (2002).
[9]CFX-ANSYS Canada Ltd. CFX-5.7: Solver Theory. Canada, (2004).
[10] H.K.Versteegand and W.Malalasekera,  An introduction to Computational Fluid Dynamics: The finite volume method, vol. Pearson Prentice Hall, (1995).
[11] H.Rouse and J.S.McNwon, Cavitation and pressure distribution head form at zero of angle of yaw, Studies in engineering bulletin 32, state university of lowa, (1948).
[12] M.Self and J.F.Ripken, steady-state cavity studies in a free-jet water tunnel, St. Anthony Falls Hydr. Lab. rep. 47, (1955).
[13] Robert T. Knapp, James W. Daily, Frederick G. Hammit, Cavitation, McGraw Hill Inc, (1970).
[14] C.Brennen, A numerical solution of axisymmtric cavity flows, Fluid Mechanics, 37, 671-688, (1968).