- Journal Home
- Volume 21 - 2024
- Volume 20 - 2023
- Volume 19 - 2022
- Volume 18 - 2021
- Volume 17 - 2020
- Volume 16 - 2019
- Volume 15 - 2018
- Volume 14 - 2017
- Volume 13 - 2016
- Volume 12 - 2015
- Volume 11 - 2014
- Volume 10 - 2013
- Volume 9 - 2012
- Volume 8 - 2011
- Volume 7 - 2010
- Volume 6 - 2009
- Volume 5 - 2008
- Volume 4 - 2007
- Volume 3 - 2006
- Volume 2 - 2005
- Volume 1 - 2004
Cited by
- BibTex
- RIS
- TXT
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.
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.