In this work, we established a novel theory for the dynamics of oscillating bubbles such as cavitation bubbles, underwater explosion bubbles, and air bubbles. For the first time, we proposed bubble dynamics equations that can simultaneously take into consideration the effects of boundaries, bubble interaction, ambient flow field, gravity, bubble migration, fluid compressibility, viscosity, and surface tension while maintaining a unified and elegant mathematical form. The present theory unifies different classical bubble equations such as the Rayleigh–Plesset equation, the Gilmore equation, and the Keller–Miksis equation. Furthermore, we validated the theory with experimental data of bubbles with a variety in scales, sources, boundaries, and ambient conditions and showed the advantages of our theory over the classical theoretical models, followed by a discussion on the applicability of the present theory based on a comparison to simulation results with different numerical methods. Finally, as a demonstration of the potential of our theory, we modeled the complex multi-cycle bubble interaction with wide ranges of energy and phase differences and gained new physical insight into inter-bubble energy transfer and coupling of bubble-induced pressure waves.

A novel theoretical model for bubble dynamics is established that simultaneously accounts for the liquid compressibility, phase transition, oscillation, migration, ambient flow field, etc. The bubble dynamics equations are presented in a unified and concise mathematical form, with clear physical meanings and extensibility. The bubble oscillation equation can be simplified to the Keller–Miksis equation by neglecting the effects of phase transition and bubble migration. The present theoretical model effectively captures the experimental results for bubbles generated in free fields, near free surfaces, adjacent to rigid walls, and in the vicinity of other bubbles. Based on the present theory, we explore the effect of the bubble content by changing the vapour proportion inside the cavitation bubble for an initial high-pressure bubble. It is found that the energy loss of the bubble shows a consistent increase with increasing Mach number and initial vapour proportion. However, the radiated pressure peak by the bubble at the collapse stage increases with decreasing Mach number and increasing vapour proportion. The energy analyses of the bubble reveal that the presence of vapour inside the bubble not only directly contributes to the energy loss of the bubble through phase transition but also intensifies the bubble collapse, which leads to greater radiation of energy into the surrounding flow field due to the fluid compressibility.

The strong bubble interactions and bursting behaviors near a free surface are studied numerically with a compressible two-phase flow solver. The interface is captured by the volume of fluid method. We investigate the effects of the dimensionless distance between the bubble and the free surface γ_f (scaled by the maximum bubble radius) ranging from 0.25 to 1.5. For the nonbursting cases, the essential evolution of the toroidal bubble is well captured, including the splitting, coalescence, and recollapse. Generally, a relatively thin spike is generated at the free surface during the first cycle of the bubble. Subsequently, a wider secondary spike at the base of the first spike is formed during the second bubble cycle, which leads to the formation of the crown-shaped spike. When γ_f is sufficiently small, the bubble bursts at the free surface and an opening cavity is generated. The pressure inside the cavity and the atmospheric pressure are not balanced immediately, and the pressure difference lasts for a while, causing the inward gas flow and the final closure of the cavity. The gas flow is found to play a vital role in the bubble bursting behaviors, which has not been well understood. By comparing with previous studies, three distinct bursting patterns are identified and discussed.

We experimentally, numerically and theoretically investigate the nonlinear interaction between a cavitation bubble and the interface of two immiscible fluids (oil and water) on multiple time scales. The underwater electric discharge method is utilized to generate a cavitation bubble near or at the interface. Both the bubble dynamics on a short time scale and the interface evolution on a much longer time scale are recorded via high-speed photography. Two mechanisms are found to contribute to the fluid mixing in our system. First, when a bubble is initiated in the oil phase or at the interface, an inertia-dominated high-speed liquid jet generated from the collapsing bubble penetrates the water–oil interface, and consequently transports fine oil droplets into the water. The critical standoff parameter for jet penetration is found to be highly dependent on the density ratio of the two fluids. Furthermore, the pinch-off of an interface jet produced long after the bubble dynamics stage is reckoned as the second mechanism, carrying water droplets into the oil bulk. The dependence of the bubble jetting behaviours and interface jet dynamics on the governing parameters is systematically studied via experiments and boundary integral simulations. Particularly, we quantitatively demonstrate the respective roles of surface tension and viscosity in interface jet dynamics. As for a bubble initiated at the interface, an extended Rayleigh–Plesset model is proposed that well predicts the asymmetric dynamics of the bubble, which accounts for a faster contraction of the bubble top and a downward liquid jet.

Vertically neutral collapse of a pulsating bubble occurs when the boundaries above or below the bubble balance the buoyancy effect over a pulsation. In this study, the vertically neutral collapse of a bubble near a vertical rigid wall below the free surface is investigated. The boundary integral method (BIM) is employed to model the bubble dynamics with an open-domain free surface. Moreover, this method is validated against several buoyant bubble experiments. Bubble dynamics in such conditions are associated with three dimensionless parameters: the bubble-free surface distance $\gamma _{{f}}$ , bubble–wall distance $\gamma _{{w}}$ and buoyancy parameter $\delta$ . We derive the Kelvin impulse of a spherical bubble and the algebraic relationship for vertically neutral collapse, which proves to be accurate for predicting vertically neutral collapse when the bubble is relatively far from the boundaries. Four patterns of the vertically neutral collapse of the bubble for different $\gamma _{{w}}$ and $\gamma _{{f}}$ are identified: (i) formally downward jet; (ii) annular collapse; (iii) horizontal jet; and (iv) weak jet. Despite the downward jet shape, the ‘formally downward jet’ is in the vertically neutral collapse state in terms of the profile of toroidal bubbles and the orientation of local high-pressure zones around the bubble at jet impact. A bulge with a high curvature above the bubble in the ‘annular collapse’ pattern is formed during bubble collapse under two local high-pressure zones at the left and right extremities of the bubble. The ‘horizontal jet’ pattern has the greatest potential to attack the wall, and the power laws of the moment of the jet impact, jet velocity and bubble displacement with respect to the theoretical Kelvin impulse are discussed. In particular, we quantitatively illustrate the role of the free surface on bubble migration towards the wall through the variational power-law exponents of the bubble displacement with respect to $\gamma _{{w}}$ .