Stability analysis of the POD reduced order method for solving the bidomain model in cardiac electrophysiology
LASSOUED, Jamila
Laboratoire de Modélisation Mathématique et Numérique dans les Sciences de l'Ingénieur [Tunis] [LR-LAMSIN-ENIT]
Laboratoire de Modélisation Mathématique et Numérique dans les Sciences de l'Ingénieur [Tunis] [LR-LAMSIN-ENIT]
MAHJOUB, Moncef
Université de Tunis El Manar [UTM]
Laboratoire de Modélisation Mathématique et Numérique dans les Sciences de l'Ingénieur [Tunis] [LR-LAMSIN-ENIT]
Leer más >
Université de Tunis El Manar [UTM]
Laboratoire de Modélisation Mathématique et Numérique dans les Sciences de l'Ingénieur [Tunis] [LR-LAMSIN-ENIT]
LASSOUED, Jamila
Laboratoire de Modélisation Mathématique et Numérique dans les Sciences de l'Ingénieur [Tunis] [LR-LAMSIN-ENIT]
Laboratoire de Modélisation Mathématique et Numérique dans les Sciences de l'Ingénieur [Tunis] [LR-LAMSIN-ENIT]
MAHJOUB, Moncef
Université de Tunis El Manar [UTM]
Laboratoire de Modélisation Mathématique et Numérique dans les Sciences de l'Ingénieur [Tunis] [LR-LAMSIN-ENIT]
< Leer menos
Université de Tunis El Manar [UTM]
Laboratoire de Modélisation Mathématique et Numérique dans les Sciences de l'Ingénieur [Tunis] [LR-LAMSIN-ENIT]
Idioma
en
Article de revue
Este ítem está publicado en
Mathematical Biosciences. 2015-12-04
Elsevier
Resumen en inglés
In this work we show the numerical stability of the Proper Orthogonal Decomposition (POD) reduced order method used in cardiac electrophysiology applications. The difficulty of proving the stability comes from the fact ...Leer más >
In this work we show the numerical stability of the Proper Orthogonal Decomposition (POD) reduced order method used in cardiac electrophysiology applications. The difficulty of proving the stability comes from the fact that we are interested in the bidomain model, which is a system of degenerate parabolic equations coupled to a system of ODEs representing the cell membrane electrical activity. The proof of the stability of this method is based an a priori estimate controlling the gap between the reduced order solution and the Galerkin finite element one. We present some numerical simulations confirming the theoretical results. We also combine the POD method with a time splitting scheme allowing a faster solution of the bidomain problem and show numerical results. Finally, we conduct numerical simulation in 2D illustrating the stability of the POD method in its sensitivity to the ionic model parameters. We also perform 3D simulation using a massively parallel code. We show the computational gain using the POD reduced order model. We also show that this method has a better scalability than the full finite element method.< Leer menos
Palabras clave en inglés
proper orthogonal decomposition
reduced order method
Bidomain equation
a priori estimates
ionic parameters
stability analysis
Mitchell-Schaeffer model
Proyecto ANR
L'Institut de Rythmologie et modélisation Cardiaque - ANR-10-IAHU-0004
Orígen
Importado de HalCentros de investigación