This paper reviews recent progress made toward modeling of cavitation and numerical simulation of cavitating water jets. Properties of existing cavitation models are discussed and a compressible mixture flow method fo...This paper reviews recent progress made toward modeling of cavitation and numerical simulation of cavitating water jets. Properties of existing cavitation models are discussed and a compressible mixture flow method for the numerical simulation of high- speed water jets accompanied by intensive cavitation is introduced. Two-phase fluids media of cavitating flow are treated as a homo- geneous bubbly mixture and the mean flow is computed by solving Reynolds-Averaged Navier-Stokes （RANS） equations for com- pressible fluid. The intensity of cavitation is evaluated by the gas volume fraction, which is governed by the compressibility of bubble-liquid mixture corresponding to the status of mean flow field. Numerical results of cavitating water jet issuing from an orifice nozzle are presented and its applicability to intensively cavitating jets is demonstrated. However, the effect of impact pressure caused by collapsing of bubbles is neglected, and effectively coupling of the present compressible mixture flow method with the dynamics of bubbles remains to be a challenge.展开更多
Cavitation is a complex multiphase flow phenomenon with an abrupt transient phase change between the liquid and the vapor, including multiscale vortical motions. The transient cavitation dynamics is closely associated...Cavitation is a complex multiphase flow phenomenon with an abrupt transient phase change between the liquid and the vapor, including multiscale vortical motions. The transient cavitation dynamics is closely associated with the evolution of the cavitation vortex structures. The present paper investigates the cavitation vortex dynamics using different vortex identification methods, including the vorticity method, the Q criterion method, the Omega method (Ω), the method and the Rortex method. The Q criterion is an eigenvalue-based criterion, and in the Ω method, the parameter is normalized, is independent of the threshold value and in most conditions Ω= 0.52 . The Rortex method is based on an eigenvector-based criterion. Numerical simulations are conducted using the implemented compressible cavitation solver in the open source software OpenFOAM for the sheet/cloud cavitating flows around a NACA66 (mod) hydrofoil fixed at a = 6°,= 1.25 and Re = 7.96 × 10^5 . The flow is characterized by the alternate interactions of the re-entrant flow and the collapse induced shock wave. Results include the vapor structures and the vortex dynamics in the unsteady sheet/cloud cavitating flows, with emphasis on the vortex structures in thecavitation region, the cavity interface, the cavity closure, the cavity wakes, and the foil wakes with the shedding cavity. The comparisons of the various methods, including that the vorticity method, the Q criterion method, the Ω method, the λ2 method and the Rortex method, show the performances of different methods in identifying the cavitation vortex structures. Generally, during the attached cavity growth stage, the Q criteria can well predict the vortex structures in the cavitation region and at the foil trailing edge in the pure liquid region, while with the Ω method and the Rortex method, the vortex structures outside the attached cavity and on the foil pressure side can also be predicted. The λ2 method can well predict the vortex structures in the cavity closure region. During the re展开更多
基金supported by the Japan Society for the Promotion of Science (JSPS), Grant-in-Aid for Scientific Research (C) (Grant No. 22560177)
文摘This paper reviews recent progress made toward modeling of cavitation and numerical simulation of cavitating water jets. Properties of existing cavitation models are discussed and a compressible mixture flow method for the numerical simulation of high- speed water jets accompanied by intensive cavitation is introduced. Two-phase fluids media of cavitating flow are treated as a homo- geneous bubbly mixture and the mean flow is computed by solving Reynolds-Averaged Navier-Stokes （RANS） equations for com- pressible fluid. The intensity of cavitation is evaluated by the gas volume fraction, which is governed by the compressibility of bubble-liquid mixture corresponding to the status of mean flow field. Numerical results of cavitating water jet issuing from an orifice nozzle are presented and its applicability to intensively cavitating jets is demonstrated. However, the effect of impact pressure caused by collapsing of bubbles is neglected, and effectively coupling of the present compressible mixture flow method with the dynamics of bubbles remains to be a challenge.
基金supported by the National Natural Science Foundation of China (Grant Nos. 51839001, 91752105).
文摘Cavitation is a complex multiphase flow phenomenon with an abrupt transient phase change between the liquid and the vapor, including multiscale vortical motions. The transient cavitation dynamics is closely associated with the evolution of the cavitation vortex structures. The present paper investigates the cavitation vortex dynamics using different vortex identification methods, including the vorticity method, the Q criterion method, the Omega method (Ω), the method and the Rortex method. The Q criterion is an eigenvalue-based criterion, and in the Ω method, the parameter is normalized, is independent of the threshold value and in most conditions Ω= 0.52 . The Rortex method is based on an eigenvector-based criterion. Numerical simulations are conducted using the implemented compressible cavitation solver in the open source software OpenFOAM for the sheet/cloud cavitating flows around a NACA66 (mod) hydrofoil fixed at a = 6°,= 1.25 and Re = 7.96 × 10^5 . The flow is characterized by the alternate interactions of the re-entrant flow and the collapse induced shock wave. Results include the vapor structures and the vortex dynamics in the unsteady sheet/cloud cavitating flows, with emphasis on the vortex structures in thecavitation region, the cavity interface, the cavity closure, the cavity wakes, and the foil wakes with the shedding cavity. The comparisons of the various methods, including that the vorticity method, the Q criterion method, the Ω method, the λ2 method and the Rortex method, show the performances of different methods in identifying the cavitation vortex structures. Generally, during the attached cavity growth stage, the Q criteria can well predict the vortex structures in the cavitation region and at the foil trailing edge in the pure liquid region, while with the Ω method and the Rortex method, the vortex structures outside the attached cavity and on the foil pressure side can also be predicted. The λ2 method can well predict the vortex structures in the cavity closure region. During the re
