University of Oxford Logo Chaste logo University of OxfordDepartment of Computer Science - Home

Cardiac Chaste: developing software for realistic heart simulations

Introduction and aims

Computational modelling of the heart is now recognised as a powerful technique in the detailed investigation of cardiac behaviour. One of the major contributions of computational approaches to cardiovascular research has been the ability to dissect various effects and to tease out important relations between parameters, which are not possible using current experimental techniques. Recently, advanced computer models of cardiac electro-mechanical activity have been developed. The electrical properties of the myocardium are generally described by the bidomain equations, a set of coupled parabolic and elliptic partial differential equations (PDEs) that represents the tissue as two separate, distinct continua - one intracellular and the other extracellular. The intracellular and the extracellular media are connected via the cell membrane, and thus the two PDEs are coupled at each point in space through a set of complex, non-linear ordinary differential equations (ODEs), which describe the ionic transport across the cell membrane. Certain modelling environments use the monodomain representation of cardiac activity, which involves solving a single parabolic PDE, by assuming either that the extracellular potentials are negligible, or that the anisotropy ratios are equal in the intracellular and the extracellular domains.

Although current cardiac solvers have achieved very good efficiency for particular applications, they all share the following important limitations: (i) they are not generic - they have been developed for specific applications based on the particular scientific interests of the developers; (ii) they have not been developed using state-of-the-art software engineering methods, and therefore have not been completely tested and validated; (iii) they have not achieved maximum efficiency to run on HPC platforms because they have not been designed for this purpose; (iv) they are not freely available to the scientific community; and (v) they do not include state-of-the-art numerical and computational techniques that could speed up the code by several orders of magnitude and would thus allow for full numerical convergence testing (in particular, current simulation tools cannot provide sufficient spatial resolution to test convergence).

The overall objective of the Chaste project is to develop a novel, generic and efficient cardiac software package, accurately tested and validated, and specifically designed to run on large scale, high performance systems such as HPCx and future generations of petascale computers. The use of advanced software engineering techniques, and close interaction with academic and industrial collaborators as well as with end users, ensures that the resulting software forms a robust, user friendly, production level platform for the investigation of important scientific questions related to cardiac function in health and disease. We expect this software to support the user community to move on to investigate the next generation of cardiac physiology research.

User functionality

Performance

Our collaborators in Fujitsu have run Chaste in HECToR (UK's high-end computing resource), obtaining the following performance for an anatomically-based rabbit ventricular mesh ~4 million nodes; see Bernabeu et al. Computers in Cardiology 2008):

Previously:

Example simulation results

Cardiac simulations can be viewed on the Chaste YouTube channel.