Thermal Oscillations of Nanobubbles

Nanobubble cavitation is advancing technologies in enhanced wastewater treatment, cancer therapy and diagnosis, and microfluidic cleaning. Current macroscale models predict that nanobubble oscillations should be isothermal, yet recent studies suggest that they are adiabatic with an associated increase in natural frequency, which becomes challenging when characterizing nanobubble sizes using ultrasound in experiments. We derive a new theoretical model that considers the nonideal nature of the nanobubble's internal gas phase and nonequilibrium effects, by employing the van der Waals (vdW) equation of state and implementing a temperature jump term at the liquid-gas interface, respectively, finding excellent agreement with molecular dynamics (MD) simulations. Our results reveal how adiabatic behavior could be erroneously interpreted when analyzing the thermal response of the gas using the commonly employed polytropic process and explain instead how nanobubble oscillations are physically closer to their isothermal limit.

Shock-induced collapse of surface nanobubbles

The collapse of cavitation bubbles often releases high-speed liquid jets capable of surface damage, with applications in drug delivery, cancer treatment, and surface cleaning. Spherical cap-shaped surface nanobubbles have previously been found to exist on immersed substrates. Despite being known nucleation sites for cavitation, their collapsing dynamics are currently unexplored. Here, we use molecular dynamics simulations to model the shock-induced collapse of different surface nanobubble sizes and contact angles. Comparisons are made with additional collapsing spherical nanobubble simulations near a substrate, to investigate the differences in their jet formation and resulting substrate pitting damage. Our main finding is that the pitting damage in the surface nanobubble simulations is greatly reduced, when compared to the spherical nanobubbles, which is primarily caused by the weaker jets formed during their collapse. Furthermore, the pit depths for surface nanobubble collapse do not depend on bubble size, unlike in the spherical nanobubble cases, but instead depend only on their contact angle. We also find a linear scaling relationship for all bubble cases between the final substrate damage and the peak pressure impulse at the impact centre, which can now be exploited to assess the relative damage in other computational studies of collapsing bubbles. We anticipate the more controlled surface-damage features produced by surface nanobubble cavitation jets will open up new applications in advanced manufacturing, medicine, and precision cleaning.

Forced oscillation dynamics of surface nanobubbles

Surface nanobubbles have potential applications in the manipulation of nanoscale and biological materials, waste-water treatment, and surface cleaning. These spherically capped bubbles of gas can exist in stable diffusive equilibrium on chemically patterned or rough hydrophobic surfaces, under supersaturated conditions. Previous studies have investigated their long-term response to pressure variations, which is governed by the surrounding liquid's local supersaturation; however, not much is known about their short-term response to rapid pressure changes, i.e., their cavitation dynamics. Here, we present molecular dynamics simulations of a surface nanobubble subjected to an external oscillating pressure field. The surface nanobubble is found to oscillate with a pinned contact line, while still retaining a mostly spherical cap shape. The amplitude-frequency response is typical of an underdamped system, with a peak amplitude near the estimated natural frequency, despite the strong viscous effects at the nanoscale. This peak is enhanced by the surface nanobubble's high internal gas pressure, a result of the Laplace pressure. We find that accurately capturing the gas pressure, bubble volume, and pinned growth mode is important for estimating the natural frequency, and we propose a simple model for the surface nanobubble frequency response, with comparisons made to other common models for a spherical bubble, a constant contact angle surface bubble, and a bubble entrapped within a cylindrical micropore. This work reveals the initial stages of growth of cavitation nanobubbles on surfaces, common in heterogeneous nucleation, where classical models based on spherical bubble growth break down.

Mechanical Stability of Surface Nanobubbles

Bubble cavitation is important in technologies such as noninvasive cancer treatment and diagnosis, surface cleaning, and waste-water treatment. The cavitation threshold is the critical external tensile pressure that induces unstable growth of the bubble. Surface nanobubbles have been previously shown experimentally to be stable down to -6 MPa, in disagreement with the Blake threshold, which is the classical cavitation model that predicts bulk bubbles with radii ∼100 nm should be unstable below -0.6 MPa. Here, we use molecular dynamics to simulate quasi-two-dimensional (2D) and three-dimensional (3D) nitrogen surface nanobubbles immersed in water, subject to a range of pressure drops until unstable growth is observed. We propose and assess new cavitation threshold models, derived from mechanical equilibrium analyses for both the quasi-2D and 3D cavitating bubbles. The discrepancies from the Blake threshold are attributed to the pinned contact line, within which the surface nanobubbles grow with constant lateral contact diameter, and consequently a reduced radius of curvature. We conclude with a critical discussion of previous experimental results on the cavitation of relatively large surface nanobubbles.