I am interested in numerical evaluation of elastic parameter inversion algorithms for MR Elastography (MRE).
(In MRE the tissue is vibrated using an external actuator and MRI is used to image the resulting microscopic 3D displacements throughout the field of view. The curl of the data is used to suppress the effects of the compressional wave and shear is thus the mode of interest)
Essentially I want to simulated (shear) wave propagation in an anisotropic (e.g. transversely isotropic or orthotropic) medium. It is often assumed that, following the application of external shear waves, the medium is harmonically oscillating in steady state. Given that simplification, can this toolbox be of use?
Instead of initial pressures could one specify initial stiffness tensors/parameters?
Thanks,
Kevin