An Unconditional Positivity-Preserving Difference Scheme for Models of Cancer Migration and Invasion

Abstract

In this paper, we consider models of cancer migration and invasion, which consist of two nonlinear parabolic equations (one of the convection–diffusion reaction type and the other of the diffusion–reaction type) and an additional nonlinear ordinary differential equation. The unknowns represent concentrations or densities that cannot be negative. Widely used approximations, such as difference schemes, can produce negative solutions because of truncation errors and can become unstable. We propose a new difference scheme that guarantees the positivity of the numerical solution for arbitrary mesh step sizes. It has explicit and fast performance even for nonlinear reaction terms that consist of sums of positive and negative functions. The numerical examples illustrate the simplicity and efficiency of the method. A numerical simulation of a model of cancer migration is also discussed.