Soft Tissue Modelling of Cardiac Fibres for Use in Coupled Mechano-Electric Simulations

Abstract The numerical solution of the coupled system of partial differential and ordinary differential equations that model the whole heart in three dimensions is a considerable computational challenge. As a consequence, it is not computationally practical—either in terms of memory or time—to repeat simulations on a finer computational mesh to ensure that convergence of the solution has been attained. In an attempt to avoid this problem while retaining mathematical rigour, we derive a one dimensional model of a cardiac fibre that takes account of elasticity properties in three structurally defined axes within the myocardial tissue. This model of a cardiac fibre is then coupled with an electrophysiological cell model and a model of cellular electromechanics to allow us to simulate the coupling of the electrical and mechanical activity of the heart. We demonstrate that currently used numerical methods for coupling electrical and mechanical activity do not work in this case, and identify appropriate numerical techniques that may be used when solving the governing equations. This allows us to perform a series of simulations that: (i) investigate the effect of some of the assumptions inherent in other models; and (ii) reproduce qualitatively some experimental observations.

[1]  M. Nash,et al.  Electromechanical model of excitable tissue to study reentrant cardiac arrhythmias. , 2004, Progress in biophysics and molecular biology.

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

[3]  A. Grimm,et al.  Deformation of the diastolic left ventricle. Nonlinear elastic effects. , 1973, Biophysical journal.

[4]  P. Hunter,et al.  New developments in a strongly coupled cardiac electromechanical model. , 2005, Europace : European pacing, arrhythmias, and cardiac electrophysiology : journal of the working groups on cardiac pacing, arrhythmias, and cardiac cellular electrophysiology of the European Society of Cardiology.

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

[6]  C. Ashley,et al.  Effect of pH, phosphate, and ADP on relaxation of myocardium after photolysis of diazo 2. , 1998, American journal of physiology. Heart and circulatory physiology.

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

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

[9]  Frederick Sachs,et al.  Cardiac mechano-electric feedback and arrhythmias : from pipette to patient , 2005 .

[10]  T. Bortfeld,et al.  Correlation between CT numbers and tissue parameters needed for Monte Carlo simulations of clinical dose distributions. , 2000, Physics in medicine and biology.

[11]  Jonathan P. Whiteley,et al.  An Efficient Numerical Technique for the Solution of the Monodomain and Bidomain Equations , 2006, IEEE Transactions on Biomedical Engineering.

[12]  C. Makridakis Finite element approximations of nonlinear elastic waves , 1993 .

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

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

[15]  D. Noble,et al.  A model for human ventricular tissue. , 2004, American journal of physiology. Heart and circulatory physiology.

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

[17]  Arun V. Holden,et al.  Computational biology of the heart , 1998, The Mathematical Gazette.

[18]  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.

[19]  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.

[20]  K. Morton,et al.  Numerical Solution of Partial Differential Equations , 1995 .

[21]  David P. Nickerson,et al.  Modelling cardiac electro-mechanics: From CellML to the whole heart , 2003 .

[22]  Andrew D. McCulloch,et al.  Effect of Laminar Orthotropic Myofiber Architecture on Regional Stress and Strain in the Canine Left Ventricle , 2000 .

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

[24]  A. Grimm,et al.  Deformation of the diastolic left ventricle--II. Nonlinear geometric effects. , 1974, Journal of biomechanics.

[25]  Guy Salama,et al.  Simultaneous maps of optical action potentials and calcium transients in guinea‐pig hearts: mechanisms underlying concordant alternans , 2000, The Journal of physiology.