Numerical Study of Fisher’s Equation by Finite Difference Schemes

Affiliation(s)

Department of Mathematics, Numerical Analysis, King Abdulaziz University, Jeddah, Saudi Arabia.

Department of Mathematics, Numerical Analysis, King Abdulaziz University, Jeddah, Saudi Arabia.

Abstract

This research paper represents a numerical approximation to three interesting equations of Fisher, which are linear, non-linear and coupled linear one dimensional reaction diffusion equations from population genetics. We studied accuracy in term of L∞ error norm by random selected grids along time levels for comparison with exact results. The test example demonstrates the accuracy, efficiency and versatility of the proposed schemes. It is shown that the numerical schemes give better solutions. Moreover, the schemes can be easily applied to a wide class of higher dimension non-linear reaction diffusion equations.

This research paper represents a numerical approximation to three interesting equations of Fisher, which are linear, non-linear and coupled linear one dimensional reaction diffusion equations from population genetics. We studied accuracy in term of L∞ error norm by random selected grids along time levels for comparison with exact results. The test example demonstrates the accuracy, efficiency and versatility of the proposed schemes. It is shown that the numerical schemes give better solutions. Moreover, the schemes can be easily applied to a wide class of higher dimension non-linear reaction diffusion equations.

Keywords

Forward in Time and Centre in Space (FTCS), Taylor’s Series, Crank Nicolson, Douglas Scheme

Forward in Time and Centre in Space (FTCS), Taylor’s Series, Crank Nicolson, Douglas Scheme

1. Introduction

Reaction diffusion equations arise as models for the densities of substances or organisms that disperse through space by Brownian motion, random walks, hydrodynamic turbulence, or similar mechanisms, and that react to each other and their surroundings in ways that affect their local densities [1] . Reaction diffusion models are in themselves deterministic, but they can be derived as limits of stochastic processes under suitable scaling. Specifically, they provide a modelling approach that allows us to translate assumptions about stochastic local movement into deterministic descriptions of global densities [1] [2] . Reaction diffusion models are spatially explicit, describe population densities, and treat space and time as continuous [1] [2] [3] . There are three major types of ecological phenomena that are supported by reaction diffusion equations: the existence of a minimal patch size necessary to support a population, the presence of travelling wave fronts corresponding to biological invasions, and the formation of spatial patterns [1] [2] [3] .

1.1. Diffusion

Diffusion is a description of movement that arises as a result of an object or organism making many short movements in random directions. The diffusive description of random motion emerges as a continuum limit of such random walks when the length $\Delta x$ of each step and the time $\Delta t$ required for each step go to zero in such a way that the ratio ${\left(\Delta x\right)}^{2}/\Delta t$ remains constant. To understand how this works it is useful to consider a simple example in one space dimension. Suppose that an organism moves along a line by moving a distance $\Delta x$ to the left with probability 1/2 or a distance $\Delta x$ to the right with probability 1/2 at each time step $\Delta t$ . Suppose that $\rho \left(x\mathrm{,}t\right)$ is the probability that the organism is at location x at time t. To arrive at that point at that time it must have been either one step to the left at time $t-\Delta t$ and then moved to the right, or one step to the right and have moved to the left. Thus, we have

$\rho \left(x,t\right)=\frac{1}{2}\rho \left(x+\Delta x,t-\Delta t\right)+\frac{1}{2}\rho \left(x-\Delta x,t-\Delta t\right)$ (1)

If we subtract $\rho \left(x\mathrm{,}t-\Delta t\right)$ from both sides and divide by $\Delta t$ in Equation (1), we obtain,

$\begin{array}{l}\frac{\rho \left(x,t\right)-\rho \left(x,t-\Delta t\right)}{\Delta t}\\ =\frac{1}{2\Delta t}\left[\rho \left(x+\Delta x,t-\Delta t\right)-2\rho \left(x,t-\Delta t\right)+\rho \left(x-\Delta x,t-\Delta t\right)\right]\end{array}$ (2)

Suppose that we now impose the diffusive scaling, ${\left(\Delta x\right)}^{2}/\Delta t=2D$ . Let us look at Equation (2),

$\begin{array}{l}\frac{\rho \left(x,t\right)-\rho \left(x,t-\Delta t\right)}{\Delta t}\\ =\frac{D}{{\left(\Delta x\right)}^{2}}\left[\rho \left(x+\Delta x,t-\Delta t\right)-2\rho \left(x,t-\Delta t\right)+\rho \left(x-\Delta x,t-\Delta t\right)\right]\end{array}$ (3)

