Fast integration method for averaging polydisperse bubble population dynamics

Ensemble-averaged polydisperse bubbly flow models require statistical moments of the evolving bubble size distribution. Under step forcing, these moments reach statistical equilibrium in finite time. However, the transitional phase before equilibrium and cases with time-dependent forcing are required to predict flow in engineering applications. Computing these moments is expensive because the integrands are highly oscillatory, even when the bubble dynamics are linear. Ensemble-averaged models compute these moments at each grid point and time step, making cost reduction important for large-scale bubbly flow simulations. Traditional methods evaluate the integrals via traditional quadrature rules. This approach requires a large number of quadrature nodes in the equilibrium bubble size, each equipped with its own advection partial differential equation (PDE), resulting in significant computational expense. We formulate a Levin collocation method to reduce this cost. Given the differential equation associated with the integrand, or moment, the method approximates it by evaluating its derivative via polynomial collocation. The differential matrix and amplitude function are well-suited to numerical differentiation via collocation, and so the computation is comparatively cheap. For an example excited polydisperse bubble population, the first moment is computed with the presented method at $10^{-3}$ relative error with 100 times fewer quadrature nodes than the trapezoidal rule. The gap increases for smaller target relative errors: the Levin method requires $10^4$ times fewer points for a relative error of $10^{-8}$. The formulated method maintains constant cost as the integrands become more oscillatory with time, making it particularly attractive for long-time simulations.