An iterative method for the computation of nonlinear, wide-angle, pulsed acoustic fields of medical diagnostic transducers

Computation of ultrasound propagation in a population of nonlinearly oscillating microbubbles including multiple scattering

In contrast-enhanced echography, the simulation of nonlinear propagation of ultrasound through a population of oscillating microbubbles imposes a computational challenge. Also, the numerical complexity increases because each scatterer has individual properties. To address these problems, the Iterative Nonlinear Contrast Source (INCS) method has been extended to include a large population of nonlinearly responding microbubbles. The original INCS method solves the Westervelt equation in a four-dimensional spatiotemporal domain by generating increasingly accurate field corrections to iteratively update the acoustic pressure. The field corrections are computed by the convolution of a nonlinear contrast source with the Green's function of the linear background medium. Because the convolution integral allows a coarse discretization, INCS can efficiently deal with large-scale problems. To include a population of microbubbles, these are considered as individual contrast point sources with their own nonlinear response. The field corrections are computed as before, but now, in each iteration, the temporal signature of each contrast point source is computed by solving the bubble's Marmottant equation. Physically, each iteration adds an order of multiple scattering. Here, the performance of the extended INCS method and the significance of multiple scattering is demonstrated through various results from different configurations.

An iterative method to evaluate one-dimensional pulsed nonlinear elastic wavefields and mixing of elastic waves in solids

Over the last 15 years, literature on nondestructive testing has shown that the generation of higher harmonics and nonlinear mixing of waves could be used to obtain the nonlinearity parameters of an elastic medium and thereby gather information about its state, e.g., aging and fatigue. To design ultrasound measurement setups based on these phenomena, efficient numerical modeling tools are needed. In this paper, the iterative nonlinear contrast source method for numerical modeling of nonlinear acoustic waves is extended to the one-dimensional elastic case. In particular, nonlinear mixing of two collinear bulk waves (one compressional, one shear) in a homogeneous, isotropic medium is considered, taking into account its third-order elastic constants ( A,  B,  and  C). The obtained results for nonlinear propagation are in good agreement with a benchmark solution based on the modified Burgers equation. The results for the resonant waves that are caused by the one-way and two-way mixing of primary waves are in quantitative agreement with the results in the literature [Chen, Tang, Zhao, Jacobs, and Qu, J. Acoust. Soc. Am. 136(5), 2389-2404 (2014)]. The contrast source approach allows the identification of the propagating and evanescent components of the scattered wavefield in the wavenumber-frequency domain, which provides physical insight into the mixing process and explains the propagation direction of the resonant wave.

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.

Linear and nonlinear ultrasound simulations using the discontinuous Galerkin method

A nodal discontinuous Galerkin (DG) code based on the nonlinear wave equation is developed to simulate transient ultrasound propagation. The DG method has high-order accuracy, geometric flexibility, low dispersion error, and excellent scalability, so DG is an ideal choice for solving this problem. A nonlinear acoustic wave equation is written in a first-order flux form and discretized using nodal DG. A dynamic sub-grid scale stabilization method for reducing Gibbs oscillations in acoustic shock waves is then established. Linear and nonlinear numerical results from a two-dimensional axisymmetric DG code are presented and compared to numerical solutions obtained from linear and Khokhlov-Zabolotskaya-Kuznetsov-based simulations in FOCUS. The numerical results indicate that these nodal DG simulations capture nonlinearity, thermoviscous absorption, and diffraction for both flat and focused pistons in homogeneous media.

Fast Nonlinear Ultrasound Propagation Simulation Using a Slowly Varying Envelope Approximation

Medical systems usually consider linear propagation of ultrasound, an approximation of reality. However, numerous studies have attempted to accurately simulate the nonlinear pressure wave distortion and to evaluate the contribution of harmonic frequencies. In such simulations, the computation time is very large, except for the method based on the angular spectrum scheme where the derivative order is reduced using the Fourier transform. However, the harmonic computation is usually limited to the second harmonic because of quasi-linear approximation. In this paper, a slowly varying envelope approximation (SVEA) is used in the Fourier domain to compute the entire nonlinear distortion induced, including high harmonics and nonlinear mixing frequencies. The simulation by SVEA is evaluated by comparison with other simulation tools. The obtained deviation and difference remain low enough to fully validate such an approximation. Moreover, the simulator is implemented on a GPU to obtain a very fast tool, where the full nonlinear distorted [Formula: see text] field is computed in less than 10 s.