From above Equation (3), the expression on the left is a difference quotient in t and also the expression on the right is a second difference in x. Taking the limit of expression in Equation (3), as $\left(\Delta x\mathrm{,}\Delta t\right)\to 0$ , while in Equation (2) remains in force yields the diffusion equation,

$\frac{\partial \rho}{\partial t}=D\frac{{\partial}^{2}\rho}{\partial {x}^{2}}.$ (4)

Mathematically this is identical to the heat equation. Note that the scaling, where D is the square of the distance $\Delta x$ moved by the organism in a time unit $\Delta t$ , produces a coefficient in front of the term ${\partial}^{2}\rho /\partial {x}^{2}$ , which is equal to 1/2 of the square of the distance moved per unit time. This interpretation of the diffusion coefficient D is valid in any number of dimensions.

1.2. Reaction

In the context of ecological models, the reaction terms in reaction diffusion equations and systems are typically the same as those that are used in non- spatial population models based on ordinary differential equations. Thus, for a single population, the reaction terms would be those that might occur in a model for a population density $\rho \left(t\right)$ of the form

$\frac{\partial \rho}{\partial t}=f\left(\rho \right),$ (5)

where $f\left(\rho \right)$ often has the form $f\left(\rho \right)=g\left(\rho \right)\rho $ . Common choices for $f\left(\rho \right)$ are $f\left(\rho \right)=R\rho $ (linear growth), $f\left(\rho \right)=R\rho \left(1-\rho /K\right)$ (logistic growth), or $f\left(\rho \right)=R\rho \left(\rho -a\right)\left(1-\rho /K\right)$ with $a\in \left(\mathrm{0,}K\right)$ (growth with Allee effect). For systems, typical reaction terms are those that occur in non-spatial models for competition, mutualism, or predator-prey interactions. Those include Lotka-Volterra models, but also more general models such as predator-prey models with a functional response. In the case of systems the stability analysis often involves the eigenvalues of matrices obtained by linearising about the equilibria. Equilibria and eigenvalues play a similar role in the analysis of reaction diffusion models, but the eigenvalues generally are associated with differential operators rather than matrices.

Here is the outline of the article. In Section 2, we mentioned literature review according to scope of the equation and numerical treatment also in section 3, we derived governing equation and its three interesting types and in section 4, methodology is explained. In the section 5, we discussed results in detail.

2. Literature Review

A well known researchers have studied such model problem, for example, Abdullaev [4] has studied the stability of symmetric travelling waves in the Cauchy problem for a more general case, also Logan has studied this problem using a perturbation method and found an approximate solution by expanding the solution in terms of a power series and in terms of some small parameters [5] , whereas numerical solution found by Gazdag and Canosa [6] [7] which exhibits consistency with partial differential equations along initial and boundary condition. Both numerical schemes intimated in [7] are totally com- plicated and source of unexpected high frequency oscillations, which must be refine at each time step. Tang and Weber [8] [9] have studied Fisher’s equation, using the Petrov-Galerkin method and Al-Khaled [10] has used the sinc col- location method to solve such model equation.

Recently, numerical solution to Fisher’s equation, have studied by many researchers, such as Wang [11] give idea of exact solution and explicit solitary wave solutions which are associated with generalized form, also Dag [12] found solution of Fisher’s equation numerically, using B-spline Galerkin method, whereas Qio and Sloan [13] built up numerical solutions of Fisher’s equation by moving mesh method, meanwhile Ting [13] studied to solve generalized Fisher’s equation by element free Galerkin method. Also modified cubic B-spline collocation method is used by Mittal and Jain [13] to study numerical solutions of non-linear Fisher’s equation. Fisher’s equation is studied numerically by Chandraker [13] , also Tomasiello studied numerical stability of differential quadrature solutions of wave problems. Korkmaz and Dag applied polynomial differential quadrature method to study numerical solutions of non-linear Burger’s equation [13] . Finite difference based methods have been applied by Kaysar [14] to solve Burger’s and Fisher’s equations numerically.

3. Governing Equation

