A Memory-Efficient Finite Volume Method for Advection-Diffusion-Reaction Systems with Non-Smooth Sources

Forskningsoutput: TidskriftsbidragArtikel i vetenskaplig tidskrift


We present a parallel matrix-free implicit nite volume scheme for the solution of unsteady three-dimensional advection-diusion-reaction equations with smooth and Dirac-Delta source terms. The scheme is formally second order in space and a Newton-Krylov method is employed for the appearing nonlinear systems in the implicit time integration. The matrix-vector product required is hardcoded without any approximations, obtaining a matrix-free method that needs little storage and is well suited for parallel implementation. We describe the matrix-free implementation of the method in detail and give numerical evidence of its second order convergence in the presence of smooth source terms. For non-smooth source terms the convergence order drops to one half. Furthermore, we demonstrate the method's applicability for the long time simulation of calcium ow in heart cells and show its parallel scaling.


  • Jonas Schäfer
  • Xuan Huang
  • Philipp Birken
  • Matthias K. Gobbert
  • Meister Andreas
Enheter & grupper

Ämnesklassifikation (UKÄ) – OBLIGATORISK

  • Matematik


Sidor (från-till)143-167
TidskriftNumerical Methods for Partial Differential Equations
Utgåva nummer1
StatusPublished - 2015
Peer review utfördJa