Accurate predictions of particle size distributions, and therefore of the underlying processes of fluid particle breakup and coalescence are of vital importance in process design, but reliable procedures are still lacking. The present paper aims at developing a modular formulation for the turbulent particle breakup process. The model is to be included in a population balance model which is formulated such as to facilitate the direct future implementation into a full multifluid CFD model. The breakup process is described without introducing adjustable parameters. The current model is a further development of an existing model by Luo and Svendsen (AlChE J. 42 (5) (1996) 1225), which has been expanded and refined, and where an inherent weakness regarding the breakup rate for small particles and small daughter particle fragments are removed. A new criterion regarding the kinetic energy density of the colliding turbulent eddy causing breakup has been introduced. This new criterion is a novel concept describing the breakup process. The details are thoroughly discussed together with possible further modifications. The results from the new model are encouraging because the breakup rate is greatly reduced when the dispersed fluid particles are reduced in size. Further, the response to changes in system variables is reasonable and the distribution of daughter sizes vary in a reasonable way for the different collision possibilities.