Verification of cardiac mechanics software: benchmark problems and solutions for testing active and passive material behaviour

Models of cardiac mechanics are increasingly used to investigate cardiac physiology. These models are characterized by a high level of complexity, including the particular anisotropic material properties of biological tissue and the actively contracting material. A large number of independent simulation codes have been developed, but a consistent way of verifying the accuracy and replicability of simulations is lacking. To aid in the verification of current and future cardiac mechanics solvers, this study provides three benchmark problems for cardiac mechanics. These benchmark problems test the ability to accurately simulate pressure-type forces that depend on the deformed objects geometry, anisotropic and spatially varying material properties similar to those seen in the left ventricle and active contractile forces. The benchmark was solved by 11 different groups to generate consensus solutions, with typical differences in higher-resolution solutions at approximately 0.5%, and consistent results between linear, quadratic and cubic finite elements as well as different approaches to simulating incompressible materials. Online tools and solutions are made available to allow these tests to be effectively used in verification of future cardiac mechanics software.

[1]  Pablo Lamata,et al.  Improving the Stability of Cardiac Mechanical Simulations , 2015, IEEE Transactions on Biomedical Engineering.

[2]  R. D. Wood,et al.  Nonlinear Continuum Mechanics for Finite Element Analysis , 1997 .

[3]  Ray W. Ogden,et al.  Nonlinear Elasticity, Anisotropy, Material Stability and Residual Stresses in Soft Tissue , 2003 .

[4]  Boyce E. Griffith,et al.  Quasi-static image-based immersed boundary-finite element model of left ventricle under diastolic loading , 2014, International journal for numerical methods in biomedical engineering.

[5]  Derek Gaston,et al.  MOOSE: A parallel computational framework for coupled systems of nonlinear equations , 2009 .

[6]  Eric A. Sobie,et al.  Regression Analysis for Constraining Free Parameters in Electrophysiological Models of Cardiac Cells , 2009, PLoS Comput. Biol..

[7]  R. A. Forster,et al.  Analytical Benchmark Test Set for Criticality Code Verification , 1999 .

[8]  Manfred Liebmann,et al.  Anatomically accurate high resolution modeling of human whole heart electromechanics: A strongly scalable algebraic multigrid solver method for nonlinear deformation , 2016, J. Comput. Phys..

[9]  Joshua J. E. Blauer,et al.  Virtual Electrophysiological Study of Atrial Fibrillation in Fibrotic Remodeling , 2015, PloS one.

[10]  J. Humphrey,et al.  On constitutive relations and finite deformations of passive cardiac tissue: I. A pseudostrain-energy function. , 1987, Journal of biomechanical engineering.

[11]  D Noble,et al.  A meta‐analysis of cardiac electrophysiology computational models , 2009, Experimental physiology.

[12]  Andrew D McCulloch,et al.  Mouse and computational models link Mlc2v dephosphorylation to altered myosin kinetics in early cardiac disease. , 2012, The Journal of clinical investigation.

[13]  Gary R. Mirams,et al.  Variability in high-throughput ion-channel screening data and consequences for cardiac safety assessment , 2013, Journal of pharmacological and toxicological methods.

[14]  Nicolas P Smith,et al.  Quantifying inter‐species differences in contractile function through biophysical modelling , 2015, The Journal of physiology.

[15]  Pras Pathmanathan,et al.  Verification of computational models of cardiac electro‐physiology , 2014, International journal for numerical methods in biomedical engineering.

[16]  Gerhard A. Holzapfel,et al.  Nonlinear Solid Mechanics: A Continuum Approach for Engineering Science , 2000 .

[17]  Wei Chen,et al.  Assessing the Reliability of Complex Models: Mathematical and Statistical Foundations of Verification, Validation, and Uncertainty Quantification , 2012 .

[18]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[19]  Nicolas P Smith,et al.  An analysis of deformation‐dependent electromechanical coupling in the mouse heart , 2012, The Journal of physiology.

[20]  A. McCulloch,et al.  Finite element stress analysis of left ventricular mechanics in the beating dog heart. , 1995, Journal of biomechanics.

[21]  A. Quarteroni,et al.  Thermodynamically consistent orthotropic activation model capturing ventricular systolic wall thickening in cardiac electromechanics , 2014 .

