A heterogenous, time harmonic, nearly incompressible transverse isotropic finite element brain simulation platform for MR elastography.

In this study, we describe numerical implementation of a heterogenous, nearly incompressible, transverse isotropic (NITI) finite element (FE) model with key advantages for use in MR elastography of fibrous soft tissue. MR elastography (MRE) estimates heterogenous property distributions from MR-measured harmonic motion fields based on assumed mechanical models of tissue response. Current MRE property estimation methods usually assume isotropic properties, which cause inconsistencies arising from model-data mismatch when anisotropy is present. In this study, we use a NITI model parameterized by a base shear modulus, shear anisotropy, tensile anisotropy, and an isotropic bulk modulus, which describes the mechanical behavior of tissues with aligned fiber structures well. Property and fiber direction heterogeneity are implemented at the level of FE Gauss points, which allows high-resolution diffusion tensor imaging (DTI) data to be incorporated easily into the model. The resulting code was validated against analytical solutions and a commercial FEM package, and is suitable for incorporation into nonlinear inversion MRE algorithms. Simulations of MRE in brain tissue with heterogeneous properties and anisotropic fiber tracts, which produced wavefields similar to experimental MRE, were generated from anatomical, DTI and MRE image data, allowing investigation of MRE inversion performance in a realistic setting where the ground truth and underlying mechanical behavior are known. Two es...

Source: Physics in Medicine and Biology - June 7, 2020 Category: Physics Authors: McGarry MDJ, Van Houten E, Guertler C, Okamoto RJ, Smith DR, Sowinski DR, Johnson CL, Bayly P, Weaver J, Paulsen KD Tags: Phys Med Biol Source Type: research