### Abstract

Modeling of the heart ventricles is one of the most challenging tasks in soft tissue mechanics because cardiac tissue is a strongly anisotropic incompressible material with an active component of stress. In most current approaches with active force, the number of degrees of freedom (DOF) is limited by the direct method of solution of linear systems of equations. We develop a new approach for high-resolution heart models with large numbers of DOF by: (1) developing a hex-dominant finite element mixed formulation and (2) developing a Krylov subspace iterative method that is able to solve the system of linearized equations for saddle-point problems with active stress. In our approach, passive cardiac tissue is modeled as a hyperelastic, incompressible material with orthotropic properties, and mixed pressure-displacement finite elements are used to enforce incompressibility. Active stress is generated by a model with force dependence on length and velocity of muscle shortening. The ventricles are coupled to a lumped circulatory model. For efficient solution of linear systems, we use Flexible GMRES with a nonlinear preconditioner based on block matrix decomposition involving the Schur complement. Three methods for approximating the inverse of the Schur complement are evaluated: inverse of the pressure mass matrix; least squares commutators; and sparse approximate inverse. The sub-matrix corresponding to the displacement variables is preconditioned by a V-cycle of hybrid geometric–algebraic multigrid followed by correction with several iterations of GMRES preconditioned by sparse approximate inverse. The overall solver is demonstrated on a high-resolution two ventricle mesh based on a human anatomy with roughly 130 K elements and 1.7 M displacement DOF. Effectiveness of the numerical method for active contraction is shown. To the best of our knowledge, this solver is the first to efficiently model ventricular contraction using an iterative linear solver for the mesh size demonstrated and therefore opens the possibility for future very high-resolution models. In addition, several relatively simple benchmark problems are designed for a verification exercise to show that the solver is functioning properly and correctly solves the underlying mathematical model. Here, the output of the newly designed solver is compared to that of the mechanics component of Chaste (‘Cancer, Heart and Soft Tissue Environment’). These benchmark tests may be used by other researchers to verify their newly developed methods and codes.

Original language | English (US) |
---|---|

Pages (from-to) | 829-849 |

Number of pages | 21 |

Journal | Biomechanics and Modeling in Mechanobiology |

Volume | 14 |

Issue number | 4 |

DOIs | |

State | Published - Aug 4 2015 |

Externally published | Yes |

### Fingerprint

### Keywords

- Cardiac mechanics
- Parallel computations
- Soft tissue mechanics
- Ventricular model

### ASJC Scopus subject areas

- Biotechnology
- Mechanical Engineering
- Modeling and Simulation

### Cite this

*Biomechanics and Modeling in Mechanobiology*,

*14*(4), 829-849. https://doi.org/10.1007/s10237-014-0639-8

**A high-resolution computational model of the deforming human heart.** / Gurev, Viatcheslav; Pathmanathan, Pras; Fattebert, Jean Luc; Wen, Hui Fang; Magerlein, John; Gray, Richard A.; Richards, David F.; Rice, J. Jeremy.

Research output: Contribution to journal › Article

*Biomechanics and Modeling in Mechanobiology*, vol. 14, no. 4, pp. 829-849. https://doi.org/10.1007/s10237-014-0639-8

}

TY - JOUR

T1 - A high-resolution computational model of the deforming human heart

AU - Gurev, Viatcheslav

AU - Pathmanathan, Pras

AU - Fattebert, Jean Luc

AU - Wen, Hui Fang

AU - Magerlein, John

AU - Gray, Richard A.

AU - Richards, David F.

AU - Rice, J. Jeremy

PY - 2015/8/4

Y1 - 2015/8/4