In 1937 Fisher [15] and Kolmogorov et al. [16] investigated independently the Fisher Kolmogorov Petrovsky Piscounov (Fisher-KPP) equation, after that it is widely known as Fisher equation. This equation has many applications in science and engineering fields [6] [17] . The researchers studied some meaningful generalization of this equation, here we considered one generalization of this equation which is called as one component reaction diffusion equation. Many reaction diffusion equations have travelling wave fronts which play an important role in the understanding of physical, chemical, and biological phenomena [18] . Reaction-diffusion systems are mathematical models which explains how the concentration of one or more substances distributed in space changes under the influence of two processes, first one is local chemical reactions in which the substances are transformed into each other and second is the diffusion which causes the substances to spread out over a surface in space. Reaction-diffusion systems are naturally applied in chemistry. However, the system can also describe the dynamical processes of non-chemical nature. In this paper, we introduce the following three major Fisher’s equations, which can be explained as.

3.1. Linear Form of Fisher’s Equation

The linear form of Fisher’s equation is as follows,

${u}_{t}=\beta {u}_{xx}+\alpha \left(1-u\left(x,t\right)\right),$ (6)

where $\beta $ is diffusive constant with value $0\le \beta \le 1$ and $\alpha $ is reactive constant with value $0\le \alpha \le 1$ . Also analytical solution to above Equation (6) is,

$u\left(x,t\right)=1-\frac{\mathrm{cosh}x}{\mathrm{cosh}1}-\frac{16}{{\text{\pi}}^{2}}{\displaystyle \underset{n=1}{\overset{\infty}{\sum}}}\frac{{\left(-1\right)}^{n}\mathrm{cos}\left(2n-1\right)\text{\pi}x/2}{\left(2n-1\right)\left[{\left(2n-1\right)}^{2}{\text{\pi}}^{2}+4\right]}{\text{e}}^{-\left(1+{\left(2n-1\right)}^{2}{\text{\pi}}^{2}/4\right)t},$ (7)

with boundary conditions are

$u\left(-1,t\right)=u\left(1,t\right)=0,$ (8)

and initial condition also,

$u\left(x,0\right)=0.$ (9)

3.2. Coupled Linear System

The coupled linear system is as follows,

$\begin{array}{l}{u}_{t}={u}_{xx}+\left(u\left(x,t\right)-v\left(x,t\right)\right)+F\left(x,t\right)\\ {v}_{t}={v}_{xx}+\left(u\left(x,t\right)+v\left(x,t\right)\right)+G\left(x,t\right)\end{array}\}$ (10)

with analytical solution,

$\begin{array}{l}u\left(x,t\right)={\text{e}}^{t}\mathrm{sin}x\\ v\left(x,t\right)={\text{e}}^{t}\mathrm{cos}x\end{array}\}$ (11)

3.3. Nonlinear Generalized Fisher’s Equation

The generalized form of nonlinear Fisher’s equation is as follows,

${u}_{t}={u}_{xx}+u\left(1-u\right)\left(u-{\alpha}_{1}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<{\alpha}_{1}<1$ (12)

with analytical solution,

$u\left(x,t\right)=\frac{1}{2}\left(1+{\alpha}_{1}\right)+\left(\frac{1}{2}-\frac{1}{2}{\alpha}_{1}\right)\mathrm{tanh}\left[\sqrt{2}\left(1-{\alpha}_{1}\right)\frac{x}{4}+\frac{\left(1-{\alpha}_{1}^{2}\right)}{4}t\right],$ (13)

the initial and boundary conditions are taken from the exact solution (13).

4. Numerical Methods

Let us apply numerical methods like Finite Difference Schemes (Forward in time and central in space (FTCS), Crank Nicolson (CN) and Douglas), to solve such Equations ((6), (10), (12)) in finite domain $\Omega =\left[0,1\right]$ . We partitioned the interval $\left[a\mathrm{,}b\right]$ into n equal parts of width h. Place a grid on the rectangle region R by drawing vertical and horizontal lines through the points with coordinates ${x}_{i}$ , where ${x}_{i}=a+ih$ for each $i=0,1,2,\cdots ,n$ also the lines $x={x}_{i}$ represent grid lines, we assume ${t}_{n}=nt,n=0,1,\cdots $ where n is the time grid step size. We denote the exact and numerical solutions at the grid point $\left({x}_{m}\mathrm{,}{t}_{n}\right)$ by ${U}_{m}^{n}$ and ${u}_{m}^{n}$ respectively.

4.1. Forward in Time and Center in Space (FTCS) Scheme

We consider forward in time and center in space (FTCS) explicit scheme by substituting the forward difference approximation for the time derivative and the central difference approximation for the space derivative in Equations ((6), (10), (12)) respectively, we get the following

${u}_{i}^{n+1}={u}_{i}^{n}+R\left({u}_{i+1}^{n}-2{u}_{i}^{n}+{u}_{i-1}^{n}\right)+Q\left(1-{u}_{i}^{n}\right)\mathrm{,}$ (14)

where $R=\frac{k\beta}{{h}^{2}}$ , and $Q=k\alpha $ . Above equation (14) represents descritezation to linear form of Fisher’s equation.

$\begin{array}{l}{u}_{i}^{n+1}={u}_{i}^{n}+R\left({u}_{i+1}^{n}-2{u}_{i}^{n}+{u}_{i-1}^{n}\right)+k{u}_{i}^{n}-k{v}_{i}^{n}+kF\left({u}_{i}^{n}\mathrm{,}{v}_{i}^{n}\mathrm{,}x\mathrm{,}t\right)\\ {v}_{i}^{n+1}={v}_{i}^{n}+R\left({v}_{i+1}^{n}-2{v}_{i}^{n}+{v}_{i-1}^{n}\right)+k{u}_{i}^{n}+k{v}_{i}^{n}+kG\left({u}_{i}^{n}\mathrm{,}{v}_{i}^{n}\mathrm{,}x\mathrm{,}t\right)\end{array}\}$ (15)

