Compared to non-cavitating flow, cavitating flow is much complex owing to the numerical difficulties caused by cavity generation and collapse. In this paper, the cavitating flow around a NACA66 hydrofoil is studied numerically with particular emphasis on understanding the cavitation structures and the shedding dynamics. Large Eddy Simulation (LES) was coupled with a homogeneous cavitation model to calculate the pressure, velocity, vapor volume fraction and vorticity around the hydrofoil. The predicted cavitation shedding dynamics behavior, including the cavity growth, break-off and collapse downstream, agrees fairly well with experiment. Some fundamental issues such as the transition of a cavitating flow structure from 2D to 3D associated with cavitation–vortex interaction are discussed using the vorticity transport equation for variable density flow. A simplified one-dimensional model for the present configuration is adopted and calibrated against the LES results to better clarify the physical mechanism for the cavitation induced pressure fluctuations. The results verify the relationship between pressure fluctuations and the cavity shedding process (e.g. the variations of the flow rate and cavity volume) and demonstrate that the cavity volume acceleration is the main source of the pressure fluctuations around the cavitating hydrofoil. This research provides a better understanding of the mechanism driving the cavitation excited pressure pulsations, which will facilitate development of engineering designs to control these vibrations

Cavitating turbulent flow around hydrofoils was simulated using the Partially-Averaged Navier–Stokes (PANS) method and a mass transfer cavitation model with the maximum density ratio (ρl/ρv,clip) effect between the liquid and the vapor. The predicted cavity length and thickness of stable cavities as well as the pressure distribution along the suction surface of a NACA66(MOD) hydrofoil compare well with experimental data when using the actual maximum density ratio (ρl/ρv,clip=43391) at room temperature. The unsteady cavitation patterns and their evolution around a Delft twisted hydrofoil were then simulated. The numerical results indicate that the cavity volume fluctuates dramatically as the cavitating flow develops with cavity growth, destabilization, and collapse. The predicted three dimensional cavity structures due to the variation of attack angle in the span-wise direction and the shedding cycle as well as its frequency agree fairly well with experimental observations. The distinct side-lobes of the attached cavity and the shedding U-shaped horse-shoe vortex are well captured. Furthermore, it is shown that the shedding horse-shoe vortex includes a primary U-shaped vapor cloud and two secondary U-shaped vapor clouds originating from the primary shedding at the cavity center and the secondary shedding at both cavity sides. The primary shedding is related to the collision of a radially-diverging re-entrant jet and the attached cavity surface, while the secondary shedding is due to the collision of side-entrant jets and the radially-diverging re-entrant jet. The local flow fields show that the interaction between the circulating flow and the shedding vapor cloud may be the main mechanism producing the cavitating horse-shoe vortex. Two side views described by iso-surfaces of the vapor volume fraction for a 10% vapor volume, and a non-dimensional Q-criterion equal to 200 are used to illustrate the formation, roll-up and transport of the shedding horse-shoe vortex. The predicted height of the shedding horse-shoe vortex increases as the vortex moves downstream. It is shown that the shape of the horse-shoe vortex for the non-dimensional Q-criterion is more complicated than that of the 10% vapor fraction iso-surface and is more consistent with the experiments. Further, though the time-averaged lift coefficient predicted by the PANS calculation is about 12% lower than the experimental value, it is better than other predictions based on RANS solvers

This study presented a multi-scale Eulerian–Lagrangian approach to simulate cavitating turbulent flow around a Clark-Y hydrofoil to study the bubble dynamics. A large eddy simulation was coupled with the volume of fluid method to capture the large vapor volumes in an Eulerian analysis. Micro-scale Lagrangian bubbles were then tracked by solving compressible the Rayleigh–Plesset equation and a bubble motion equation. A Gaussian kernel function was used to model the interactions between the flow field and the vapor bubbles in a coupled two-way algorithm. The predictions give satisfactory agreement with experimental data for the bubble size oscillations, bubble motion, and cavity shedding characteristics. Further investigations analyzed the influence of various parameters on the transformation between the Euler and Lagrange models. The numerical results provide detailed information about the influence of the cavitating turbulent flow on the bubble behavior, especially how the reentrant jet significantly affects the bubble generation and motion. The calculations also capture the bubble size oscillations caused by the surrounding liquid pressure variations and how these generate very high local pressures near the surface. The results show that the pressure wave released as a bubble is compressed reaches 10⁷ Pa, which may cause cavitation erosion of the hydrofoil surface. This research provides a promising method to better investigate the bubble motion characteristics in macroscopic flows and demonstrates that the cavitation erosion caused by bubble size oscillations is significant and deserves attention.