Vapour bubbles produced by long-pulsed laser often have complex non-spherical shapes that reflect some characteristics of the laser beam. The transition between two commonly observed shapes, namely, a rounded pear-like shape and an elongated conical shape, is studied using a new computational model that combines compressible multiphase fluid dynamics with laser radiation and phase transition. Two laboratory experiments are simulated, in which a holmium:YAG or thulium fibre laser is used to generate bubbles of different shapes. In both cases, the predicted bubble nucleation and morphology agree reasonably well with the experimental observation. The full-field results of laser irradiance, temperature, velocity and pressure are analysed to explain bubble dynamics and energy transmission. It is found that due to the lasting energy input, the vapour bubble's dynamics is driven not only by advection, but also by the continued vaporisation at its surface. Vaporisation lasts less than $1~{\rm \mu}$s in the case of the pear-shaped bubble, compared with over $50~{\rm \mu}$s for the elongated bubble. It is thus hypothesised that the bubble's morphology is determined by competition. When the speed of advection is higher than that of vaporisation, the bubble tends to grow spherically. Otherwise, it elongates along the laser beam direction. To test this hypothesis, the two speeds are defined analytically using a model problem, then estimated for the experiments using simulation results. The results support the hypothesis. They also suggest that when the laser's power is fixed, a higher laser absorption coefficient and a narrower beam facilitate bubble elongation.

The collapse of a vapour bubble over a material surface has been widely studied over the past few decades, but a comprehensive and quantitative analysis of the cavitation dynamics and its effects on solid materials at the mesoscale (nanometre up to micrometre), which would be of particular interest in applications exploiting cavitation power, is still lacking. In this work, we adopt a diffuse interface model to describe the microbubble dynamics, and a dynamic plasticity model for the solid. The former is particularly suited to studying the rich phenomenology characterising bubble collapse at the mesoscale, which comprises transitions to supercritical conditions, emission and propagation of shock waves, generation of liquid microjets and topological transitions, whereas the latter is used to characterise the permanent plastic deformation caused by the bubble collapse, and has been augmented to consider inertial effects, to assess whether or not an interaction between elastic and plastic waves may influence the resulting deformation. Results concerning the collapse of a microbubble at different liquid overpressures and initial standoff ratios are discussed, and the elastoplastic wave propagation in the solid, together with plastic deformation, is studied for different cases, depending on elastic and plastic material parameters.

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.

We investigate the heterogeneous cavitation phenomenon in water when a spherical surface is abruptly separated from a nearby flat substrate, at a distance of approximately 10 nm. By tracking the surface separation using Newton ring positions, and capturing the bubble evolution with a high-speed camera on a microscope, we compare our experimental findings with hydrodynamic predictions at low Reynolds numbers. Upon upward movement of the spherical surface, the resulting bubble develops branched fingers through the Saffman–Taylor instability. Simultaneously, negative liquid pressures in the range $\sim$10 atm are observed. These large tension values occasionally lead to secondary nucleation events. The bubble sizes satisfy a predicted Family–Vicsek scaling law where the bubble area is proportional to the inverse bubble lifetime. The fact that creeping flow cavitation bubbles are more short-lived the larger they are separates them from bubbles that are governed by inertial dynamics.

In this paper, we numerically study the mechanism of the oscillatory flow dynamics associated with the tip vortex cavitation (TVC) over an elliptical hydrofoil section. Using our recently developed three-dimensional variational multiphase flow solver, we investigate the TVC phenomenon at Reynolds number $Re = 8.95 \times 10^5$ via dynamic subgrid-scale modelling and the homogeneous mixture theory. To begin, we examine the grid resolution requirements and introduce a length scale considering both the tip vortex strength and the core radius. This length scale is then employed to non-dimensionalize the spatial resolution in the tip vortex region, the results of which serve as a basis for estimation of the required mesh resolution in large eddy simulations of TVC. We next perform simulations to analyse the dynamical modes of tip vortex cavity oscillation at different cavitation numbers, and compare them with the semi-analytical solution. The breathing mode of cavity surface oscillation is extracted from the spatial-temporal evolution of the cavity's effective radius. The temporally averaged effective radius demonstrates that the columnar cavity experiences a growth region followed by decay as it progresses away from the tip. Further examination of the characteristics of local breathing mode oscillations in the growth and decay regions indicates the alteration of the cavity's oscillatory behaviour as it travels from the growth region to the decay region, with the oscillations within the growth region being characterized by lower frequencies. For representative cavitation numbers $\sigma \in [1.2,2.6]$, we find that pressure fluctuations exhibit a shift of the spectrum towards lower frequencies as the cavitation number decreases, similar to its influence on breathing mode oscillations. The results indicate the existence of correlations between the breathing mode oscillations and the pressure fluctuations. While the low-frequency pressure fluctuations are found to be correlated with the growth region, the breathing mode oscillations within the decay region are related to higher-frequency pressure fluctuations. The proposed mechanism can play an important role in developing mitigation strategies for TVC, which can reduce the underwater radiated noise by marine propellers.