above Equation (15) represents descritezation to coupled linear system.

${u}_{i}^{n+1}={u}_{i}^{n}+R\left({u}_{i+1}^{n}-2{u}_{i}^{n}+{u}_{i-1}^{n}\right)+k{u}_{i}^{n}\left(1-{u}_{i}^{n}\right)\left({u}_{i}^{n}-{\alpha}_{1}\right)\mathrm{,}$ (16)

above Equation (16) represents descritezation to nonlinear generalized form of Fisher’s equation.

Since the one dimensional Fisher’s equation or system is well posed, make sure the spacing h for spatial and k for time of the finite difference grid are made sufficiently small [6] [17] [19] . The FTCS scheme, from Equations (14)-(16), is classified as explicit because the value of ${u}_{i}^{n+1}$ at the $\left(n+1\right)th$ time level may be calculated directly from known value of ${u}_{i}^{n}$ at previous time levels. It is a two level method because values of $u\left(x\mathrm{,}t\right)$ at only two levels of time are involved in the approximating finite difference equation [6] [17] [19] . There is no best method for obtaining approximating difference formula, the only requirement is that the formula, having been obtained, must pass certain tests of accuracy, consistency, stability and convergence [20] [21] [22] . By Von Neumann stability

analysis, the FTCS scheme is always conditionally stable, which is $0<R\le \frac{1}{2}$ .

4.2. Crank Nicolson Implicit Scheme

Let us apply implicit finite difference scheme, which is Crank Nicolson. This method uses central finite difference approximation for both time and space derivatives at the point $\left({x}_{m}\mathrm{,}{t}_{n}\right)$ [6] . For diffusion equations (and many other equations), it can be shown that, the Crank Nicolson method is unconditionally stable [23] [24] [25] . However, the approximate solutions can still contain (decaying) spurious oscillations if the ratio of time step k times to the square of space step ${h}^{2}$ , is large (typically larger than 1/2 per Von-Neumann stability analysis) [21] [26] [27] . For this reason, whenever large time steps or high spatial resolution is necessary, the less accurate backward Euler method is often used, which is both stable and immune to oscillations [24] [25] [27] [28] [29] [30] . In this method, we consider the Equations ((6), (10), (12)) respectively, in the following way,

$\begin{array}{c}\frac{{u}_{i}^{n+1}-{u}_{i}^{n}}{k}=\frac{k\beta}{2{h}^{2}}{\delta}_{x}^{2}\left[{u}_{i}^{n+1}+{u}_{i}^{n}\right]+\alpha \left(1-\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-{R}_{1}{u}_{i+1}^{n+1}+\left(1+2{R}_{1}+0.5{Q}_{1}\right){u}_{i}^{n+1}-{R}_{1}{u}_{i-1}^{n+1}\\ ={R}_{1}{u}_{i+1}^{n}+\left(1-2{R}_{1}-0.5{Q}_{1}\right){u}_{i}^{n}+{R}_{1}{u}_{i-1}^{n}+{Q}_{1}\end{array}$ (17)

where ${R}_{1}=\frac{k\beta}{2{h}^{2}}$ and ${Q}_{1}=k\alpha $ . Above Equation (17) represents descritezation

using Crank Nicolson to linear form of Fisher’s equation. Now let us look at coupled linear system, in the following way,

$\begin{array}{l}-{R}_{1}{u}_{i+1}^{n+1}+\left(1+2{R}_{1}-0.5k\right){u}_{i}^{n+1}+0.5k{v}_{i}^{n+1}-{R}_{1}{u}_{i-1}^{n+1}\\ ={R}_{1}{u}_{i+1}^{n}+\left(1-2{R}_{1}+0.5k\right){u}_{i}^{n}-0.5k{v}_{i}^{n}+kF\left(\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\mathrm{,}\frac{1}{2}\left({v}_{i}^{n+1}+{v}_{i}^{n}\right)\mathrm{,}x\mathrm{,}t\right)\\ -{R}_{1}{v}_{i+1}^{n+1}+\left(1+2{R}_{1}-0.5k\right){v}_{i}^{n+1}-0.5k{u}_{i}^{n+1}-{R}_{1}{v}_{i-1}^{n+1}\\ ={R}_{1}{v}_{i+1}^{n}+\left(1-2{R}_{1}+0.5k\right){v}_{i}^{n}+0.5k{u}_{i}^{n}+kG\left(\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\mathrm{,}\frac{1}{2}\left({v}_{i}^{n+1}+{v}_{i}^{n}\right)\mathrm{,}x\mathrm{,}t\right)\end{array}\}$ (18)