Design of HIFU Transducers for Generating Specified Nonlinear Ultrasound Fields

Various clinical applications of high-intensity focused ultrasound have different requirements for the pressure levels and degree of nonlinear waveform distortion at the focus. The goal of this paper is to determine transducer design parameters that produce either a specified shock amplitude in the focal waveform or specified peak pressures while still maintaining quasi-linear conditions at the focus. Multiparametric nonlinear modeling based on the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation with an equivalent source boundary condition was employed. Peak pressures, shock amplitudes at the focus, and corresponding source outputs were determined for different transducer geometries and levels of nonlinear distortion. The results are presented in terms of the parameters of an equivalent single-element spherically shaped transducer. The accuracy of the method and its applicability to cases of strongly focused transducers were validated by comparing the KZK modeling data with measurements and nonlinear full diffraction simulations for a single-element source and arrays with 7 and 256 elements. The results provide look-up data for evaluating nonlinear distortions at the focus of existing therapeutic systems as well as for guiding the design of new transducers that generate specified nonlinear fields.

Focused ultrasound actuation of shape memory polymers; acoustic-thermoelastic modeling and testing

Controlled drug delivery (CDD) technologies have received extensive attention recently.

Superharmonic Imaging for Medical Ultrasound: a Review

Ultrasound with harmonics has emerged as an exceptional alternative to competitively low resolution fundamental ultrasound imaging. The use of second harmonic is already a trend now but higher harmonics are also being seen as even better option due to its improved resolution. The resolution improved with frequency but achieves penetration of reduced energy. The cumulative addition of higher harmonics during propagation yields higher harmonics giving better resolution with adequate penetration. This paper summarizes the progress of such similar decade old harmonic ultrasound imaging technique i.e., superharmonic imaging (SHI) geared towards medical field. It comprises of harmonics higher than second harmonic preferably up to 5th harmonic. We conclude that SHI can be an advanced ultrasound imaging with comprehensive high resolution and adequate penetration depth on sole and coded modes.

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.

First-arrival traveltime sound speed inversion with a priori information

PURPOSE

A first-arrival travel-time sound speed algorithm presented by Tarantola [Inverse Problem Theory and Methods for Model Parameter Estimation (SIAM, Philadelphia, PA, 2005)] is adapted to the medical ultrasonics setting. Through specification of a covariance matrix for the object model, the algorithm allows for natural inclusion of physical a priori information of the object. The algorithm's ability to accurately and robustly reconstruct a complex sound speed distribution is demonstrated on simulation and experimental data using a limited aperture.

METHODS

The algorithm is first demonstrated generally in simulation with a numerical breast phantom imaged in different geometries. As this work is motivated by the authors' limited aperture dual sided ultrasound breast imaging system, experimental data are acquired with a Verasonics system with dual, 128 element, linear L7-4 arrays. The transducers are automatically calibrated for usage in the eikonal forward model.A priori information such as knowledge of correlated regions within the object is obtained via segmentation of B-mode images generated from synthetic aperture imaging.

RESULTS

As one illustration of the algorithm's facility for inclusion ofa priori information, physically grounded regularization is demonstrated in simulation. The algorithm's practicality is then demonstrated through experimental realization in limited aperture cases. Reconstructions of sound speed distributions of various complexity are improved through inclusion of a priori information. The sound speed maps are generally reconstructed with accuracy within a few m/s.

CONCLUSIONS

This paper demonstrates the ability to form sound speed images using two opposed commercial linear arrays to mimic ultrasound image acquisition in the compressed mammographic geometry. The ability to create reasonably good speed of sound images in the compressed mammographic geometry allows images to be readily coregistered to tomosynthesis image volumes for breast cancer detection and characterization studies.

Ultrasound coefficient of nonlinearity imaging

