Numerical solution of Fisher's equation

Abstract

The accurate space derivative (ASD) method for the numerical treatment of nonlinear partial differential equations has been applied to the solution of Fisher's equation, a nonlinear diffusion equation describing the rate of advance of a new advantageous gene, and which is also related to certain water waves and plasma shock waves described by the Korteweg-de-Vries-Burgers equation. The numerical experiments performed indicate how from a variety of initial conditions, (including a step function, and a wave with local perturbation) the concentration of advantageous gene evolves into the travelling wave of minimal speed. For an initial superspeed wave this evolution depends on the cutting off of the right-hand tail of the wave, which is physically plausible; this condition is necessary for the convergence of the ASD method. Detailed comparisons with an analytic solution for the travelling waves illustrate the striking accuracy of the ASD method for other than very small values of the concentration.