Critical corrections to models of nonlinear power dissipation of ultrasonically excited bubbles

Sound attenuation in high mach number oscillating bubble media

We study theoretically and numerically sound attenuation in bubble-containing media when the bubbles are freely oscillating at high Mach numbers. This paper expands one of the main forms of bubble-related acoustic damping factors by extending the previous theories to higher Mach numbers, further improves the theories of nonlinear sound propagation in bubble-containing media. A nonlinear sound propagation model incorporating second-order liquid compression terms is developed, expressing the sound velocity and density in the medium as a function of the driving pressure, and taking into account the higher-order liquid compression effects on sound propagation. The correctness of the proposed model is verified by comparing with a linear model and a nonlinear model containing only low-order Mach number terms. When the bubble oscillates at a high Mach number, radiation damping, which is directly related to Mach number, becomes the main damping component affecting sound attenuation. The higher the driving amplitude, the stronger the nonlinear effect, and the greater the impact of high-order liquid compression effects on the sound attenuation, the more necessary it is to use the proposed model to calculate the sound attenuation. For high Mach numbers, varying the bubble radius and bubble number density, respectively, the difference between the proposed model and the model containing only low-order Mach number terms in capturing the pressure-dependent attenuation is calculated. Due to stronger radiation damping in smaller bubbles, the effect of compressibility becomes more important. The smaller the bubble radius, the greater the half-quality factor of the curve related to the difference in attenuation calculated by the two models, the more necessary it is to calculate the pressure-dependent attenuation using the proposed model. Here, the half-quality factor is defined as the corresponding frequency bandwidth when the curve falls from the maximum value to 22 times. Without considering the coupling effect between bubbles, the half-quality factor of the curve is not affected by the bubble number density.

Probing the pressure dependence of sound speed and attenuation in bubbly media: Experimental observations, a theoretical model and numerical calculations

The problem of attenuation and sound speed of bubbly media has remained partially unsolved. Comprehensive data regarding pressure-dependent changes of the attenuation and sound speed of a bubbly medium are not available. Our theoretical understanding of the problem is limited to linear or semi-linear theoretical models, which are not accurate in the regime of large amplitude bubble oscillations. Here, by controlling the size of the lipid coated bubbles (mean diameter of ≈5.4μm), we report the first time observation and characterization of the simultaneous pressure dependence of sound speed and attenuation in bubbly water below, at and above microbubbles resonance (frequency range between 1-3 MHz). With increasing acoustic pressure (between 12.5-100 kPa), the frequency of the peak attenuation and sound speed decreases while maximum and minimum amplitudes of the sound speed increase. We propose a nonlinear model for the estimation of the pressure dependent sound speed and attenuation with good agreement with the experiments. The model calculations are validated by comparing with the linear and semi-linear models predictions. One of the major challenges of the previously developed models is the significant overestimation of the attenuation at the bubble resonance at higher void fractions (e.g. 0.005). We addressed this problem by incorporating bubble-bubble interactions and comparing the results to experiments. Influence of the bubble-bubble interactions increases with increasing pressure. Within the examined exposure parameters, we numerically show that, even for low void fractions (e.g. 5.1×10-6) with increasing pressure the sound speed may become 4 times higher than the sound speed in the non-bubbly medium.

Numerical investigation of the translational motion of bubbles: The comparison of capabilities of the time-resolved and the time-averaged methods

In the present study, the accuracies of two different numerical approaches used to model the translational motion of acoustic cavitational bubble in a standing acoustic field are compared. The less accurate but less computational demanding approach is to decouple the equation of translational motion from the radial oscillation, and solve it by calculating the time-averaged forces exerted on the bubble for one acoustic cycle. The second approach is to solve the coupled ordinary differential equations directly, which provides more accurate results with higher computational effort. The investigations are carried out in the parameter space of the driving frequency, pressure amplitude and equilibrium radius. Results showed that both models are capable to reveal stable equilibrium positions; however, in the case of oscillatory solutions, the difference in terms of translational frequency may be more than three fold, and the amplitude of translational motion obtained by the time-averaged method is roughly 1.5 times higher compared to the time-resolved solution at particular sets of parameters. This observation implies that where the transient behaviour is important, the time-resolved approach is the proper choice for reliable results.

Characterization of the interaction of nanobubble ultrasound contrast agents with human blood components