[22]  Pablo Lamata,et al.  Analyses of the Redistribution of Work following Cardiac Resynchronisation Therapy in a Patient Specific Model , 2012, PloS one.

[23]  Viatcheslav Gurev,et al.  A high-resolution computational model of the deforming human heart , 2015, Biomechanics and modeling in mechanobiology.

[24]  R. Rivlin Large elastic deformations of isotropic materials. III. Some simple problems in cyclindrical polar co-ordinates , 1948, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[25]  Eric Kerfoot,et al.  Verification of cardiac tissue electrophysiology simulators using an N-version benchmark , 2011, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[26]  O. Dössel,et al.  Simulation of the contraction of the ventricles in a human heart model including atria and pericardium , 2014, Biomechanics and modeling in mechanobiology.

[27]  P. Hunter,et al.  Laminar structure of the heart: a mathematical model. , 1997, The American journal of physiology.

[28]  Edward W Hsu,et al.  Quantitative comparison of myocardial fiber structure between mice, rabbit, and sheep using diffusion tensor cardiovascular magnetic resonance , 2011, Journal of cardiovascular magnetic resonance : official journal of the Society for Cardiovascular Magnetic Resonance.

[29]  Alfio Quarteroni,et al.  An orthotropic active{strain model for the myocardium mechanics and its numerical approximation , 2014 .

[30]  Robert H. Anderson,et al.  The three‐dimensional arrangement of the myocytes in the ventricular walls , 2009, Clinical anatomy.

[31]  Sander Land,et al.  Computational modeling of Takotsubo cardiomyopathy: effect of spatially varying β-adrenergic stimulation in the rat left ventricle , 2014, American journal of physiology. Heart and circulatory physiology.

[32]  Heye Zhang,et al.  OpenCMISS: a multi-physics & multi-scale computational infrastructure for the VPH/Physiome project. , 2011, Progress in biophysics and molecular biology.

[33]  Alexander G. Fletcher,et al.  Chaste: An Open Source C++ Library for Computational Physiology and Biology , 2013, PLoS Comput. Biol..

[34]  Pras Pathmanathan,et al.  Ensuring reliability of safety-critical clinical applications of computational cardiac models , 2013, Front. Physiol..

[35]  Steven Niederer,et al.  Efficient Computational Methods for Strongly Coupled Cardiac Electromechanics , 2012, IEEE Transactions on Biomedical Engineering.

[36]  P J Hunter,et al.  A three-dimensional finite element method for large elastic deformations of ventricular myocardium: II--Prolate spheroidal coordinates. , 1996, Journal of biomechanical engineering.

[37]  J. C. Simo,et al.  Quasi-incompressible finite elasticity in principal stretches. Continuum basis and numerical algorithms , 1991 .

[38]  A. McCulloch,et al.  Modelling cardiac mechanical properties in three dimensions , 2001, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[39]  Dongbin Xiu,et al.  Parametric uncertainty analysis of pulse wave propagation in a model of a human arterial network , 2007, J. Comput. Phys..

[40]  Gerhard A Holzapfel,et al.  Constitutive modelling of passive myocardium: a structurally based framework for material characterization , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[41]  Kevin Burrage,et al.  Population of Computational Rabbit-Specific Ventricular Action Potential Models for Investigating Sources of Variability in Cellular Repolarisation , 2014, PloS one.

[42]  Paul R. Heyliger,et al.  Nonlinear continuum mechanics for finite element analysis (2nd edn). Javier Bonet and Richard D. Wood, Cambridge University Press, Cambridge, 2008. No. of pages: 318. ISBN: 978-0-521-83870-2 , 2008 .

[43]  P. Hunter,et al.  Computational mechanics of the heart : From tissue structure to ventricular function , 2000 .

[44]  A. Cristiano I. Malossi,et al.  Numerical Comparison and Calibration of Geometrical Multiscale Models for the Simulation of Arterial Flows , 2013 .

[45]  P. Hunter,et al.  Computational Mechanics of the Heart , 2000 .

[46]  Timothy G. Trucano,et al.  Verification and Validation in Computational Fluid Dynamics , 2002 .