A 2D finite element model for shear wave propagation in biological soft tissues: Application to magnetic resonance elastography

In silico analysis reveals the prospects of renal anisotropy in improving chronic kidney disease detection using ultrasound shear wave elastography

Renal anisotropy is a complex property of the kidney and often poses a challenge in obtaining consistent measurements when using shear wave elastography to detect chronic kidney disease. To circumvent the challenge posed by renal anisotropy in clinical settings, a dimensionless biomarker termed the 'anisotropic ratio' was introduced to establish a correlation between changes in degree of renal anisotropy and progression of chronic kidney disease through an in silico perspective. To achieve this, an efficient model reduction approach was developed to model the anisotropic property of kidneys. Good agreement between the numerical and experimental data were obtained, as percentage errors of less than 5.5% were reported when compared against experimental phantom measurement from the literature. To demonstrate the applicability of the model to clinical measurements, the anisotropic ratio of sheep kidneys was quantified, with both numerical and derived experimental results reporting a value of .667. Analysis of the anisotropic ratio with progression of chronic kidney disease demonstrated that patients with normal kidneys would have a lower anisotropic ratio of .872 as opposed to patients suffering from renal impairment, in which the anisotropic ratio may increase to .904, as determined from this study. The findings demonstrate the potential of the anisotropic ratio in improving the detection of chronic kidney disease using shear wave elastography.

Numerical simulation of wave propagation through interfaces using the extended finite element method for magnetic resonance elastography

Magnetic resonance elastography (MRE) is an elasticity imaging technique for quantitatively assessing the stiffness of human tissues. In MRE, finite element method (FEM) is widely used for modeling wave propagation and stiffness reconstruction. However, in front of inclusions with complex interfaces, FEM can become burdensome in terms of the model partition and computationally expensive. In this work, we implement a formulation of FEM, known as the eXtended finite element method (XFEM), which is a method used for modeling discontinuity like crack and heterogeneity. Using a level-set method, it makes the interface independent of the mesh, thus relieving the meshing efforts. We investigate this method in two studies: wave propagation across an oblique linear interface and stiffness reconstruction of a random-shape inclusion. In the first study, numerical results by XFEM and FEM models revealing the wave conversion rules at linear interface are presented and successfully compared to the theoretical predictions. The second study, investigated in a pseudo-practical application, demonstrates further the applicability of XFEM in MRE and the convenience, accuracy, and speed of XFEM with respect to FEM. XFEM can be regarded as a promising alternative to FEM for inclusion modeling in MRE.