Statistical equilibrium of bubble oscillations in dilute bubbly flows

The problem of predicting the moments of the distribution of bubble radius in bubbly flows is considered. The particular case where bubble oscillations occur due to a rapid (impulsive or step change) change in pressure is analyzed, and it is mathematically shown that in this case, inviscid bubble oscillations reach a stationary statistical equilibrium, whereby phase cancellations among bubbles with different sizes lead to time-invariant values of the statistics. It is also shown that at statistical equilibrium, moments of the bubble radius may be computed using the period-averaged bubble radius in place of the instantaneous one. For sufficiently broad distributions of bubble equilibrium (or initial) radius, it is demonstrated that bubble statistics reach equilibrium on a time scale that is fast compared to physical damping of bubble oscillations due to viscosity, heat transfer, and liquid compressibility. The period-averaged bubble radius may then be used to predict the slow changes in the moments caused by the damping. A benefit is that period averaging gives a much smoother integrand, and accurate statistics can be obtained by tracking as few as five bubbles from the broad distribution. The period-averaged formula may therefore prove useful in reducing computational effort in models of dilute bubbly flow wherein bubbles are forced by shock waves or other rapid pressure changes, for which, at present, the strong effects caused by a distribution in bubble size can only be accurately predicted by tracking thousands of bubbles. Some challenges associated with extending the results to more general (nonimpulsive) forcing and strong two-way coupled bubbly flows are briefly discussed.