Numerical Simulation of High Intensity Focused Ultrasound Temperature Distribution for Transcranial Brain Therapy

To provide a reference for the HIFU clinical therapeutic planning, the temperature distribution and lesion volume are analyzed by the numerical simulation. The adopted numerical simulation is based on a transcranial ultrasound therapy model, including an 8 annular-element curved phased array transducer. The acoustic pressure and temperature elevation are calculated by using the approximation of Westervelt Formula and the Pennes Heat Transfer Equation. In addition, the Time Reversal theory and eliminating hot spot technique are combined to optimize the temperature distribution. With different input powers and exposure times, the lesion volume is evaluated based on temperature threshold theory. The lesion region could be restored at the expected location by the time reversal theory. Although the lesion volume reduces after eliminating the peak temperature in the skull and more input power and exposure time is required, the injury of normal tissue around skull could be reduced during the HIFU therapy. The prediction of thermal deposition in the skull and the lesion region could provide a reference for clinical therapeutic dose.