In this paper, phase correction and amplitude compensation are introduced to a previously developed mixed domain method (MDM), which is only accurate for modeling wave propagation in weakly heterogeneous media. Multiple reflections are also incorporated with the one-way model to improve the accuracy. The resulting model is denoted as the modified mixed-domain method (MMDM) and is numerically evaluated for its accuracy and efficiency using two distinct cases: a layered medium and a human skull. It is found that the MMDM is significantly more accurate than the MDM for strongly heterogeneous media, especially when the phase aberrating layer is roughly perpendicular to the acoustic beam. Additionally, convergence study suggests that the secondorder reflection is sufficient for wave modeling in lossy biological media. The method developed in this work could be used to facilitate therapeutic ultrasound for treating brain-related diseases and disorders.
2
I. IntroductionNumerical modeling of acoustic wave propagation in heterogeneous media is of great importance for medical ultrasound. In therapeutic ultrasound applications, for example, numerical simulations can be used to study the phase aberration in MR-guided focused ultrasound surgery[1], [2] and to improve the treatment outcome. For diagnostic ultrasound, numerical modeling has been used as an important tool for image reconstruction[3], [4], [5] as well as to understand the sources of image degradation in ultrasound imaging[6]. A myriad of wave propagation algorithms that take medium heterogeneities into account have been developed. A vast majority of these algorithms operate in the time-domain. Treeby et al.[7] developed a k-space time-domain (KSTD) method using the coupled nonlinear wave equation. Jing et al.[8] alternatively developed the KSTD from the Westervelt equation. Pinton et al.[9] studied a heterogeneous nonlinear attenuating full-wave model based on the finite-difference timedomain (FDTD) method. Frequency-domain methods have also been investigated. For example, Clement and Hynynen[10] combined the Angular Spectrum Approach (ASA) with ray theory to describe the propagation of ultrasound through randomly oriented, dissipative, layered media. Vyas and Christensen[11] modified the conventional ASA method to model linear wave propagation in inhomogeneous media. Most recently, a mixed domain method (MDM) for modeling linear/nonlinear wave propagation in dissipative, weakly heterogeneous media has been presented[12], [13]. A detailed summary of modern ultrasound modeling algorithms can be found in a review paper[14]. Although there are many existing ultrasound numerical models, none can currently achieve efficient yet sufficiently accurate simulations for linear/nonlinear acoustic wave propagation in large-scale, strongly heterogeneous media. Driven by this motivation, this paper aims to establish 3 and validate an accuracy-efficiency balanced numerical model for simulating acoustic wave propagation in strongly heterogeneous media. Within the realm of biome...