Hostname: page-component-586b7cd67f-2brh9 Total loading time: 0 Render date: 2024-11-22T11:37:18.493Z Has data issue: false hasContentIssue false

Jetting enhancement from wall-proximal cavitation bubbles by a distant wall

Published online by Cambridge University Press:  16 May 2024

Qingyun Zeng
Affiliation:
College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, PR China Nanhai Institute of Harbin Engineering University, Sanya 572024, PR China
A-Man Zhang*
Affiliation:
College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, PR China Nanhai Institute of Harbin Engineering University, Sanya 572024, PR China
Beng Hau Tan
Affiliation:
Shenzhen Institute for Advanced Study, University of Electronic Science and Technology of China, Shenzhen 518109, PR China
Hongjie An
Affiliation:
Queensland Micro and Nanotechnology Centre, Griffith University, Nathan, QLD 4111, Australia
Claus-Dieter Ohl*
Affiliation:
Department for Soft Matter, Otto-von-Guericke-University Magdeburg, Magdeburg 39016, Germany
*
Email addresses for correspondence: [email protected], [email protected]
Email addresses for correspondence: [email protected], [email protected]

Abstract

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.

Type
JFM Rapids
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2024. Published by Cambridge University Press.

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.

Figure 1. (a) Schematic diagram of the problem and parameters. The dashed circle indicates the nucleated bubble. The solid hat-shaped curve represents a bubble at its maximum expansion with a spherical equivalent radius $R_{max}$. The distant wall is located at a distance $H$ from the proximal wall. The liquid film between the bubble's lower wall and the boundary has a thickness of $h$. (b) Schematic diagrams at the moments prior to (solid line) and immediately after (dashed line) the formation of the micro-jet. Here $R_s$ and $R_r$ indicate the radius of the bubble's apex and the bubble's radial extension, respectively. Curved arrows sketch the reverse stagnation-point flow during the collapse. (c) Comparison of experiments and simulation on bubble radii as a function of time for two selected cases with different $\eta =2.23$ and $\eta =2.99$ but similar stand-off distance $\gamma \approx 0.48$; areas coded by light colour indicate regions of experimental uncertainties. Black dashed line is the numerical bubble dynamics for the unconfined case, $\eta \rightarrow \infty$.

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:

(2.1)\begin{equation} \frac{\partial \alpha}{\partial t}+\boldsymbol{\nabla}\boldsymbol{\cdot}(\alpha \boldsymbol{u})+\boldsymbol{\nabla}\boldsymbol{\cdot}(\alpha(1-\alpha) \boldsymbol{U_r})=\alpha(1-\alpha)\left(\frac{\psi_g}{\rho_g} -\frac{\psi_l}{\rho_l} \right) \frac{{\rm D}p}{{\rm D}t}+\alpha\boldsymbol{\nabla}\boldsymbol{\cdot} \boldsymbol{u} , \end{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:

(2.2)$$\begin{gather} \frac{\partial \rho}{\partial t}+\boldsymbol{\nabla}\boldsymbol{\cdot}(\rho \boldsymbol{u})=0, \end{gather}$$
(2.3)$$\begin{gather}\frac{\partial \rho \boldsymbol{u}}{\partial t}+\boldsymbol{\nabla}\boldsymbol{\cdot}(\rho \boldsymbol{uu})={-}\boldsymbol{\nabla} p +\boldsymbol{\nabla}\boldsymbol{\cdot}\boldsymbol{\tau} + \boldsymbol{f_\delta} , \end{gather}$$

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.

Figure 2. Comparison between experiments (background images) and simulation (red contours). Experiment, $R_{max}=566\pm 18\ \mathrm {\mu }{\rm m}$, $d=273\pm 18\ \mathrm {\mu }{\rm m}$, $H=1260\pm 18\ \mathrm {\mu }{\rm m}$, $U_{jet}=95.6\pm 34.1\ {\rm m}\ {\rm s}^{-1}$, $\gamma \approx 0.48$, $\eta \approx 2.23$; simulation, $R_{max}=570.8\ \mathrm {\mu }{\rm m}$, $d=270\ \mathrm {\mu }{\rm m}$, $H=1260\ \mathrm {\mu }{\rm m}$, $U_{jet}=109.2\ {\rm m}\ {\rm s}^{-1}$, $\gamma =0.47$, $\eta =2.21$, $p_0=946$ bar. The blue contours are computed bubble dynamics for a larger $\eta =2.95$ but the same $\gamma =0.47$, $U_{jet}=71.9\ {\rm m}\ {\rm s}^{-1}$. Numbers in the frames are time in microseconds.

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}$.

