Adaptive numerical simulation of intracellular calcium dynamics using domain decomposition methods

Authors:
Ch. Nagaiah;S. Rüdiger;G. Warnecke;M. Falcke
Affiliations:
Institute of Mathematics and Scientific Computing, University of Graz, A-8010, Austria and Institute of Analysis and Numerics, Otto-von-Guericke University, Universitätsplatz 2, 39106 Magdebu ...;Institut für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15, 12489 Berlin, Germany and Hahn-Meitner Institut, Glienicker Str. 100, 14109 Berlin, Germany;Institute of Analysis and Numerics, Otto-von-Guericke University, Universitätsplatz 2, 39106 Magdeburg, Germany;Hahn-Meitner Institut, Glienicker Str. 100, 14109 Berlin, Germany
Venue:
Applied Numerical Mathematics
Year:
2008

Citing 3
Cited 1

A PI stepsize control for the numerical solution of ordinary differential equations

BIT
BI-CGSTAB: a fast and smoothly converging variant of BI-CG for the solution of nonsymmetric linear systems

SIAM Journal on Scientific and Statistical Computing
Matrix-free W-methods using a multiple Arnoldi iteration

NUMDIFF-7 Selected papers of the seventh conference on Numerical treatment of differential equations

Spectral method and high-order finite differences for the nonlinear cable equation

Neural Computation

Quantified Score

Hi-index	0.01

Visualization

Abstract

In this paper we present adaptive numerical simulations of intracellular calcium dynamics using domain decomposition methods. The modeling of diffusion, binding and membrane transport of calcium ions in cells leads to a system of reaction-diffusion equations. We describe the modulation of cytosolic and ER calcium concentration close to the membrane of the cell organelle. This leads to a two-dimensional model. Local grid refinement is adjusted to the strongly local effects of the calcium dynamics. A finite element method is used for spatial discretization and a Rosenbrock solver for time integration which is third-order accurate and very suitable for nonlinear parabolic problems. We present sequential and parallelized numerical results. The latter are based on domain decomposition methods.