Year: 2011
Author: M. Hanslien, R. Artebrant, A. Tveito, G. T. Lines, X. Cai
International Journal of Numerical Analysis and Modeling, Vol. 8 (2011), Iss. 3 : pp. 427–442
Abstract
We propose a second-order accurate method for computing the
solutions to the Aliev-Panfilov model of cardiac excitation. This two-variable
reaction-diffusion system is due to its simplicity a popular choice for modeling
important problems in electrocardiology; e.g. cardiac arrhythmias. The
solutions might be very complicated in structure, and hence highly resolved
numerical simulations are called for to capture the fine details. Usually the forward
Euler time-integrator is applied in these computations; it is very simple
to implement and can be effective for coarse grids. For fine-scale simulations,
however, the forward Euler method suffers from a severe time-step restriction,
rendering it less efficient for simulations where high resolution and accuracy
are important.
We analyze the stability of the proposed second-order method and the forward
Euler scheme when applied to the Aliev-Panfilov model. Compared to the Euler
method the suggested scheme has a much weaker time-step restriction, and
promises to be more efficient for computations on finer meshes.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/2011-IJNAM-694
International Journal of Numerical Analysis and Modeling, Vol. 8 (2011), Iss. 3 : pp. 427–442
Published online: 2011-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 16
Keywords: reaction-diffusion system implict Runge-Kutta electrocardiology.