OpenCMISS: a multi-physics & multi-scale computational infrastructure for the VPH/Physiome project.

The VPH/Physiome Project is developing the model encoding standards CellML (cellml.org) and FieldML (fieldml.org) as well as web-accessible model repositories based on these standards (models.physiome.org). Freely available open source computational modelling software is also being developed to solve the partial differential equations described by the models and to visualise results. The OpenCMISS code (opencmiss.org), described here, has been developed by the authors over the last six years to replace the CMISS code that has supported a number of organ system Physiome projects. OpenCMISS is designed to encompass multiple sets of physical equations and to link subcellular and tissue-level biophysical processes into organ-level processes. In the Heart Physiome project, for example, the large deformation mechanics of the myocardial wall need to be coupled to both ventricular flow and embedded coronary flow, and the reaction-diffusion equations that govern the propagation of electrical waves through myocardial tissue need to be coupled with equations that describe the ion channel currents that flow through the cardiac cell membranes. In this paper we discuss the design principles and distributed memory architecture behind the OpenCMISS code. We also discuss the design of the interfaces that link the sets of physical equations across common boundaries (such as fluid-structure coupling), or between spatial fields over the same domain (such as coupled electromechanics), and the concepts behind CellML and FieldML that are embodied in the OpenCMISS data structures. We show how all of these provide a flexible infrastructure for combining models developed across the VPH/Physiome community.

[1]  William L. Briggs,et al.  A multigrid tutorial , 1987 .

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

[3]  Geoffrey McLennan,et al.  CT-based geometry analysis and finite element models of the human and ovine bronchial tree. , 2004, Journal of applied physiology.

[4]  Andrea Toselli,et al.  Domain decomposition methods : algorithms and theory , 2005 .

[5]  D. Rypl,et al.  Parallel explicit finite element dynamics with nonlocal constitutive models , 2001 .

[6]  P. Hunter,et al.  Integration from proteins to organs: the Physiome Project , 2003, Nature Reviews Molecular Cell Biology.

[7]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[8]  I. G. Currie Fundamental mechanics of fluids , 1974 .

[9]  Nicolas Le Novère,et al.  SED-ML - An XML Format for the Implementation of the MIASE Guidelines , 2008, CMSB.

[10]  J. Keener An eikonal-curvature equation for action potential propagation in myocardium , 1991, Journal of mathematical biology.

[11]  G. Richard Christie,et al.  Modelling and visualising the heart , 2002 .

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

[13]  Tarek P. Mathew,et al.  Domain Decomposition Methods for the Numerical Solution of Partial Differential Equations , 2008, Lecture Notes in Computational Science and Engineering.

[14]  D. Noble,et al.  Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes. , 1998, The Canadian journal of cardiology.

[15]  N P Smith,et al.  Coupling multi-physics models to cardiac mechanics. , 2011, Progress in biophysics and molecular biology.

[16]  P. J. Hunter,et al.  Generation of an Anatomically Based Three-Dimensional Model of the Conducting Airways , 2000, Annals of Biomedical Engineering.

[17]  Hervé Delingette,et al.  Proceedings of the 5th International Conference on Functional Imaging and Modeling of the Heart , 2009 .

[18]  Peter Hunter,et al.  A strategy for integrative computational physiology. , 2005, Physiology.

[19]  David P. Nickerson,et al.  An overview of the CellML API and its implementation , 2010, BMC Bioinformatics.

[20]  P. Hunter,et al.  The analysis of cardiac function: a continuum approach. , 1988, Progress in biophysics and molecular biology.

[21]  Peter J Hunter,et al.  Modeling total heart function. , 2003, Annual review of biomedical engineering.

[22]  P. Hunter,et al.  A quantitative analysis of cardiac myocyte relaxation: a simulation study. , 2006, Biophysical journal.

[23]  Barry F. Smith,et al.  Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations , 1996 .

[24]  Andrew J. Pullan,et al.  A Finite Element Method for an Eikonal Equation Model of Myocardial Excitation Wavefront Propagation , 2002, SIAM J. Appl. Math..

[25]  D Gavaghan,et al.  Rabbit-specific ventricular model of cardiac electrophysiological function including specialized conduction system. , 2011, Progress in biophysics and molecular biology.

[26]  Andrew J. Pullan,et al.  An Anatomically Based Model of Transient Coronary Blood Flow in the Heart , 2002, SIAM J. Appl. Math..

[27]  P J Hunter,et al.  The IUPS Physiome Project: a framework for computational physiology. , 2004, Progress in biophysics and molecular biology.