Above equation (18) represents descritezation using Crank Nicolson to coupled linear system. Now let us look at generalized nonlinear Fisher’s equation using Crank Nicolson,

$\begin{array}{l}-{R}_{1}{u}_{i+1}^{n+1}+\left(1+2{R}_{1}\right){u}_{i}^{n+1}-{R}_{1}{u}_{i-1}^{n+1}\\ ={R}_{1}{u}_{i+1}^{n}+\left(1-2{R}_{1}\right){u}_{i}^{n}+\frac{k}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\left(1-0.5\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\right)\left(0.5\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)-{\alpha}_{1}\right)\end{array}$ (19)

4.3. Fourth Order Accurate Implicit Scheme

Let us apply another implicit scheme to Equations ((6), (12)) in an order respectively.

$\begin{array}{l}\frac{{u}_{i}^{n+1}-{u}_{i}^{n}}{k}=\frac{\beta}{2{h}^{2}}{\left[1+{\delta}_{x}^{2}\right]}^{-1}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)+\alpha \left(1-\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\right)\\ \left[1+{\delta}_{x}^{2}\right]\left(\frac{{u}_{i}^{n+1}-{u}_{i}^{n}}{k}\right)=\frac{\beta}{2{h}^{2}}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)+\left[1+{\delta}_{x}^{2}\right]\alpha \left(1-\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\right)\\ \left({u}_{i}^{n+1}-{u}_{i}^{n}\right)+{\delta}_{x}^{2}\left({u}_{i}^{n+1}-{u}_{i}^{n}\right)={R}_{1}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)+\left[1+{\delta}_{x}^{2}\right]\alpha \left(1-\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\right)\end{array}\}$ (20)

Now Douglas scheme to nonlinear generalized Fisher’s equation,

$\begin{array}{l}\frac{{u}_{i}^{n+1}-{u}_{i}^{n}}{k}=\frac{\beta}{2{h}^{2}}{\left[1+{\delta}_{x}^{2}\right]}^{-1}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)+\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\left(1-\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\right)\left(\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)-{\alpha}_{1}\right)\\ \left[1+{\delta}_{x}^{2}\right]\left(\frac{{u}_{i}^{n+1}-{u}_{i}^{n}}{k}\right)=\frac{\beta}{2{h}^{2}}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)+\left[1+{\delta}_{x}^{2}\right]\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\left(1-\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\right)\left(\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)-{\alpha}_{1}\right)\\ \left[1+{\delta}_{x}^{2}\right]\left({u}_{i}^{n+1}-{u}_{i}^{n}\right)={R}_{1}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)+\left[1+{\delta}_{x}^{2}\right]\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\left(1-\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)\right)\left(\frac{1}{2}\left({u}_{i}^{n+1}+{u}_{i}^{n}\right)-{\alpha}_{1}\right)\end{array}\}$ (21)

5. Error Norms

The aim of the accuracy is assessed by some redefined norms, associated with the consistency of the finite difference schemes, such scaled measurement to error defined in term of norms specially ${L}_{\infty}$ , which is outlined below:

${L}_{\infty}={\mathrm{max}}_{i}\left|{u}_{i}^{\text{Exact}}-{u}_{i}^{\text{Approximation}}\right|$ (22)

6. Results

