In this paper we report progress toward developing advanced cavitation models with Eulerian multifluid method. The bubble number density and the interfacial area equations are introduced into the general framework of multifluid method for multiphase cavitating flows in order to account for the variable size nature of the bubble distribution. The physical models for bubble breakup and coalescence are based on the work by Ishii group’s work in recent years. Simulation results of a cavitating flow are compared with the corresponding experimental data, which include the bubble size distribution, bubble volume fraction and bubble number density. The ability of predicting bubble distribution characteristics is particularly useful as an input for cavitation erosion analysis.Copyright © 2005 by ASME