Numerical and Experimental Evaluation of High-Intensity Focused Ultrasound-Induced Lesions in Liver Tissue Ex Vivo

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.

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.

Multiphysics Analysis of Ultrasonic Shock Wave Lithotripsy and Side Effects on Surrounding Tissues

Background:

Today, the most common method for kidney stone therapy is extracorporeal shock wave lithotripsy. Current research is a numerical simulation of kidney stone fragmentation via ultrasonic shock waves. Most numerical studies in lithotripsy have been carried out using the elasticity or energy method and neglected the dissipation phenomenon. In the current study, it is solved by not only the linear acoustics equation, but also the Westervelt acoustics equation which nonlinearity and dissipation are involved.

Objective:

This study is to compare two methods for simulation of shock wave lithotripsy, clarifying the effect of shock wave profiles and stones’ material, and investigating side effects on surrounding tissues

Material and Methods:

Computational study is done using COMSOL Multiphysics, commercial software based on the finite element method. Nonlinear governing equations of acoustics, elasticity and bioheat-transfer are coupled and solved.

Results:

A decrease in the rise time of shock wave leads to increase the produced acoustic pressure and enlarge focus region. The shock wave damages kidney tissues in both linear and nonlinear simulation but the damage due to high temperature is very negligible compared to the High Intensity Focused Ultrasound (HIFU).

Conclusion:

Disaffiliation of wave nonlinearity causes a high incompatibility with reality. Stone’s material is an important factor, affecting the fragmentation

Validation of hybrid angular spectrum acoustic and thermal modelling in phantoms

In focused ultrasound (FUS) thermal ablation of diseased tissue, acoustic beam and thermal simulations enable treatment planning and optimization. In this study, a treatment-planning methodology that uses the hybrid angular spectrum (HAS) method and the Pennes' bioheat equation (PBHE) is experimentally validated in homogeneous tissue-mimicking phantoms. Simulated three-dimensional temperature profiles are compared to volumetric MR thermometry imaging (MRTI) of FUS sonications in the phantoms, whose acoustic and thermal properties are independently measured. Additionally, Monte Carlo (MC) uncertainty analysis is performed to quantify the effect of tissue property uncertainties on simulation results. The mean error between simulated and experimental spatiotemporal peak temperature rise was +0.33°C (+6.9%). Despite this error, the experimental temperature rise fell within the expected uncertainty of the simulation, as determined by the MC analysis. The average errors of the simulated transverse and longitudinal full width half maximum (FWHM) of the profiles were -1.9% and 7.5%, respectively. A linear regression and local sensitivity analysis revealed that simulated temperature amplitude is more sensitive to uncertainties in simulation inputs than in the profile width and shape. Acoustic power, acoustic attenuation and thermal conductivity had the greatest impact on peak temperature rise uncertainty; thermal conductivity and volumetric heat capacity had the greatest impact on FWHM uncertainty. This study validates that using the HAS and PBHE method can adequately predict temperature profiles from single sonications in homogeneous media. Further, it informs the need to accurately measure or predict patient-specific properties for improved treatment planning of ablative FUS surgeries.