This is a MATLAB implementation of a 2D pseudo-spectral,
time domain solution of the acoustic wave equation described in the literature
(reference below). It uses pseudo-spectral methods to calculate spatial
derivatives and a staggered Adams-Bashforth method to integrate forward in
time. A perfectly-matched boundary layer is applied at the edges of the
calculation domain. This model can include nonlinear propagation and frequency
dependent attenuation. This model was developed for research and instruction
in the acoustics of medical ultrasound.
|
Author: |
|
|
Organization: |
Philips Medical Systems, Ultrasound R&D (code developed at Duke University) |
|
Date posted: |
November 1, 2000 |
|
Code type / platform: |
|
|
Use restrictions: |
None. (Acknowledgement requested.) |
|
Download format(s): |
Zip archive: |
The core MATLAB function solve_sps2d.m is provided along with several examples in which this function is called to simulate various acoustic conditions, as described in the READ_ME_sps2d.txt file.
The author welcomes bug reports and suggestions for refinements.
Wojcik G. L., B. Fornberg, R. Waag, L. Carcione, J. Mould, L. Nikodym, and T. Driscoll, "Pseudospectral methods for large-scale bioacoustic models," Proc. IEEE-UFFC Ultrasonics Symposium, pp. 1501-1506, 1997.
G. L. Wojcik, J. Mould, S. Ayter and L. M. Carcione, "A study of second harmonic generation by focused medical transducer pulses," Proc. IEEE-UFFC Ultrasonics Symposium, pp. 1583-1588, 1998.
Nachman, A. I., J. F. Smith, and R. C. Waag, "An equation for acoustic propagation in inhomogeneous media with relaxation losses," JASA, 88, pp. 1584-1595, 1990.
Yuan, X., D. Borup, J. W. Wiskin, M. Berggren, R. Eidens, and S. Johnson, "Formulation and validation of Berenger's PML absorbing boundary for the FDTD simulation of acoustic scattering," IEEE Trans. UFFC, 44, 816-822, 1997.
Yuan, X., D. Borup, J. W. Wiskin, M. Berggren, and S. Johnson, "Simulation of acoustic wave propagation in dispersive media with relaxation losses by using FDTD method with PML absorbing boundary conditions," IEEE Trans. UFFC, 46, 14-23, 1999.
Ghrist, M., B. Fornberg, and T. A. Driscoll, "Staggered time integrators for wave equations," SIAM Journal on Numerical Analysis, 38(3), pp. 718-741, 2000.
Simulation of broadband focused ultrasound pulse emitted from a phased
array:
Simulation of a planar ultrasound wavefront impinging on a circular inclusion of material having a slower sound speed:
