Scattering and Diffraction of Elastodynamic Waves in a Concentric Cylindrical Phantom for MR Elastography

Analytical solution for converging elliptic shear wave in a bounded transverse isotropic viscoelastic material with nonhomogeneous outer boundary

Dynamic elastography methods-based on optical, ultrasonic, or magnetic resonance imaging-are being developed for quantitatively mapping the shear viscoelastic properties of biological tissues, which are often altered by disease and injury. These diagnostic imaging methods involve analysis of shear wave motion in order to estimate or reconstruct the tissue's shear viscoelastic properties. Most reconstruction methods to date have assumed isotropic tissue properties. However, application to tissues like skeletal muscle and brain white matter with aligned fibrous structure resulting in local transverse isotropic mechanical properties would benefit from analysis that takes into consideration anisotropy. A theoretical approach is developed for the elliptic shear wave pattern observed in transverse isotropic materials subjected to axisymmetric excitation creating radially converging shear waves normal to the fiber axis. This approach, utilizing Mathieu functions, is enabled via a transformation to an elliptic coordinate system with isotropic properties and a ratio of minor and major axes matching the ratio of shear wavelengths perpendicular and parallel to the plane of isotropy in the transverse isotropic material. The approach is validated via numerical finite element analysis case studies. This strategy of coordinate transformation to equivalent isotropic systems could aid in analysis of other anisotropic tissue structures.

Mechanical analysis of an axially symmetric cylindrical phantom with a spherical heterogeneity for MR elastography

Cylindrical homogenous phantoms for magnetic resonance (MR) elastography in biomedical research provide one way to validate an imaging systems performance, but the simplified geometry and boundary conditions can cloak complexity arising at tissue interfaces. In an effort to develop a more realistic gel tissue phantom for MRE, we have constructed a heterogenous gel phantom (a sphere centrally embedded in a cylinder). The actuation comes from the phantom container, with the mechanical waves propagating toward the center, focusing the energy and thus allowing for the visualization of high-frequency waves that would otherwise be damped. The phantom was imaged and its stiffness determined using a 9.4 T horizontal MRI with a custom build piezo-elastic MRE actuator. The phantom was vibrated at three frequencies, 250, 500, and 750 Hz. The resulting shear wave images were first used to reconstruct material stiffness maps for thin (1 mm) axial slices at each frequency, from which the complex shear moduli μ were estimated, and then compared with forward modeling using a recently developed theoretical model which took μ as inputs. The overall accuracy of the measurement process was assessed by comparing theory with experiment for selected values of the shear modulus (real and imaginary parts). Close agreement is shown between the experimentally obtained and theoretically predicted wave fields.