Nanoscale ultrasound contrast agents, or nanobubbles, are being explored in preclinical applications ranging from vascular and cardiac imaging to targeted drug delivery in cancer. These sub-micron particles are approximately 10x smaller than clinically available microbubbles. This allows them to effectively traverse compromised physiological barriers and circulate for extended periods of time. While various aspects of nanobubble behavior have been previously examined, their behavior in human whole blood has not yet been explored. Accordingly, herein we examined, for the first time, the short and long-term effects of blood components on nanobubble acoustic response. We observed differences in the kinetics of backscatter from nanobubble suspensions in whole blood compared to bubbles in phosphate buffered saline (PBS), plasma, or red blood cell solutions (RBCs). Specifically, after introducing nanobubbles to fresh human whole blood, signal enhancement, or the magnitude of nonlinear ultrasound signal, gradually increased by 22.8 ± 13.1% throughout our experiment, with peak intensity reached within 145 s. In contrast, nanobubbles in PBS had a stable signal with negligible change in intensity (−1.7 ± 3.2%) over 8 min. Under the same conditions, microbubbles made with the same lipid formulation showed a −56.8 ± 6.1% decrease in enhancement in whole blood. Subsequent confocal, fluorescent, and scanning electron microscopy analysis revealed attachment of the nanobubbles to the surface of RBCs, suggesting that direct interactions, or hitchhiking, of nanobubbles on RBCs in the presence of plasma may be a possible mechanism for the observed effects. This phenomenon could be key to extending nanobubble circulation time and has broad implications in drug delivery, where RBC interaction with nanoparticles could be exploited to improve delivery efficiency.

•

Show distinct signal enhancement curve of nanobubbles (NB) in human whole blood

•

Compare signal enhancement over time in blood with NBs and microbubbles (MBs)

•

Examine the effect of blood components on NB and MB movement over time

•

Demonstrate, with microscopy, that NBs localize to the RBC surface in whole blood

Comparison of frequency domain and time domain methods for the numerical simulation of contactless ultrasonic cavitation

The use of a top-mounted electromagnetic induction coil has been demonstrated as a contactless alternative to traditional ultrasonic treatment (UST) techniques that use an immersed mechanical sonotrode for the treatment of metals in the liquid state. This method offers similar benefits to existing UST approaches, including degassing, grain refinement, and dispersion of nanoparticles, while also preventing contact contamination due to erosion of the sonotrode. Contactless treatment potentially extends UST to high temperature or reactive melts. Generally, the method relies on acoustic resonance to reach pressure levels suitable for inertial cavitation and as a result the active cavitation volume tends to lie deep in the melt rather than in the small volume surrounding the immersed sonotrode probe. Consequently, (i) with suitable tuning of the coil supply frequency for resonance, the treatment volume can be made arbitrarily large, (ii) the problem of shielding and pressure wave attenuation suffered by the immersed sonotrode is avoided. However, relying on acoustic resonance presents problems: (i) the emergence of bubbles alters the speed of sound, resonance is momentarily lost, and cavitation becomes intermittent, (ii) as sound waves travel through and reflect on all the materials surrounding the melt, the sound characteristics of the crucible and supporting structures need to be carefully considered. The physics of cavitation coupled with this intermittent behaviour poses a challenge to sonotrode modelling orthodoxy, a problem we are trying to address in this publication. Two alternative approaches will be discussed, one of which is in the time domain and one in the frequency domain, which couple the solution of a bubble dynamics solver with that of an acoustics solver, to give an accurate prediction of the acoustic pressure generated by the induction coil. The time domain solver uses a novel algorithm to improve simulation time, by detecting an imminent bubble collapse and prescribing its subsequent behaviour, rather than directly solving a region that would normally require extremely small time steps. This way, it is shown to predict intermittent cavitation. The frequency domain solver for the first time couples the nonlinear Helmholtz model used for studying cavitation, with a background source term for the contribution of Lorentz forces. It predicts comparable RMS pressures to the time domain solver, but not the intermittent behaviour due to the underlying harmonic assumption. As further validation, the frequency domain method is also used to compare the generated acoustic pressure with that of traditional UST using a mechanical sonotrode.

Weakly nonlinear propagation of focused ultrasound in bubbly liquids with a thermal effect: Derivation of two cases of Khokolov-Zabolotskaya-Kuznetsoz equations

A physico-mathematical model composed of a single equation that consistently describes nonlinear focused ultrasound, bubble oscillations, and temperature fluctuations is theoretically proposed for microbubble-enhanced medical applications. The Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation that has been widely used as a simplified model for nonlinear propagation of focused ultrasound in pure liquid is extended to that in liquid containing many spherical microbubbles, by applying the method of multiple scales to the volumetric averaged basic equations for bubbly liquids. As a result, for two-dimensional and three-dimensional cases, KZK equations composed of the linear combination of nonlinear, dissipation, dispersion, and focusing terms are derived. Especially, the dissipation term depends on three factors, i.e., interfacial liquid viscosity, liquid compressibility, and thermal conductivity of gas inside bubbles; the thermal conduction is evaluated by using four types of temperature gradient models. Finally, we numerically solve the derived KZK equation and show a moderate temperature rise appropriate to medical applications.

A single oscillating bubble in liquids with high Mach number

The oscillation characteristics of a single bubble and its induced radiation pressure and the dissipated power are essential for a wide range of applications. For bubble oscillations with high Mach number, the influence of the liquid compressibility is significantly strong and should be fully considered. In the present paper, the bubble wall motion equation with the second-order Mach number is employed for investigating a free oscillating bubble in the liquid with numerical and experimental verifications. For the purpose of comparisons, the revised Keller-Miksis equation up to the first-order Mach number is solved with the same conditions (e.g. the initial conditions and the ambient pressure). Through our simulations, comparing with the predictions by the first-order equation, we find that: (1) The bubble radius, the bubble wall radial velocity and the bubble wall radial acceleration predicted by the second-order equation with high Mach number are significantly different respectively, and the dimensionless differences increase with the increase of the Mach number. (2) The valid prediction range of the second-order equation is much larger. (3) The dissipated power predicted by the second-order equation with high Mach number is smaller.

