A Source Term Approach for Generation of One-way Acoustic Waves in the Euler and Navier-Stokes equations

MFC: An open-source high-order multi-component, multi-phase, and multi-scale compressible flow solver

MFC is an open-source tool for solving multi-component, multi-phase, and bubbly compressible flows. It is capable of efficiently solving a wide range of flows, including droplet atomization, shock-bubble interaction, and bubble dynamics. We present the 5- and 6-equation thermodynamically-consistent diffuse-interface models we use to handle such flows, which are coupled to high-order interface-capturing methods, HLL-type Riemann solvers, and TVD time-integration schemes that are capable of simulating unsteady flows with strong shocks. The numerical methods are implemented in a flexible, modular framework that is amenable to future development. The methods we employ are validated via comparisons to experimental results for shock-bubble, shock-droplet, and shock-water-cylinder interaction problems and verified to be free of spurious oscillations for material-interface advection and gas-liquid Riemann problems. For smooth solutions, such as the advection of an isentropic vortex, the methods are verified to be high-order accurate. Illustrative examples involving shock-bubble-vessel-wall and acoustic-bubble-net interactions are used to demonstrate the full capabilities of MFC.

One-Way Wave Equation Derived from Impedance Theorem

The wave equations for longitudinal and transverse waves being used in seismic calculations are based on the classical force/moment balance. Mathematically, these equations are 2nd order partial differential equations (PDE) and contain two solutions with a forward and a backward propagating wave, therefore also called “Two-way wave equation”. In order to solve this inherent ambiguity many auxiliary equations were developed being summarized under “One-way wave equation”. In this article the impedance theorem is interpreted as a wave equation with a unique solution. This 1st order PDE is mathematically more convenient than the 2nd order PDE. Furthermore the 1st order wave equation being valid for three-dimensional wave propagation in an inhomogeneous continuum is derived.

Bubble cloud dynamics in an ultrasound field

The dynamics of bubble clouds induced by high-intensity focused ultrasound are investigated in a regime where the cloud size is similar to the ultrasound wavelength. High-speed images show that the cloud is asymmetrical; the bubbles nearest the source grow to a larger radius than the distal ones. Similar structures of bubble clouds are observed in numerical simulations that mimic the laboratory experiment. To elucidate the structure, a parametric study is conducted for plane ultrasound waves with various amplitudes and diffuse clouds with different initial void fractions. Based on an analysis of the kinetic energy of liquid induced by bubble oscillations, a new scaling parameter is introduced to characterize the dynamics. The new parameter generalizes the cloud interaction parameter originally introduced by d'Agostino & Brennen (1989). The dynamic interaction parameter controls the energy localization and consequent anisotropy of the cloud. Moreover, the amplitude of the far-field, bubble-scattered acoustics is likewise correlated with the proposed parameter. Findings of the present study not only shed light on the physics of cloud cavitation, but may also be of use to quantification of the effects of cavitation on outcomes of ultrasound therapies including HIFU-based lithotripsy.

Modeling and numerical simulation of the bubble cloud dynamics in an ultrasound field for burst wave lithotripsy

Modeling and numerical simulation of bubble clouds induced by intense ultrasound waves are conducted to quantify the effect of cloud cavitation on burst wave lithotripsy, a proposed non-invasive alternative to shock wave lithotripsy that uses pulses of ultrasound with an amplitude of O(1) MPa and a frequency of O(100) kHz. A unidirectional acoustic source model and an Eulerian-Lagrangian method are developed for simulation of ultrasound generation from a multi-element array transducer and cavitation bubbles, respectively. Parametric simulations of the spherical bubble cloud dynamics reveal a new scaling parameter that dictates both the structure of the bubble cloud and the amplitude of the far-field, bubble-scattered acoustics. The simulation further shows that a thin layer of bubble clouds nucleated near a kidney stone model can shield up to 90% of the incoming wave energy, indicating a potential loss of efficacy during the treatment due to cavitation. Strong correlations are identified between the far-field, bubble-scattered acoustics and the magnitude of the shielding, which could be used for ultrasound monitoring of cavitation during treatments. The simulations are validated by companion experiments in vitro.

Energy shielding by cavitation bubble clouds in burst wave lithotripsy

Combined laboratory experiment and numerical simulation are conducted on bubble clouds nucleated on the surface of a model kidney stone to quantify the energy shielding of the stone caused by cavitation during burst wave lithotripsy (BWL). In the experiment, the bubble clouds are visualized and bubble-scattered acoustics are measured. In the simulation, a compressible, multi-component flow solver is used to capture complex interactions among cavitation bubbles, the stone, and the burst wave. Quantitative agreement is confirmed between results of the experiment and the simulation. In the simulation, a significant shielding of incident wave energy by the bubble clouds is quantified. The magnitude of shielding can reach up to 90% of the energy of the incoming burst wave that otherwise would be transmitted into the stone, suggesting a potential loss of efficacy of stone comminution. There is a strong correlation between the magnitude of the energy shielding and the amplitude of the bubble-scattered acoustics, independent of the initial size and the void fraction of the bubble cloud within a range addressed in the simulation. This correlation could provide for real-time monitoring of cavitation activity in BWL.

Eulerian-Lagrangian method for simulation of cloud cavitation

We present a coupled Eulerian-Lagrangian method to simulate cloud cavitation in a compressible liquid. The method is designed to capture the strong, volumetric oscillations of each bubble and the bubble-scattered acoustics. The dynamics of the bubbly mixture is formulated using volume-averaged equations of motion. The continuous phase is discretized on an Eulerian grid and integrated using a high-order, finite-volume weighted essentially non-oscillatory (WENO) scheme, while the gas phase is modeled as spherical, Lagrangian point-bubbles at the sub-grid scale, each of whose radial evolution is tracked by solving the Keller-Miksis equation. The volume of bubbles is mapped onto the Eulerian grid as the void fraction by using a regularization (smearing) kernel. In the most general case, where the bubble distribution is arbitrary, three-dimensional Cartesian grids are used for spatial discretization. In order to reduce the computational cost for problems possessing translational or rotational homogeneities, we spatially average the governing equations along the direction of symmetry and discretize the continuous phase on two-dimensional or axi-symmetric grids, respectively. We specify a regularization kernel that maps the three-dimensional distribution of bubbles onto the field of an averaged two-dimensional or axi-symmetric void fraction. A closure is developed to model the pressure fluctuations at the sub-grid scale as synthetic noise. For the examples considered here, modeling the sub-grid pressure fluctuations as white noise agrees a priori with computed distributions from three-dimensional simulations, and suffices, a posteriori, to accurately reproduce the statistics of the bubble dynamics. The numerical method and its verification are described by considering test cases of the dynamics of a single bubble and cloud cavitaiton induced by ultrasound fields.