Temperature elevation by HIFU in ex vivo porcine muscle: MRI measurement and simulation study.

PURPOSE: High-intensity focused ultrasound is a rapidly developing medical technology with a large number of potential clinical applications. Computational model can play a pivotal role in the planning and optimization of the treatment based on the patient's image. Nonlinear propagation effects can significantly affect the temperature elevation and should be taken into account. In order to investigate the importance of nonlinear propagation effects, nonlinear Westervelt equation was solved. Weak nonlinear propagation effects were studied. The purpose of this study was to investigate the correlation between the predicted and measured temperature elevations and lesion in a porcine muscle. METHODS: The investigated single-element transducer has a focal length of 12 cm, an aperture of 8 cm, and frequency of 1.08 MHz. Porcine muscle was heated for 30 s by focused ultrasound transducer with an acoustic power in the range of 24-56 W. The theoretical model consists of nonlinear Westervelt equation with relaxation effects being taken into account and Pennes bioheat equation. RESULTS: Excellent agreement between the measured and simulated temperature rises was found. For peak temperatures above 85-90 °C "preboiling" or cavitation activity appears and lesion distortion starts, causing small discrepancy between the measured and simulated temperature rises. From the measurements and simulations, it was shown that distortion of the lesion was caused by the "preboiling" activity. CONCLUSIONS: The present study demonstrated that for peak temperatures below 85-90 °C numerical simulation results are in excellent agreement with the experimental data in three dimensions. Both temperature rise and lesion size can be well predicted. Due to nonlinear effect the temperature in the focal region can be increased compared with the linear case. The current magnetic resonance imaging (MRI) resolution is not sufficient. Due to the inevitable averaging the measured temperature can be 10-30 °C lower than the peak temperature. Computational fluid dynamics can provide additional important information that is lost using a state of the art MRI device.

Prediction of strong-shock structure using the bimodal distribution function.

A modified Mott-Smith method for predicting the one-dimensional shock wave solution at very high Mach numbers is constructed by developing a system of fluid dynamics equations. The predicted shock solutions in a gas of Maxwell molecules, a hard-sphere gas, and in argon using the newly proposed formalism are compared with the experimental data, direct-simulation Monte Carlo (DSMC) solution, and other solutions computed from some existing theories for Mach numbers M<50. In the limit of an infinitely large Mach number, the predicted shock profiles are also compared with the DSMC solution. The density, temperature and heat flux profiles calculated at different Mach numbers have been shown to have good agreement with the experimental and DSMC solutions.

Prediction of shock structure using the bimodal distribution function.

A modification of Mott-Smith method for predicting the one-dimensional shock wave solution is presented. Mott-Smith distribution function is used to construct the system of moment equations to study the steady-state structure of shock wave in a gas of Maxwell molecules and in argon. The predicted shock solutions using the newly proposed formalism are compared to the experimental data, direct-simulation Monte Carlo (DSMC) solution, and the solutions predicted by other existing theories for Mach numbers M<11 . The density, temperature, heat flux profiles, and shock thickness calculated at different Mach numbers have been shown to have good agreement with the experimental and DSMC solutions. In addition, the predicted shock thickness is in good agreement with the DSMC simulation result at low Mach numbers.