In this paper, we sought to investigate unsteady cavitating turbulent flows around a twisted three-dimensional NACA 16012 hydrofoil using the large eddy simulation (LES) and volume of fluid (VOF) methods. The one equation eddy viscosity model (OEEVM) was used to calculate the sub-grid scale (SGS) stress tensor. The compressive velocity VOF technique and the Kunz cavitation model were also employed. Numerical simulation was performed using the interPhaseChangeFoam solver within the OpenFOAM ramework. Cavitation simulation was performed at three cavitation numbers within the cloud cavitation regime; that is to say, cavitation numbers: 0.95, 1.15, and 1.35. Afterwards, the details of flow predictions including the pressure, velocity, streamlines, volume fraction of water phase, and turbulent kinetic energy were reported. Cavity behaviours, including growth, shedding, and collapse of the cavity, were considered in detail. Three-dimensional cavity structures such as primary and secondary shedding and the shedding of the U-shaped horseshoe vortex were reported as well. The present work illustrated the side-entrant jets and the radially diverging re-entrant jet corresponding to the three-dimensional effect of the twisted wing. The cavity pattern and the shedding cycle frequency agreed well with the available experimental observations.

Keywords