New semi-analytical method for fast transcranial ultrasonic field simulation.

Objective. To optimize and ensure the safety of ultrasound brain therapy, personalized transcranial ultrasound simulations are very useful. They allow to predict the pressure field, depending on the patient skull and probe position. Most transcranial ultrasound simulations are based on numerical methods which have a long computation time and a high memory usage. The goal of this study is to develop a new semi-analytical field computation method that combines realism and computation speed. Approach. Instead of the classic ray tracing, the ultrasonic paths are computed by time of flight minimization. Then the pressure field is computed using the pencil method. This method requires a smooth and homogeneous skull model. The simulation algorithm, so-called SplineBeam, was numerically validated, by comparison with existing solvers, and experimentally validated by comparison with hydrophone measured pressure fields through an ex vivo human skull. Main results. SplineBeam simulated pressure fields were close to the experimentally measured ones, with a focus position difference of the order of the positioning error and a maximum pressure difference lower than 6.02%. In addition, for those configurations, SplineBeam computation time was lower than another simulation software, k-Wave's, by two orders of magnitude, thanks to its capacity to compute the field only at the focal spot. Significance. These results show the potential of this new method to compute fast and realistic transcranial pressure fields. The combination of this two assets makes it a promising tool for real time transcranial pressure field prediction during ultrasound brain therapy interventions.