1. Introduction
The collapse of a wall-proximal cavitation bubble results in a micro-jet that pierces the bubble, morphs it into a torus, and the jet flow impacts onto the boundary, and then spreads along the boundary (Plesset & Chapman Reference Plesset and Chapman1971; Lauterborn & Bolle Reference Lauterborn and Bolle1975; Zhang et al. Reference Zhang, Li, Cui, Li and Liu2023). This rich fluid mechanics occurs in a brief period of time, including ring vortices formation and spreading along or away from the wall (Reuter et al. Reference Reuter, Gonzalez-Avila, Mettin and Ohl2017), creation of large wall-normal and tangential shear stresses through the jet impact and its spreading (Dijkink & Ohl Reference Dijkink and Ohl2008; Zeng, An & Ohl Reference Zeng, An and Ohl2022). Depending on the application, these flows may be unwanted due to the destructive effect or exploited for their benefit, e.g. for ultrasonic cleaning, drug delivery, histotripsy, lithotripsy and cavitation peening.
Recent experiment (Reuter & Ohl Reference Reuter and Ohl2021) and simulation (Lechner et al. Reference Lechner, Lauterborn, Koch and Mettin2020) demonstrated that a micro-jet could be either a needle-type or a regular jet depending on the stand-off distance $\gamma$. Under atmospheric pressure, simulations and experiments find a needle jet velocity of the order $1000\ {\rm m}\ {\rm s}^{-1}$, yet it is extremely thin and develops only for $\gamma < 0.2$. In contrast, the regular jet has a speed of the order of $80\ {\rm m}\ {\rm s}^{-1}$ that increases slowly with $\gamma$. The micro-jet shape and velocity are dependent on the particular flow conditions, where acoustic (Rosselló et al. Reference Rosselló, Lauterborn, Koch, Wilken, Kurz and Mettin2018) and electric fields (Taleghani et al. Reference Taleghani, Khodadadi, Maddahian and Mokhtari-Dizaji2023), bubble interaction (Han et al. Reference Han, Köhler, Jungnickel, Mettin, Lauterborn and Vogel2015), and the presence of more than one boundary (Zeng, Gonzalez-Avila & Ohl Reference Zeng, Gonzalez-Avila and Ohl2020) have been reported to affect the velocity.
A recent study highlighted that an additional distant wall leads to more intricate bubble dynamics, unveiling three distinct jetting phenomena: transferred, double and directed jets (Zeng et al. Reference Zeng, Gonzalez-Avila and Ohl2020). The conditions for the formation of different jets were studied in detail, combining high-speed photography and Computational Fluid Dynamics (CFD) simulations, taking into account the stand-off distance $\gamma$ from the proximal wall and the normalized gap height $\eta =H/R_{max}$ as crucial parameters. This finding opens avenues for manipulating jets, as the direction and strength of the jet can be tuned for specific applications based on the identified parameters.
Here we focus on the directed jet regime (Zeng et al. Reference Zeng, Gonzalez-Avila and Ohl2020), where the bubble neither splits (Ogasawara et al. Reference Ogasawara, Tsubota, Seki, Shigaki and Takahira2015; Seki et al. Reference Seki, Fujita, Ogasawara and Takahira2015) nor translates (Gonzalez-Avila et al. Reference Gonzalez-Avila, van Blokland, Zeng and Ohl2020) to the distant wall. The aim is to quantitatively elucidate the micro-jet enhancement by the distant wall. We combine high-speed photography at up to 500 000 frames per second (f.p.s.) and compressible volume of fluid (VoF) simulations with high spatio-temporal resolution to investigate the details of the flow and derive an analytic expression for the jet velocity.
2. Experiments and simulation
Figures 1(a) and 1(b) sketch the geometry and introduce the coordinates and parameters of the problem. The experimental set-up is similar to Zeng et al. (Reference Zeng, Gonzalez-Avila and Ohl2020). In the experiment, a laser pulse from a $Q$-switched Nd: YAG laser (Litron Nano S, 532 nm wavelength, 6 ns pulse duration) is focused into the deionized water at a distance of $d$ to a no-slip proximal wall to nucleate the cavitation bubble. The distant wall is located at a distance of $H$ from the proximal wall. Both walls are made of glass. A high-speed camera (Shimadzu Hypervision HPV-2) connected to a macro lens (Nikor) and an adapter ring (Nikon, model PB-6) records the bubble dynamics. A weakly diffused light pulse from a strobe light (Sunpak 3075G, 4.4 ms flash duration) illuminates the scene. The camera operates at a frame rate of 500 000 f.p.s. and a resolution of $18\ \mathrm {\mu }\mathrm {m}$ per pixel. A delay generator (Quantum, 9520 series) controls the synchronization of the camera, the pulsed laser and the strobe. Experiments are conducted under atmospheric pressure, i.e. $p_{\infty }=101\,325$ Pa.
The bubble is created with a focused laser pulse that induces a dielectric breakdown in the liquid. Once formed, the laser-induced cavitation bubble expands and collapses violently, and continues to oscillate a few times. Yet most of the potential energy of the bubble is transferred into the liquid during the first cycle, i.e. when the micro-jet is formed. During this period, the cavitation bubble content can be treated as a gas bubble with a polytropic equation of state (EOS) and the liquid is modelled as a Tait liquid to account for compressibility (Koch et al. Reference Koch, Lechner, Reuter, Köhler, Mettin and Lauterborn2016).
We define a volume fraction of the liquid phase $\alpha$, whose evolution is solved by the transport equation:
where $t$ is time and $\boldsymbol {u}$ is the flow velocity. The compressibilities $\psi ={\rm D}\rho /{\rm D}p$ of both phases are modelled based on their EOS with fluid density $\rho$ and the flow pressure $p$. Subscripts $g$ and $l$ denote the gas and liquid phases, respectively. Here $\boldsymbol {U_r}$ is an artificial compressible term obtained from the velocity difference between two phases, which enables tuning of the sharpness of the interface (Rusche Reference Rusche2003).
Here, the heat and mass transfer are neglected (Koch et al. Reference Koch, Lechner, Reuter, Köhler, Mettin and Lauterborn2016). Thus, the flow is governed by the equations of continuity and momentum:
where $\boldsymbol {\tau }=\mu [\boldsymbol {\nabla }\boldsymbol {u}+\boldsymbol {\nabla } \boldsymbol {u}^T- \frac {2}{3}(\boldsymbol {\nabla }\boldsymbol {\cdot }\boldsymbol {u}){\boldsymbol {I}}]$ is the viscous stress tensor calculated based on the identity tensor ${\boldsymbol {I}}$ and the dynamic viscosity $\mu$. The surface tension term $\boldsymbol {f}_{\!\boldsymbol {\delta }}$ is accounted for using the continuous-surface-force (CSF) method (Brackbill, Kothe & Zemach Reference Brackbill, Kothe and Zemach1992).
Due to the axisymmetric nature of the problem, the domain is simplified as a wedge. To achieve a well-resolved flow and meanwhile keep the computational efficiency, we used a locally and recursively refined mesh and kept the cell size to $1.5\ \mathrm {\mu }{\rm m}$ in the area where the bubble is located. To resolve the boundary-layer flow, the near-wall region is further refined to a cell size of 50 nm. Both walls are modelled as no-slip walls. The ambient pressure $p_\infty = 101\,325\ {\rm Pa}$ is applied at the exit of the gap. See Zeng et al. (Reference Zeng, Gonzalez-Avila, Ten Voorde and Ohl2018, Reference Zeng, Gonzalez-Avila and Ohl2020); Reese et al. (Reference Reese, Schadel, Reuter and Ohl2022) for further details on the numerical implementation.
3. Bubble dynamics and micro-jet formation
In the simulation, the explosive growth of a laser-induced bubble is achieved by initiating a high-pressure gas pocket with an initial radius $R_{0}=50\ \mathrm {\mu }{\rm m}$ at $t=0$. The initial pressure of the bubble $p_0$ is the only free parameter and is adjusted to match the experimental period of growth and collapse of the bubble. Figure 1(c) compares the experimental and numerical results for a volume-equivalent spherical bubble radius $R$ for two selected $\eta$. This demonstrates favourable agreement of the volume dynamics in size and timing. The results also indicate that the bubble acquires a longer oscillation period as the distant wall moves closer to the proximal wall, i.e. for a more confined bubble.
Figure 2 compares the experimental and numerical bubble dynamics. The shadow images are experimental recordings, while the red contours superimposed are computed bubble interfaces at the same time instance. As well known, we find for small stand-off distance $\gamma =0.47$ that the lower bubble wall flattens considerably and the bubble obtains a nearly hemispherical shape that later changes into a cone before collapsing as a toroidal bubble with a micro-jet impacting on the proximal wall. The bubble dynamics is rather weakly affected during the expansion and early shrinkage by the gap height $\eta$, see red and blue contours. Yet during collapse, the bubble develops a narrower apex during the collapse when the distance reduces to $\eta =2.21$ (see $t=140\ \mathrm {\mu }{\rm s}$). This particular shape helps to build up a larger liquid pressure that accelerates the top bubble interface into a faster micro-jet. The computed velocity of the micro-jet $U_{jet}=109.2\ {\rm m}\ {\rm s}^{-1}$ agrees well with the experimental measurement $U_{jet}=95.6\pm 24.3\ {\rm m}\ {\rm s}^{-1}$. The jet velocity increases significantly from $U_{jet}=71.9\ {\rm m}\ {\rm s}^{-1}$ for $\eta =2.95$ to $U_{jet}=109.2\ {\rm m}\ {\rm s}^{-1}$ for $\eta =2.23$, which provides us with a robust method to enhance the micro-jet velocity.
Next, we concentrate on the jet formation process to develop a quantitative description of the micro-jet dynamics. Figure 3 compares the jet formation process for three selected gap heights of $H=5$, 1.4 and 1.15 mm, respectively. Although the bubbles acquire similar outlines at their maximum expansion (see figure 2), the bubble's top surface behaves very differently during the collapse for various $\eta$. The main reason is that the pressure distribution varies with $\eta$ at the maximum expansion. Since the flow between the bubble's apex and distant wall resembles a stagnation-point flow during the bubble's expansion, we analyse this flow with a potential flow coupled to a boundary-layer flow and to obtain the pressure difference between the bubble's top surface and the distant wall: $\Delta p=p_{wall}-p_b\sim k(H-R_{y,m})^{2}+\mu f(y)$ (Schlichting & Gersten Reference Schlichting and Gersten2016), where $p_{wall}$ and $p_{b}$ are the pressures at the distant wall and bubble's top surface, $f(y)$ is a function accounting for the viscous boundary layer, $R_{y,m}$ the distance between the top surface and proximal wall that only slightly changes with $\eta$, $H-R_{y,m}$ is the distance from the bubble's top surface to the distant wall, and $k$ is the coefficient. This pressure difference induces a pressure gradient that drives the bubble's top surface to move downwards once the collapse starts, which increases with $H$. Let's compare two cases $\eta =8.42$ and $\eta =2.00$. Figure 4 shows that the top surface of a bubble with a large distance $\eta =8.42$ moves earlier and reaches the bubble's bottom surface sooner, yet it has a slower final velocity of $64.3\ {\rm m}\ {\rm s}^{-1}$. Here $y_{top}$ and $U_{top}$ are the vertical postion and the vertical velocity of the bubble's top surface, respectively. For a smaller distance $\eta =2.00$, the top surface moves slower at the beginning but is accelerated dramatically at approximately $t=151\ \mathrm {\mu }{\rm s}$ and eventually develops a fast micro-jet with a speed of $280\ {\rm m}\ {\rm s}^{-1}$.
One may ask, why does the bubble with a small $\eta$ move slower initially, but at the end develop a faster micro-jet? The mechanism is as follows. During the collapse, since the top surface is exposed to a smaller pressure gradient, the bubble shrinks faster in the annular direction, subsequently, this annular flow shapes the top surface into a conical shape, whose curvature $1/R_s$ increases with decreasing $\eta$ (see figure 3). The annular flow meets on the axis of symmetry and is redirected downwards, which turns into a reverse stagnation-point flow that flows towards the conical shape (see sketch in figure 1(b) and local flow field in figure 3). As a result, hydrodynamic pressure builds up above the top surface. Once this pressure exceeds a threshold value, the flow starts to push the conical shape downwards. The convex interface then goes through a nearly flat shape and eventually develops as a jet at the axis. A straightforward comparison in figure 3 indicates that the bubble with a higher curvature acquires a larger pressure and therefore develops a faster micro-jet.
4. Micro-jet velocity
Next we determine the jet velocity $U_{jet}$ to quantify the enhancement causing by the distant wall. Figure 5(a) compares $U_{jet}$ as a function of $\eta$ for 13 selected values of $\gamma$ between 0.46 and 1.45. Again, we see favourable agreement of $U_{jet}$ between experimental measurements and computed results for the large range of $\eta$ studied. Our results show a similar trend for all $\gamma$ that $U_{jet}$ decreases dramatically with increasing $\eta$ and eventually reaches a constant value. To determine an expression for the jet velocity, we start with an approximate acceleration on the bubble's top surface by simplifying the Navier–Stokes equation:
where $\delta$ is the distance between the bubble's top surface and the high pressure $p$ location. Note that once the jet is formed, the acceleration drops dramatically as $\delta$ increases. Therefore, the acceleration is short-lived (see figures 3 and 4), and we expect the highest acceleration to occur immediately before the apex develops a concave shape. Below, we give an estimate for the driven pressure $p_m$ at this instance of time. Assuming the short-lived acceleration is increasing linearly, $U_{jet}$ can be approximated as
considering that the pressure $p_m \geqslant p_b$ and $\Delta t$ is the duration of the pressure pulse. It is convenient to decompose the pressure into two parts $p_m \propto p_\infty +p_h$, where $p_\infty$ and $p_h$ are the pressure induced by the proximal and distant walls, respectively. For a certain $\gamma$, $p_\infty$ generates a jet velocity $U_{\infty }$ for $\eta \to \infty$ which is well studied in the literature; here we obtain an approximation of $U_{\infty }=87.94\gamma ^{0.5}$ for an ambient pressure of 1 atm. From figure 3 we observe that a relatively strong reverse stagnation-point flow impacting bubble's apex is formed due to the distant wall. Before jet formation the bubble's apex is hardly moved along the $y$ direction, therefore applying potential flow theory $p_m$ can be simplified to a hydrodynamic pressure caused by the flow with $p_h\sim \rho U_s^2$, where $U_s$ is the additional flow velocity due to the distant wall. Next, we estimate $U_s$. The maximum boundary-layer thickness of the distant wall can be approximated to be of the order of $\sim \sqrt {\mu t_c/\rho }\approx 14\ \mathrm {\mu }{\rm m}$ with the collapse time $t_{c}$, therefore it is reasonable to apply the potential flow approximation where the velocity potential $\phi$ satisfies the Laplace equation $\nabla ^{2} \phi =0$. The solution for an axisymmetric stagnation-point flow is $U_r=ar$ and $U_y=-k(H-y)$, where a and k are the coefficients to be determined. For the cylindrical flow regime between the bubble's top surface and the distant wall, the flow is mostly downwards, therefore we can simplify the kinetic energy of this cylindrical flow as
where $R_{y}$ is the vertical position of the bubble's top surface (see figure 1b), $R_s$ is the radius of the bubble's apex, $R_t$ the equivalent radius at time $t$, $f$ is the vertical portion of the kinetic energy $(4/3){\rm \pi} \Delta p(R_{max}^3-R_t^3)$ with $\Delta p=p_{\infty }-p_b\approx p_{\infty }$ (Lamb Reference Lamb1924); the right part of (4.3) also assumes the energy is proportional to area in the radial direction before the jet is formed, therefore a factor ${R_s^2}/{R_r^2}$ is applied, where $R_r$ is the bubble's radial extension. At $r=0$, $U_r=0$ and $U_y=-k(H-R_y)$, thus together with (4.3) we find
Combining (4.2) and (4.4) and accounting for the time of acceleration $\Delta t\sim R_s (\rho /\Delta p)^{1/2}$, the jet velocity becomes
Considering $U_\infty \sim \sqrt {\Delta p/\rho }$, and the fact that $R_s$ and $\delta$, $R_{max}$ and $R_r$ are of the same order of magnitude, respectively, we expect the jet velocity can be expressed as
with two proportionality factors $a$ and $b$ that we expect $a\sim O(0.1)$ and $b\sim O(1)$. Here, to relate $R_y$ with the more easily accessible parameter $\lambda =R_{y,m}/R_{max}$, we introduce the proportionality factor $b$ that accounts for the changing between $R_{y,m}$ and $R_y$ and the effect of local curvature of the apex. We obtain $(a, b) = (1/3, 1.2)$ for all data (see figure 5b), and therefore the jet velocity under atmospheric pressure can be predicted as
5. Viscous effect
Viscosity is known to affect the boundary layer and induce viscous drag during the expanding and collapsing of a cavitation bubble. As a result, the jet formation and velocity of a single wall-proximal bubble are highly influenced, too. Yet with an additional distant wall, we observe distinct fluid mechanics. Figure 6 compares the jet velocities for three selected stand-off distances $\gamma =0.55$, 0.97 and 1.37, which are plotted with different colours. For each $\gamma$, we also compare the results of five viscosities indicated by different symbols. Overall, we observe a rapid reduction of the jet velocity as $\eta$ increases, which indicates that the distant wall still holds a significant impact on $U_{jet}$ even for liquids with high viscosities.
Looking closer into the data, we find that for small stand-off distance $\gamma =0.55$, the data for all viscosities collapse onto the theoretical curve $U_{jet}=87.94\gamma ^{0.5}(1+(1/3)(\eta -(1+\gamma )^{1.2})^{-2})$. Here $\lambda \approx 1+\gamma$ (Reuter, Zeng & Ohl Reference Reuter, Zeng and Ohl2022). The result suggests that the jet formation is still dominated by the geometry, i.e. by the proximal and distant walls. Within this particular regime, a thin liquid film occurs between the bubble's lower surface and the proximal wall $h$, which ensures the bubble's lower wall is trapped. Except for the lower interface, the viscous drag that acts on all collapsing parts of the bubble's surface is proportional to $\mu$; as a result the relative motion remains, and the distant wall holds the same role in building up the high pressure to drive the bubble's apex, thus (4.7) remains appropriate.
When $\gamma$ increases to 0.97, our result shows that for small $\eta <3$ the jet velocities still conform well to (4.7), implying the dominance of the geometry remains. Yet for large $\eta >3$ the results display a more scattered distribution and deviate from the theoretical prediction. This is due to the increase in the liquid-film thickness $h$ between the bubble's lower surface and the proximal wall (Zeng et al. Reference Zeng, An and Ohl2022), resulting in a greater motion of the lower-bubble interface during the collapse. Here, a thicker viscous boundary layer on the proximal wall results in a more complex vertical motion of the bubble. This may explain the scatter of the data for $U_{jet}$ below $150\ {\rm m}\ {\rm s}^{-1}$.
For a large stand-off distance $\gamma =1.37$, viscosity not only slows down the bubble's motion but also helps to form complex structures. Unexpectedly large curvatures at the bubble apex are known to be formed due to a combined effect of viscous drag and the collapsing flow, which results in the occurrence of higher velocity jets in a more viscous liquid (Zeng et al. Reference Zeng, An and Ohl2022). Here we measure higher velocities for liquids with viscosities $\mu =0.02\ \mathrm {Pa}\ \mathrm {s}$, $\mu =0.04\ \mathrm {Pa}\ \mathrm {s}$ and $\mu =0.06\ \mathrm {Pa}\ \mathrm {s}$ as compared with water. As viscosity increases to $\mu =0.08\ \mathrm {Pa}\ \mathrm {s}$, the large viscous drag forces the bubble surface to a convex and smooth shape, thereby decreasing the jet velocity.
6. Conclusion
In summary, a distant wall is found to enhance the jetting from a collapsing wall-proximal cavitation bubble. Combining high-speed photography and compressible VoF simulations, we find that the distant wall influences the micro-jet by converting the collapsing flow into a stagnation-point flow that drives the bubble's top surface. This flow, which changes dramatically for small $\eta$ but becomes constant for relatively large $\eta$, builds up a high pressure that drives the bubble's top surface to collapse into micro-jet. We derive a formula for the jet velocity of $U_{jet}=87.94\gamma ^{0.5}(1+(1/3)(\eta -\lambda ^{1.2})^{-2})$, which correctly predicts our experimental and simulation results. Although viscosity affects the bubble dynamics in a rather complex way, the jet velocity for small $\gamma$ is found to still follow the above-derived equation, where geometric effects dominate the collapsing flow.
Acknowledgements
The anonymous referees are acknowledged for their very helpful suggestions to improve the manuscript.
Funding
This research was partly funded by the National Natural Science Foundation of China under grant nos. 51925904 and 12302313, the Fundamental Research Funds for the Central Universities of Ministry of Education of China under grant no. 3072023CFJ0105, the Finance Science and Technology Project of Hainan Province under grant no. ZDKJ2021020, the Australian Research Council Discovery Project under grant no. DP230100556 and the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under grant no. DE 3156/6-1.
Declaration of interests
The authors report no conflict of interest.