The development and optimization of medical ultrasound transducers and imaging modalities require a computational method that accurately predicts the nonlinear acoustic pressure field. A prospective method should provide the wide-angle, pulsed field emitted by an arbitrary planar source distribution and propagating in a three-dimensional, large scale domain holding a nonlinear acoustic medium. In this paper, a method is presented that is free of any assumed wavefield directionality. The nonlinear acoustic wave equation is solved by treating the nonlinear term as a contrast source. This formulation leads to an iterative scheme that involves the repetitive solution of a linear wave problem through Green's function method. It is shown that accurate field predictions may be obtained within a few iterations. Moreover, by employing a dedicated numerical convolution technique, the method allows for a discretization down to two points per wavelength or period of the highest frequency of interest. The performance of the method is evaluated through a number of nonlinear field predictions for pulsed transducers with various geometries. The results demonstrate the directional independence of the method. Moreover, comparison with results from several existing methods shows that the method accurately predicts the nonlinear field for weak to moderate nonlinearity.

The full-wave computation of transient acoustic fields with sizes in the order of 100ϫ 100ϫ 100 wavelengths by 100 periods requires a numerical method that is extremely efficient in terms of storage and computation. Iterative integral equation methods offer a good performance on these points, provided that the recurring spatiotemporal convolutions are computed with a coarse sampling and relatively few computational operations. This paper describes a method for the numerical evaluation of very large-scale, four-dimensional convolutions that employs a fast Fourier transformation and that uses a sampling rate close to or at the limit of two points per wavelength and per period. To achieve this, the functions involved are systematically filtered, windowed, and zero-padded with respect to all relevant coordinates prior to sampling. The method is developed in the context of the Neumann iterative solution of the acoustic contrast source problem for an inhomogeneous medium. The implementation of the method on a parallel computer is discussed. The obtained numerical results have a relative root mean square error of a few percent when sampling at two points per wavelength and per period. Further, the results prove that the method enables the computation of transient fields in the order of the indicated size.

No abstract

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.

No abstract