[28]  I. E. Vignon-Clementel,et al.  A poroelastic model valid in large strains with applications to perfusion in cardiac modeling , 2010 .

[29]  David Nordsletten,et al.  A non-conforming monolithic finite element method for problems of coupled mechanics , 2010, J. Comput. Phys..

[30]  T. Mathew,et al.  Domain Decomposition Methods for the Numerical Solution of Partial Differential Equations (Lecture Notes in Computational Science and Engineering) , 2008 .

[31]  A. Garfinkel,et al.  An advanced algorithm for solving partial differential equation in cardiac conduction , 1999, IEEE Transactions on Biomedical Engineering.

[32]  Alexander G. Fletcher,et al.  Chaste: A test-driven approach to software development for biological modelling , 2009, Comput. Phys. Commun..

[33]  A. Tveito,et al.  An operator splitting method for solving the bidomain equations coupled to a volume conductor model for the torso. , 2005, Mathematical biosciences.

[34]  YingLiang Ma,et al.  Anisotropic Wave Propagation and Apparent Conductivity Estimation in a Fast Electrophysiological Model: Application to XMR Interventional Imaging , 2007, MICCAI.

[35]  Michael J. Holst,et al.  Adaptive Numerical Treatment of Elliptic Systems on Manifolds , 2001, Adv. Comput. Math..

[36]  Giorgio Fotia,et al.  Simulation of cardiac electrophysiology on next-generation high-performance computers , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[37]  J. Lambert Numerical Methods for Ordinary Differential Equations , 1991 .

[38]  James P. Keener,et al.  Mathematical physiology , 1998 .

[39]  Alan Garfinkel,et al.  Multi-scale modeling in biology: how to bridge the gaps between scales? , 2011, Progress in biophysics and molecular biology.

[40]  Benjamin S. Kirk,et al.  Library for Parallel Adaptive Mesh Refinement / Coarsening Simulations , 2006 .

[41]  Carole A. Goble,et al.  Taverna: a tool for building and running workflows of services , 2006, Nucleic Acids Res..

[42]  J. Oden Finite Elements of Nonlinear Continua , 1971 .

[43]  W. Bangerth,et al.  deal.II—A general-purpose object-oriented finite element library , 2007, TOMS.

[44]  Alejandro F. Frangi,et al.  GIMIAS: An Open Source Framework for Efficient Development of Research Tools and Clinical Prototypes , 2009, FIMH.

[45]  P J Hunter,et al.  The Noble cardiac ventricular electrophysiology models in CellML. , 2006, Progress in biophysics and molecular biology.

[46]  J. Sethian,et al.  Fast methods for the Eikonal and related Hamilton- Jacobi equations on unstructured meshes. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[47]  S. Tentoni,et al.  Mathematical modeling of the excitation process in myocardial tissue: influence of fiber rotation on wavefront propagation and potential field. , 1990, Mathematical biosciences.

[48]  P. Hunter,et al.  Modelling the mechanical properties of cardiac muscle. , 1998, Progress in biophysics and molecular biology.

[49]  A. J. Pullan,et al.  Geometric modeling of the human torso using cubic hermite elements , 2007, Annals of Biomedical Engineering.

[50]  Nicholas J. Higham,et al.  INVERSE PROBLEMS NEWSLETTER , 1991 .

[51]  Andrew J. Pullan,et al.  Mathematically Modelling the Electrical Activity of the Heart: From Cell to Body Surface and Back Again , 2005 .

[52]  Bořek Patzák,et al.  Design of object oriented finite element code , 2001 .

[53]  L. E. Malvern Introduction to the mechanics of a continuous medium , 1969 .

[54]  G Plank,et al.  Computational tools for modeling electrical activity in cardiac tissue. , 2003, Journal of electrocardiology.

[55]  Edmund J. Crampin,et al.  Multiscale computational modelling of the heart , 2004, Acta Numerica.

[56]  P. J. Hunter,et al.  Generation of an Anatomically Based Geometric Coronary Model , 2004, Annals of Biomedical Engineering.

[57]  G. Batchelor,et al.  An Introduction to Fluid Dynamics , 1968 .

[58]  Peter J. Hunter,et al.  The CellML Model Repository , 2008, Bioinform..

[59]  Peter J. Hunter,et al.  FieldML: concepts and implementation , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[60]  Peter J. Hunter,et al.  An Overview of CellML 1.1, a Biological Model Description Language , 2003, Simul..