Imaging the acoustical coefficient of nonlinearity, β, is of interest in several healthcare interventional applications. It is an important feature that can be used for discriminating tissues. In this paper, we propose a nonlinearity characterization method with the goal of locally estimating the coefficient of nonlinearity. The proposed method is based on a 1-D solution of the nonlinear lossy Westerfelt equation, thereby deriving a local relation between β and the pressure wave field. Based on several assumptions, a β imaging method is then presented that is based on the ratio between the harmonic and fundamental fields, thereby reducing the effect of spatial amplitude variations of the speckle pattern. By testing the method on simulated ultrasound pressure fields and an in vitro B-mode ultrasound acquisition, we show that the designed algorithm is able to estimate the coefficient of nonlinearity, and that the tissue types of interest are well discriminable. The proposed imaging method provides a new approach to β estimation, not requiring a special measurement setup or transducer, that seems particularly promising for in vivo imaging.

Simulations and measurements of 3-D ultrasonic fields radiated by phased-array transducers using the westervelt equation

The purpose of this work is to validate, by comparing numerical and experimental results, the ability of the Westervelt equation to predict the behavior of ultrasound beams generated by phased-array transducers. To this end, the full Westervelt equation is solved numerically and the results obtained are compared with experimental measurements. The numerical implementation of the Westervelt equation is performed using the explicit finite-difference time-domain method on a three-dimensional Cartesian grid. The validation of the developed numerical code is first carried out by using experimental data obtained for two different focused circular transducers in the regimes of small-amplitude and finite-amplitude excitations. Then, the comparison of simulated and measured ultrasonic fields is extended to the case of a modified 32-element array transducer. It is shown that the developed code is capable of correctly predicting the behavior of the main lobe and the grating lobes in the cases of zero and nonzero steering angles for both the fundamental and the second-harmonic components.

An improved wave-vector frequency-domain method for nonlinear wave modeling

In this paper, a recently developed wave-vector frequency-domain method for nonlinear wave modeling is improved and verified by numerical simulations and underwater experiments. Higher order numeric schemes are proposed that significantly increase the modeling accuracy, thereby allowing for a larger step size and shorter computation time. The improved algorithms replace the left-point Riemann sum in the original algorithm by the trapezoidal or Simpson's integration. Plane waves and a phased array were first studied to numerically validate the model. It is shown that the left-point Riemann sum, trapezoidal, and Simpson's integration have first-, second-, and third-order global accuracy, respectively. A highly focused therapeutic transducer was then used for experimental verifications. Short high-intensity pulses were generated. 2-D scans were conducted at a prefocal plane, which were later used as the input to the numerical model to predict the acoustic field at other planes. Good agreement is observed between simulations and experiments.

Modeling nonlinear wave propagation on nonuniform grids using a mapped k-space pseudospectral method

Simulating the propagation of nonlinear ultrasound waves is computationally difficult because of the dense grids needed to capture high-frequency harmonics. Here, a mapped k-space pseudospectral method is presented which allows the use of nonuniform grid spacings. This enables grid points to be clustered around steep regions of the wave field. Compared with using a uniform grid, this significantly reduces the total number of grid points needed for accurate simulations. Two methods for selecting a suitable nonuniform grid mapping are discussed.

Unmitigated numerical solution to the diffraction term in the parabolic nonlinear ultrasound wave equation

Various numerical algorithms have been developed to solve the Khokhlov-Kuznetsov-Zabolotskaya (KZK) parabolic nonlinear wave equation. In this work, a generalized time-domain numerical algorithm is proposed to solve the diffraction term of the KZK equation. This algorithm solves the transverse Laplacian operator of the KZK equation in three-dimensional (3D) Cartesian coordinates using a finite-difference method based on the five-point implicit backward finite difference and the five-point Crank-Nicolson finite difference discretization techniques. This leads to a more uniform discretization of the Laplacian operator which in turn results in fewer calculation gridding nodes without compromising accuracy in the diffraction term. In addition, a new empirical algorithm based on the LU decomposition technique is proposed to solve the system of linear equations obtained from this discretization. The proposed empirical algorithm improves the calculation speed and memory usage, while the order of computational complexity remains linear in calculation of the diffraction term in the KZK equation. For evaluating the accuracy of the proposed algorithm, two previously published algorithms are used as comparison references: the conventional 2D Texas code and its generalization for 3D geometries. The results show that the accuracy/efficiency performance of the proposed algorithm is comparable with the established time-domain methods.

Computation of nonlinear ultrasound fields using a linearized contrast source method

