A 0D model of the Heart

This model is a simplified "0D" version of the 3D electromechanical cardiac model described in [Marchesseau, S. et al.] (1) and developed at Inria Sophia Antipolis - Asclepios Team (2).

It describes the mechanical behavior of a single heart ventricle, based on the Bestel-Clement-Sorine model of myocardial tissue as described in [Chapelle et al.] (3) and the Mooney-Rivlin (4) strain energy for passive elasticity. Hemodynamics are represented through global values of pressures and flows in the ventricular chambers. The mechanical equations are coupled to the 3-element Windkessel model (5) for the after-load. Finally, cardiac pre-load is modeled by a (prescribed) varying atrial pressure with a "pressure bump" at the beginning of the cardiac cycle, to account for atrial contraction.

Some features related to the dynamics of a heart valve are currently relying on an unorthodox CellML (6) implementation (HOD.cellml). Thus, to run this model in OpenCOR (7), specific ODE solver settings are needed:

  • Name: Euler (forward); and
  • Step: 1e-05 second.

For the simulation itself, we recommend using the following settings, which are used to compute 3 heartbeats:

  • Starting point: 0 second;
  • Ending point: 2.4 second; and
  • Point interval: 0.001 second.

Finally, for the output, we recommend creating 3 graph panels to, respectively, plot the following graphs:

  • Volume curve: V_mL against the variable of integration;
  • Ventricular, atrial and aortic pressure curves: P_v, P_at and P_ar against the variable of integration; and
  • Pressure-volume curve: P_v against V_mL.

This exposure also comes with both a SED-ML (7) file (HOD.sedml) and a COMBINE archive (8) file (HOD.omex). They contain exactly this configuration and can be readily opened and run using OpenCOR 0.5 (available at http://www.opencor.ws/downloads).

Contact: Roch Molléro at rochmollero@hotmail.com.

Source
Derived from workspace A 0D model of the Heart at changeset b4f1199958a2.
Collaboration
To begin collaborating on this work, please use your git client and issue this command: