Přístupnostní navigace
E-application
Search Search Close
Publication detail
TREEBY, B. JAROŠ, J. ROHRBACH, D. COX, B.
Original Title
Modelling Elastic Wave Propagation Using the k-Wave MATLAB Toolbox
Type
conference paper
Language
English
Original Abstract
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.
Keywords
elastic wave, k-Wave, Matlab
Authors
TREEBY, B.; JAROŠ, J.; ROHRBACH, D.; COX, B.
RIV year
2014
Released
31. 8. 2014
Publisher
Institute of Electrical and Electronics Engineers
Location
Chicago, IL
ISBN
978-1-4799-7049-0
Book
IEEE International Ultrasonics Symposium, IUS
Pages from
146
Pages to
149
Pages count
4
URL
http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=6931789&tag=1
BibTex
@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" }