A numerical scheme for modeling wavefront propagation on a monolayer of arbitrary geometry.
The majority of models of wavefront propagation in cardiac tissue have assumed relatively simple geometries. Extensions to complicated three-dimensional (3-D) representations are computationally challenging due to issues related both to problem size and to the correct implementation of flux conservation. In this paper, we present a generalized finite difference scheme (GDFS) to simulate the reaction-diffusion system on a 3-D monolayer of arbitrary shape. GDFS is a vertex-centered variant of the finite-volume method that ensures local flux conservation. Owing to an effectively lower dimensionality, the overall computation time is reduced compared to full 3-D models at the same spatial resolution. We present the theoretical background to compute both the wavefront conduction and local electrograms using a matrix formulation. The same matrix is used for both these quantities. We then give some results of simulation for simple monolayers and complex monolayers resembling a human atria.
Signal Processing Institute, EPFL, CH-1015 Lausanne, Switzerland. steeve.zoror@lis.inpg.fr
Steeve
S
Zozor
Olivier
O
Blanc
Vincent
V
Jacquemet
Nathalie
N
Virag
Jean-Marc
JM
Vesin
Etienne
E
Pruvot
Lukas
L
Kappenberger
Craig
C
Henriquez
D000465:Algorithms; D018780:Body Surface Potential Mapping / Q000379:methods; D003198:Computer Simulation; D004562:Electrocardiography / Q000379:methods*; D004574:Electromagnetic Fields; D006321:Heart / Q000502:physiology*; D006329:Heart Conduction System / Q000502:physiology*; D006801:Humans; D008564:Membrane Potentials / Q000502:physiology; D008955:Models, Cardiovascular; D032383:Myocytes, Cardiac / Q000502:physiology