On a boundedness-preserving semi-linear discretization of a two-dimensional nonlinear diffusion-reaction model

Abstract

Departing from a method to approximate the solutions of a two-dimensional generalization of the well-known Fisher's equation from population dynamics, we extend this computational technique to calculate the solutions of a FitzHugh-Nagumo model and derive conditions under which its positive and bounded analytic solutions are estimated consistently by positive and bounded numerical approximations. The constraints are relatively flexible, and they are provided exclusively in terms of the model coefficients and the computational parameters. The proofs are established with the help of the theory of M-matrices, using the facts that such matrices are non-singular, and that the entries of their inverses are positive numbers. Some numerical experiments are performed in order to show that our method is capable of preserving the positivity and the boundedness of the numerical solutions. The simulations evince a good agreement between the numerical estimations and the corresponding exact solutions derived in this work. © 2012 Copyright Taylor and Francis Group, LLC.