The bubbly shock-driven partial cavitation in an axisymmetric venturi is studied with time-resolved two-dimensional X-ray densitometry. The bubbly shock waves are characterised using the vapour fraction and pressure changes across it, propagation velocity, and Mach number. The sharp changes in vapour fraction measured with X-ray densitometry, combined with high-frequency dynamic pressure measurements, reveal that the interaction of the pressure wave with the vapour cavity dictates the shedding dynamics. At the lowest cavitation number ($\sigma \sim 0.47$), the condensation shock front is the predominant shedding mechanism. However, as $\sigma$ increases ($\sigma \sim 0.78$), we observe an upstream travelling pressure discontinuity that changes into a condensation shock as it approaches the venturi throat. This coincides with the increasing strength of the bubbly shock wave as it propagates upstream, manifested by the increasing velocity of the shock front and the pressure rise across it. Consequently, the Mach number of the shock front increases and surpasses the critical value 1, favouring condensation shocks. Further, at higher $\sigma$ (${\sim }0.84\unicode{x2013}0.9$), both the re-entrant jet and pressure wave can cause cavity detachment. However, at such $\sigma$, the pressure wave likely remains subsonic. Hence cavity condensation is not favoured readily. This leads to the re-entrant jet causing the cavity detachment at higher $\sigma$. The shock front is accelerated as it propagates upstream through the variable cross-section of the venturi. This enhances its strength, favouring cavity condensation and eventual shedding. These observations explain the existence of shock fronts in an axisymmetric venturi for a large range of $\sigma$.

An additional distant wall is known to highly alter the jetting scenarios of wall-proximal bubbles. Here, we combine high-speed photography and axisymmetric volume of fluid (VoF) simulations to quantitatively describe its role in enhancing the micro-jet dynamics within the directed jet regime (Zeng et al., J. Fluid Mech., vol. 896, 2020, A28). Upon a favourable agreement on the bubble and micro-jet dynamics, both experimental and simulation results indicate that the micro-jet velocity increases dramatically as $\eta$ decreases, where $\eta =H/R_{max}$ is the distance between two walls $H$ normalized by the maximum bubble radius $R_{max}$. The mechanism is related to the collapsing flow, which is constrained by the distant wall into a reverse stagnation-point flow that builds up pressure near the bubble's top surface and accelerates it into micro-jets. We further derive an equation expressing the micro-jet velocity $U_{jet}=87.94\gamma ^{0.5}(1+(1/3)(\eta -\lambda ^{1.2})^{-2})$, where ${\gamma =d/R_{max}}$ is the stand-off distance to the proximal wall with $d$ the distance between the initial bubble centre and the wall, $\lambda =R_{y,m}/R_{max}$ with $R_{y,m}$ the distance between the top surface and the proximal wall at the bubble's maximum expansion. Viscosity has a minimal impact on the jet velocity for small $\gamma$, where the pressure buildup is predominantly influenced by geometry.

Pulmonary tuberculosis (PTB) elimination efforts must consider the global growth of the ageing population. Here we used TB surveillance data from Texas, United States (2008–2020; total n = 10656) to identify unique characteristics and outcomes in older adults (OA, ≥65 years) with PTB, compared to young adults (YA, 18–39 years) or middle-aged adults (40–64 years). We found that the proportion of OA with PTB increased from 15% in 2008 to 24% in 2020 (trend p < 0.05). Diabetes was highly prevalent in OA (32%) but not associated with adverse outcomes. Death was 13-fold higher in OA compared to YA and was 7% at the time of diagnosis which suggests diagnostic delays. However, once TB was suspected, we found no differences in culture, smear, or nucleic acid detection of mycobacteria (although less lung cavitations) in OA. During treatment, OA had less drug-resistant TB, few adverse reactions and adhered with TB treatment. We recommend training healthcare workers to ‘think TB’ in OA, for prompt treatment initiation to diminish deaths. Furthermore, OA should be added as a priority group to the latent TB treatment guidelines by the World Health Organization, to prevent TB disease in this highly vulnerable group.