Figure 3. Jet formation for three selected cases with a constant $\gamma =0.55$ but various $\eta$. (a) $H=5$ mm, $\eta =8.42$; (b) $H=1.4$ mm, $\eta =2.42$; (c) $H=1.1$5 mm, $\eta =2.00$. All images focus on a region with a height of 1 mm and width of 0.5 mm. Red lines represent the bubble interfaces. The pressure and velocity magnitude are colour-coded on the left and right parts of each panel, respectively. Arrows indicate only the direction of the flow. Note the range of pressure and velocity for the colours are different for the three cases.

Figure 4. Comparison of the motion of the bubble's top surface for six selected $\eta$ but constant $\gamma \approx 0.55$. (a) Position of top surfaces vs time. (b) Velocities of the top surfaces.

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:

(4.1)\begin{equation} \frac{{\rm D}U}{{\rm D}t}=\frac{1}{\rho}\frac{p-p_b}{\delta} , \end{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

(4.2)\begin{equation} U_{jet}\approx \frac{p_m\Delta t}{2\rho\delta}, \end{equation}

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

(4.3)\begin{equation} \int_{R_y}^{H} \rho{\rm \pi} R_s^2({-}k(H-y)) \,{{\rm d}y} \approx \frac{4}{3}{\rm \pi}\Delta p(R_{max}^3-R_t^3)\frac{R_s^2}{R_r^2} f, \end{equation}

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

(4.4)\begin{equation} p_h\sim \rho U_s^2\approx \rho k^2(H-R_y)^2\sim \frac{4}{9}{\rm \pi}\Delta p\left(\frac{R_{max}^3-R_t^3}{R_r^2}\right)^2(H-R_y)^{{-}2}. \end{equation}

Figure 5. (a) Micro-jet velocities as a function of $\eta$. Solid lines are computed results for 13 selected values of $\gamma$ between 0.46 and 1.45; black filled squares, diamonds and circles represent the experimental jet velocities for $\gamma \approx 0.47$, 1.11 and 1.43, respectively. A longitude scale for $\eta$ is used to show better comparisons between different $\gamma$. (b) All data collapse onto an orange dashed line $U_{jet}^{*}=1+1/3(\eta -\lambda ^{1.2})^{-2}$, $U_{jet}^{*}=U_{jet}/U_{\infty }$ with $U_{\infty }$ the jet velocity for $\eta \to \infty$ for a certain $\gamma$. Opacity from dark to light indicates increasing $\eta$. The inset depicts $U_{\infty }$ as a function of $\gamma$; the orange solid line is $U_{\infty }=87.94\gamma ^{0.5}$ for an ambient pressure of 1 atm.

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

(4.5)\begin{equation} U_{jet}\sim \frac{2}{9} \left(\frac{H-R_y}{R_{max}}\right)^{{-}2} \frac{R_s}{\delta} \left(\frac{R_{max}^3-R_t^3}{R_{max}R_r^2}\right)^2 \sqrt{\frac{\Delta p}{\rho}} +U_{\infty}. \end{equation}

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

(4.6)\begin{equation} U_{jet}^*=U_{jet}/U_{\infty}=1+a(\eta-\lambda^b)^{{-}2}, \end{equation}

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

(4.7)\begin{equation} U_{jet}=87.94\gamma^{0.5}(1+(1/3)(\eta-\lambda^{1.2})^{{-}2}). \end{equation}

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.

Figure 6. Comparison of $U_{jet}$ as a function of $\eta$ for bubbles with three selected $\gamma$ and liquids of five different viscosities. Dashed lines are calculated using (4.7).

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.

References

Brackbill, J.U., Kothe, D.B. & Zemach, C. 1992 A continuum method for modeling surface tension. J. Comput. Phys. 100 (2), 335354.CrossRefGoogle Scholar
Dijkink, R. & Ohl, C.-D. 2008 Measurement of cavitation induced wall shear stress. Appl. Phys. Lett. 93 (25), 254107.CrossRefGoogle Scholar
Gonzalez-Avila, S.R., van Blokland, A.C., Zeng, Q. & Ohl, C.-D. 2020 Jetting and shear stress enhancement from cavitation bubbles collapsing in a narrow gap. J. Fluid Mech. 884, A23.CrossRefGoogle Scholar
Han, B., Köhler, K., Jungnickel, K., Mettin, R., Lauterborn, W. & Vogel, A. 2015 Dynamics of laser-induced bubble pairs. J. Fluid Mech. 771, 706742.CrossRefGoogle Scholar
Koch, M., Lechner, C., Reuter, F., Köhler, K., Mettin, R. & Lauterborn, W. 2016 Numerical modeling of laser generated cavitation bubbles with the finite volume and volume of fluid method, using openfoam. Comput. Fluids 126, 7190.CrossRefGoogle Scholar
Lamb, H. 1924 Hydrodynamics. University Press.Google Scholar
Lauterborn, W. & Bolle, H. 1975 Experimental investigations of cavitation-bubble collapse in the neighbourhood of a solid boundary. J. Fluid Mech. 72 (2), 391399.CrossRefGoogle Scholar
Lechner, C., Lauterborn, W., Koch, M. & Mettin, R. 2020 Jet formation from bubbles near a solid boundary in a compressible liquid: numerical study of distance dependence. Phys. Rev. Fluids 5 (9), 093604.CrossRefGoogle Scholar
Ogasawara, T., Tsubota, N., Seki, H., Shigaki, Y. & Takahira, H. 2015 Experimental and numerical investigations of the bubble collapse at the center between rigid walls. In Journal of Physics: Conference Series (ed. M. Farhat & A. Müller), vol. 656, p. 012031. IOP Publishing.CrossRefGoogle Scholar
Plesset, M.S. & Chapman, R.B. 1971 Collapse of an initially spherical vapour cavity in the neighbourhood of a solid boundary. J. Fluid Mech. 47 (2), 283290.CrossRefGoogle Scholar
Reese, H., Schadel, R., Reuter, F. & Ohl, C.-D. 2022 Microscopic pumping of viscous liquids with single cavitation bubbles. J. Fluid Mech. 944, A17.CrossRefGoogle Scholar
Reuter, F., Gonzalez-Avila, S.R., Mettin, R. & Ohl, C.-D. 2017 Flow fields and vortex dynamics of bubbles collapsing near a solid boundary. Phys. Rev. Fluids 2, 064202.CrossRefGoogle Scholar
Reuter, F. & Ohl, C.-D. 2021 Supersonic needle-jet generation with single cavitation bubbles. Appl. Phys. Lett. 118 (13), 134103.CrossRefGoogle Scholar
Reuter, F., Zeng, Q. & Ohl, C.-D. 2022 The Rayleigh prolongation factor at small bubble to wall stand-off distances. J. Fluid Mech. 944, A11.CrossRefGoogle Scholar
Rosselló, J.M., Lauterborn, W., Koch, M., Wilken, T., Kurz, T. & Mettin, R. 2018 Acoustically induced bubble jets. Phys. Fluids 30 (12), 122004.CrossRefGoogle Scholar
Rusche, H. 2003 Computational fluid dynamics of dispersed two-phase flows at high phase fractions. PhD thesis, Imperial College London.Google Scholar
Schlichting, H. & Gersten, K. 2016 Boundary-Layer Theory. Springer.Google Scholar
Seki, H., Fujita, K., Ogasawara, T. & Takahira, H. 2015 Numerical simulations of growth and collapse of a bubble between two parallel walls. In Fluids Engineering Division Summer Meeting, vol. 57230, p. V02AT05A004. ASME.CrossRefGoogle Scholar
Taleghani, M.H., Khodadadi, S., Maddahian, R. & Mokhtari-Dizaji, M. 2023 Enhancing the bubble collapse energy using the electrohydrodynamic force. Phys. Fluids 35 (5), 053316.CrossRefGoogle Scholar
Zeng, Q., An, H. & Ohl, C.-D. 2022 Wall shear stress from jetting cavitation bubbles: influence of the stand-off distance and liquid viscosity. J. Fluid Mech. 932, A14.CrossRefGoogle Scholar
Zeng, Q., Gonzalez-Avila, S.R. & Ohl, C.-D. 2020 Splitting and jetting of cavitation bubbles in thin gaps. J. Fluid Mech. 896, A28.CrossRefGoogle Scholar
Zeng, Q., Gonzalez-Avila, S.R., Ten Voorde, S. & Ohl, C.-D. 2018 Jetting of viscous droplets from cavitation-induced Rayleigh–Taylor instability. J. Fluid Mech. 846, 916943.CrossRefGoogle Scholar
Zhang, A., Li, S., Cui, P., Li, S. & Liu, Y. 2023 A unified theory for bubble dynamics. Phys. Fluids 35 (3), 033323.Google Scholar
Figure 0

Figure 1. (a) Schematic diagram of the problem and parameters. The dashed circle indicates the nucleated bubble. The solid hat-shaped curve represents a bubble at its maximum expansion with a spherical equivalent radius $R_{max}$. The distant wall is located at a distance $H$ from the proximal wall. The liquid film between the bubble's lower wall and the boundary has a thickness of $h$. (b) Schematic diagrams at the moments prior to (solid line) and immediately after (dashed line) the formation of the micro-jet. Here $R_s$ and $R_r$ indicate the radius of the bubble's apex and the bubble's radial extension, respectively. Curved arrows sketch the reverse stagnation-point flow during the collapse. (c) Comparison of experiments and simulation on bubble radii as a function of time for two selected cases with different $\eta =2.23$ and $\eta =2.99$ but similar stand-off distance $\gamma \approx 0.48$; areas coded by light colour indicate regions of experimental uncertainties. Black dashed line is the numerical bubble dynamics for the unconfined case, $\eta \rightarrow \infty$.

Figure 1

Figure 2. Comparison between experiments (background images) and simulation (red contours). Experiment, $R_{max}=566\pm 18\ \mathrm {\mu }{\rm m}$, $d=273\pm 18\ \mathrm {\mu }{\rm m}$, $H=1260\pm 18\ \mathrm {\mu }{\rm m}$, $U_{jet}=95.6\pm 34.1\ {\rm m}\ {\rm s}^{-1}$, $\gamma \approx 0.48$, $\eta \approx 2.23$; simulation, $R_{max}=570.8\ \mathrm {\mu }{\rm m}$, $d=270\ \mathrm {\mu }{\rm m}$, $H=1260\ \mathrm {\mu }{\rm m}$, $U_{jet}=109.2\ {\rm m}\ {\rm s}^{-1}$, $\gamma =0.47$, $\eta =2.21$, $p_0=946$ bar. The blue contours are computed bubble dynamics for a larger $\eta =2.95$ but the same $\gamma =0.47$, $U_{jet}=71.9\ {\rm m}\ {\rm s}^{-1}$. Numbers in the frames are time in microseconds.

Figure 2

Figure 3. Jet formation for three selected cases with a constant $\gamma =0.55$ but various $\eta$. (a) $H=5$ mm, $\eta =8.42$; (b) $H=1.4$ mm, $\eta =2.42$; (c) $H=1.1$5 mm, $\eta =2.00$. All images focus on a region with a height of 1 mm and width of 0.5 mm. Red lines represent the bubble interfaces. The pressure and velocity magnitude are colour-coded on the left and right parts of each panel, respectively. Arrows indicate only the direction of the flow. Note the range of pressure and velocity for the colours are different for the three cases.

Figure 3

Figure 4. Comparison of the motion of the bubble's top surface for six selected $\eta$ but constant $\gamma \approx 0.55$. (a) Position of top surfaces vs time. (b) Velocities of the top surfaces.

Figure 4

Figure 5. (a) Micro-jet velocities as a function of $\eta$. Solid lines are computed results for 13 selected values of $\gamma$ between 0.46 and 1.45; black filled squares, diamonds and circles represent the experimental jet velocities for $\gamma \approx 0.47$, 1.11 and 1.43, respectively. A longitude scale for $\eta$ is used to show better comparisons between different $\gamma$. (b) All data collapse onto an orange dashed line $U_{jet}^{*}=1+1/3(\eta -\lambda ^{1.2})^{-2}$, $U_{jet}^{*}=U_{jet}/U_{\infty }$ with $U_{\infty }$ the jet velocity for $\eta \to \infty$ for a certain $\gamma$. Opacity from dark to light indicates increasing $\eta$. The inset depicts $U_{\infty }$ as a function of $\gamma$; the orange solid line is $U_{\infty }=87.94\gamma ^{0.5}$ for an ambient pressure of 1 atm.

Figure 5

Figure 6. Comparison of $U_{jet}$ as a function of $\eta$ for bubbles with three selected $\gamma$ and liquids of five different viscosities. Dashed lines are calculated using (4.7).