Modelling Elastic Wave Propagation Using the k-Wave MATLAB Toolbox

被引：108

作者：

Treeby, Bradley E. ^{[1
]}

Jaros, Jiri ^{[2
]}

Rohrbac, Daniel ^{[3
]}

Cox, B. T. ^{[1
]}

机构：

[1] UCL, Dept Med Phys & Biomed Engn, London, England

[2] Brno Univ Technol, Fac Informat Technol, CS-61090 Brno, Czech Republic

[3] Riverside Res, Lizzi Ctr Biomed Engn, New York, NY USA

来源：

2014 IEEE INTERNATIONAL ULTRASONICS SYMPOSIUM (IUS) | 2014年

基金：

英国工程与自然科学研究理事会;

关键词：

ULTRASOUND PROPAGATION; HETEROGENEOUS MEDIA; SIMULATION;

D O I：

10.1109/ULTSYM.2014.0037

中图分类号：

O42 [声学];

学科分类号：

070206 ; 082403 ;

摘要：

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.

引用

页码：146 / 149

页数：4

共 11 条

[1] Three-dimensional simulation of ultrasound propagation through trabecular bone structures measured by synchrotron microtomography [J].