The coherent dynamics of bubble clusters are of fundamental and industrial importance, and are elusive due to the complex interactions of disordered bubble oscillations. Here we introduce and demonstrate a method for decomposition of the Lagrangian time series of bubble dynamics data by combining theory and principal component analysis. The decomposition extracts coherent features of bubble oscillations based on their energy, in a way similar to proper orthogonal decomposition of Eulerian flow field data. This method is applied to a dataset of spherical clusters under harmonic excitation at different amplitudes, with various nuclei density and polydispersity parameters. Results indicate that the underlying correlated mode of oscillations is isolated in a single dominant feature in cavitating regimes, independent of the nuclei's parameters. A systematic data analysis procedure further suggests that this feature is globally controlled by the dynamic cloud interaction parameter of Maeda and Colonius (J. Fluid Mech., vol. 862, 2019, pp. 1105–1134) that quantifies the mean-field interactions, regardless of initial polydispersity or nonlinearity. The method provides a simplified and comprehensive representation of complex bubble dynamics as well as a new path to reduced-order modelling of cavitation and nucleation.

In this article, we present direct numerical simulation results for the expansion of spherical cap bubbles attached to a rigid wall due to a sudden drop in the ambient pressure. The critical pressure drop beyond which the bubble growth becomes unstable is found to match well with the predictions from classical theory of heterogeneous nucleation imposing a quasi-static bubble evolution. When the pressure drop is significantly higher than the critical value, a liquid microlayer appears between the bubble and the wall. In this regime, the interface outside the microlayer grows at an asymptotic velocity that can be predicted from the Rayleigh–Plesset equation, while the contact line evolves with another asymptotic velocity that scales with a visco-capillary velocity that obeys the Cox–Voinov law. In general, three distinctive regions can be distinguished: the region very close to the contact line where dynamics is governed by visco-capillary effects, an intermediate region controlled by inertio-viscous effects away from the contact line yet inside the viscous boundary layer, and the region outside the boundary layer dominated by inertial effects. The microlayer forms in a regime where the capillary effects are confined in a region much smaller than the viscous boundary layer thickness. In this regime, the global capillary number takes values much larger then the critical capillary number for bubble nucleation, and the microlayer height is controlled by viscous effects and not surface tension.

Cavitation inception originates from nuclei in a liquid. This paper proposes a Gibbs free energy approach that provides a smooth transition from homogeneous to heterogeneous nucleation when gas is present. The impact of gas content on nucleation is explored. It is found that the gas content stabilises nuclei, a phenomenon not present in pure liquid–vapour systems. This reduces the energy barrier over that required to nucleate a vapour bubble. Different gas saturation levels are studied. Gas content can significantly reduce the energy barrier required for nucleation, and under certain circumstances eliminate it. An analytic solution for the critical radius and activation energy is obtained that accounts for gas content. The classical Blake radius is recovered as a limiting case. The hysteresis between incipience and desinence is explained using the asymmetry observed in the critical radii. The solution is used to obtain the initial bubble radius, given a critical pressure condition in cavitation susceptibility meter experiments. The relationship between initial bubble diameter and critical pressure is described by an analytic solution that accounts for gas content. A model for the derivative of the cumulative nuclei histogram with respect to bubble diameter is proposed. An analytic expression is obtained that shows good agreement with decades worth of experimental data compiled by Khoo et al. (Exp. Fluids, vol. 61, issue 2, 2020, pp. 1–20) from ocean to water tunnels. The expression recovers the $-4$ power law that is observed experimentally.

This study gives insights into the interfacial instabilities of a ventilation cavity by injecting gas vertically into the horizontal liquid crossflow through both numerical and experimental investigations. We identified four distinct regimes of the ventilation cavity based on their topological characteristics: (I) discrete bubble, (II) continuous cavity, (III) bifurcated cavity, and (IV) bubble plume. The boundaries for these regimes are delineated within the parameter space of crossflow velocity and jet speed. A comprehensive analysis of the flow characteristics associated with each regime is presented, encompassing the phase mixing properties, the dominant frequency of pulsation, and the time-averaged profile of the cavity. This study conducted a detailed investigation of the periodic pulsation at the leading-edge interface of the cavity, also known as the ‘puffing phenomenon’. The results of local spectral analysis and dynamic mode decomposition indicate that the high-frequency instability in the near-field region exhibits the most significant growth rate. In contrast, the low-frequency mode with the largest amplitude spans a broader region from the orifice to the cavity branches. A conceptual model has been proposed to elucidate the mechanism behind the pulsation phenomenon observed along the cavity interface: the pulsation results from the alternate intrusion of the crossflow and the cavity recovery at the leading edge, being governed mainly by the periodic oscillating imbalance between the static pressure of gas near the orifice and the stagnation pressure of crossflow at the leading edge.

