Accounting for sliding motion in fast numerical simulations of abdominal HIFU applications for targets under respiratory motion
© Schwenke et al; licensee BioMed Central Ltd. 2015
Published: 30 June 2015
Non-invasive ablation and hyperthermia treatments of abdominal organs using High Intensity Focused Ultrasound (HIFU) still impose severe challenges. Especially for liver and kidney treatments, the target motion and accessibility over the entire respiratory cycle make the application of HIFU difficult. To allow for a safe, effective, and efficient treatment, detailed planning and monitoring of the intervention is needed. Current general-purpose physics simulation software cannot fulfill the real-time requirements of the therapy-monitoring application. Our goal is to improve this by numerical simulation to allow for fast and accurate treatment planning and real-time detection of motion-induced risks during HIFU-therapy.
We developed a method to numerically simulate HIFU using patient-specific motion information in real time. The bio-heat equation describing the temperature distribution in the patient is mathematically transformed to a static reference anatomy. The numerical solution can then be performed in the reference anatomy on a static computational domain. Patient-specific respiratory motion information is provided by an abdominal motion model that can be fed with sparse tracking data in the clinical setting. The respiratory motion is modelled allowing sliding motion of the inner organs along the abdominal wall.
Motion information is fed into the simulation by an unstructured moving point cloud. On top of this point cloud a tetrahedral mesh is build to interpolate the motion information. The mesh is split at the abdominal wall into individually moving parts to allow the sliding motion of the liver along the abdominal wall. The numerical simulation thereby resolves the sliding boundary and allows the prediction of the temperature rise also in the beam path between transducer and target. Fast US beam steering is simulated using ultrasound pressure field pre-computations for key states over the respiratory cycle. Any remaining location mismatch between closest pre-computation and actual target location is compensated using an additional transformation. The numerical simulation method is integrated into a HIFU-therapy planning tool. Based on image data the target volume, anatomical and risk structures are defined. Treatment planning is then performed guided by direct feedback from the numerical simulation including respiratory motion of the domain. Using the feedback from the simulation, the user can manually optimize transducer and target location and HIFU parameters.
Results and conclusions
The research leading to these results has received funding from the European Community’s Seventh Framework Programme FP7/2007-2013 under grant agreement n° 611889.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.