Journal
MATHEMATICAL BIOSCIENCES
Volume 208, Issue 2, Pages 375-392Publisher
ELSEVIER SCIENCE INC
DOI: 10.1016/j.mbs.2006.09.006
Keywords
stable scheme; maximum principle; defibrillation; Luo-Rudy 1
Categories
Ask authors/readers for more resources
Numerical simulations of defibrillation using the Bidomain model coupled to a model of membrane kinetics represent a serious numerical challenge. This is because very high voltages close to defibrillation electrodes demand that extreme time step restrictions be placed on standard numerical schemes, e.g. the forward Euler scheme. A common solution to this problem is to modify the cell model by simple if-tests applied to several equations and rate functions. These changes are motivated by numerical problems rather than physiology, and should therefore be avoided whenever possible. The purpose of this paper is to present a numerical scheme that handles the original model without modifications and which is unconditionally stable for the Luo-Rudy phase I model. This also shows that the cell model is mathematically well-behaved, even in the presence of very high voltages. Our theoretical results are illustrated by numerical computations. (C) 2006 Elsevier Inc. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available