On the use of an absorption layer for the angular spectrum approach (L)

Nonlinear Acoustic Holography With Adaptive Sampling

Accurate and efficient numerical simulation of highly nonlinear ultrasound propagation is essential for a wide range of therapeutic and physical ultrasound applications. However, due to large domain sizes and the generation of higher harmonics, such simulations are computationally challenging, particularly in 3-D problems with shock waves. Current numerical methods are based on computationally inefficient uniform meshes that resolve the highest harmonics across the entire spatial domain. To address this challenge, we present an adaptive numerical algorithm for computationally efficient nonlinear acoustic holography. At each propagation step, the algorithm monitors the harmonic content of the acoustic signal and adjusts its discretization parameters accordingly. This enables efficient local resolution of higher harmonics in areas of high nonlinearity while avoiding unnecessary resolution elsewhere. Furthermore, the algorithm actively adapts to the signal's nonlinearity level, eliminating the need for prior reference simulations or information about the spatial distribution of the harmonic content of the acoustic field. The proposed algorithm incorporates an upsampling process in the frequency domain to accommodate the generation of higher harmonics in forward propagation and a downsampling process when higher harmonics are decimated in backward propagation. The efficiency of the algorithm was evaluated for highly nonlinear 3-D problems, demonstrating a significant reduction in computational cost with a nearly 50-fold speedup over a uniform mesh implementation. Our findings enable a more rapid and efficient approach to modeling nonlinear high-intensity focused ultrasound (HIFU) wave propagation.

Corrections to "mSOUND: An Open Source Toolbox for Modeling Acoustic Wave Propagation in Heterogeneous Media"

In the above article [1]>, the authors regret that there were some mistakes pertaining to (1)-(3) and (5).

mSOUND: An Open Source Toolbox for Modeling Acoustic Wave Propagation in Heterogeneous Media

mSOUND is an open-source toolbox written in MATLAB. This toolbox is intended for modeling linear/ nonlinear acoustic wave propagation in media (primarily biological tissues) with arbitrary heterogeneities, in which, the speed of sound, density, attenuation coefficient, power-law exponent, and nonlinear coefficient are all spatially varying functions. The computational model is an iterative one-way model based on a mixed domain method. In this article, a general guideline is given along with three representative examples to illustrate how to set up simulations using mSOUND. The first example uses the transient mixed-domain method (TMDM) forward projection to compute the transient acoustic field for a given source defined on a plane. The second example uses the frequency-specific mixed-domain method (FSMDM) forward projection to rapidly obtain the pressure distribution directly at the frequencies of interest, assuming linear or weakly nonlinear wave propagation. The third example demonstrates how to use TMDM backward projection to reconstruct the initial acoustic pressure field to facilitate photoacoustic tomography (PAT). mSOUND (https://m-sound.github.io/mSOUND/home) is designed to be complementary to existing ultrasound modeling toolboxes and is expected to be useful for a wide range of applications in medical ultrasound including treatment planning, PAT, transducer design, and characterization.

A modified mixed domain method for modeling acoustic wave propagation in strongly heterogeneous media

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 four distinct cases. It is found that the MMDM is significantly more accurate than the MDM for strongly heterogeneous media, especially when the phase aberrating layer is approximately perpendicular to the acoustic beam. Additionally, a convergence study suggests that the second-order reflection could be sufficient for cases involving high contrast inhomogeneities and large loss values (e.g., skulls). The method developed in this work could facilitate therapeutic ultrasound for treating brain-related diseases and disorders.

Simulation of the Second-Harmonic Ultrasound Field in Heterogeneous Soft Tissue Using a Mixed-Domain Method

A mixed-domain method (MDM) dubbed frequency-specific MDM (FSMDM) is introduced for the simulation of the second-harmonic ultrasound field in weakly heterogeneous media. The governing equation for the second harmonics is derived based on the quasi-linear theory. The speed of sound, nonlinear coefficient, and attenuation coefficient are all spatially varying functions in the equation. The fundamental frequency pressure field is first solved by the FSMDM and it is subsequently used as the source term for the second-harmonics equation. This equation can be again solved by the FSMDM to rapidly obtain the second-harmonic pressure field. Five 2-D cases, including one with a realistic human tissue map, are studied to systematically verify the proposed method. Results from the previously developed transient MDM are used as the benchmark solutions. Comparisons show that the two methods give similar results for all cases. More importantly, the FSMDM has a crucial advantage over the transient MDM in that it can be two orders of magnitude faster.

Numerical Modeling of Ultrasound Propagation in Weakly Heterogeneous Media Using a Mixed-Domain Method

A mixed-domain method (MDM) is presented in this paper for modeling one-way linear/nonlinear wave propagation in biological tissue with arbitrary heterogeneities, in which sound speed, density, attenuation coefficients, and nonlinear coefficients are all spatial varying functions. The present method is based on solving an integral equation derived from a Westervelt-like equation. One-dimensional problems are first studied to verify the MDM and to reveal its limitations. It is shown that this method is accurate for cases with small variation of sound speed. A 2-D case is further studied with focused ultrasound beams to validate the application of the method in the medical field. Results from the MATLAB toolbox k-Wave are used as the benchmark. Normalized root-mean-square (rms) error estimated at the focus of the transducer is 0.0133 when the coarsest mesh (1/3 of the wavelength) is used in the MDM. Fundamental and second-harmonic fields throughout the considered computational domains are compared and good agreement is observed. Overall, this paper demonstrates that the MDM is a computationally efficient and accurate method when used to model wave propagation in biological tissue with relatively weak heterogeneities.

Modeling of wave propagation for medical ultrasound: a review

Numerical modeling of medical ultrasound has advanced tremendously in the past two decades. This opens up a great number of opportunities for medical ultrasound and associated technologies. Numerous new governing equations and algorithms have emerged and been applied to studying various medical ultrasound applications, including ultrasound imaging, photo-acoustic imaging, and therapeutic ultrasound. In addition, thanks to the rapid development of computers, modeling acoustic wave propagation in three-dimensional, large-scale domains has become a reality. This article will provide an indepth literature and technical review of recent progress on numerical modeling of medical ultrasound. Future challenges will also be discussed.