Workshop on Virtual Reality Interaction and Physical Simulation Vriphys (2011) Focused Ultrasound - Efficient Gpu Simulation Methods for Therapy Planning

Over the past years, high intensity focused ultrasound therapy has beco me a promising therapeutic alternative for non-invasive tumor treatment. The basic idea of this interventional approa ch is to apply focused ultrasound waves to the tumor tissue such that the cells are heated and hence destroyed. Sinc e it is quite difficult to assess the quality of this non-invasive therapy, there is a dire need for computer support in planning, conduction, and monitoring of such treatments. In this work, we propose efficient simulation techniques for focused ultras ound waves as well as their heat dissemination using current graphics hardware as a numerical co-proc essor. We achieve speed-ups between 10 and 700 for the single simulation steps compared to an optimized CPU solution, over all resulting in a significant performance gain over previous approaches for simulation of focused u ltrasound.

[1]  Jens H. Krüger,et al.  GPGPU: general purpose computation on graphics hardware , 2004, SIGGRAPH '04.

[2]  H. O'neil Theory of Focusing Radiators , 1949 .

[3]  V. Frenkel,et al.  Preliminary optimization of non-destructive high intensity focused ultrasound exposures for hyperthermia applications , 2009, 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[4]  P. Meaney,et al.  The intensity dependence of lesion position shift during focused ultrasound surgery. , 2000, Ultrasound in medicine & biology.

[5]  P.J. White,et al.  A full-wave Helmholtz model for continuous-wave ultrasound transmission , 2005, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[6]  Xiaozheng Zeng,et al.  Optimal simulations of ultrasonic fields produced by large thermal therapy arrays using the angular spectrum approach. , 2009, The Journal of the Acoustical Society of America.

[7]  G T Clement,et al.  A non-invasive method for focusing ultrasound through the human skull. , 2002, Physics in medicine and biology.

[8]  Urvi Vyas,et al.  Ultrasound beam propagation using the hybrid angular spectrum method , 2008, 2008 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[9]  Nassir Navab,et al.  Visualization and GPU-accelerated simulation of medical ultrasound from CT images , 2009, Comput. Methods Programs Biomed..

[10]  Murat Manguoglu,et al.  Parallel scalable PDE-constrained optimization: antenna identification in hyperthermia cancer treatment planning , 2009, Computer Science - Research and Development.

[11]  Ruobing Li OpenCL Fast Fourier Transform , 2013 .

[12]  R. Cleveland,et al.  Numerical simulations of heating patterns and tissue temperature response due to high-intensity focused ultrasound , 2000, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[13]  Olivier Salvado,et al.  Ultrasound goes GPU: real-time simulation using CUDA , 2009, Medical Imaging.

[14]  Hao-Li Liu,et al.  Focal beam distortion and treatment planning for transrib focused ultrasound thermal therapy: a feasibility study using a two-dimensional ultrasound phased array. , 2010, Medical physics.

[15]  R. McGough,et al.  Evaluation of the angular spectrum approach for simulations of near-field pressures. , 2008, The Journal of the Acoustical Society of America.

[16]  F. Jolesz MRI-guided focused ultrasound surgery. , 2007, Annual review of medicine.

[17]  Jeffrey C. Bamber,et al.  Physical principles of medical ultrasonics , 2004 .

[18]  Jens H. Krüger,et al.  A Survey of General‐Purpose Computation on Graphics Hardware , 2007, Eurographics.

[19]  John D. Owens,et al.  GPU Computing , 2008, Proceedings of the IEEE.

[20]  Naga K. Govindaraju,et al.  A Survey of General‐Purpose Computation on Graphics Hardware , 2007 .

[21]  John D. Owens,et al.  General Purpose Computation on Graphics Hardware , 2005, IEEE Visualization.

[22]  M Tanter,et al.  Experimental demonstration of noninvasive transskull adaptive focusing based on prior computed tomography scans. , 2003, The Journal of the Acoustical Society of America.

[23]  J. Kulpa,et al.  Time-frequency analysis using NVIDIA compute unified device architecture (CUDA) , 2009, Symposium on Photonics Applications in Astronomy, Communications, Industry, and High-Energy Physics Experiments (WILGA).

[24]  Dong Zhang,et al.  Influence of ribs on the nonlinear sound field of therapeutic ultrasound. , 2007, Ultrasound in medicine & biology.