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

Weakly nonlinear focused ultrasound in viscoelastic media containing multiple bubbles

To facilitate practical medical applications such as cancer treatment utilizing focused ultrasound and bubbles, a mathematical model that can describe the soft viscoelasticity of human body, the nonlinear propagation of focused ultrasound, and the nonlinear oscillations of multiple bubbles is theoretically derived and numerically solved. The Zener viscoelastic model and Keller-Miksis bubble equation, which have been used for analyses of single or few bubbles in viscoelastic liquid, are used to model the liquid containing multiple bubbles. From the theoretical analysis based on the perturbation expansion with the multiple-scales method, the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation, which has been used as a mathematical model of weakly nonlinear propagation in single phase liquid, is extended to viscoelastic liquid containing multiple bubbles. The results show that liquid elasticity decreases the magnitudes of the nonlinearity, dissipation, and dispersion of ultrasound and increases the phase velocity of the ultrasound and linear natural frequency of the bubble oscillation. From the numerical calculation of resultant KZK equation, the spatial distribution of the liquid pressure fluctuation for the focused ultrasound is obtained for cases in which the liquid is water or liver tissue. In addition, frequency analysis is carried out using the fast Fourier transform, and the generation of higher harmonic components is compared for water and liver tissue. The elasticity supresses the generation of higher harmonic components and promotes the remnant of the fundamental frequency components. This indicates that the elasticity of liquid suppresses shock wave formation in practical applications.

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.

Production and dispersion of free radicals from transient cavitation Bubbles: An integrated numerical scheme and applications

As an advanced oxidation process with a wide range of applications, sonochemistry relies on acoustic cavitation to induce free radicals for degrading chemical contaminants. The complete process includes two critical steps: the radical production inside the cavitation bubble, and the ensuing dispersion of these radicals into the bulk solution. To grasp the physicochemical details in this process, we developed an integrated numerical scheme with the ability to quantitatively describe the radical production-dispersion behavior. It employs coupled simulations of bubble dynamics, intracavity chemical reactions, and diffusion-reaction-dominated mass transport in aqueous solutions. Applying this method to the typical case of argon and oxygen bubbles, the production mechanism for the main radicals is revealed. Moreover, the temporal-spatial distribution of the radicals in the liquid phase is presented. The results demonstrate that the enhanced radical production observed in oxygen bubbles can be traced to the initiation reaction O₂ + H₂O → OH+HO₂, which requires relatively low activation energy. In the outside liquid region, the dispersion of radicals is limited by robust recombination reactions. The simulated penetration depth of OH is around 0.2 μm and agrees with reported experimental measurements. The proposed numerical approach can be employed to better capture the radical activity and is instrumental in optimizing the engineering application of sonochemistry.