Based on energy principies, we propose a statistical model to describe the bubble size probability density function of the daughter bubbles resulting from the shattering of a mother bubble of size Do immersed in a fully developed turbulent water flow. The model shows that the bubble size p.d.f. depends not only on Do, but also on the valué of the dissipation rate of turbulent kinetic energy of the underlying turbulenee of the water, e. The phenomenological model is simple, yet it prediets detailed experimental measurements of the transient bubble size p.d.f.s performed over a range of bubble sizes and dissipation rates e in a very consistent manner. The agreement between the model and the experiments is particularly good for low and modérate bubble turbulent Weber numbers, We t = pAu 2 (D 0 )D 0 /a where the assumption of the binary breakup is shown to be consistent with the experimental observations. At larger valúes of We t , it was found that the most probable number of daughter bubbles increases and the assumption of tertiary breakup is shown to lead to a better fit of the experimental measurements.