Numerical computations have been performed using the uniform grid. We used FTCS, Crank Nicolson and Douglas finite difference schemes to analyse numerical behaviour of simple linear Fisher;s equation, one dimensional linear coupled system and non-linear Fisher’s equation respectively. First we look at the linear Fisher’s equation by finite difference schemes as in Table 1, we used FTCS explicit scheme with some variations in grid size and h (space step) is changed according to the grid sizes. this table explains the second order accuracy in term of ${L}_{\infty}$ norm, of the explicit numerical scheme also Table 2 explains results for FTCS with different time steps (k). As we know that FTCS scheme is always conditionally stable and we can see from Table 2. In Table 3, we used implicit Crank Nicolson finite difference scheme, which shows that results with $k=0.0001,t=1,\text{interval}=\left[-1,1\right]$ and different grid sizes with h changes accordingly also Table 4 explains the method with different k (time steps). Table 5 shows results using Douglas scheme with $k=0.0001,t=1,\text{interval}=\left[-1,1\right]$ and different grid sizes with h changes accordingly also Table 6 shows results using Douglas scheme with $\text{Grid}=71\times 71,t=1,\text{interval}=\left[-1,1\right]$ and different k. In Table 7 we represent results for linear Fisher’s equation with comparison of two implicit schemes and Douglas improves and encourages our solution.

Secondly, we look at the coupled linear system by finite difference schemes as in Table 8 and Table 9, we used FTCS explicit scheme with some variations in grid size and h (space step) is changed according to the grid sizes. these tables explain the second order accuracy in term of ${L}_{\infty}$ norm and also classical simple error as $\text{Error}=\left|U\left(x,t\right)-u\left(x,t\right)\right|$ , by using both explicit and implicit schemes.

Lastly, we look at the generalized Fisher’s equation by finite difference schemes as in Table 10, we used FTCS explicit scheme with some variations in grid size and h (space step) is changed according to the grid sizes. this table explains the second order accuracy in term of ${L}_{\infty}$ norm, of the explicit numerical scheme. Table 11 shows results for Crank Nicolson to generalized non-linear Fisher’s equation along Table 12 shows Douglas results at different grid sizes.

Table 1. This table shows results using FTCS explicit scheme with $k=0.0001,t=1,$ $\text{interval}=\left[-1,1\right]$ and different grid sizes with h changes accordingly.

Table 2. This table shows results using FTCS explicit scheme with $\text{Grid}=71\times 71,t=1,$ Math_85# and different k.

Table 3. This table shows results using Crank Nicolson scheme with $k=0.0001,t=1,$ Math_89# and different grid sizes with h changes accordingly.

Table 4. This table shows results using Crank Nicolson scheme with $\text{Grid}=71\times 71,t=1,$ Math_96# and different k.

Table 5. This table shows results using Douglas scheme with $k=0.0001,t=1,$ Math_102# and different grid sizes with h changes accordingly.

Table 6. This table shows results using Douglas scheme with $\text{Grid}=71\times 71,t=1,$ Math_109# and different k.

Table 7. This table shows comparison between Crank Nicolson and Douglas schemes with $k=0.0001,t=1$ and $\text{interval}=\left[-\mathrm{1,1}\right].$

Table 8. This table shows results using FTCS explicit scheme with $k=0.0001,t=0.1$ and $\text{interval}=\left[-\mathrm{3,3}\right]$ .

Table 9. This table shows results using Crank Nicolson implicit scheme with $k=0.00001,t=0.001$ and $\text{interval}=\left[-\mathrm{10,10}\right]$ .

Table 10. This Table shows results using FTCS explicit scheme with $k=0.0001,t=0.1,$ Math_165#, and different grid sizes with h changes accordingly.

Table 11. This Table shows results using Crank Nicolson scheme with $k=0.0001,t=0.1,$ Math_168# and different grid sizes with h changes accordingly.

Table 12. This Table shows results using Douglas scheme with $t=0.1$ and $\text{interval}=$ Math_176#.

To analyse the graphic representation to linear Fisher’s equation, we have Figure 1, by FTCS along Figures 2-4 by Crank Nicolson and by Douglas, Figure 5 enhance our knowledge. Figure 6 & Figure 7 show results for coupled linear system by FTCS and Crank Nicolson. Also Figure 8 & Figure 9 show results for non-linear generalized Fisher’s equation by Crank Nicolson and Douglas schemes respectively.

7. Conclusion

