The modified Rayleigh equation is integrated to obtain the bubble wall velocities.
The coupling strength during the evaporation and condensation are different.
The second-order derivative is treated by recording the initial void fraction.
This cavitation model provides high accuracy and strong numerical stability.