Nonlinear ultrasound is important in medical diagnostics because imaging of the higher harmonics improves resolution and reduces scattering artifacts. Second harmonic imaging is currently standard, and higher harmonic imaging is under investigation. The efficient development of novel imaging modalities and equipment requires accurate simulations of nonlinear wave fields in large volumes of realistic (lossy, inhomogeneous) media. The Iterative Nonlinear Contrast Source (INCS) method has been developed to deal with spatiotemporal domains measuring hundreds of wavelengths and periods. This full wave method considers the nonlinear term of the Westervelt equation as a nonlinear contrast source, and solves the equivalent integral equation via the Neumann iterative solution. Recently, the method has been extended with a contrast source that accounts for spatially varying attenuation. The current paper addresses the problem that the Neumann iterative solution converges badly for strong contrast sources. The remedy is linearization of the nonlinear contrast source, combined with application of more advanced methods for solving the resulting integral equation. Numerical results show that linearization in combination with a Bi-Conjugate Gradient Stabilized method allows the INCS method to deal with fairly strong, inhomogeneous attenuation, while the error due to the linearization can be eliminated by restarting the iterative scheme.

Single pulse frequency compounding protocol for superharmonic imaging

Second harmonic imaging is currently accepted as the standard in commercial echographic systems. A new imaging technique, coined as superharmonic imaging (SHI), combines the third till the fifth harmonics, arising during nonlinear sound propagation. It could further enhance the resolution and quality of echographic images. To meet the bandwidth requirement for SHI a dedicated phased array has been developed: a low frequency subarray, intended for transmission, interleaved with a high frequency subarray, used in reception. As the bandwidth of the elements is limited, the spectral gaps in between the harmonics cause multiple (ghost) reflection artifacts. A dual-pulse frequency compounding method aims at suppressing those artifacts at a price of a reduced frame rate. In this study we explore a possibility of performing frequency compounding within a single transmission. The traditional frequency compounding method suppresses the ripples by consecutively emitting two short Gaussian bursts with a slightly different center frequency. In the newly proposed method, the transmit aperture is divided into two parts: the first half is used to send a pulse at the lower center frequency, while the other half simultaneously transmits at a slightly higher center frequency. The suitability of the protocol for medical imaging applications in terms of the steering capabilities was performed in a simulation study with INCS and the hydrophone measurements. Moreover, an experimental study was carried out to find the optimal parameters for the clinical imaging protocol. The latter was subsequently used to obtain the images of a tissue mimicking phantom containing strongly reflecting wires. Additionally, the images of a human heart in the parasternal projection were acquired. The scanning aperture with the developed protocol amounts to approximately 90°, which is sufficient to capture the cardiac structures in the standard anatomical projections. The theoretically estimated and experimentally measured grating lobe levels are equal to -28.3 dB and -35.9 dB, respectively. A considerable improvement in the axial resolution of the SHI component (0.73 mm) at -6 dB in comparison with the third harmonic (2.23 mm) was observed. A similar comparison in terms of the lateral resolution slightly favored the superharmonic component by 0.2 mm. Additionally, the images of the tissue mimicking phantom exhibited the absence of the multiple reflection artifacts. The in-vivo acquisition allows one to clearly observe the dynamic of the mitral valve leaflets. The new method is equally effective in eliminating the ripple artifacts associated with SHI as the dual-pulse technique, while the full frame rate is maintained.

Application of the split-step Padé approach to nonlinear field predictions

We herein propose a new theoretical approach for analyzing the nonlinear propagation of directive sound beams emitted from a planar piston source with a circular aperture. The proposed approach relies on the split-step Padé approximation, which is an efficient method for obtaining wide-angle one-way wave equations, especially in underwater acoustics. Despite including only two Padé terms in the expansion, the theory was applicable to a beam angle of up to ±40° relative to the main propagation direction, the angle of which is approximately twice that of the Khokhlov-Zabolotskaya-Kuznetsov equation, which is based on parabolic approximation. In order to demonstrate the effectiveness of the newly proposed theoretical approach, we performed an experiment using an airborne ultrasonic emitter with a circular aperture of 7.5cm in radius. We drove the emitter powerfully at a 36-kHz and 40-kHz bi-frequency signal and measured the beam patterns of the primary and secondary waves, such as parametric sounds within wide propagation angles. Excellent agreement between measured data and the corresponding numerical simulations supports the validity of the proposed model equations and the computational methods for their numerical solutions.

