Abstract: Based on the boundary data immersion method (BDIM), a numerical solver for multiphase ventilated cavitation flow around complex axisymmetric bodies with cavitators by combining implicit large-eddy simulation (ILES) with non-conformal Cartesian grids was developed. The numerical method’s accuracy and effectiveness in simulating ventilated cavitation around complex axisymmetric geometries were validated through experimental comparisons. Numerical simulations were then conducted to analyze ventilated cavitation under different ventilation rates. The results demonstrate that ventilation rate significantly affects the evolution of cavity morphology and vortex structures. Increasing ventilation rate enhances gravitational effects on cavity morphology, reduces near-wall velocity gradients of the axisymmetric body, and accelerates the breakup and shedding of annular vortex structures. Furthermore, variations in ventilation rate lead to corresponding changes in both magnitude and distribution of enstrophy and turbulent kinetic energy transport, while maintaining a dynamic balance between energy production and dissipation throughout the process.