In this paper, the solution to linear form of the Fishers equation, coupled linear system and generalized Fisher’s equation is successfully approximated by a various numerical finite difference schemes. Two of them are implicit in nature such as Crank Nicolson and Douglas and one is explicit FTCS schemes. We have to pay attention to parameter $R$ , which can stabilize the results as we can see from figures and tables. For instant, Von-Neumann’s method of stability analysis can not be used other than locally, since it only applies to linear finite

Figure 1. Results for different time level, with fixed grid size and time step, using FTCS scheme.

Figure 2. Results for different time level, with fixed grid size and time step, using Crank Nicolson scheme after zoom out.

Figure 3. Results for different grid sizes with time step, and find common region of interaction, using Crank Nicolson scheme.

Figure 4. Results for Douglas scheme and error profile.

Figure 5. Results for Douglas scheme and error profile with common region.

Figure 6. Results for FTCS scheme on Coupled Linear System.

Figure 7. Results for Douglas scheme and error profile with common region.

Figure 8. Results for Crank Nicolson scheme and error profile on Non-linear Generalized Fisher’s equation.

difference schemes. In many cases, numerical experimentation, such as solving the finite difference schemes using progressively smaller grid spacing and examining the behaviour of the sequence of the values of $u\left(x\mathrm{,}t\right)$ obtained at given points, is the suitable method available with which to assess the numerical model. The various methods of obtaining a finite difference numerical model corresponding to a particular mathematical model may result in either explicit or implicit finite difference schemes. Explicit schemes are conditionally stable and implicit schemes are unconditionally stable. Two implicit schemes are also applied to improve accuracy, stability restrictions and consistency in solution. It can be observed that the computed results show excellent agreement with the analytical solution. Our main purpose of this research is to improve accuracy in result. Accuracy in results is glanced from figures and tables.

Figure 9. Results for Douglas scheme and error profile on Non-linear Generalized Fisher’s equation.

Acknowledgements

Bader Saad Alshammari and Prof. Daoud Mashat are very thankful to Dr Muhammad Faheem Afzaal, Department of Chemical Engineering, Imperial College London and Vineet K. Srivastava, Scientist, ISTRAC/ISRO, Bangalore, India for thoughtful remarks. This research was supported by Department of Mathematics, division of Numerical Analysis, King Abdulaziz University, Jeddah, Saudi Arabia.

Conflict of Interest

There is no conflict of interest in this research paper.

Cite this paper

Alshammari, B. and Mashat, D. (2017) Numerical Study of Fisher’s Equation by Finite Difference Schemes.*Applied Mathematics*, **8**, 1100-1116. doi: 10.4236/am.2017.88083.

Alshammari, B. and Mashat, D. (2017) Numerical Study of Fisher’s Equation by Finite Difference Schemes.

References

[1] Kristiansen, K.K. (2008) Reaction Diffusion Models in Mathematical Biology. Technical University of Denmark, the Department of Mathematics, Lyngby.

[2] Cosner, C.C. (2014) Reaction Diffusion Equations and Ecological Modelling. Department of Mathematics, University of Miami, Coral Gables.

[3] Li, J.L. and Chen, Y.T. (2008) Computational Partial Differential Equations Using MATLAB. University of Nevada, Las Vegas, Paradise, 79-120.

[4] Abdullaev, U.G. (1994) Stability of Symmetric Travelling Waves in the Cauchy Problem for the KPP Equation. Differential Equations, 30, 377-386.

[5] Logan, D.J. (1984) An Introduction to Nonlinear Partial Differential Equations. John Wiley & Sons, New York.

[6] Canosa, J.C. (1973) On a Nonlinear Diffusion Equation Describing Population Growth. IBM Journal of Research Development, 17, 307-313.
https://doi.org/10.1147/rd.174.0307

[7] Gazdag, J.G. and Canosa, J.C. (1974) Numerical Solutions of Fisher’s Equation. Journal of Applied Probability, 11, 445-457.
https://doi.org/10.1017/S0021900200096236

[8] Tang, S.T., Qin, S.Q. and Weber, R.O. (1991) Numerical Solution of a Non Linear Reaction Diffusion Equation. Chinese J. of Appl. Math, and Mech., 12, 703-709.

[9] Tang, S.T. and Weber, R.O. (1991) Numerical Study of Fisher’s Equation by a Petrov-Galerkin Finite Element Method. Journal of the Australian Mathematical Society Series B, 33, 27-38. https://doi.org/10.1017/S0334270000008602

[10] Khaled, K.A. (2001) Numerical Study of Fisher’s Diffusion Reaction Equation by the Sinc Collocation Method. Journal of Computational and Applied Mathematics, 137, 245-255. https://doi.org/10.1016/S0377-0427(01)00356-9

