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  . 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   . Reaction diffusion models are spatially explicit, describe population densities, and treat space and time as continuous    . 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    .
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 of each step and the time required for each step go to zero in such a way that the ratio 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 to the left with probability 1/2 or a distance to the right with probability 1/2 at each time step . Suppose that 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 and then moved to the right, or one step to the right and have moved to the left. Thus, we have
If we subtract from both sides and divide by in Equation (1), we obtain,
Suppose that we now impose the diffusive scaling, . Let us look at Equation (2),
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 , while in Equation (2) remains in force yields the diffusion equation,
Mathematically this is identical to the heat equation. Note that the scaling, where D is the square of the distance moved by the organism in a time unit , produces a coefficient in front of the term , 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.
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 of the form
where often has the form . Common choices for are (linear growth), (logistic growth), or with (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  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  , whereas numerical solution found by Gazdag and Canosa   which exhibits consistency with partial differential equations along initial and boundary condition. Both numerical schemes intimated in  are totally com- plicated and source of unexpected high frequency oscillations, which must be refine at each time step. Tang and Weber   have studied Fisher’s equation, using the Petrov-Galerkin method and Al-Khaled  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  give idea of exact solution and explicit solitary wave solutions which are associated with generalized form, also Dag  found solution of Fisher’s equation numerically, using B-spline Galerkin method, whereas Qio and Sloan  built up numerical solutions of Fisher’s equation by moving mesh method, meanwhile Ting  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  to study numerical solutions of non-linear Fisher’s equation. Fisher’s equation is studied numerically by Chandraker  , 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  . Finite difference based methods have been applied by Kaysar  to solve Burger’s and Fisher’s equations numerically.
3. Governing Equation
In 1937 Fisher  and Kolmogorov et al.  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   . 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  . 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,
where is diffusive constant with value and is reactive constant with value . Also analytical solution to above Equation (6) is,
with boundary conditions are
and initial condition also,
3.2. Coupled Linear System
The coupled linear system is as follows,
with analytical solution,
3.3. Nonlinear Generalized Fisher’s Equation
The generalized form of nonlinear Fisher’s equation is as follows,
with analytical solution,
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 . We partitioned the interval 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 , where for each also the lines represent grid lines, we assume where n is the time grid step size. We denote the exact and numerical solutions at the grid point by and 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
where , and . Above equation (14) represents descritezation to linear form of Fisher’s equation.
above Equation (15) represents descritezation to coupled linear system.
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    . The FTCS scheme, from Equations (14)-(16), is classified as explicit because the value of at the time level may be calculated directly from known value of at previous time levels. It is a two level method because values of at only two levels of time are involved in the approximating finite difference equation    . 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    . By Von Neumann stability
analysis, the FTCS scheme is always conditionally stable, which is .
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  . For diffusion equations (and many other equations), it can be shown that, the Crank Nicolson method is unconditionally stable    . However, the approximate solutions can still contain (decaying) spurious oscillations if the ratio of time step k times to the square of space step , is large (typically larger than 1/2 per Von-Neumann stability analysis)    . 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       . In this method, we consider the Equations ((6), (10), (12)) respectively, in the following way,
where and . 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,
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,
4.3. Fourth Order Accurate Implicit Scheme
Let us apply another implicit scheme to Equations ((6), (12)) in an order respectively.
Now Douglas scheme to nonlinear generalized Fisher’s equation,
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 , which is outlined below:
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 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 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 and different grid sizes with h changes accordingly also Table 6 shows results using Douglas scheme with 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 norm and also classical simple error as , 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 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 and different grid sizes with h changes accordingly.
Table 2. This table shows results using FTCS explicit scheme with Math_85# and different k.
Table 3. This table shows results using Crank Nicolson scheme with Math_89# and different grid sizes with h changes accordingly.
Table 4. This table shows results using Crank Nicolson scheme with Math_96# and different k.
Table 5. This table shows results using Douglas scheme with Math_102# and different grid sizes with h changes accordingly.
Table 6. This table shows results using Douglas scheme with Math_109# and different k.
Table 7. This table shows comparison between Crank Nicolson and Douglas schemes with and
Table 8. This table shows results using FTCS explicit scheme with and .
Table 9. This table shows results using Crank Nicolson implicit scheme with and .
Table 10. This Table shows results using FTCS explicit scheme with Math_165#, and different grid sizes with h changes accordingly.
Table 11. This Table shows results using Crank Nicolson scheme with Math_168# and different grid sizes with h changes accordingly.
Table 12. This Table shows results using Douglas scheme with and 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.
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 , 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 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.
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.
 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
 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
 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.
 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.
 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
 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.
 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.
 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
 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
 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.
 Williams, S.W. and Chow, P.L. (1978) Nonlinear Reaction-Diffusion Models for Interacting Populations. Journal of Mathematical Analysis and Applications, 62, 157-169.