WIAS Preprint No. 1827, (2013)

Convergence of an implicit Voronoi finite volume method for reaction-diffusion problems



Authors

  • Fiebach, André
  • Glitzky, Annegret
    ORCID: 0000-0003-1995-5491
  • Linke, Alexander
    ORCID: 0000-0002-0165-2698

2010 Mathematics Subject Classification

  • 35K57 35R05 65M08, 65M12, 80A30

Keywords

  • reaction-diffusion systems, heterostructures, finite volume method, convergence, long-term simulation

DOI

10.20347/WIAS.PREPRINT.1827

Abstract

We investigate the convergence of an implicit Voronoi finite volume method for reaction- diffusion problems including nonlinear diffusion in two space dimensions. The model allows to handle heterogeneous materials and uses the chemical potentials of the involved species as primary variables. The numerical scheme uses boundary conforming Delaunay meshes and preserves positivity and the dissipative property of the continuous system. Starting from a result on the global stability of the scheme (uniform, mesh-independent global upper and lower bounds), we prove strong convergence of the chemical activities and their gradients to a weak solution of the continuous problem. In order to illustrate the preservation of qualitative properties by the numerical scheme, we present a long-term simulation of the Michaelis-Menten-Henri system. Especially, we investigate the decay properties of the relative free energy and the evolution of the dissipation rate over several magnitudes of time, and obtain experimental orders of convergence for these quantities.

Appeared in

  • Numer. Methods Partial Differential Equations, 32 (2016), pp. 141--174.

Download Documents