Title :
Accuracy and speed considerations of a new numerical method for modeling cardiac excitation
Author :
Allexandre, D. ; Otani, NF
Author_Institution :
Dept. of Biomed. Eng., Case Western Reserve Univ., Cleveland, OH, USA
Abstract :
A fast algorithm called the modified backward Euler (mBE) method, an explicit, variable timestep method, has been developed. Despite its explicit nature, the algorithm is unconditionally stable, allowing the use of timesteps several orders of magnitude larger than the timestep limit imposed by forward Euler methods. A guard cell algorithm also allows the use of timesteps larger than those imposed by the Courant Friedrich Levy (CFL) stability condition. The algorithm was applied to a two-dimensional tissue using the advanced Luo-Rudy ventricular cell model. Characteristic properties of action potential (AP) propagation including velocity, AP amplitude and duration, maximum membrane voltage and V˙max were computed versus computational speed and compared to those obtained from a forward Euler fixed timestep (FEFT) method. Overall, for a given accuracy in velocity of 1% to 3%, the computational speed was increased by a factor of 30 to 40 over the FEFT, while the other errors were kept below 1%.
Keywords :
bioelectric potentials; biological tissues; biomembranes; cardiology; medical computing; numerical analysis; 2D tissue; action potential propagation; advanced Luo-Rudy ventricular cell model; cardiac excitation modeling; computational speed; explicit variable timestep method; fast algorithm; forward Euler fixed timestep method; guard cell algorithm; maximum membrane voltage; modified backward Euler method; numerical method; timesteps; velocity; Biomedical engineering; Biomembranes; Cardiac tissue; Computational modeling; Equations; Iron; Noise measurement; Numerical models; Numerical stability; Voltage;
Conference_Titel :
Computers in Cardiology, 2002
Print_ISBN :
0-7803-7735-4
DOI :
10.1109/CIC.2002.1166694