Numerical simulations of magnetic resonance elastography using finite element analysis with a linear heterogeneous viscoelastic model

Magnetic resonance elastography (MRE) is a technique to identify the viscoelastic moduli of biological tissues by solving the inverse problem from the displacement field of viscoelastic wave propagation in a tissue measured by MRI. Because finite element analysis (FEA) of MRE evaluates not only the viscoelastic model for a tissue but also the efficiency of the inversion algorithm, we developed FEA for MRE using commercial software called ANSYS, the Zener model for displacement field of a wave inside tissue, and an inversion algorithm called the modified integral method. The profile of the simulated displacement field by FEA agrees well with the experimental data measured by MRE for gel phantoms. Similarly, the value of storage modulus (i.e., stiffness) recovered using the modified integral method with the simulation data is consistent with the value given in FEA. Furthermore, applying the suggested FEA to a human liver demonstrates the effectiveness of the present simulation scheme.