Development of a Numerical Model of High-Intensity Focused Ultrasound Treatment in Mobile and Elastic Organs: Application to a Beating Heart.

High-intensity focused ultrasound (HIFU) is a promising method used to treat cardiac arrhythmias, as it can induce lesions at a distance throughout myocardium thickness. Numerical modeling is commonly used for ultrasound probe development and optimization of HIFU treatment strategies. This study was aimed at describing a numerical method to simulate HIFU thermal ablation in elastic and mobile heart models. The ultrasound pressure field is computed on a 3-D orthonormal grid using the Rayleigh integral method, and the attenuation is calculated step by step between cells. The temperature distribution is obtained by resolution of the bioheat transfer equation on a 3-D non-orthogonally structured curvilinear grid using the finite-volume method. The simulation method is applied on two regions of the heart (atrioventricular node and ventricular apex) to compare the thermal effects of HIFU ablation depending on deformation, motion type and amplitude. The atrioventricular node requires longer sonication than the ventricular apex to reach the same lesion volume. Motion considerably influences treatment duration, lesion shape and distribution in cardiac HIFU treatment. These results emphasize the importance of considering local motion and deformation in numerical studies to define efficient and accurate treatment strategies.