We study the existence and regularity of minimizers of the neo-Hookean energy in the closure of classes of deformations without cavitation. The exclusion of cavitation is imposed in the form of the divergence identities, which is equivalent to the well-known condition (INV) with $\operatorname{Det} = \operatorname{det}$. We show that the neo-Hookean energy admits minimizers in classes of maps that are one-to-one a.e. with positive Jacobians, provided that these maps are the weak limits of sequences of maps that satisfy the divergence identities. In particular, these classes include the weak closure of diffeomorphisms and the weak closure of homeomorphisms satisfying Lusin’s condition N. Moreover, if the minimizers satisfy condition (INV), then their inverses have Sobolev regularity. This extends a recent result by Doležalová, Hencl, and Molchanova by showing that the minimizers they obtained enjoy extra regularity properties and that the existence of minimizers can still be obtained even when their coercivity assumption is relaxed.

Understanding the mechanism of hydrodynamic cloud cavitation is crucial to reducing noise, vibration and wear. Recent studies have clarified the physics of two distinct formation mechanisms of cloud cavitation. Ganesh et al. (J. Fluid Mech., vol. 802, 2016, pp. 37–78) identified the propagation of bubbly shockwaves as a cloud detachment mechanism. Pelz et al. (J. Fluid Mech., vol. 817, 2017, pp. 439–454) explained the influence of Reynolds number and cavitation number on asymptotic growth of the cavity sheet and its periodic shedding caused by re-entrant flow. In this paper the two mechanisms are set in relation to each other. For this, we show firstly that the transition from re-entrant flow to shockwave-driven cloud cavitation is given by a kinematic condition, namely the asymptotic sheet length equates to the chord length, $\hat {a}=L$. For $\hat {a}>L$ shockwave-driven cloud cavitation dominates. For $\hat {a}< L$ re-entrant flow-driven cloud cavitation dominates. As the cavitation number decreases, the closure region of the cavity sheet reaches the trailing edge of the hydrofoil, identifying the trailing edge as a trigger for condensation shockwaves, particularly as re-entrant flow-driven cavitation diminishes. Since the sheet length is an implicit function of the cavitation number, the kinematic condition $\hat {a}/L=1$ results in a critical cavitation number ${\sigma _\mathrm {II,III}}$ that is calculated analytically and validated by experiments. Secondly, we derive the relationship between the Strouhal number and the asymptotic sheet length for re-entrant flow-driven cloud cavitation. The model presented here is thoroughly validated by experiments.

Despite the extensive research on bubble collapse near rigid walls, the bubble collapse dynamics in the presence of shear flow near a rigid wall is poorly understood. We conduct direct simulations of the Navier–Stokes equations to explore the bubble dynamics and pressures during bubble collapse near a rigid, flat wall under linear shear flow conditions. We examine the dependence of the bubble collapse morphology and wall pressures on the initial bubble location and shear rate. We find that shear distorts the bubble, generating two re-entrant jets – one developing from the side opposite to the mean flow and the other from the far end toward the wall. Upon impact of the jet on the opposite side of the bubble, water-hammer shocks are produced, which propagate outward and interact with the convoluted bubble shape. The shock stretches the bubble towards the wall, resulting in a closer impact location for the jet originating from the far end compared with the case with no shear flow. The water-hammer pressure location can be approximated as the theoretical distance travelled by a particle initialised at the bubble centre with the corresponding constant shear flow velocity. The maximum wall pressures can thus be predicted by considering the distance between the far jet impingement location and the wall along the wall-normal direction. As the shear rate is increased, the maximum wall pressure increases, although only marginally. We determine the critical initial stand-off distance from the wall at which the bubble morphology is shear dominated, i.e. characterised by converging re-entrant jets.

This study investigates experimentally the pressure fluctuations of liquids in a column under short-time acceleration. It demonstrates that the Strouhal number $St=L/(c\,\Delta t)$, where $L$, $c$ and $\Delta t$ are the liquid column length, speed of sound, and acceleration duration, respectively, provides a measure of the pressure fluctuations for intermediate $St$ values. On the one hand, the incompressible fluid theory implies that the magnitude of the averaged pressure fluctuation $\bar {P}$ becomes negligible for $St\ll 1$. On the other hand, the water hammer theory predicts that the pressure tends to $\rho cu_0$ (where $u_0$ is the change in the liquid velocity) for $St\geq O(1)$. For intermediate $St$ values, there is no consensus on the value of $\bar {P}$. In our experiments, $L$, $c$ and $\Delta t$ are varied so that $0.02 \leq St \leq 2.2$. The results suggest that the incompressible fluid theory holds only up to $St\sim 0.2$, and that $St$ governs the pressure fluctuations under different experimental conditions for higher $St$ values. The data relating to a hydrogel also tend to collapse to a unified trend. The inception of cavitation in the liquid starts at $St\sim 0.2$ for various $\Delta t$, indicating that the liquid pressure goes lower than the liquid vapour pressure. To understand this mechanism, we employ a one-dimensional wave propagation model with a pressure wavefront of finite thickness that scales with $\Delta t$. The model provides a reasonable description of the experimental results as a function of $St$.

The collapse of an initially spherical cavitation bubble near a free surface leads to the formation of two jets: a downward jet into the liquid, and an upward jet penetrating the free surface. In this study, we examine the surprising interaction of a bubble trapped in a stable cavitating vortex ring approaching a free surface. As a result, a single fast and tall liquid jet forms. We find that this jet is observed only above critical Froude numbers ($Fr$) and Weber numbers ($We$) when ${Fr}^2 (1.6-2.73/{We}) > 1$, illustrating the importance of inertia, gravity and surface tension in accelerating this novel jet and thereby reaching heights several hundred times the radius of the vortex ring. Our experimental results are supported by numerical simulations, revealing that the underlying mechanism driving the vortex ring acceleration is the disruption of the equilibrium of high-pressure regions at the front and rear of the vortex ring caused by the free surface. Quantitative analysis based on the energy relationships elucidates that the velocity ratio between the maximum velocity of the free-surface jet and the translational velocity of the vortex ring is relatively stable yet is attenuated by surface tension when the jet is mild.

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact, partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

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.