A k-space method for moderately nonlinear wave propagation

A k-space method for moderately nonlinear wave propagation in absorptive media is presented. The Westervelt equation is first transferred into k-space via Fourier transformation, and is solved by a modified wave-vector time-domain scheme. The present approach is not limited to forward propagation or parabolic approximation. One- and two-dimensional problems are investigated to verify the method by comparing results to analytic solutions and finite-difference time-domain (FDTD) method. It is found that to obtain accurate results in homogeneous media, the grid size can be as little as two points per wavelength, and for a moderately nonlinear problem, the Courant-Friedrichs-Lewy number can be as large as 0.4. Through comparisons with the conventional FDTD method, the k-space method for nonlinear wave propagation is shown here to be computationally more efficient and accurate. The k-space method is then employed to study three-dimensional nonlinear wave propagation through the skull, which shows that a relatively accurate focusing can be achieved in the brain at a high frequency by sending a low frequency from the transducer. Finally, implementations of the k-space method using a single graphics processing unit shows that it required about one-seventh the computation time of a single-core CPU calculation.

Modeling nonlinear ultrasound propagation in heterogeneous media with power law absorption using a k-space pseudospectral method

The simulation of nonlinear ultrasound propagation through tissue realistic media has a wide range of practical applications. However, this is a computationally difficult problem due to the large size of the computational domain compared to the acoustic wavelength. Here, the k-space pseudospectral method is used to reduce the number of grid points required per wavelength for accurate simulations. The model is based on coupled first-order acoustic equations valid for nonlinear wave propagation in heterogeneous media with power law absorption. These are derived from the equations of fluid mechanics and include a pressure-density relation that incorporates the effects of nonlinearity, power law absorption, and medium heterogeneities. The additional terms accounting for convective nonlinearity and power law absorption are expressed as spatial gradients making them efficient to numerically encode. The governing equations are then discretized using a k-space pseudospectral technique in which the spatial gradients are computed using the Fourier-collocation method. This increases the accuracy of the gradient calculation and thus relaxes the requirement for dense computational grids compared to conventional finite difference methods. The accuracy and utility of the developed model is demonstrated via several numerical experiments, including the 3D simulation of the beam pattern from a clinical ultrasound probe.

Comparison of fundamental, second harmonic, and superharmonic imaging: a simulation study

In medical ultrasound, fundamental imaging (FI) uses the reflected echoes from the same spectral band as that of the emitted pulse. The transmission frequency determines the trade-off between penetration depth and spatial resolution. Tissue harmonic imaging (THI) employs the second harmonic of the emitted frequency band to construct images. Recently, superharmonic imaging (SHI) has been introduced, which uses the third to the fifth (super) harmonics. The harmonic level is determined by two competing phenomena: nonlinear propagation and frequency dependent attenuation. Thus, the transmission frequency yielding the optimal trade-off between the spatial resolution and the penetration depth differs for THI and SHI. This paper quantitatively compares the concepts of fundamental, second harmonic, and superharmonic echocardiography at their optimal transmission frequencies. Forward propagation is modeled using a 3D-KZK implementation and the iterative nonlinear contrast source (INCS) method. Backpropagation is assumed to be linear. Results show that the fundamental lateral beamwidth is the narrowest at focus, while the superharmonic one is narrower outside the focus. The lateral superharmonic roll-off exceeds the fundamental and second harmonic roll-off. Also, the axial resolution of SHI exceeds that of FI and THI. The far-field pulse-echo superharmonic pressure is lower than that of the fundamental and second harmonic. SHI appears suited for echocardiography and is expected to improve its image quality at the cost of a slight reduction in depth-of-field.

Verification of the Westervelt equation for focused transducers

This study investigates the validity of the Westervelt equation for focused transducers. The angular spectrum method is employed to analyze the second-harmonic acoustic field under the weakly nonlinear approximation. Although it is well known that the Westervelt equation is accurate for the case of quasi-plane waves, the present work demonstrates accurate solution for the highly focused case of a spherically-curved ultrasound transducer having an aperture angle of 80°. It is further found that the solution error is inversely dependent on the nonlinearity coefficient.

A contrast source method for nonlinear acoustic wave fields in media with spatially inhomogeneous attenuation