N2 - Modeling of the heart ventricles is one of the most challenging tasks in soft tissue mechanics because cardiac tissue is a strongly anisotropic incompressible material with an active component of stress. In most current approaches with active force, the number of degrees of freedom (DOF) is limited by the direct method of solution of linear systems of equations. We develop a new approach for high-resolution heart models with large numbers of DOF by: (1) developing a hex-dominant finite element mixed formulation and (2) developing a Krylov subspace iterative method that is able to solve the system of linearized equations for saddle-point problems with active stress. In our approach, passive cardiac tissue is modeled as a hyperelastic, incompressible material with orthotropic properties, and mixed pressure-displacement finite elements are used to enforce incompressibility. Active stress is generated by a model with force dependence on length and velocity of muscle shortening. The ventricles are coupled to a lumped circulatory model. For efficient solution of linear systems, we use Flexible GMRES with a nonlinear preconditioner based on block matrix decomposition involving the Schur complement. Three methods for approximating the inverse of the Schur complement are evaluated: inverse of the pressure mass matrix; least squares commutators; and sparse approximate inverse. The sub-matrix corresponding to the displacement variables is preconditioned by a V-cycle of hybrid geometric–algebraic multigrid followed by correction with several iterations of GMRES preconditioned by sparse approximate inverse. The overall solver is demonstrated on a high-resolution two ventricle mesh based on a human anatomy with roughly 130 K elements and 1.7 M displacement DOF. Effectiveness of the numerical method for active contraction is shown. To the best of our knowledge, this solver is the first to efficiently model ventricular contraction using an iterative linear solver for the mesh size demonstrated and therefore opens the possibility for future very high-resolution models. In addition, several relatively simple benchmark problems are designed for a verification exercise to show that the solver is functioning properly and correctly solves the underlying mathematical model. Here, the output of the newly designed solver is compared to that of the mechanics component of Chaste (‘Cancer, Heart and Soft Tissue Environment’). These benchmark tests may be used by other researchers to verify their newly developed methods and codes.

AB - Modeling of the heart ventricles is one of the most challenging tasks in soft tissue mechanics because cardiac tissue is a strongly anisotropic incompressible material with an active component of stress. In most current approaches with active force, the number of degrees of freedom (DOF) is limited by the direct method of solution of linear systems of equations. We develop a new approach for high-resolution heart models with large numbers of DOF by: (1) developing a hex-dominant finite element mixed formulation and (2) developing a Krylov subspace iterative method that is able to solve the system of linearized equations for saddle-point problems with active stress. In our approach, passive cardiac tissue is modeled as a hyperelastic, incompressible material with orthotropic properties, and mixed pressure-displacement finite elements are used to enforce incompressibility. Active stress is generated by a model with force dependence on length and velocity of muscle shortening. The ventricles are coupled to a lumped circulatory model. For efficient solution of linear systems, we use Flexible GMRES with a nonlinear preconditioner based on block matrix decomposition involving the Schur complement. Three methods for approximating the inverse of the Schur complement are evaluated: inverse of the pressure mass matrix; least squares commutators; and sparse approximate inverse. The sub-matrix corresponding to the displacement variables is preconditioned by a V-cycle of hybrid geometric–algebraic multigrid followed by correction with several iterations of GMRES preconditioned by sparse approximate inverse. The overall solver is demonstrated on a high-resolution two ventricle mesh based on a human anatomy with roughly 130 K elements and 1.7 M displacement DOF. Effectiveness of the numerical method for active contraction is shown. To the best of our knowledge, this solver is the first to efficiently model ventricular contraction using an iterative linear solver for the mesh size demonstrated and therefore opens the possibility for future very high-resolution models. In addition, several relatively simple benchmark problems are designed for a verification exercise to show that the solver is functioning properly and correctly solves the underlying mathematical model. Here, the output of the newly designed solver is compared to that of the mechanics component of Chaste (‘Cancer, Heart and Soft Tissue Environment’). These benchmark tests may be used by other researchers to verify their newly developed methods and codes.

KW - Cardiac mechanics

KW - Parallel computations

KW - Soft tissue mechanics

KW - Ventricular model

UR - http://www.scopus.com/inward/record.url?scp=84939418056&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84939418056&partnerID=8YFLogxK

U2 - 10.1007/s10237-014-0639-8

DO - 10.1007/s10237-014-0639-8

M3 - Article

C2 - 25567753

AN - SCOPUS:84939418056

VL - 14

SP - 829

EP - 849

JO - Biomechanics and Modeling in Mechanobiology

JF - Biomechanics and Modeling in Mechanobiology

SN - 1617-7959

IS - 4

ER -