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

Forskningsoutput: TidskriftsbidragArtikel i vetenskaplig tidskrift

Abstract

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.

Detaljer

Författare
  • Jonas Schäfer
  • Xuan Huang
  • Philipp Birken
  • Matthias K. Gobbert
  • Meister Andreas
Enheter & grupper
Forskningsområden

Ämnesklassifikation (UKÄ) – OBLIGATORISK

  • Matematik

Nyckelord

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