Journal
MATHEMATICAL MODELS & METHODS IN APPLIED SCIENCES
Volume 14, Issue 6, Pages 883-911Publisher
WORLD SCIENTIFIC PUBL CO PTE LTD
DOI: 10.1142/S0218202504003489
Keywords
reaction-diffusion equations; bidomain model; finite elements; parallel solver
Categories
Ask authors/readers for more resources
In this work, a parallel three-dimensional solver for numerical simulations in computational electrocardiology is introduced and studied. The solver is based on the anisotropic Bidomain cardiac model, consisting of a system of two degenerate parabolic reaction-diffusion equations describing the intra and extracellular potentials of the myocardial tissue. This model includes intramural fiber rotation and anisotropic conductivity coefficients that can be fully orthotropic or axially symmetric around the fiber direction. The solver also includes the simpler anisotropic Monodomain model, consisting of only one reaction-diffusion equation. These cardiac models are coupled with a membrane model for the ionic currents, consisting of a system of ordinary differential equations that can vary from the simple FitzHugh-Nagumo (FHN) model to the more complex phase-I Luo-Rudy model (LR1). The solver employs structured isoparametric Q(1) finite elements in space and a semi-implicit adaptive method in time. Parallelization and portability are based on the PETSc parallel library. Large-scale computations with up to O(10(7)) unknowns have been run on parallel computers, simulating excitation and repolarization phenomena in three-dimensional domains.
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