Mapped Chebyshev pseudo-spectral method for simulating the shear wave propagation in the plane of symmetry of a transversely isotropic viscoelastic medium

Bo Qiang, John C. Brigham, Robert J. McGough, James F Greenleaf, Matthew W Urban

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

Shear wave elastography is a versatile technique that is being applied to many organs. However, in tissues that exhibit anisotropic material properties, special care must be taken to estimate shear wave propagation accurately and efficiently. A two-dimensional simulation method is implemented to simulate the shear wave propagation in the plane of symmetry in transversely isotropic viscoelastic media. The method uses a mapped Chebyshev pseudo-spectral method to calculate the spatial derivatives and an Adams–Bashforth–Moulton integrator with variable step sizes for time marching. The boundaries of the two-dimensional domain are surrounded by perfectly matched layers to approximate an infinite domain and minimize reflection errors. In an earlier work, we proposed a solution for estimating the apparent shear wave elasticity and viscosity of the spatial group velocity as a function of rotation angle through a low-frequency approximation by a Taylor expansion. With the solver implemented in MATLAB, the simulated results in this paper match well with the theory. Compared to the finite element method simulations we used before, the pseudo-spectral solver consumes less memory and is faster and achieves better accuracy.

Original languageEnglish (US)
Pages (from-to)1-13
Number of pages13
JournalMedical and Biological Engineering and Computing
DOIs
StateAccepted/In press - May 25 2016

Fingerprint

Shear waves
Wave propagation
MATLAB
Elasticity
Materials properties
Viscosity
Tissue
Derivatives
Finite element method
Data storage equipment

Keywords

  • Perfectly matched layer
  • Pseudo-spectral method
  • Shear wave
  • Transversely isotropic
  • Viscoelastic

ASJC Scopus subject areas

  • Biomedical Engineering
  • Computer Science Applications

Cite this

@article{40d4672b7f0246a2a18efcf407dde76c,
title = "Mapped Chebyshev pseudo-spectral method for simulating the shear wave propagation in the plane of symmetry of a transversely isotropic viscoelastic medium",
abstract = "Shear wave elastography is a versatile technique that is being applied to many organs. However, in tissues that exhibit anisotropic material properties, special care must be taken to estimate shear wave propagation accurately and efficiently. A two-dimensional simulation method is implemented to simulate the shear wave propagation in the plane of symmetry in transversely isotropic viscoelastic media. The method uses a mapped Chebyshev pseudo-spectral method to calculate the spatial derivatives and an Adams–Bashforth–Moulton integrator with variable step sizes for time marching. The boundaries of the two-dimensional domain are surrounded by perfectly matched layers to approximate an infinite domain and minimize reflection errors. In an earlier work, we proposed a solution for estimating the apparent shear wave elasticity and viscosity of the spatial group velocity as a function of rotation angle through a low-frequency approximation by a Taylor expansion. With the solver implemented in MATLAB, the simulated results in this paper match well with the theory. Compared to the finite element method simulations we used before, the pseudo-spectral solver consumes less memory and is faster and achieves better accuracy.",
keywords = "Perfectly matched layer, Pseudo-spectral method, Shear wave, Transversely isotropic, Viscoelastic",
author = "Bo Qiang and Brigham, {John C.} and McGough, {Robert J.} and Greenleaf, {James F} and Urban, {Matthew W}",
year = "2016",
month = "5",
day = "25",
doi = "10.1007/s11517-016-1522-9",
language = "English (US)",
pages = "1--13",
journal = "Medical and Biological Engineering and Computing",
issn = "0140-0118",
publisher = "Springer Verlag",

}

TY - JOUR

T1 - Mapped Chebyshev pseudo-spectral method for simulating the shear wave propagation in the plane of symmetry of a transversely isotropic viscoelastic medium

AU - Qiang, Bo

AU - Brigham, John C.

AU - McGough, Robert J.

AU - Greenleaf, James F

AU - Urban, Matthew W

PY - 2016/5/25

Y1 - 2016/5/25

N2 - Shear wave elastography is a versatile technique that is being applied to many organs. However, in tissues that exhibit anisotropic material properties, special care must be taken to estimate shear wave propagation accurately and efficiently. A two-dimensional simulation method is implemented to simulate the shear wave propagation in the plane of symmetry in transversely isotropic viscoelastic media. The method uses a mapped Chebyshev pseudo-spectral method to calculate the spatial derivatives and an Adams–Bashforth–Moulton integrator with variable step sizes for time marching. The boundaries of the two-dimensional domain are surrounded by perfectly matched layers to approximate an infinite domain and minimize reflection errors. In an earlier work, we proposed a solution for estimating the apparent shear wave elasticity and viscosity of the spatial group velocity as a function of rotation angle through a low-frequency approximation by a Taylor expansion. With the solver implemented in MATLAB, the simulated results in this paper match well with the theory. Compared to the finite element method simulations we used before, the pseudo-spectral solver consumes less memory and is faster and achieves better accuracy.

AB - Shear wave elastography is a versatile technique that is being applied to many organs. However, in tissues that exhibit anisotropic material properties, special care must be taken to estimate shear wave propagation accurately and efficiently. A two-dimensional simulation method is implemented to simulate the shear wave propagation in the plane of symmetry in transversely isotropic viscoelastic media. The method uses a mapped Chebyshev pseudo-spectral method to calculate the spatial derivatives and an Adams–Bashforth–Moulton integrator with variable step sizes for time marching. The boundaries of the two-dimensional domain are surrounded by perfectly matched layers to approximate an infinite domain and minimize reflection errors. In an earlier work, we proposed a solution for estimating the apparent shear wave elasticity and viscosity of the spatial group velocity as a function of rotation angle through a low-frequency approximation by a Taylor expansion. With the solver implemented in MATLAB, the simulated results in this paper match well with the theory. Compared to the finite element method simulations we used before, the pseudo-spectral solver consumes less memory and is faster and achieves better accuracy.

KW - Perfectly matched layer

KW - Pseudo-spectral method

KW - Shear wave

KW - Transversely isotropic

KW - Viscoelastic

UR - http://www.scopus.com/inward/record.url?scp=84971653769&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84971653769&partnerID=8YFLogxK

U2 - 10.1007/s11517-016-1522-9

DO - 10.1007/s11517-016-1522-9

M3 - Article

SP - 1

EP - 13

JO - Medical and Biological Engineering and Computing

JF - Medical and Biological Engineering and Computing

SN - 0140-0118

ER -