Experimental data reveals that attenuation is an important phenomenon in medical ultrasound. Attenuation is particularly important for medical applications based on nonlinear acoustics, since higher harmonics experience higher attenuation than the fundamental. Here, a method is presented to accurately solve the wave equation for nonlinear acoustic media with spatially inhomogeneous attenuation. Losses are modeled by a spatially dependent compliance relaxation function, which is included in the Westervelt equation. Introduction of absorption in the form of a causal relaxation function automatically results in the appearance of dispersion. The appearance of inhomogeneities implies the presence of a spatially inhomogeneous contrast source in the presented full-wave method leading to inclusion of forward and backward scattering. The contrast source problem is solved iteratively using a Neumann scheme, similar to the iterative nonlinear contrast source (INCS) method. The presented method is directionally independent and capable of dealing with weakly to moderately nonlinear, large scale, three-dimensional wave fields occurring in diagnostic ultrasound. Convergence of the method has been investigated and results for homogeneous, lossy, linear media show full agreement with the exact results. Moreover, the performance of the method is demonstrated through simulations involving steered and unsteered beams in nonlinear media with spatially homogeneous and inhomogeneous attenuation.

Optimization of a phased-array transducer for multiple harmonic imaging in medical applications: frequency and topology

Second-harmonic imaging is currently one of the standards in commercial echographic systems for diagnosis, because of its high spatial resolution and low sensitivity to clutter and near-field artifacts. The use of nonlinear phenomena mirrors is a great set of solutions to improve echographic image resolution. To further enhance the resolution and image quality, the combination of the 3rd to 5th harmonics--dubbed the superharmonics--could be used. However, this requires a bandwidth exceeding that of conventional transducers. A promising solution features a phased-array design with interleaved low- and high-frequency elements for transmission and reception, respectively. Because the amplitude of the backscattered higher harmonics at the transducer surface is relatively low, it is highly desirable to increase the sensitivity in reception. Therefore, we investigated the optimization of the number of elements in the receiving aperture as well as their arrangement (topology). A variety of configurations was considered, including one transmit element for each receive element (1/2) up to one transmit for 7 receive elements (1/8). The topologies are assessed based on the ratio of the harmonic peak pressures in the main and grating lobes. Further, the higher harmonic level is maximized by optimization of the center frequency of the transmitted pulse. The achievable SNR for a specific application is a compromise between the frequency-dependent attenuation and nonlinearity at a required penetration depth. To calculate the SNR of the complete imaging chain, we use an approach analogous to the sonar equation used in underwater acoustics. The generated harmonic pressure fields caused by nonlinear wave propagation were modeled with the iterative nonlinear contrast source (INCS) method, the KZK, or the Burger's equation. The optimal topology for superharmonic imaging was an interleaved design with 1 transmit element per 6 receive elements. It improves the SNR by ~5 dB compared with the interleaved (1/2) design reported in literature. The optimal transmit frequency for superharmonic echocardiography was found to be 1.0 to 1.2 MHz. For superharmonic abdominal imaging this frequency was found to be 1.7 to 1.9 MHz. For 2nd-harmonic echocardiography, the optimal transmit frequency of 1.8 MHz reported in the literature was corroborated with our simulation results.

Evaluation of a wave-vector-frequency-domain method for nonlinear wave propagation

A wave-vector-frequency-domain method is presented to describe one-directional forward or backward acoustic wave propagation in a nonlinear homogeneous medium. Starting from a frequency-domain representation of the second-order nonlinear acoustic wave equation, an implicit solution for the nonlinear term is proposed by employing the Green's function. Its approximation, which is more suitable for numerical implementation, is used. An error study is carried out to test the efficiency of the model by comparing the results with the Fubini solution. It is shown that the error grows as the propagation distance and step-size increase. However, for the specific case tested, even at a step size as large as one wavelength, sufficient accuracy for plane-wave propagation is observed. A two-dimensional steered transducer problem is explored to verify the nonlinear acoustic field directional independence of the model. A three-dimensional single-element transducer problem is solved to verify the forward model by comparing it with an existing nonlinear wave propagation code. Finally, backward-projection behavior is examined. The sound field over a plane in an absorptive medium is backward projected to the source and compared with the initial field, where good agreement is observed.