GPU accelerated numerical investigation of the spherical stability of an acoustic cavitation bubble excited by dual-frequency

The spherical stability of an acoustic cavitation bubble under dual-frequency excitation is investigated numerically. The radial dynamics is described by the Keller-Miksis equation, which is a second-order ordinary differential equation. The surface dynamics is modelled by a set of linear ordinary differential equation according to Hao and Prosperetti (1999), which takes into account the effect of vorticity by boundary layer approximation. Due to the large amount of investigated parameter combinations, the numerical computations were carried out on graphics processing units. The results showed that for bubble size between R_E=2μm and 4μm, the combination of a low and a high frequency, and the combination of two close but not equal frequencies are important to prevent the bubble losing its shape stability, while reaching the chemical threshold (R_max/R_E=3) (Kalmár et al., 2020). The phase shift between harmonic components of dual-frequency excitation has no effect on the shape stability.

Bursting Microbubbles: How Nanobubble Contrast Agents Can Enable the Future of Medical Ultrasound Molecular Imaging and Image-Guided Therapy

The field of medical ultrasound has undergone a significant evolution since the development of microbubbles as contrast agents. However, due to their size, microbubbles remain in the vasculature, and therefore have limited clinical applications. Building a better - and smaller - bubble can expand the applications of contrast-enhanced ultrasound by allowing bubbles to extravasate from blood vessels - creating new opportunities. In this review, we summarize recent research on the formulation and use of NBs as imaging agents and as therapeutic vehicles. We discuss the ongoing debates in the field and reluctance to accepting NBs as an acoustically active construct and a potentially impactful clinical tool that can help shape the future of medical ultrasound. We hope that the overview of key experimental and theoretical findings in the NB field presented in this paper provides a fundamental framework that will help clarify NB-ultrasound interactions and inspire engagement in the field.

On the threshold of 1/2 order subharmonic emissions in the oscillations of ultrasonically excited bubbles

The pressure threshold for 1/2 order subharmonic (SH) emissions and period doubling during the oscillations of ultrasonically excited bubbles is thought to be minimum when the bubble is sonicated with twice its resonance frequency (f_r). This estimate is based on studies that simplified or neglected the effects of thermal damping. In this work, the nonlinear dynamics of ultrasonically excited bubbles is investigated accounting for the thermal dissipation. Results are visualized using bifurcation diagrams as a function of pressure. Here we show that, and depending on the gas, the pressure threshold for 1/2 order SHs can be minimum at a frequency between 0.5f_r≤f≤0.6f_r. In this frequency range, the generation of 1/2 order SHs are due to the occurrence of 5/2 order ultra-harmonic resonance. The stability of such oscillations is size dependent. For an air bubble immersed in water, only bubbles bigger than 1 μm in diameter are able to emit non-destructive SHs in these frequency ranges.

Nonlinear power loss in the oscillations of coated and uncoated bubbles: Role of thermal, radiation and encapsulating shell damping at various excitation pressures

This study presents the fundamental equations governing the pressure dependent disipation mechanisms in the oscillations of coated bubbles. A simple generalized model (GM) for coated bubbles accounting for the effect of compressibility of the liquid is presented. The GM was then coupled with nonlinear ODEs that account for the thermal effects. Starting with mass and momentum conservation equations for a bubbly liquid and using the GM, nonlinear pressure dependent terms were derived for power dissipation due to thermal damping (Td), radiation damping (Rd) and dissipation due to the viscosity of liquid (Ld) and coating (Cd). The pressure dependence of the dissipation mechanisms of the coated bubble have been analyzed. The dissipated energies were solved for uncoated and coated 2-20 μm in bubbles over a frequency range of 0.25f_r-2.5f_r (f_r is the bubble resonance) and for various acoustic pressures (1 kPa-300 kPa). Thermal effects were examined for air and C3F8 gas cores. In the case of air bubbles, as pressure increases, the linear thermal model looses accuracy and accurate modeling requires inclusion of the full thermal model. However, for coated C3F8 bubbles of diameter 1-8 μm, which are typically used in medical ultrasound, thermal effects maybe neglected even at higher pressures. For uncoated bubbles, when pressure increases, the contributions of Rd grow faster and become the dominant damping mechanism for pressure dependent resonance frequencies (e.g. fundamental and super harmonic resonances). For coated bubbles, Cd is the strongest damping mechanism. As pressure increases, Rd contributes more to damping compared to Ld and Td. For coated bubbles, the often neglected compressibility of the liquid has a strong effect on the oscillations and should be incorporated in models. We show that the scattering to damping ratio (STDR), a measure of the effectiveness of the bubble as contrast agent, is pressure dependent and can be maximized for specific frequency ranges and pressures.