We consider discretizations for reaction-diffusion systems with nonlinear
diffusion in two space dimensions. The applied model allows to handle heterogeneous
materials and uses the chemical potentials of the involved species as primary variables.
We propose an implicit Voronoi finite volume discretization on regular Delaunay
meshes that allows to prove uniform, mesh-independent global upper and lower L1
bounds for the chemical potentials. These bounds provide the main step for a convergence
analysis for the full discretized nonlinear evolution problem. The fundamental
ideas are energy estimates, a discrete Moser iteration and the use of discrete
Gagliardo-Nirenberg inequalities. For the proof of the Gagliardo-Nirenberg inequalities
we exploit that the discrete Voronoi finite volume gradient norm in 2d coincides
with the gradient norm of continuous piecewise linear finite elements.
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.