As the 1st step of the development of supportive software to design a cryogenic high-efficiency pump allowed cavitation, a bubble size distribution model has been developed by the author for this software. In this model, bubble growth/decay and bubble advection are solved for each class of the bubble size, strictly mass, when there are various mass bubbles in the same calculation region. The above calculations are treated as Eulerian approach with respect to the bubble mass. The numerical results based on this model have agreed with the experimental results as a whole, however, some inconsistency still remained. It is suspected that the model of the bubble growth/decay causes the difference between the numerical and experimental results because heat transfer around the bubble was approximately computed by an analytical solution of unsteady heat transfer based on the elapsed-time from the bubble nucleation. In this paper, a new bubble size distribution model was developed, in which the bubble growth/decay calculation employ a new method combining two rigorous methods, namely, a Rayleigh-Plesset equation for bubble oscillation, and a heat conduction equation in a thermal boundary layer around the bubble to evaluate mass rate of evaporation/condensation.