The composition of the gaseous phase of cavitation bubbles and its role on the collapse remains to date poorly understood. In this work, experiments of single cavitation bubbles in aqueous ammonia serve as a novel approach to investigate the effect of the vapor contained in a bubble on its collapse. We find that the higher vapor pressure of more concentrated aqueous ammonia acts as a resistance to the collapse, reducing the total energy dissipation. In line with visual observation, acoustic measurements, and luminescence recordings, it is also observed that higher vapor pressures contribute to a more spherical collapse, likely hindering the growth of interface instabilities by decreasing the collapse velocities and accelerations. Remarkably, we evidence a strong difference between the effective damping and the energy of the shock emission, suggesting that the latter is not the dominant dissipation mechanism at collapse as predicted from classical correction models accounting for slightly compressible liquids. Furthermore, our results suggest that the vapor inside collapsing bubbles gets compressed, consistently with previous studies performed in the context of single bubble sonoluminescence, addressing the question about the ability of vapors to readily condense during a bubble collapse in similar regimes. These findings provide insight into the identification of the influence of the bubble content and the energy exchanges of the bubble with its surrounding media, eventually paving the way to a more efficient use of cavitation in engineering and biomedical applications.
I. INTRODUCTION
The collapse of cavitation bubbles often leads to shock waves emission, light radiation, and rebound bubbles. The occurrence of these phenomena suggests the presence of a gaseous phase within the bubble, which is highly compressed during the collapse. However, its nature and influence on the bubble dynamics is to date still a subject of debate. Furthermore, although it is widely accepted that the non-condensable gas within the bubble undergoes adiabatic compression, the role of condensable vapors remains vague as equilibrium conditions are at best only satisfied at the bubble interface but not in the bubble interior.1–3 This problem was already introduced in the last century by Plesset,4 who speculated about the inability of vapor to change phase at the same rate of the bubble shrinkage. Successively, several sophisticated numerical models have been proposed to capture the influence of phase change on the bubble dynamics. Fujikawa and Akamatsu,5 and later Akhatov et al.,6 developed numerical models that include liquid compressibility, and heat and mass transfer. They investigated the incidence of non-equilibrium processes at the bubble wall due to thermal inertia of condensing vapor, concluding that this could lead to the occurrence of supercritical conditions at the final stage of the collapse. They therefore highlighted the possibility of the vapor to behave as a non-condensable gas in the case where the volume reduction rate of the bubble was much higher than the condensation rate. In addition, Akhatov et al.6 introduced a sticking coefficient of water vapor, which played as a tuning parameter to fit experimental data and predicted the condensed vapor at the bubble–liquid interface. Their work was further endorsed by Szeri et al.,7 who studied the heat and mass transfer during cavitation bubbles collapse, concluding that the latter occurs so fast that thermal diffusion and phase change effects are nearly obviated, as the vapor condensation rate is much slower than the bubble volume reduction rate. More recently, Magaletti et al.8 conducted a similar numerical investigation considering phase change, occurrence of supercritical conditions, thermal conduction, and liquid compressibility effects, reporting the disappearance and reappearance of the liquid–vapor interface during the final stage of the collapse because of a transition to supercritical conditions of the vapor. They also concluded that, in agreement with Fujikawa and Akamatsu,5 purely vapor bubble may be able to emit shock waves at collapse. Lately, Liang et al.9 studied the transition from nonlinear to linear oscillations of collapsing cavitation bubbles. They developed a novel approach based on the Gilmore10 model, with which they could fit the progressive condensation of water vapor during nonlinear bubble oscillations from experimental data by means of a tuning parameter, obtaining in turn the partial pressure of condensable vapor and non-condensable gas within the bubble. The results were in good agreement with the ones of Akhatov et al.6 Following this approach, Wen et al.11 were able to track the bubble dynamics of millimeter-sized spherical cavitation bubbles up to the fourth oscillation. The relevance of phase change on the bubble motion has been also discussed theoretically. Already for linear oscillations, it is possible to distinguish regimes where the vapor is trapped inside the bubble while keeping equilibrium conditions at the interface.12,13 For strongly nonlinear oscillations, Fuster et al.3 have shown that this asymptotic limit is reached for large values of the accommodation coefficient, where the net flux across the interface is eventually dominated by diffusion effects.7 Thus, the relevance of non-equilibrium conditions at the interface and the consequences of it on the bubble motion remain an open problem that needs of careful experimental investigations.
Irrespective from the fact that equilibrium conditions are sustained at the interface or not, the influence of phase change on the dynamics of bubbles has been experimentally confirmed in several works. In relation to single bubble sonoluminescence, Vazquez and Putterman14 observed an increased collapse cushioning and decreased light emission at increasing water temperature. The latter observation was further confirmed by Toegel et al.15 and Hopkins et al.,16 who highlighted the importance of the partial pressure of vapor trapped within the bubble during the collapse on the intensity of light emission. Later, Tinguely17 and Phan et al.18 investigated the dynamics of laser-induced cavitation bubbles in water at different temperatures, showing that the higher the water temperature, hence the vapor pressure, the larger the rebound bubble. To explain these effects, numerical models based on the slightly compressible versions of the Rayleigh–Plesset equation proposed by Keller and Miksis19 and Gilmore10 use the bubble internal pressure to fit the bubble radius evolution.9,20,21 However, due to technological challenges involved in measuring the bubble contents at the sub-millimeter and sub-millisecond scale within the bubbles, direct probing of the inner bubble pressure remains uncertain22 and may hinder some limitations of the model when using experimental data to fit the evolution of the bubble radius. Among others, some limitations of these models are that numerically describe the dynamics of perfectly spherical cavitation bubbles neglecting effects such as chemical reactions, phase change, or strongly nonlinear effects related to liquid and gas compressibility. In addition to the aforementioned uncertainties regarding the modeling of phase change processes, the influence of non-spherical deformations deserves particular attention as its effect has been indeed observed in various studies. Brennen23 investigated the bubble fission process due to bubble shape instabilities, concluding that the energy dissipated by the mixing and turbulence due to bubble fission may be preponderant compared to the conventional viscous and acoustic damping. Delale and Tunç24 developed a numerical model also accounting for deviations from sphericity, confirming the results of Brennen.23 Moreover, Supponen et al.25 experimentally investigated deformed cavitation bubbles and reported that bubbles experience weak jetting phenomena even with reduced anisotropy. Bubble shape perturbations have been also largely investigated in relation to single bubble sonoluminescence, where the development of hydrodynamic instabilities at the bubble–liquid interface has been shown to determine the stability diagrams in which light emission is observed.26
In this work, we investigate the influence of the vapor content during the collapse of laser-generated single cavitation bubbles in aqueous ammonia by systematically varying the ammonia mass fraction in solution. The latter two-component solutions have similar densities but different values, which we use to investigate the influence of the bubble internal composition. Compared to single bubble sonoluminescence experiments, the main difference is that it is possible to investigate the influence of the vapor content in regimes in conditions where the bubble oscillation is not stable gaining further insight about the role of phase change on extremely violent transient collapses. In addition, because the system's temperature is kept constant, we avoid some problems related to thermal expansion-related misalignment of optical components, sensitivity variation of measuring instruments, and change of laser energy absorbed by the liquid at bubble generation. These experiments allow us to investigate the influence of vapor content on various variables, including (i) the rebound size and collapse time, (ii) microscopic shape of rebound bubble, (iii) luminescence, (iv) radiated shock at collapse, and (v) liquid pressure build up prior to final collapse. Our findings show that all five vary significantly with the NH3 concentration in solution, evidencing the significant role of the latter on the bubble collapse and supporting the notion of vapor compression. Furthermore, our measurements provide evidence of the effect of the bubble contents on the acoustic emission at bubble collapse.
II. MATERIALS AND METHODS
A. Laser-induced cavitation bubbles
The experimental apparatus relies on a laser-based technique for the generation of single cavitation bubbles27–30 schematically shown in Fig. 1. A bubble arises from the plasma generated by a 9-ns Nd:YAG laser pulse (Quantel CFR 400, 532 nm) focused into a point at the center of an extended volume of liquid (aqueous ammonia in this work). The laser beam, redirected with a set of high-intensity mirrors, is enlarged tenfold with a beam expander and focused with an off-axis parabolic mirror, immersed into the liquid in the test chamber. The anisotropy parameter ζ for the generated bubbles was kept below the topological limit between spherical and toroidal collapse ( ), such that any reentrant jet did not pierce the bubble at collapse (a detailed description of ζ is given by Obreschkow et al.31 and Supponen et al.25). To this end, we generated bubbles with a maximum radius mm. This size was achieved by adjusting the laser beam energy with a neutral-density filter. Aqueous ammonia was contained in a transparent gas-tight box (cubic shape, 18 cm edge length), which prevented ammonia leaks and hence concentration variations. A total of six ammonia mass fractions were exploited in the experiments, ranging from 0 (pure water) to 0.05. The pressure and the temperature of the liquid at rest were kept constant and equal to the atmospheric pressure ( kPa) and ambient temperature ( ), respectively.
A high-speed camera (Shimadzu HPV-X2) filming at up to 10 million frames per second, backlighted by a collimated LED, recorded shadowgrams of the bubble. The camera is equipped with a 105 mm objective (Nikon AF-S Micro 1:2.8 GED), and a 2× teleconverter (Nikon AF-S TC-20E). A safety filter (532 nm high-pass filter) is mounted in front of the objective to prevent reflected laser beams to accidentally reach the camera sensor. The camera records 256 frames per film at a fixed resolution of 400 × 250 pixels at all frame rates below 10 × 106. At the latter, the resolution is halved. The results of bubble dynamics presented in this work were obtained from high-speed recordings at 500 000 fps, whereas luminescence was recorded at 10 × 106 fps. The instantaneous radius R of the bubble was obtained with automated image processing from high-speed recordings by retrieving the equivalent radius of the bubble axial cross sectional area A assuming spherical symmetry: . The bubble cross sectional area is represented by the black area in the bubble shadowgram recorded by the camera.
B. Aqueous ammonia solutions
Vapor pressure data of aqueous ammonia solutions are reported in Fig. 2(a). The density ρ of aqueous ammonia with a mass fraction from 0 to 0.05 is taken from Green and Perry32 as 998.2, 993.9, 989.5, 985.3, 981.1, and 977.0 kg m−3, respectively. The speed of sound c of aqueous ammonia is retrieved from high-speed movies and is 1479, 1487, 1494, 1510, 1516, and 1535 m s−1 for the same solutions cited before, respectively. The critical point of ammonia, as reported in Fig. 2(b), is defined at 132.5 °C and 11.28 MPa, whereas the critical point of water is defined at 374 °C and 22.06 MPa.32 At experimental conditions, ammonia and water coexist with the liquid phase in vapor form and can therefore be liquefied by compression only. As reported by Narita et al.,40 who investigated the dissolved oxygen in aqueous ammonia, the concentration of dissolved gas decreases with increasing ammonia concentration.
(a) Total vapor pressure of aqueous ammonia solutions as a function of the mass fraction of ammonia (and relative water mass fraction ) at 21 °C from Green and Perry.32 The plot also displays the partial vapor pressure of the solution's components. (b) Phase diagram of pure water and pure ammonia from Green and Perry.32 The solid lines show the liquid–vapor phase boundary for both substances. The black circles indicate the triple point, whereas the black squares indicate the critical point. L, V, and G stand for liquid, vapor, and gas, respectively. The triangle indicates the experimental conditions.
(a) Total vapor pressure of aqueous ammonia solutions as a function of the mass fraction of ammonia (and relative water mass fraction ) at 21 °C from Green and Perry.32 The plot also displays the partial vapor pressure of the solution's components. (b) Phase diagram of pure water and pure ammonia from Green and Perry.32 The solid lines show the liquid–vapor phase boundary for both substances. The black circles indicate the triple point, whereas the black squares indicate the critical point. L, V, and G stand for liquid, vapor, and gas, respectively. The triangle indicates the experimental conditions.
C. Analytical models for bubble motion
A detailed procedure for fitting the bubble dynamics from experimental data is described in the Appendix.
III. RESULTS AND DISCUSSION
A. Bubble dynamics in aqueous ammonia
Figure 3(a) shows the radial evolution in time in normalized coordinates, vs , of single cavitation bubbles in aqueous ammonia, with R0 the radius of the bubble at its maximum expansion, and the Rayleigh43 collapse time for an empty bubble collapsing without dissipation effects. The plot distinctly illustrates the dependence of the maximum radius of the rebound bubble R1 on the concentration of the ammonia solution, which increases from 19% of R0 for the least concentrated solution to 57% of R0 for the most concentrated solution. On the other hand, R0 remained nearly constant at 1.5 mm (±6% due to laser power oscillation, see Table I) independently of the ammonia concentration. This pointed out that the influence of the ammonia content on the first bubble oscillation is negligible, suggesting that the physics of the plasma at generation is not altered by the presence of ammonia in solution and that mass transfer effects do not affect the growth of the bubble. Nevertheless, as shown in the nested plot in Fig. 3(a), and in Fig. 3(b), the bubble lifetime was gradually prolonged up to approximately 6% with increasing ammonia concentration.
(a) Radial evolution over time in normalized coordinates of single laser-induced cavitation bubbles generated in six different aqueous ammonia solutions ( from 0 to 0.05). The points are discrete data averaged over ten different experimental measurements. The color code indicates . The shaded area shows the standard deviation of the experimental measurements. The black dashed line is the solution of the Rayleigh43 model for an empty collapsing bubble. The nested plot is a magnification of the final stage of the collapse. (b) Raw hydrophone signal from acoustic measurements over time in normalized coordinate of bubbles collapsing in six different aqueous ammonia. The time t = 0 is set at the instant when the generation shock is recorded.
(a) Radial evolution over time in normalized coordinates of single laser-induced cavitation bubbles generated in six different aqueous ammonia solutions ( from 0 to 0.05). The points are discrete data averaged over ten different experimental measurements. The color code indicates . The shaded area shows the standard deviation of the experimental measurements. The black dashed line is the solution of the Rayleigh43 model for an empty collapsing bubble. The nested plot is a magnification of the final stage of the collapse. (b) Raw hydrophone signal from acoustic measurements over time in normalized coordinate of bubbles collapsing in six different aqueous ammonia. The time t = 0 is set at the instant when the generation shock is recorded.
Average radius of the bubbles at maximum expansion R0, and average maximum radius of the first rebound R1 for each ammonia concentration exploited. indicates the standard deviation of the measurements.
(-) . | R0 (mm) . | . | R1 (mm) . | . |
---|---|---|---|---|
0 | 1.53 | 6.8 | 0.29 | 9.7 |
0.01 | 1.62 | 6.5 | 0.45 | 12.9 |
0.02 | 1.55 | 6.5 | 0.60 | 7.3 |
0.03 | 1.53 | 13.1 | 0.69 | 11.8 |
0.04 | 1.51 | 10.2 | 0.79 | 8.2 |
0.05 | 1.50 | 9.3 | 0.85 | 8.7 |
(-) . | R0 (mm) . | . | R1 (mm) . | . |
---|---|---|---|---|
0 | 1.53 | 6.8 | 0.29 | 9.7 |
0.01 | 1.62 | 6.5 | 0.45 | 12.9 |
0.02 | 1.55 | 6.5 | 0.60 | 7.3 |
0.03 | 1.53 | 13.1 | 0.69 | 11.8 |
0.04 | 1.51 | 10.2 | 0.79 | 8.2 |
0.05 | 1.50 | 9.3 | 0.85 | 8.7 |
Furthermore, we noticed additional details strictly related to ammonia content in solution. For instance, high-speed shadowgraphs, displayed in Fig. 4(a), show that the rebound bubble in pure water has a rather irregular shape with various defects at its interface. The development of these deformations seems to be prevented by an augmentation of the ammonia content in solution as schematically represented in Fig. 4(b). The bubble sketched on the left-hand side has a lower vapor content, leading to a larger growth of reentrant jets and surface perturbations. On the contrary, greater vapor pressure in the bubble on the right-hand side opposes to the liquid contraction, leading to a more spherical collapse.
(a) Shadowgrams of the first rebound bubble for different ammonia mass fraction in water . The pictures display the rebound bubble taken at about the same expansion stage. The white line indicates 0.5-mm scale. (b) Illustration of the possible effect of vapor at the final stage of the bubble collapse. (c) Snapshots of the bubble in the final stage of the collapse (filmed at 10 × 106 frames per second with an exposure time of 50 ns) in aqueous solutions with and 0.01. The white line indicates the 0.5-mm scale. (d) Hydrophone signal over time recorded at collapse of two bubbles in aqueous ammonia at of 0 and 0.05. The peak of the signals is centered at t = 0.
(a) Shadowgrams of the first rebound bubble for different ammonia mass fraction in water . The pictures display the rebound bubble taken at about the same expansion stage. The white line indicates 0.5-mm scale. (b) Illustration of the possible effect of vapor at the final stage of the bubble collapse. (c) Snapshots of the bubble in the final stage of the collapse (filmed at 10 × 106 frames per second with an exposure time of 50 ns) in aqueous solutions with and 0.01. The white line indicates the 0.5-mm scale. (d) Hydrophone signal over time recorded at collapse of two bubbles in aqueous ammonia at of 0 and 0.05. The peak of the signals is centered at t = 0.
Observation of snapshots of the last instant of the collapsing bubble displayed in Fig. 4(c) reveals a brighter hot-spot for bubbles collapsing in pure water than in aqueous ammonia at , indicating a larger energy density achieved due to a smaller minimum radius at collapse. Although light emission is no longer detected at larger ammonia concentrations, the displayed snapshots of collapsing bubbles in the most diluted solutions show multiple irregular hot-spots.
The hydrophone signals U over time, displayed in Fig. 4(d), also pointed out the dependence of the emitted shock waves on the ammonia content in solution. The signal of the shock wave emitted upon bubble collapse in pure water ( ) exhibits multiple peaks, while it smooths out at larger ammonia concentrations ( ). In addition, we surprisingly observed that the pressure build up recorded by the hydrophone before the main shock is more pronounced at larger ammonia concentration.
B. Effects of the gaseous phase on the bubble dynamics
We noticed several details about cavitation bubbles in aqueous ammonia at increasing concentration, namely, (i) longer collapse time, followed by larger rebound bubbles, (ii) increased regularity of the spherical microscopic shape of rebound bubble, (iii) reduced light emission at collapse, (iv) enhanced regularity of the radiated shock's shape at collapse, and (v) larger liquid pressure build up prior to final collapse. These observations, strictly related to the bubble contents, evidence an increased bubble internal pressure during the collapse, likely due to an increased vapor pressure of the two-component aqueous ammonia solution. In particular, the longer collapse time, followed by larger rebounds, is likely due to the cushioning effect of the compressed gaseous phase, whose pressure increases with increasing ammonia content.14,16,17
Bubble shape's irregularities, probably due to a loss of sphericity during the last moments of its collapse, decrease when the ammonia concentration is increased. Although we are unable to observe the very final stage of the collapse, a faster increase in the bubble pressure at collapse reduces the retraction of the liquid, increasing the minimum radius at collapse, and inhibits the crumpling of the bubble, thus the growth of deviations from sphericity, as illustrated in Fig. 4(b). Such deviations from sphericity are likely brought on by the development of reentrant non-piercing jets and/or hydrodynamic instabilities at the bubble–liquid interface such as Rayleigh–Taylor instabilities.20,31,44–47 Optical systems, as the one employed for this work, are well known to generate bubbles with surface perturbations.20,48 Our reasoning about an augmented bubble pressure at collapse is also backed by visualizations of the light emitted at bubble collapse.
Observations of light emission are consistent with the works of Toegel et al.,15 Vazquez and Putterman,14 Moss et al.,49 and Hopkins et al.,16 who also observed higher light intensity with lower water vapor pressure within the bubble.
Finally, the analysis of the shock waves emitted at bubble collapse support the idea of a faster increase in the bubble internal pressure at collapse due to the increased vapor pressure in more concentrated ammonia solutions. The indented shock wave signal recorded at bubble collapse in water, displayed in Fig. 4(d), probably arises due to successive collapses of different bubble fragments, whereas the smoother signal for the bubble in concentrated aqueous ammonia is synonym with more uniform bubble collapse. Moreover, the more pronounced pressure build up in the liquid surrounding the bubbles in more concentrated aqueous ammonia is likely attributed to a greater deceleration of the liquid surrounding the bubble as the internal pressure of the latter increases.33,43
All else equal, i.e., , R0, and ζ, a larger bubble internal pressure would only occur if a larger amount of matter in the gaseous form is compressed during the collapse. In this regard, although we cannot exclude that the non-condensable gas within the bubbles may have an effect on their dynamics,6 we argue that our observations are not a consequence of the varying partial pressure of non-condensable components. Specifically, non-condensable gases within laser-induced cavitation bubbles may result from vaporization of dissolved non-condensable gas at bubble generation, chemical reactions from plasma recombination, and diffusion from the liquid. For bubbles similar to those presented in this work, diffusion effects were shown to be negligible with respect to the other two gas sources.4,6,26,48 Furthermore, laser-generated gas due to chemical reactions is assumed to be proportional to the initial plasma energy density, hence to the bubble potential energy ,9,20 which was nearly constant. We could also neglect any laser-generated gas due to the different liquid composition, as the dissociation energy of water and ammonia is similar (498 kJ/mol for the first O–H bond in water and 435 kJ/mol for the NH2–H bond of ammonia).50,51 In turn, the constant bubble volume at maximum expansion highlights the independence of the plasma nature on the ammonia content. Conversely, this would likely reflect on the bubble size, absorbed energy, and/or energy partitioning at bubble generation. This actually corroborates the effectiveness of employing aqueous ammonia as a successful tool for the investigation of cavitation bubbles. Finally, vaporization due to plasma generation, proportional to the dissolved gas saturation, may give rise to variations in the partial pressure of non-condensable gas. As mentioned before, Narita et al.40 reported a decreased dissolved oxygen in aqueous ammonia with increasing ammonia concentration. However, an increased ammonia content would translate into a decreased vaporized non-condensable gas, which contravene our observations of increased bubble pressure.
Altogether, despite we are unable to directly probe the bubble contents and their behavior during the bubble collapse, our observations provide evidence that the vapor contained within a collapsing bubble is compressed. As claimed in several numerical works present in the literature,1,3,5–8,12,13,52 irrespective if equilibrium conditions are meet at the interface or not, condensation kinetics of vapors within the bubble are slower than the volume reduction rate. Similar works have also been published recently.53–56 Vapors are then unable to fully condense during the collapse due to the thermal lag at the bubble interface, likely transitioning to the non-condensable gaseous state and providing the bubble with stronger means to resist and delay its collapse. The development of a thermal boundary layer near the bubble wall was introduced earlier by Fujikawa and Akamatsu,5 Yasui,57 Storey and Szeri,1 and Akhatov et al.6 in their numerical models. In light of these theories, as the vapor pressure promptly rises with increasing ammonia concentration in the liquid [Fig. 2(a)], a greater amount of vapor remains trapped into the bubble at collapse due to the higher heat dissipation required for condensations at the bubble wall.
It is to be noted that in our experiments, the compressed gaseous phase is a two-component mixture of ammonia vapor and water vapor. As shown in Fig. 2(a), the water vapor partial pressure decreases in more concentrated solutions, whereas the ammonia vapor partial pressure quickly increases. Therefore, during the bubble collapse at increasingly concentrated aqueous ammonia solutions, the ratio of ammonia to water vapor that gets compressed rapidly changes. Moreover, as the binodal curve of ammonia locates well above that of water [Fig. 2(b)], ammonia condenses more slowly and under higher pressure than water vapor. It is however interesting to note that supercritical conditions may be reached more easily for ammonia. Furthermore, in our multicomponent system, a mass-transfer boundary layer may also concurrently limit the condensation process along with heat transfer, adding an additional level of complexity to the problem.
In addition to supporting compression of vapors, it is to be noted that our results seem to be in contrast with theories focused on the idea that luminescence is due to shock-initiated thermal emission within the collapsing bubbles.49,57 In particular, Evans58 found that small deformations of a nearly spherical converging shock wave increase as the shock converges, eventually limiting the process of light emission. This would be the case of our bubbles in more diluted solutions, which were observed with a more deformed surface. Nevertheless, they emitted more light at collapse than the most spherical bubbles in more concentrated aqueous ammonia. Our results are thus more in line with the theory of luminescence phenomena due to adiabatic compression of the bubble contents, also supported by Brenner et al.26
Finally, it is interesting and surprising to observe that the presence of ammonia did not influence the maximum radius of the bubbles. This highlights the purely inertial nature of the bubble growth. Furthermore, as the rebound ensuing the collapse is largely affected by the ammonia content, our results evidence that a higher concentration of ammonia in the water build the bubble contents up, likely increasing the vapor pressure within the bubble, highlighting the effectiveness of this technique to study the effect of the vapor pressure on cavitation bubbles.
C. Energy dissipation at bubble collapse
The energy balance in Eq. (2) revealed that a mere 1% of the initial potential energy of the bubble was reconverted into potential energy at the rebound in pure water, resulting in a substantial normalized damping of the bubble oscillation of approximately 99%, where . Conversely, in aqueous solutions with increasing ammonia content, gradually decreases down to approximately 80% at .
(a) Energy of shock waves emitted at bubble generation, normalized by the potential energy of the bubble at maximum expansion, as a function of . (b) Energy of shock waves emitted at bubble collapse, normalized by the potential energy of the bubble at maximum expansion (experiments left y-axis, and simulations right y-axis), as a function of . In both (a) and (b), points with error bars indicate the average normalized energy of the shock waves at bubble collapse obtained from 10 experiments in different aqueous ammonia with their relative standard deviation. Triangles and squares in (b) refer to the energy of shock wave at bubble collapse normalized by the potential energy of the bubble at maximum expansion obtained from simulations with the Keller and Miksis19 model and the Gilmore10 model, respectively, for bubbles at different initial non-condensable gas pressure .
(a) Energy of shock waves emitted at bubble generation, normalized by the potential energy of the bubble at maximum expansion, as a function of . (b) Energy of shock waves emitted at bubble collapse, normalized by the potential energy of the bubble at maximum expansion (experiments left y-axis, and simulations right y-axis), as a function of . In both (a) and (b), points with error bars indicate the average normalized energy of the shock waves at bubble collapse obtained from 10 experiments in different aqueous ammonia with their relative standard deviation. Triangles and squares in (b) refer to the energy of shock wave at bubble collapse normalized by the potential energy of the bubble at maximum expansion obtained from simulations with the Keller and Miksis19 model and the Gilmore10 model, respectively, for bubbles at different initial non-condensable gas pressure .
This result seems to contradict the predictions of analytical models based on the Keller and Miksis19 and Gilmore10 models, as shown in Fig. 5(b) (right y-axis), where the value of the energy emitted in the shock wave at bubble collapse , where , for different values of initial non-condensable gas is per construction comparable to (see the Appendix). Interestingly, the trend of is not reproduced by the reported experimental results.
D. Mechanisms of energy dissipation
The trend of shown in Fig. 5 indicates the presence of two regimes associated with bubble collapse in aqueous ammonia at various concentrations. At low ammonia contents ( ), results indicate that the controlling mechanism of oscillation damping is different from that of liquid compressibility, as monotonically decreases as a function of , contrary to what we expected. Conversely, at large ammonia concentrations ( ), the liquid compressibility becomes the controlling mechanism of oscillation damping as is monotonically increasing with . This experimental observation seems contrary to the predictions provided by simple models based on the Keller and Miksis19 and Gilmore10 equation, where the main damping mechanism is controlled by liquid compressibility implying that the energy of the emitted wave and the damping of the bubble should be directly correlated. The reasons why simplified models based on Keller and Miksis19 and Gilmore10 equation are not able to explain experimental observations should be topic of more careful future investigations. One possible explanation in light of diverse theories developed in the early 2000s is that deviations from sphericity are energy-consuming, which eventually dissipate energy into the liquid as heat.23,24,59 The bubble fragmentation and later coalescence phenomena have related extremely small length scales at which viscous dissipation may be magnified by several orders of magnitude compared to the spherical case. Although other mechanisms related to strongly nonlinear effects in the liquid or the presence of chemical reactions may also contribute to the effective damping of the bubble and will need detailed investigations, the results are in line with previous observations, as deviations from sphericity are the largest in the most diluted solutions, where we expect a smaller minimum radius at bubble collapse. In this case, the damping effect due to the liquid compressibility seems to be overtaken by this second mechanisms of energy dissipation. Our reasoning is further endorsed by the work of Baghdassarian et al.48 and Supponen et al.,20 who reported that laser-induced bubbles are naturally formed with surface perturbations.
Finally, although this could be a mere coincidence, it is also interesting to note that the curve peaks at approximately the same ammonia concentration whereby the ammonia vapor pressure overtakes the one of water (Fig. 2). In practice, our results suggest that the perfectly spherical bubble collapse hypothesized in the Keller and Miksis19 and the Gilmore10 models is inaccurate under certain conditions. As these models account for dissipation mainly due to liquid compressibility, the use of these models may hinder the relevance of other dissipation effects during the collapse of the bubble which are not taken into account. This may in turn compromise the accuracy of the bubble dynamics fitting from experimental data.
IV. CONCLUSION
In this study, we provide experimental evidence about the effect of vapor content on transient cavitation bubbles:
-
Consistently with the predictions for single bubble sonoluminescence experiments, the vapor is trapped inside the bubble and compressed during the violent collapse. Through an increase in the vapor pressure within the bubble by adding ammonia to water, we observe a striking role of the vapor on the amplitude of the rebound bubble. We infer that the progressively decreasing bubble deformation as a function of the ammonia concentration is due to a larger resistance of the bubble content to compression, which limits the growth of deformations and/or hydrodynamic instabilities as a consequence of the smaller values of the acceleration during the last stages of the collapse.
-
The results are endorsed by the luminescence phenomena, whose brighter emission in cavitation bubbles generated in pure water, thus with a smaller minimum radius, is synonym of larger energy density.
-
The liquid pressure build up prior to the final collapse and the shape of the shock wave emitted successively also indicate that the bubble internal pressure increases faster during the bubble collapse at the instant of collapse for increasing ammonia in solutions. Our reasoning behind the compression of vapors is that vapor–liquid phase transition are restrained, likely by mass and heat transfer.
Moreover, we demonstrate that:
-
While the initial growth and subsequent compression of a laser-induced bubble is nearly insensitive to the partial vapor pressure of the host fluid, the damping mechanisms controlling the amplitude of the rebound are greatly influenced by the vapor content. Liquid compressibility effects alone cannot explain the oscillation damping of collapsing bubbles. Some phenomena that could explain the experimental observations include phase change effects during the last stages of the bubble collapse, chemical dissociation reactions, and bubble fragmentation due to Rayleigh–Taylor instabilities, which may eventually lead to the appearance of extremely small scale agitation that would ultimately increase the viscous dissipation.
-
These other dissipation phenomena may play a preponderant role particularly in bubbles with low internal pressure at collapse, hence diluted solutions, prevalent in most of real-life applications.
It would be however interesting to obtain absolute pressure measurement obtaining the calibration constant G as a function of the ammonia content in solution. Furthermore, it would be of great relevance to study the shock waves at bubble collapse in more details, ideally investigating experimentally the shock waves behavior in the near-field.
Our findings may help improving the development of more sophisticated numerical models for improved prediction of the cavitation process. We therefore point out the urge to include more precise phase transition prediction in numerical models and call into question the disregard of aspherical perturbed collapsing bubbles. These findings also contribute to our general understanding of cavitation bubbles, with several opportunities for practical applications, including understanding and controlling cavitation erosion processes, regulating the collapse intensity for controlled sonochemical reactions, and improving cavitation-based biomedical technologies. Moreover, they also provide novel information of cavitation in aqueous ammonia, whose occurrence is to be reckoned. Alongside an already-extensive use in chemical plants, ammonia as carbon-free fuel and green hydrogen carrier has gained significant attention over the last decade.60,61
ACKNOWLEDGMENTS
This work was supported by the MSCA-ITN of the EU Horizon 2020 Research and Innovation program (Grant Agreement No. 813766), and the Swiss National Science Foundation (Grant No. 179018).
AUTHOR DECLARATIONS
Conflict of Interest
The authors have no conflicts to disclose.
Author Contributions
Davide Bernardo Preso: Conceptualization (equal); Formal analysis (equal); Investigation (lead); Methodology (equal); Validation (equal); Visualization (equal); Writing – original draft (lead). Daniel Fuster: Conceptualization (equal); Formal analysis (equal); Methodology (equal); Supervision (equal); Writing – review & editing (equal). Armand Baptiste Sieber: Investigation (supporting); Methodology (equal); Visualization (equal); Writing – review & editing (equal). Danail Obreschkow: Formal analysis (equal); Methodology (equal); Validation (equal); Writing – review & editing (equal). Mohamed Farhat: Conceptualization (equal); Formal analysis (equal); Funding acquisition (lead); Project administration (equal); Supervision (equal); Writing – review & editing (equal).
DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.
APPENDIX: ANALYTICAL MODELS FOR BUBBLE DYNAMICS
The Rayleigh–Plesset,4 the Keller and Miksis,19 and the Gilmore10 models are numerically solved with a fourth-order Runge–Kutta method with adaptive time stepping. With the last two, the internal pressure of a collapsing bubble is estimated by fitting the experimental data up to the second oscillation. The model is constrained to find the value of , the free parameter, with which the model best fits the maximum amplitude of the rebound bubble after collapse. The parameters used for simulations with the Keller and Miksis19 and the Gilmore10 model in this work are reported in Table II. An example of fitting with the two models is reported in Fig. 6. In the Keller and Miksis19 and the Gilmore10 models, we neglected viscous dissipation effects for simplicity.
(at 20 °C) (Pa) . | (kg m−3) . | c (m s−1) . | γ . | (kPa) . |
---|---|---|---|---|
2502.8 | 998.2 | 1479 | 1.35 | 97.3 |
(at 20 °C) (Pa) . | (kg m−3) . | c (m s−1) . | γ . | (kPa) . |
---|---|---|---|---|
2502.8 | 998.2 | 1479 | 1.35 | 97.3 |
Radial evolution in time, in normalized coordinates, of cavitation bubbles obtained from experiments and numerical models. In both cases, R0 is 1.5 mm. The used in the simulations is 50, 48, and 42 Pa, for the Rayleigh–Plesset,4 Keller and Miksis,19 and Gilmore10 models, respectively.