Cavitation is described by formation and collapse of the bubbles in a liquid when the ambient pressure decreases. Formed bubbles grow and collapse by change of pressure, and when they collapse, shockwave by high pressure is generated. In general, bubble behavior can be described by Rayleigh-Plesset equation under adiabatic or isothermal condition and hence, phase shift by the pressure change in a bubble cannot be considered in the equation. In our study, a numerical model is developed from the mathematical model considering the phase shift from the previous study. In the developed numerical model, size of single spherical bubble is calculated by the change of mass calculated from the change of the ambient pressure in a liquid. The developed numerical model is verified by a case of liquid flow in a narrow channel.