Fourth-order Time-stepping Scheme in Simulation of Ultrasound Propagation

Modeling of ultrasound waves propagation in hard biological materials such as bones and skull has a rapidly growing area of applications, e.g. brain cancer treatment planing, deep brain neurostimulation and neuromodulation, and opening blood brain barriers. Recently, we have developed a novel numerical model of elastic wave propagation based on the Kelvin-Voigt model accounting for linear elastic wave proration in heterogeneous absorption media. So far, the model has been using leap-frog time stepping scheme to integrate forward in time. In this poster, we present an improved model incorporating the 4th-order backward difference method and effect of this time-stepping scheme on the accuracy, runtime and memory consumption of the application