In this paper, we present a theoretical, experimental and numerical study of the dynamics of cavitation bubbles inside a droplet suspended in another host fluid. On the theoretical side, we provided a modified Rayleigh collapse time and natural frequency for spherical bubbles in our particular context, characterized by the density ratio between the two liquids and the bubble-to-droplet size ratio. Regarding the experimental aspect, experiments were carried out for laser-induced cavitation bubbles inside oil-in-water (O/W) or water-in-oil (W/O) droplets. Two distinct fluid-mixing mechanisms were unveiled in the two systems, respectively. In the case of O/W droplets, a liquid jet emerges around the end of the bubble collapse phase, effectively penetrating the droplet interface. We offer a detailed analysis of the criteria governing jet penetration, involving the standoff parameter and impact velocity of the bubble jet on the droplet surface. Conversely, in the scenario involving W/O droplets, the bubble traverses the droplet interior, inducing global motion and eventually leading to droplet pinch-off when the local Weber number exceeds a critical value. This phenomenon is elucidated through the equilibrium between interfacial and kinetic energies. Lastly, our boundary integral model faithfully reproduces the essential physics of the non-spherical bubble dynamics observed in the experiments. We conduct a parametric study spanning a wide parameter space to investigate bubble–droplet interactions. The insights from this study could serve as a valuable reference for practical applications in the field of ultrasonic emulsification, pharmacy, etc.

Jetting of collapsing bubbles is a key aspect in cavitation-driven fluid–solid interactions as it shapes the bubble dynamics and additionally due to its direct interaction with the wall. We study experimentally and numerically the near-wall collapse and jetting of a single bubble seeded into the stagnation flow of a wall jet, i.e. a jet that impinges perpendicular onto a solid wall. High-speed imaging shows rich and rather distinct bubble dynamics for different wall jet flow velocities and bubble-to-wall stand-off distances. The simulations use a volume-of-fluid method that allows us to numerically determine the microscopic and transient pressures and shear stresses on the wall. It is shown that a wall jet at moderate flow velocities of a few metres per second already shapes the bubble ellipsoidally inducing a planar and convergent jet flow. The distinct bubble dynamics allow us to tailor the wall interaction. In particular, the shear stresses can be increased by orders of magnitude without increasing impact pressures the same way. Interestingly, at small seeding stand-offs, the bubble during the final collapse stage can lift off the wall and migrate against the flow direction of the wall jet such that the violent collapse occurs away from the wall.

In this paper, we look for minimizers of the energy functional for isotropic compressible elasticity taking into consideration the effect of a gravitational field induced by the body itself. We consider two types of problems: the displacement problem in which the outer boundary of the body is subjected to a Dirichlet-type boundary condition, and the one with zero traction on the boundary but with an internal pressure function. For a spherically symmetric body occupying the unit ball $\mathcal {B}\in \mathbb {R}^3$, the minimization is done within the class of radially symmetric deformations. We give conditions for the existence of such minimizers, for satisfaction of the Euler–Lagrange equations, and show that for large displacements or large internal pressures, the minimizer must develop a cavity at the centre. We discuss a numerical scheme for approximating the minimizers for the displacement problem, together with some simulations that show the dependence of the cavity radius and minimum energy on the displacement and mass density of the body.

We investigate experimentally and theoretically the interactions between a cavitation bubble and a hemispherical pendant oil droplet immersed in water. In experiments, the cavitation bubble is generated by a focused laser pulse right below the pendant droplet with well-controlled bubble–wall distances and bubble–droplet size ratios. By high-speed imaging, four typical interactions are observed, namely: oil droplet rupture; water droplet entrapment; oil droplet large deformation; and oil droplet mild deformation. The bubble jetting at the end of collapse and the migration of the bubble centroid are particularly different in each bubble–droplet interaction. We propose theoretical models based on the method of images for calculating the Kelvin impulse and the anisotropy parameter which quantitatively reflects the migration of the bubble centroid at the end of the collapse. Finally, we explain that a combination of the Weber number and the anisotropy parameter determines the regimes of the bubble–droplet interactions.

The maximum size and lifetime of an acoustically nucleated cavitation bubble is inversely proportional to the driving frequency and has achieved a limit of approximately 10 MHz. Smaller cavitation bubbles that are critical to microscopic applications require shorter lifetimes that correspond to higher oscillation frequencies. Here, we demonstrate that acoustic cavitation in the 100 MHz range and beyond can be achieved through wave propagation in a solid rather than in a liquid. The cavitation bubble is nucleated at a nano-sized fracture on a glass substrate, and its expansion is driven by a leaky Rayleigh wave, while the inertial collapse is induced by a trailing shock wave. As both waves travel at different velocities, the time interval between these two events is a function of the distance to the source. In this way, we demonstrate experimentally control of the lifetime of the bubbles in a range between 6 and 80 ns, corresponding to oscillation frequencies between 13 and 166 MHz. Our results agree with finite-volume fluid–structure interaction simulations.