[11] Wang, X.Y. (1988) Exact and Explicit Solitary Wave Solutions for the Generalized Fishers Equation. Physics Letters A, 131, 227-279.
https://doi.org/10.1016/0375-9601(88)90027-8

[12] Dag, I.D., Sahin, A.S. and Korkmaz, A.K. (2010) Numerical Investigation of the Solution of Fishers Equation via the B-Spline Galerkin Method. Numerical Methods for Partial Differential Equations, 26, 1483-1503.

[13] Mittal, R.C. and Rajni, R.R. (2016) A Study of One Dimensional Nonlinear Diffusion Equations by Bernstein Polynomial Based Differential Quadrature Method. Journal of Mathematical Chemistry, 2016, 1-23.

[14] Rahman, K.R., Helil, N.H. and Yimin, A.R. (2010) Some New Semi-Implicit Finite Difference Schemes for Numerical Solution of Burgers Equation. International Conference on Computer Application and System Modeling, 14, 451-455.
https://doi.org/10.1109/ICCASM.2010.5622119

[15] Fisher, R.A. (1936) The Wave of Advance of Advantageous Genes. Annals of Eugenics, 7, 355-369. https://doi.org/10.1111/j.1469-1809.1937.tb02153.x

[16] Kolmogorov, A.K., Petrovsky, N.P. and Piscounov, S.P. (1937) Etude de I Equations de la Diffusion avec Croissance de la Quantitate de Matiere et Son Application a un Probolome Biologique. Bulletin University Moscow, 1, 1-25.

[17] Newman, W.I. (1980) Some Exact Solutions to a Non-Linear Diffusion Problem in Population Genetics and Combustion. Journal of Theoretical Biology, 85, 325-334.

https://doi.org/10.1016/0022-5193(80)90024-7

[18] Tuckwell, H.C. (1998) Introduction to Theoretical Neurobiology. Cambridge University Press, Cambridge.

[19] Arnold, R.A., Showalter, K. and Tyson, J.J. (1987) Propagation of Chemical Reactions in Space. Journal of Chemical Education, 64, 740-744.
https://doi.org/10.1021/ed064p740

[20] Argyris, J.A., Haase, M.H. and Heinrich, J.C. (1991) Finite Approximation to Two- Dimensional Sine Gordon Equations. Computer Methods in Applied Mechanics and Engineering, 86, 1-26. https://doi.org/10.1016/0045-7825(91)90136-T

[21] Grimshaw, R.G. and Tang, S.T. (1990) The Rotation-Modified Kadomtsev Petviashvili Equation: An Analytical and Numerical Study. Studies in Applied Mathematics, 83, 223-248. https://doi.org/10.1002/sapm1990833223

[22] Ames, W.F. (1969) Finite Difference Methods for Partial Differential Equations. Academic Press, New York.

[23] Ablowitz, M.J. and Zeppetella, A.Z. (1979) Explicit Solutions of Fisher’s Equation for a Special Wave Speed. Bulletin of Mathematical Biology, 41, 835-840.

https://doi.org/10.1007/BF02462380

[24] Kanti, P.K. and Lajja, V.L. (2011) A Note on Crank-Nicolson Scheme for Burgers Equation. Applied Mathematics, 2, 888-899.

[25] Srivastava, V.K. and Tamsir, M.T. (2012) Crank Nicolson Semi Implicit Approach for Numerical Solution of Two Dimensional Coupled Nonlinear Burgers Equations. International Journal of Applied Mechanics and Engineering, 17, 571-581.

[26] Aronson, D.G. and Weinberger, H.F. (1978) Multidimensional Non-Linear Diffusion Arising in Population Genetics. Advance in Mathematics, 30, 33-76.
https://doi.org/10.1016/0001-8708(78)90130-5

[27] Mazumder, S.M. (2015) Numerical Methods for Partial Differential Equations: Finite Difference and Finite Volume Methods. Academic Press, New York.

[28] Williams, S.W. and Chow, P.L. (1978) Nonlinear Reaction-Diffusion Models for Interacting Populations. Journal of Mathematical Analysis and Applications, 62, 157-169.

https://doi.org/10.1016/0022-247X(78)90227-5

[29] Whitham, G.B. (1974) Linear and Nonlinear Waves. John Wiley Sons, New York.

[30] Babuska, I.B. (1968) Numerical Stability in Mathematical Analysis. IFIP Congress, Amsterdam, 11-23.