Numerical simulation of the unsteady behaviour of cavitating flows

TLDR

The authors propose a two‑dimensional numerical model, including a modified k‑ε turbulence model, to simulate unsteady cavitating flows. The model solves Reynolds‑averaged Navier–Stokes equations for a variable‑density liquid–vapour mixture using a barotropic state law, a SIMPLE‑derived pressure‑correction finite‑volume scheme, and a cavitation‑aware modification applied to two Venturi sections. The modified k‑ε model eliminates non‑physical stable cavities and accurately reproduces experimental unsteady cavitation behaviour in 2‑D Venturi geometries. © 2003 John Wiley & Sons, Ltd.

Abstract

Abstract A 2D numerical model is proposed to simulate unsteady cavitating flows. The Reynolds‐averaged Navier–Stokes equations are solved for the mixture of liquid and vapour, which is considered as a single fluid with variable density. The vapourization and condensation processes are controlled by a barotropic state law that relates the fluid density to the pressure variations. The numerical resolution is a pressure‐correction method derived from the SIMPLE algorithm, with a finite volume discretization. The standard scheme is slightly modified to take into account the cavitation phenomenon. That numerical model is used to calculate unsteady cavitating flows in two Venturi type sections. The choice of the turbulence model is discussed, and the standard RNG k–εmodel is found to lead to non‐physical stable cavities. A modified k–εmodel is proposed to improve the simulation. The influence of numerical and physical parameters is presented, and the numerical results are compared to previous experimental observations and measurements. The proposed model seems to describe the unsteady cavitation behaviour in 2D geometries well. Copyright © 2003 John Wiley & Sons, Ltd.