The three-dimensional dynamic behaviour of a tip vortex generated by a NACA $66_2$-415 hydrofoil is investigated under wetted flow and cavitating conditions using time-resolved tomographic particle image velocimetry. Two main cavitation modes are studied, namely the breathing and double-helical modes. The time-averaged flow field consists of a system of two streamwise vortices for all three conditions. The shape of the tip vortex resembles that of the cavitation mode, instead of a circular cylinder. Multi-scale vortical structures are captured in the instantaneous flow field. The surface oscillation of the cavity contributes to the growth of Kelvin–Helmholtz instability over the tip vortex, leading to the onset of hairpin vortices. Stronger spanwise interaction between the tip vortex and flow separation of the hydrofoil is produced by cavitation, further intensifying the perturbation growth. The proper orthogonal decomposition analysis gives insight into the relationship between cavity oscillation and vortex instability. Two major types of unstable modes of the tip vortex are obtained, leading to serpentine centreline displacement and elliptical deformation motion. The selection of the dominant unstable mode is associated with cavity surface oscillation. For the wetted flow condition, the displacement mode dominates the growth of vortex perturbation, while for breathing mode cavitation, the most energetic unstable mode changes into an elliptical deformation pattern, the disturbance energy of which is negligible in the wetted flow condition. Consistency is found between the peak frequency of the deformation mode and the cavity resonance frequency, indicating the contribution of cavity oscillation to the disturbance growth and breakdown of the tip vortex.

Cavitation inception in the turbulent shear layer developing behind a backward-facing step occurs at multiple points along quasi-streamwise vortices (QSVs), at a rate that increases with the Reynolds number (Re). This study investigates the evolution of the unsteady pressure field and the distribution of nuclei within and around the QSVs. The time-resolved volumetric velocity in the non-cavitating flow is measured using tomographic particle tracking, and the pressure is determined by spatial integration of the material acceleration. Analysis in Eulerian and Lagrangian reference frames reveals that the pressure is lower, and its minima last longer within the QSVs compared with the surrounding flow. The intermittent low pressure regions, whose sizes and shapes are consistent with those of the cavities, are likely to be preceded by axial vortex stretching and followed by contraction. Such phenomena have been observed before in simulations of stretched vortex elements. For the same axial straining, the pressure minima last longer with increasing Re, a trend elucidated in terms of viscous diffusion of the stretched vortex core. The impact of nuclei availability is studied under ‘natural’ and controlled seeding. Owing to differences in the saturation level of non-condensable gas, the microbubble concentration in the shear layer decreases with increasing Re, in contrast to the rate of cavitation events. Minor differences in entrainment rate into the shear layer also do not explain the substantial Re effects on cavitation inception. Hence, the Re scaling of inception appears to be dominated by trends of the pressure field.

Ventilated cavitating flows are investigated via direct numerical simulations, using a coupled level set and volume of fluid method to capture the interface between the air and water phases. A ventilated disk cavitator is used to create the cavity and is modelled by a sharp-interface immersed boundary method. The simulation data provide a comprehensive description of the two-phase flow and the air leakage and vortex shedding processes in the cavitating flow. The mean velocity of the air phase suggests the existence of three characteristic flow structures, namely the shear layer (SL), recirculating area (RA) and jet layer (JL). The turbulent kinetic energy (TKE) is concentrated in the JL in the closure region, and streamwise turbulent fluctuations dominate transverse fluctuations in both SL and JL. Budget analyses of the TKE show that the production term causes the TKE to increase in the SL due to the high velocity gradients, and decrease in the JL due to streamwise stretching effects. Air leakage and vortex shedding occur periodically in the closure region, and the one-to-one correspondence between these two processes is confirmed by the velocity and volume fluid spectra results, and the autocorrelation function of the air volume fraction. Moreover, the coherent flow structures are analysed using the spectral proper orthogonal decomposition method. We identify several fine coherent structures, including $SL_{KH}$ induced by the Kelvin–Helmholtz instability, $SL_{out}$ associated with large-scale vortex shedding, $SL_{in}$ associated with small-scale vortex shedding, and $SL_{r}$ associated with upstream turbulent convection. The present study complements previous research by providing detailed descriptions of the turbulent motions associated with the violent mixing of air and water, and the complex interactions between different characteristic structures in cavitating flows.