Systematic Investigation of the Explicit, Dynamically Consistent Methods for Fisher’s Equation

We systematically investigate the performance of numerical methods to solve Fisher’s equation, which contains a linear diffusion term and a nonlinear logistic term. The usual explicit finite difference algorithms are only conditionally stable for this equation, and they can yield concentrations below zero or above one, even if they are stable. Here, we collect the stable and explicit algorithms, most of which we invented recently. All of them are unconditionally dynamically consistent for Fisher’s equation; thus, the concentration remains in the unit interval for arbitrary parameters. We perform tests in the cases of 1D and 2D systems to explore how the errors depend on the coefficient of the nonlinear term, the stiffness ratio, and the anisotropy of the system. We also measure running times and recommend which algorithms should be used in specific circumstances.