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.
Similar content being viewed by others
References
Alexandrescu A, Bueno-Orovio A, Salgueiro JR, Perez-Garcia VM (2009) Mapped Chebyshev pseudospectral method for the study of multiple scale phenomena. Comput Phys Commun 180:912–919. doi:10.1016/j.cpc.2008.12.018
Amador C, Urban MW, Chen S, Greenleaf JF (2011) Shearwave dispersion ultrasound vibrometry (SDUV) on swine kidney. IEEE Trans Ultrason Ferroelectr Freq Control 58:2608–2619
Aristizabal S, Amador C, Qiang B, Kinnick RR, Nenadic IZ, Greenleaf JF, Urban MW (2014) Shear wave vibrometry evaluation in transverse isotropic tissue mimicking phantoms and skeletal muscle. Phys Med Biol 59:7735–7752. doi:10.1088/0031-9155/59/24/7735
Bammer R, Acar B, Moseley ME (2003) In vivo MR tractography using diffusion imaging. Eur J Radiol 45:223–234
Basser PJ, Mattiello J, LeBihan D (1994) MR diffusion tensor spectroscopy and imaging. Biophys J 66:259–267. doi:10.1016/S0006-3495(94)80775-1
Basser PJ, Pierpaoli C (1996) Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. J Magn Reson, Ser B 111:209–219. doi:10.1006/jmrb.1996.0086
Bercoff J, Tanter M, Fink M (2004) Supersonic shear imaging: a new technique for soft tissue elasticity mapping. IEEE Trans Ultrason Ferroelectr Freq Control 51:396–409
Boyd JP (2001) Chebyshev and Fourier spectral methods, 2nd edn. Dover Publications, Mineola
Brum J, Bernal M, Gennisson JL, Tanter M (2014) In vivo evaluation of the elastic anisotropy of the human Achilles tendon using shear wave dispersion analysis. Phys Med Biol 59:505–523
Budzik JF, Balbi V, Verclytte S, Pansini V, Le Thuc V, Cotten A (2014) Diffusion tensor imaging in musculoskeletal disorders. Radiographics 34:E56–E72. doi:10.1148/rg.343125062
Carcione JM (1993) Seismic modeling in viscoelastic media. Geophysics 58:110–120. doi:10.1190/1.1443340
Carcione JM (2007) Wave fields in real media: wave propagation in anisotropic, anelastic, porous and electromagnetic media. Handbook of geophysical exploration Seismic exploration, vol 38, 2nd edn. Elsevier, Amsterdam
Carcione JM, Kosloff D, Kosloff R (1988) Wave propagation simulation in a linear viscoelastic medium. Geophys J 95:597–611
Carcione JM, Poletto F, Gei D (2004) 3-D wave simulation in anelastic media using the Kelvin–Voigt constitutive equation. J Comput Phys 196:282–297. doi:10.1016/j.jcp.2003.10.024
Catheline S, Gennisson JL, Delon G, Fink M, Sinkus R, Abouelkaram S, Culioli J (2004) Measurement of viscoelastic properties of homogeneous soft solid using transient elastography: an inverse problem approach. J Acoust Soc Am 116:3734–3741
Chatelin S, Bernal M, Deffieux T, Papadacci C, Flaud P, Nahas A, Boccara C, Gennisson JL, Tanter M, Pernot M (2014) Anisotropic polyvinyl alcohol hydrogel phantom for shear wave elastography in fibrous biological soft tissue: a multimodality characterization. Phys Med Biol 59:6923–6940. doi:10.1088/0031-9155/59/22/6923
Chui C, Kobayashi E, Chen X, Hisada T, Sakuma I (2007) Transversely isotropic properties of porcine liver tissue: experiments and constitutive modelling. Med Biol Eng Comput 45:99–106
Corovic S, Zupanic A, Kranjc S, Al Sakere B, Leroy-Willig A, Mir LM, Miklavcic D (2010) The influence of skeletal muscle anisotropy on electroporation: in vivo study and numerical modeling. Med Biol Eng Comput 48:637–648
Eby SF, Song P, Chen S, Chen Q, Greenleaf JF, An KN (2013) Validation of shear wave elastography in skeletal muscle. J Biomech 46:2381–2387. doi:10.1016/j.jbiomech.2013.07.033
Galban CJ, Maderwald S, Stock F, Ladd ME (2007) Age-related changes in skeletal muscle as detected by diffusion tensor magnetic resonance Imaging. J Gerontol Ser A Biol Sci Med Sci 62:453–458
Garmirian LP, Chin AB, Rutkove SB (2009) Discriminating neurogenic from myopathic disease via measurement of muscle anisotropy. Muscle Nerve 39:16–24. doi:10.1002/mus.21115
Gates F, McCammond D, Zingg W, Kunov H (1980) In vivo stiffness properties of the canine diaphragm muscle. Med Biol Eng Comput 18:625–632
Gennisson JL, Catheline S, Chaffai S, Fink M (2003) Transient elastography in anisotropic medium: application to the measurement of slow and fast shear wave speeds in muscles. J Acoust Soc Am 114:536–541
Gennisson JL, Deffieux T, Mace E, Montaldo G, Fink M, Tanter M (2010) Viscoelastic and anisotropic mechanical properties of in vivo muscle tissue assessed by supersonic shear imaging. Ultrasound Med Biol 36:789–801. doi:10.1016/j.ultrasmedbio.2010.02.013
Gennisson JL, Grenier N, Combe C, Tanter M (2012) Supersonic shear wave elastography of in vivo pig kidney: influence of blood pressure, urinary pressure and tissue anisotropy. Ultrasound Med Biol 38:1559–1567. doi:10.1016/j.ultrasmedbio.2012.04.013
Guillard H, Male JM, Peyret R (1992) Adaptive spectral methods with application to mixing layer computations. J Comput Phys 102:114–127. doi:10.1016/S0021-9991(05)80010-5
Khalil C, Budzik JF, Kermarrec E, Balbi V, Le Thuc V, Cotten A (2010) Tractography of peripheral nerves and skeletal muscles. Eur J Radiol 76:391–397. doi:10.1016/j.ejrad.2010.03.012
Komatitsch D, Martin R (2007) An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation. Geophysics Sm72:Sm155–Sm167. doi:10.1190/1.2757586
Kosloff D, Talezer H (1993) A modified Chebyshev pseudospectral method with an O(N-1) time step restriction. J Comput Phys 104:457–469. doi:10.1006/jcph.1993.1044
Le Bihan D, Mangin JF, Poupon C, Clark CA, Pappata S, Molko N, Chabriat H (2001) Diffusion tensor imaging: concepts and applications. J Magn Reson Imaging 13:534–546
Lee W-N, Pernot M, Couade M, Messas E, Bruneval P, Bel A, Hagege AA, Fink M, Tanter M (2012) Mapping myocardial fiber orientation using echocardiography-based shear wave imaging. IEEE Trans Med Imaging 31:554–562
Lee WN, Larrat B, Pernot M, Tanter M (2012) Ultrasound elastic tensor imaging: comparison with MR diffusion tensor imaging in the myocardium. Phys Med Biol 57:5075–5095. doi:10.1088/0031-9155/57/16/5075
Li Y, Bou Matar O (2010) Convolutional perfectly matched layer for elastic second-order wave equation. J Acoust Soc Am 127:1318–1327. doi:10.1121/1.3290999
Martin R, Komatitsch D (2009) An unsplit convolutional perfectly matched layer technique improved at grazing incidence for the viscoelastic wave equation. Geophys J Int 179:333–344. doi:10.1111/j.1365-246X.2009.04278.x
Martin R, Komatitsch D, Gedney SD, Bruthiaux E (2010) A high-order time and space formulation of the unsplit perfectly matched layer for the seismic wave equation using auxiliary differential equations (ADE-PML). Comp Model Eng 56:17–41
Masutani Y, Aoki S, Abe O, Hayashi N, Otomo K (2003) MR diffusion tensor imaging: recent advance and new techniques for diffusion tensor visualization. Eur J Radiol 46:53–66
MATLAB manual (2013). Mathworks Inc.
McGarry MDJ, Van Houten EEW (2008) Use of a Rayleigh damping model in elastography. Med Biol Eng Comput 46:759–766
Moseley ME, Cohen Y, Kucharczyk J, Mintorovitch J, Asgari HS, Wendland MF, Tsuruda J, Norman D (1990) Diffusion-weighted MR imaging of anisotropic water diffusion in cat central nervous system. Radiology 176:439–445. doi:10.1148/radiology.176.2.2367658
Muthupillai R, Ehman RL (1996) Magnetic resonance elastography. Nat Med 2:601–603
Ng KT, Yan R (2003) Three-dimensional pseudospectral modelling of cardiac propagation in an inhomogeneous anisotropic tissue. Med Biol Eng Comput 41:618–624
Palmeri ML, Sharma AC, Bouchard RR, Nightingale RW, Nightingale KR (2005) A finite-element method model of soft tissue response to impulsive acoustic radiation force. IEEE Trans Ultrason Ferroelectr Freq Control 52:1699–1712
Palmeri ML, Wang MH, Dahl JJ, Frinkley KD, Nightingale KR (2008) Quantifying hepatic shear modulus in vivo using acoustic radiation force. Ultrasound Med Biol 34:546–558. doi:10.1016/j.ultrasmedbio.2007.10.009
Pierpaoli C, Basser PJ (1996) Toward a quantitative assessment of diffusion anisotropy. Magn Reson Med 36:893–906. doi:10.1002/mrm.1910360612
Qi J, Olsen NJ, Price RR, Winston JA, Park JH (2008) Diffusion-weighted imaging of inflammatory myopathies: polymyositis and dermatomyositis. J Magn Reson Imaging 27:212–217. doi:10.1002/jmri.21209
Qiang B, Brigham JC, Aristizabal S, Greenleaf JF, Zhang X, Urban MW (2015) Modeling transversely isotropic, viscoelastic, incompressible tissue-like materials with application in ultrasound shear wave elastography. Phys Med Biol 60:1289–1306. doi:10.1088/0031-9155/60/3/1289
Rudenko OV, Sarvazyan AP (2014) Wave anisotropy of shear viscosity and elasticity. Acoust Phys 60:710–718. doi:10.1134/s1063771014060141
Sansalone M, Carino NJ, Hsu NN (1987) A finite-element study of transient wave-propagation in plates. J Res Nat Bur Stand 92:267–278. doi:10.6028/Jres.092.025
Saotome T, Sekino M, Eto F, Ueno S (2006) Evaluation of diffusional anisotropy and microscopic structure in skeletal muscles using magnetic resonance. Magn Reson Imaging 24:19–25. doi:10.1016/j.mri.2005.09.009
Sarvazyan AP, Rudenko OV, Swanson SD, Fowlkes JB, Emelianov SY (1998) Shear wave elasticity imaging: a new ultrasonic technology of medical diagnostics. Ultrasound Med Biol 24:1419–1435. doi:10.1016/S0301-5629(98)00110-0
Trefethen LN (2000) Spectral methods in MATLAB. Software, environments, tools. Society for Industrial and Applied Mathematics, Philadelphia, PA
Urban MW, Chen S, Fatemi M (2012) A review of shearwave dispersion ultrasound vibrometry (SDUV) and its applications. Curr Med Imaging Rev 8:27–36
Virieux J, Calandra H, Plessix RE (2011) A review of the spectral, pseudo-spectral, finite-difference and finite-element modelling techniques for geophysical imaging. Geophys Prospect 59:794–813. doi:10.1111/j.1365-2478.2011.00967.x
Wang M, Byram B, Palmeri M, Rouze N, Nightingale K (2013) Imaging transverse isotropic properties of muscle by monitoring acoustic radiation force induced shear waves using a 2-D matrix ultrasound array. IEEE Trans Med Imaging 32:1671–1684. doi:10.1109/TMI.2013.2262948
Zhan Z, Ng KT (2000) Two-dimensional Chebyshev pseudospectral modelling of cardiac propagation. Med Biol Eng Comput 38:311–318
Zhang J, Zhang G, Morrison B, Mori S, Sheikh KA (2008) Magnetic resonance imaging of mouse skeletal muscle to measure denervation atrophy. Exp Neurol 212:448–457. doi:10.1016/j.expneurol.2008.04.033
Zhao H, Song P, Urban MW, Greenleaf JF, Chen S (2012) Shear wave speed measurement using an unfocused ultrasound beam. Ultrasound Med Biol 38:1646–1655
Acknowledgments
This work was supported by Grant R01DK092255 from the National Institute of Diabetes and digestive and Kidney Diseases (NIDDK) and National Institutes of Health (NIH). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIDDK and NIH.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Qiang, B., Brigham, J.C., McGough, R.J. et al. Mapped Chebyshev pseudo-spectral method for simulating the shear wave propagation in the plane of symmetry of a transversely isotropic viscoelastic medium. Med Biol Eng Comput 55, 389–401 (2017). https://doi.org/10.1007/s11517-016-1522-9
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11517-016-1522-9