Publication Details
Modelling Elastic Wave Propagation Using the k-Wave MATLAB Toolbox
elastic wave, k-Wave, Matlab
A new model for simulating elastic wave propagation using the open-source k-Wave MATLAB Toolbox is described. The model is based on two coupled first-order equations describing the stress and particle velocity within an isotropic medium. For absorbing media, the Kelvin-Voigt model of viscoelasticity is used. The equations are discretised in 2D and 3D using an efficient time-stepping pseudospectral scheme. This uses the Fourier collocation spectral method to compute spatial derivatives and a leapfrog finite-difference scheme to integrate forwards in time. A multi-axial perfectly matched layer (M-PML) is implemented to allow free-field simulations using a finite-sized computational grid. Acceleration using a graphics processing unit (GPU) is supported via the MATLAB Parallel Computing Toolbox. An overview of the simulation functions and their theoretical and numerical foundations is described.
@inproceedings{BUT111641,
author="Bradley {Treeby} and Jiří {Jaroš} and Daniel {Rohrbach} and Ben {Cox}",
title="Modelling Elastic Wave Propagation Using the k-Wave MATLAB Toolbox",
booktitle="IEEE International Ultrasonics Symposium, IUS",
year="2014",
pages="146--149",
publisher="Institute of Electrical and Electronics Engineers",
address="Chicago, IL",
doi="10.1109/ULTSYM.2014.0037",
isbn="978-1-4799-7049-0",
url="http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=6931789&tag=1"
}