Efficient solution of ordinary differential equations modeling electrical activity in cardiac cells
J Sundnes1, G T Lines, A Tveito
1Department of Informatics, University of Oslo, P.O. Box 1080, Blindern, N-0316 Oslo, Norway. sundnes@ifi.uio.no
Abstract:
The contraction of the heart is preceded and caused by a cellular electro-chemical reaction, causing an electrical field to be generated. Performing realistic computer simulations of this process involves solving a set of partial differential equations, as well as a large number of ordinary differential equations (ODEs) characterizing the reactive behavior of the cardiac tissue. Experiments have shown that the solution of the ODEs contribute significantly to the total work of a simulation, and there is thus a strong need to utilize efficient solution methods for this part of the problem. This paper presents how an efficient implicit Runge-Kutta method may be adapted to solve a complicated cardiac cell model consisting of 31 ODEs, and how this solver may be coupled to a set of PDE solvers to provide complete simulations of the electrical activity.
Related Concept Videos
Specialized Characteristics of Cardiac Muscles
Cardiac muscle cells are smaller than skeletal muscles, averaging 10–20 mm in diameter and 50–100 mm in length. However, they have large energy demands for continuous contraction and relaxation. This energy is almost exclusively derived from aerobic metabolism of energy reserves in...
Electrophysiology of Normal Cardiac Rhythm
Bode Plots Construction
The Electrical Double Layer
Linear Approximation in Frequency Domain
In contrast, nonlinear systems do not inherently possess these properties. However, for small deviations around an operating point, a nonlinear system can often be approximated as linear.
Cardiac Action Potential
The cardiac action potential process involves a series of phases characterized by the movement of ions across the cardiac cell membranes, leading to the depolarization and repolarization of the cardiac myocytes.
Ionic Basis of Cardiac Action Potentials


