A memory-efficient finite element method for systems of reaction-diffusion equations with non-smooth forcing

Alexander L. Hanhart, Matthias K. Gobbert, Leighton T Izu

Research output: Contribution to journalArticle

17 Citations (Scopus)

Abstract

The release of calcium ions in a human heart cell is modeled by a system of reaction-diffusion equations, which describe the interaction of the chemical species and the effects of various cell processes on them. The release is modeled by a forcing term in the calcium equation that involves a superposition of many Dirac delta functions in space; such a nonsmooth right-hand side leads to divergence for many numerical methods. The calcium ions enter the cell at a large number of regularly spaced points throughout the cell; to resolve those points adequately for a cell with realistic three-dimensional dimensions, an extremely fine spatial mesh is needed. A finite element method is developed that addresses the two crucial issues for this and similar applications: Convergence of the method is demonstrated in extension of the classical theory that does not apply to nonsmooth forcing functions like the Dirac delta function; and the memory usage of the method is optimal and thus allows for extremely fine three-dimensional meshes with many millions of degrees of freedom, already on a serial computer. Additionally, a coarse-grained parallel implementation of the algorithm allows for the solution on meshes with yet finer resolution than possible in serial.

Original languageEnglish (US)
Pages (from-to)431-458
Number of pages28
JournalJournal of Computational and Applied Mathematics
Volume169
Issue number2
DOIs
StatePublished - Aug 15 2004
Externally publishedYes

Fingerprint

Reaction-diffusion Equations
Forcing
Delta functions
Calcium
Computer systems
Finite Element Method
Finite element method
Data storage equipment
Cell
Dirac delta function
Mesh
Ions
Numerical methods
Three-dimensional
Forcing Term
Parallel Implementation
Superposition
Resolve
Divergence
Degree of freedom

Keywords

  • Cluster computing
  • Galerkin method
  • Matrix-free iterative method
  • Non-smooth data
  • Reaction-diffusion equation

ASJC Scopus subject areas

  • Applied Mathematics
  • Computational Mathematics
  • Numerical Analysis

Cite this

A memory-efficient finite element method for systems of reaction-diffusion equations with non-smooth forcing. / Hanhart, Alexander L.; Gobbert, Matthias K.; Izu, Leighton T.

In: Journal of Computational and Applied Mathematics, Vol. 169, No. 2, 15.08.2004, p. 431-458.

Research output: Contribution to journalArticle

@article{a1b29b1307464b9dad7524006237d9d9,
title = "A memory-efficient finite element method for systems of reaction-diffusion equations with non-smooth forcing",
abstract = "The release of calcium ions in a human heart cell is modeled by a system of reaction-diffusion equations, which describe the interaction of the chemical species and the effects of various cell processes on them. The release is modeled by a forcing term in the calcium equation that involves a superposition of many Dirac delta functions in space; such a nonsmooth right-hand side leads to divergence for many numerical methods. The calcium ions enter the cell at a large number of regularly spaced points throughout the cell; to resolve those points adequately for a cell with realistic three-dimensional dimensions, an extremely fine spatial mesh is needed. A finite element method is developed that addresses the two crucial issues for this and similar applications: Convergence of the method is demonstrated in extension of the classical theory that does not apply to nonsmooth forcing functions like the Dirac delta function; and the memory usage of the method is optimal and thus allows for extremely fine three-dimensional meshes with many millions of degrees of freedom, already on a serial computer. Additionally, a coarse-grained parallel implementation of the algorithm allows for the solution on meshes with yet finer resolution than possible in serial.",
keywords = "Cluster computing, Galerkin method, Matrix-free iterative method, Non-smooth data, Reaction-diffusion equation",
author = "Hanhart, {Alexander L.} and Gobbert, {Matthias K.} and Izu, {Leighton T}",
year = "2004",
month = "8",
day = "15",
doi = "10.1016/j.cam.2003.12.035",
language = "English (US)",
volume = "169",
pages = "431--458",
journal = "Journal of Computational and Applied Mathematics",
issn = "0377-0427",
publisher = "Elsevier",
number = "2",

}

TY - JOUR

T1 - A memory-efficient finite element method for systems of reaction-diffusion equations with non-smooth forcing

AU - Hanhart, Alexander L.

AU - Gobbert, Matthias K.

AU - Izu, Leighton T

PY - 2004/8/15

Y1 - 2004/8/15

N2 - The release of calcium ions in a human heart cell is modeled by a system of reaction-diffusion equations, which describe the interaction of the chemical species and the effects of various cell processes on them. The release is modeled by a forcing term in the calcium equation that involves a superposition of many Dirac delta functions in space; such a nonsmooth right-hand side leads to divergence for many numerical methods. The calcium ions enter the cell at a large number of regularly spaced points throughout the cell; to resolve those points adequately for a cell with realistic three-dimensional dimensions, an extremely fine spatial mesh is needed. A finite element method is developed that addresses the two crucial issues for this and similar applications: Convergence of the method is demonstrated in extension of the classical theory that does not apply to nonsmooth forcing functions like the Dirac delta function; and the memory usage of the method is optimal and thus allows for extremely fine three-dimensional meshes with many millions of degrees of freedom, already on a serial computer. Additionally, a coarse-grained parallel implementation of the algorithm allows for the solution on meshes with yet finer resolution than possible in serial.

AB - The release of calcium ions in a human heart cell is modeled by a system of reaction-diffusion equations, which describe the interaction of the chemical species and the effects of various cell processes on them. The release is modeled by a forcing term in the calcium equation that involves a superposition of many Dirac delta functions in space; such a nonsmooth right-hand side leads to divergence for many numerical methods. The calcium ions enter the cell at a large number of regularly spaced points throughout the cell; to resolve those points adequately for a cell with realistic three-dimensional dimensions, an extremely fine spatial mesh is needed. A finite element method is developed that addresses the two crucial issues for this and similar applications: Convergence of the method is demonstrated in extension of the classical theory that does not apply to nonsmooth forcing functions like the Dirac delta function; and the memory usage of the method is optimal and thus allows for extremely fine three-dimensional meshes with many millions of degrees of freedom, already on a serial computer. Additionally, a coarse-grained parallel implementation of the algorithm allows for the solution on meshes with yet finer resolution than possible in serial.

KW - Cluster computing

KW - Galerkin method

KW - Matrix-free iterative method

KW - Non-smooth data

KW - Reaction-diffusion equation

UR - http://www.scopus.com/inward/record.url?scp=4444290686&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=4444290686&partnerID=8YFLogxK

U2 - 10.1016/j.cam.2003.12.035

DO - 10.1016/j.cam.2003.12.035

M3 - Article

AN - SCOPUS:4444290686

VL - 169

SP - 431

EP - 458

JO - Journal of Computational and Applied Mathematics

JF - Journal of Computational and Applied Mathematics

SN - 0377-0427

IS - 2

ER -