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.

This paper is concerned with the dynamics of large bubbles subject to various strengths of buoyancy effects, which are associated with applications for underwater explosion. The bubble is produced by electric discharge in a low-pressure tank to enhance the buoyancy effects. Experiments are carried out for a bubble in an infinite field, below a free surface and above a rigid boundary. The effects of buoyancy are reflected by the dimensionless parameter ${\it\delta}=\sqrt{{\it\rho}gR_{m}/(p_{amb}-p_{v})}$ , where $R_{m}$ , $p_{amb}$ , $p_{v}$ , ${\it\rho}$ and $g$ are the maximum bubble radius, ambient pressure, saturated vapour pressure, density of water and the acceleration of gravity respectively. A systematic study of buoyancy effects is carried out for a wide range of ${\it\delta}$ from 0.034 to 0.95. A series of new phenomena and new features is observed. The bubbles recorded are transparent, and thus we are able to display and study the jet formation, development and impact on the opposite bubble surface as well as the subsequent collapsing and rebounding of the ring bubble. Qualitative analyses are carried out for the bubble migration, jet velocity and jet initiation time, etc. for different values of ${\it\delta}$ . When a bubble oscillates below a free surface or above a rigid boundary, the Bjerknes force due to the free surface (or rigid boundary) and the buoyancy are in opposite directions. Three situations are studied for each of the two configurations: (i) the Bjerknes force being dominant, (ii) the buoyancy force being dominant and (iii) the two forces being approximately balanced. For case (iii), we further consider two subcases, where both the balanced Bjerknes and buoyancy forces are weak or strong. When the Bjerknes and buoyancy forces are approximately balanced over the pulsation, some representative bubble behaviours are observed: the bubble near free surface is found to split into two parts jetting away from each other for small ${\it\delta}$ , or involutes from both top and bottom for large ${\it\delta}$ . A bubble above a rigid wall is found to be subject to contraction from the lateral part leading to bubble splitting. New criteria are established based on experimental results for neutral collapses where there is no dominant jetting along one direction, which correlate well with the criteria of Blake et al. (J. Fluid Mech., vol. 170, 1986, pp. 479–497; J. Fluid Mech., vol. 181, 1987, pp. 197–212) but agree better with the experimental and computational results.

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.

This work focuses on using the power of a collapsing bubble in ice breaking. We experimentally validated the possibility and investigated the mechanism of ice breaking with a single collapsing bubble, where the bubble was generated by underwater electric discharge and collapsed at various distances under ice plates with different thicknesses. Characteristics of the ice fracturing, bubble jets and shock waves emitted during the collapse of the bubble were captured. The pattern of the ice fracturing is related to the ice thickness and the bubble–ice distance. Fractures develop from the top of the ice plate, i.e. the ice–air interface, and this is attributed to the tension caused by the reflection of the shock waves at the interface. Such fracturing is lessened when the thickness of the ice plate or the bubble–ice distance increases. Fractures may also form from the bottom of the ice plate upon the shock wave incidence when the bubble–ice distance is sufficiently small. The ice plate motion and its effect on the bubble behaviour were analysed. The ice plate motion results in higher jet speed and greater elongation of the bubble shape along the vertical direction. It also causes the bubble initiated close to the ice plate to split and emit multiple shock waves at the end of the collapse. The findings suggest that collapsing bubbles can be used as a brand new way of ice breaking.

High-speed water entry is a transient hydrodynamic process that is accompanied by strongly compressible flow, free surface splash, cavity evolution and other nonlinear hydrodynamic phenomena. To address these problems, a novel fluid–structure interaction (FSI) scheme based on the immersed boundary method is proposed which is suitable for strongly compressible multiphase flows. In this scheme, considering the multiphase interfaces at the immersed boundary, an improved immersed boundary method for effectively suppressing the non-physical force oscillation is proposed. Additionally, a quaternion-based six degrees of freedom motion system is used to describe rigid body motion, and the multiphase flow Eulerian finite element method is applied as the fluid solver. Using analytical solutions, experimental data and literature data, the accuracy and robustness of the FSI scheme are validated. Finally, the high-speed water entry of the slender body with different noses is investigated, and the hydrodynamic loads including the axial and normal drag forces and the bending moment are extensively discussed. The hydrodynamic load and motion trajectory are determined by the nose configuration. The tail slamming phenomenon is the primary focus, and it is revealed that its formation is primarily related to the pitch moment formed at the stage of crossing the free surface. Tail slamming also causes violent impact loads, especially bending moments, which may cause slender projectiles to break off. Finally, to combine the features of the flat and hemispherical noses, the water entry of the projectile with a truncated hemispherical nose is simulated 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.