4D simulation of nonlinear pressure field propagation on GPU with the angular spectrum method

The simulation of nonlinear propagation of ultrasound wave in biological tissue is a very time consuming operation. Different simulators, based on finite difference or angular spectrum methods have been reported in the literature and the second one provide faster simulations by considering separately the different harmonics. In this paper we proposed to use a generalized angular spectrum method (GASM) to quickly compute the first three harmonics of the nonlinear field. Moreover, the possible inhomogeneity of the medium is considered in the mathematical background. The GASM has been implemented on a graphic processing unit (GPU) in order to improve the computation time. The obtained results show a good agreement between the GASM and another finite difference nonlinear simulator. In terms of computation time, the GPU implementation allows a speed-up factor comprised between 3.5 and 13 depending on the GPU uses.

[1]  O. Basset,et al.  Fundamental and second-harmonic ultrasound field computation of inhomogeneous nonlinear medium with a generalized angular spectrum method , 2011, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[2]  Gabriel Kiss,et al.  GPU volume rendering in 3D echocardiography: Real-time pre-processing and ray-casting , 2010, 2010 IEEE International Ultrasonics Symposium.

[3]  J. Jensen,et al.  Simulation of second harmonic ultrasound fields , 2010, 2010 IEEE International Ultrasonics Symposium.

[4]  Nassir Navab,et al.  Fast Ultrasound Image Simulation Using the Westervelt Equation , 2010, MICCAI.

[5]  C. J. Martin,et al.  Using GPUs for beamforming acceleration on SAFT imaging , 2009, 2009 IEEE International Ultrasonics Symposium.

[6]  Li-Wen Chang,et al.  GPU-based color Doppler ultrasound processing , 2009, 2009 IEEE International Ultrasonics Symposium.

[7]  T. Varslot,et al.  Computer simulation of forward wave propagation in soft tissue , 2005, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[8]  Steven G. Johnson,et al.  The Design and Implementation of FFTW3 , 2005, Proceedings of the IEEE.

[9]  M. Averkiou,et al.  A new imaging technique based on the nonlinear properties of tissues , 1997, 1997 IEEE Ultrasonics Symposium Proceedings. An International Symposium (Cat. No.97CH36118).

[10]  R. Cleveland,et al.  Time‐domain modeling of finite‐amplitude sound in relaxing fluids , 1995 .

[11]  S. Aanonsen,et al.  Distortion and harmonic generation in the nearfield of a finite amplitude sound beam , 1984 .

[12]  Xiang Yan,et al.  Angular Spectrum Decomposition Analysis of Second Harmonic Ultrasound Propagation and its Relation to Tissue Harmonic Imaging , 2006 .

[13]  M. Hamilton,et al.  Time‐domain modeling of pulsed finite‐amplitude sound beams , 1995 .

[14]  T. Szabo.  Generalized Fourier transform diffraction theory for parabolically anisotropic media , 1978 .