A stabilized linear finite element method for anisotropic poroelastodynamics with application to cardiac perfusion

We propose a variational multiscale method stabilization of a linear finite element method for nonlinear poroelasticity. Our approach is suitable for the implicit time integration of poroelastic formulations in which the solid skeleton is anisotropic and incompressible. A detailed numerical methodology is presented for a monolithic formulation that includes both structural dynamics and Darcy flow. Our implementation of this methodology is verified using several hyperelastic and poroelastic benchmark cases, and excellent agreement is obtained with the literature. Grid convergence studies for both anisotropic hyperelastodynamics and poroelastodynamics demonstrate that the method is second-order accurate. The capabilities of our approach are demonstrated using a model of the left ventricle (LV) of the heart derived from human imaging data. Simulations using this model indicate that the anisotropicity of the myocardium has a substantial influence on the pore pressure. Furthermore, the temporal variations of the various components of the pore pressure (hydrostatic pressure and pressure resulting from changes in the volume of the pore fluid) are correlated with the variation of the added mass and dynamics of the LV, with maximum pore pressure being obtained at peak systole. The order of magnitude and the temporal variation of the pore pressure are in good agreement with the literature.