A compressible multi-scale model to simulate cavitating flows

We propose a compressible multi-scale model that (i) captures the dynamics of both large vapour cavities (resolved vapour) and micro-bubbles (unresolved vapour), and (ii) accounts for medium compressibility. The vapour mass, momentum and energy in the compressible homogeneous mixture equations are explicitly decomposed into constituent resolved and unresolved components that are independently treated. The homogeneous mixture of liquid and resolved vapour is tracked as a continuum in an Eulerian sense. The unresolved vapour terms are expressed in terms of subgrid bubble velocities and radii that are tracked in a Lagrangian sense using a novel kR-RP equation' (k, constant multiple; R, bubble size; RP, Rayleigh-Plesset). The kR-RP equation is formally derived in terms of the pressure at a finite distance (kR) from the bubble while accounting for the effects of neighbouring bubbles; p(kR) may therefore be either a near-field or far-field pressure. The equation exactly recovers the classical Rayleigh-Plesset and Keller-Miksis equations in the limits that k and c (speed of sound) become very large. Also, the results are independent of k for a single bubble for all k, and for multiple bubbles when kR < d (where d denotes separation distance). Numerical results show this robustness of the model to the choice of k, which can be different for each bubble. The multi-scale model is validated for the collapse of a single resolved/unresolved bubble. Its ability to capture inter-bubble interactions is demonstrated for multiple bubbles exposed to an acoustic pulse. The model is then applied to a problem where resolved and unresolved bubbles co-exist. Finally, it is validated using a cluster of 1200 bubbles exposed to a strong acoustic pulse. The results show the impact of the bubble cluster on the transmitted and reflected waves and the shielding effect where bubbles at the edge of the cluster shield the interior bubbles by dampening the incident acoustic wave.