A new cavitating model by using bubble size disla'ibution based on mass of bubbles is proposed. Liquid phase istreated with Eulerian framework as a mixture containing minute cavitating bubbles. Vapor phase consists ofvarious sizes of minute vapor bubbles, which is distributed to classes based on their mass. The change of bubblenumber density for each class was solved by considering the change of bubble mass due to phase change as wellas generation of new bubbles due to heterogeneous nucleation. In this method the mass of bubbles is treated as anindependent variable, in other word, a new coordinate, and dependant variables are solved in Eulerian frameworkfor spatial coordinates and bubble-mass coordinate. The present method is applied to a cavitating flow in aconvergent-divergent nozzle, and the two-phase flow with bubble size distribution and phase change wassuccessfully predicted.