Numerical Study of One Dimensional Fishers KPP Equation with Finite Difference Schemes

In this paper, we originate results with finite difference schemes to approximate the solution of the classical Fisher Kolmogorov Petrovsky Piscounov (KPP) equation from population dynamics. Fisher’s equation describes a balance between linear diffusion and nonlinear reaction. Numerical example illustrates the efficiency of the proposed schemes, also the Neumann stability analysis reveals that our schemes are indeed stable under certain choices of the model and numerical parameters. Numerical comparisons with analytical solution are also discussed. Numerical results show that Crank Nicolson and Richardson extrapolation are very efficient and reliably numerical schemes for solving one dimension fisher’s KPP equation.


Introduction
Fisher gives introduction to nonlinear evolution equation to inquisitive, the proliferation of an beneficial gene in a population dynamics [1].Fisher's equation also specify the logistic diffusion process [1].It has the form ( ) where 0 β > is a diffusion constant with 0 α > is the linear growth rate [1].The reaction diffusion Equation (1) also express a model equation for the evolution of a neutron population in a nuclear reactor [2] and also appears in chemical engineering applications [2].This equation accommodates the effects of linear diffusion along xx u and nonlinear local multiplication or reaction along ( ) [4].It has become one of the most important nonlinear equations and occurs in many biological and chemical processes [4] [5].Recently, many researcher are working on this type of model to understand growth rate and diffusion aspect, for example, Abdullaev has studied the stability of symmetric travelling waves in the Cauchy problem for a more general case than Equation (1) [6].Also Logan has studied this problem using a perturbation method and settled up with an approximate solution by expanding the solution in terms of a power series and in terms of some small parameter [7].The finite difference schemes and auxiliary conditions of the numerical model must be consistent with the partial differential equations and initial and boundary conditions of the mathematical model [8] [9].The numerical model is consistent if the truncation error, that is the discrepancy between the finite difference approximation and the continuous derivatives, tends to zero as the grid spacing get smaller and smaller.Secondly the solution to finite difference schemes must converge to the solution of the partial differential equations as the grid spacing gets smaller and smaller [9].Thus we can say that, difference between the exact solution and approximated solution must vanish as the grid spacing tends to zero.
In this paper, we started the solution of Fishers equation with various finite difference schemes.We discussed Forward in time and Central in space (FTCS) and Lax-Wendroff schemes, which are explicit and conditionally stable.Also FTCS is first order accurate in time and second order accurate in space and Lax-Wendroff is second order accurate in both space and time, which shows improvement in accuracy in later method.Crank-Nicolson scheme is unconditionally stable and implicit with second order accuracy in both space and time [9].This applies computational stability for any size of the time increment [9].However, the size of t ∆ is still limited by accuracy required in the solution.
Usage of very large values of time step results in poor approximation to solution because of unacceptably large truncation errors produced [10].We applied Richardson extrapolation method to improve accuracy with no issue of stability.
We came to know this method is highly accurate and easy to implement with no mess of computations [10] [11] and we can see from our results that this method is excellent agreement to the exact solution.

Exact Solution
Let us consider Equation (1), within domain ( ) To derive the exact solution of the given Equation (1), we have the following exact solution [10] [11] to above equation, ( )

Numerical Methods
We consider the numerical solution of the non-linear Equation (2) in a finite domain Ω .The first step is to choose integers n to define step sizes U respectively.We consider four finite different schemes as we mentioned in keywords in abstract.

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 Equation (1), which leads to the following, ( ) ( ) , final form of above equation is, ( ) ( ) Since the one dimensional diffusion reaction equation is nonlinear and wellposed [10] [11], Lax's equivalence theorem indicates that consistency and stability of the FTCS finite difference approximation is necessary and sufficient for FD solution to converge to diffusion reaction equation [12].Once convergence has been proved, the solution to the given partial differential equation can be obtained to any desired degree of accuracy [12].Make sure the spacing h for spatial and k for time of the finite difference grid are made sufficiently small.The FTCS scheme, from Equation (3), is classified as explicit because the value of time level may be calculated directly from known value of n i u at previous time levels.It is a two level method because values of u at only two levels of time are involved in the approximating finite difference equation [12].
Accuracy and Stability to FTCS Scheme: Accuracy: To find accuracy of the FTCS scheme for Fisher-KPP equation, we apply Taylor's series on each term of the equation.Let us consider the Taylor's series in the following way, ( , , Apply Taylor s series 2 6 Take Equation ( 1) into account, and apply this scheme with Taylor's series on each term, updated equation is as follows.

Eq u x t k u x t u x h t u x t u x h t ku x t u x t
Now principle part of the truncation error is along with Equation (1).So first part of the above equation goes to zero if we consider Equation (1).
Which shows that this scheme is first order accurate in time and 2nd order accurate in space, such as ( )

Stability:
We want to study under what condition the error can be magnified.Many methods can be used to study this issue.We consider only Von-Neumann stability analysis to explain this method on FTCS scheme.Consider the scheme in the following way, ( ) .
According to Von-Neumann stability analysis, let us consider the solution as: .
The Von-neumann stability condition is Note that, Apply above terms to Equation (1), we get the following ( ) (7) According to Equation ( 6), Take left hand side of above equation along Equation ( 7), ( ) Take right hand side of above equation along Equation ( 7), ( ) ( ) , according to Von-Neumann stability analysis on both sides as left and right, which shows that FTCS scheme is conditionally stable for Fisher-KPP equation.

Lax Wendroff Scheme
A numerical technique proposed in 1960 by P.D. Lax and B. Wendroff [13] [14] [15] for solving partial differential equations and system numerically.In spite of the impressive developments on numerical methods for partial differential equations from 1970s onwards, in which the Lax Wendroff method has played a historic role, there are presently (1998) substantial research activities aimed at further improvements of methods [15].Lax Wendroff's method is also explicit method but needs improvement in accuracy in time.This method is an example of explicit time integration where the function that defines governing equation is evaluated at the current time [15].Purpose of this method to achieve good enough accuracy in time.This method can be explained in two steps which is as follows, ( ) ( ) ( ) Accuracy and Stability to Lax-Wendroff Scheme: Accuracy: To find accuracy of the Lax-Wendroff scheme for Fisher-KPP equation, we apply Taylor's series on each term of the discritized scheme.Let us consider the following,

(
) ( ) , , Apply Taylor s series 2 6 above mentioned scheme with Taylor's series, we have Apply above terms to Equation (1), we get the following According to above result, the Von-Neumann stability criterion 1 k e α ≤ is satisfied as long as 1 α ≤ , or equivalently, as long as the CFL condition is satisfied.In this respect, the dissipative properties of the Lax-Friedrichs scheme are not completely lost in the Lax-Wendroff scheme but are much less severe [15].The Lax-Wendroff scheme is a two-level scheme, but can be recast in a one-level form by means of algebraic manipulations [15].This is clear from expressions where quantities at time-levels n and n + 1 only appear [15].

Crank Nicolson Scheme
Let us apply an implicit method to improve the accuracy.The implicit method we consider, is Crank Nicolson.By substituting the value of n m u by average value of mesh points, the Crank Nicolson method is a finite difference method used for numerically solving linear and nonlinear partial differential equations [16] [17] [18].It is implicit in time and can be written in form of an implicit Runge Kutta method, and it is numerically stable [16] [17] [18].The method was developed by John Crank and Phyllis Nicolson in the mid 20th century [18].For diffusion equations (and many other equations), it can be shown that, the Crank Nicolson method is unconditionally stable [19].However, the approximate solutions can still contain (decaying) spurious oscillations if the ratio of time step k times to the square of space step 2 h , 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 [20].In this method, consider the followings,

Accuracy to Crank Nicolson Scheme:
To find accuracy of the CN scheme for Fisher-KPP equation,we apply Taylor's series on each term of the discritized scheme.Let us consider the following,

(
) ( ) . 0.5 2 0.5 0.5 0.5 0.25 2 2 Now principle part of the truncation error along with Equation (1).So first and second part of the above equation goes to zero if we consider Equation (1).
Which show that this scheme is second order accurate in time and space, such as , O k h .

Stability:
According to Von-Neumann stability analysis, we have linear form of the Equation (1) is, ( ) According to Von-Neumann stability analysis as ( 5), (6), after some simplifications, we reached at, ( ) ( ) Hence the Von-Neumann stability condition is satisfied for all positive values of 1 R .This shows that CN scheme is unconditionally stable.

Richardson Extrapolation Technique
Techniques for obtaining more accurate solutions by using finite difference methods which have already been discussed, including reducing the truncation error by using central differences, using methods which minimize numerical damping and dispersion as well as using more accurate approximation to derivative boundary conditions [21] [22].Other methods include, 1) reducing the size of grid spacing, 2) using higher order difference approximation, and 3) reducing the discretisation error by extrapolation (Richardson).Richardson extrapolation method, which was introduced by Richardson (1910) and later on called "deferred approach to the limit", may lead to considerable improvement of numerical results which solving the partial differential equation system by finite difference method.Applications of this method to equations solved on rectangular regions in space are described in Richardson and Gaunt (1927) and Salvadori (1951) [22].The reduction of the error depends upon there being a reliable estimate of the discretisation error as a function of the grid spacing h .Let φ be the exact solution of the partial differential equation and 1 , φ 2 φ the finite difference solutions at a particular point in the solution domain computed using grid spacing 1 2 , h h respectively [22] [23].Richardsons extrapolation formula is as follows, where ( ) j N h is the approximations which is even more accurate then the previous O h .This gives more accurate and precise result which lead to convergence [23].This method improvement the results upto fourth order accuracy in space.Keep in mind, accuracy in time still upto second order.Also this method hold the stable results.Clearly, this is a useful way of increasing the accuracy of the result at a given grid point without having to excessively increase the number of calculations required.This method may be of uncertain value near the boundaries which are not straight, such type of problems can be discussed with mixed boundaries with interpolation technique [23].

Results
Numerical computations have been performed by using the uniform grid [23].For the test problem, the approximated and exact values of ( ) , u x t and ( ) , U x t have been given in Figure 1 with some fixed parameters such as 1, and 0.0001 k α = = using FTCS.Figure 2 shows the results taken at different time to maintain stability.Figure 3 shows results by using Crank Nicolson scheme; comparison has been done with existence of exact solution.Figure 4 represents results from Richardson Extrapolation.These values has been taken at some typical grid point (TGP) as we presents in Table 1 and Table 2. Some critical values can be seen from the following data Table 1 and Table 2 at different location along x-axis.In Table 1, we compare our results, calculated from Crank-Nicolson scheme with exact solution.Error can be seen from last column.

Discussion
In this paper, the solution of the Fishers equation is successfully approximated by a various numerical finite difference schemes.Two of them are explicit in nature such as FTCS and Lax-Wendroff.We have to pay attention to parameter 1 R ,  which can stabilize the results as we can see from Figure 1 and Figure 2.Such schemes on nonlinear one dimensional PDE, consistency, stability and convergence are more difficult to prove [24] [25] [26].For instant, Von-Neumann's    method of stability analysis can not be used other than locally, since it is only applied to linear finite difference schemes.In many cases, numerical experimentation, such as solving the finite difference schemes by using progressively smaller grid spacing and examining the behaviour of the sequence of the values of ( ) , u x t obtained at given points, is the suitable method available with which to assess the numerical model [26].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 [26].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 exact solution.Our main purpose of this research to improve accuracy in result of Fisher KPP equation.At the end, we introduce Richardson extrapolation method; we observe improvement in accuracy and this method also stable and show excellent agreement with the exact solution.Accuracy in results are glanced from figures and tables.
of the truncation error along with Equation (1).So first part of the above equation goes to zero if we consider Equation (1).


this scheme is second order accurate in time and space, such as consider Von-Neumann stability analysis to explain on Lax-Wendroff scheme.Consider the scheme in the following way, According to Von-Neumann stability analysis as (5), (6) Note that, (

2 N
In our case we use 2 j = to get the approximation formula ( ) h with truncation error ( ) 4

Figure 1 .
Figure 1. Figure shows results with FTCS, by fixing some parameters to be in range of stability.Due to explicit nature of the scheme, we choose h and k in such a way that 1 r ≤ .

Figure 2 .
Figure 2. Figure shows results with Lax-Wendroff, by fixing some parameters to stabilize.Due to explicit nature of the scheme, we choose h and k in such a way that 1 r ≤ with acceleration in computations.Results at different time levels, also compared with exact solution.

Figure 3 .
Figure 3. Shows results with implicit Crank Nicolson, different α and time level along comparison with exact solution.

Figure 4 .
Figure 4. Shows results with Richardson Extrapolation, different time levels along comparison with exact solution.

Table 1 .
Estimates the results at different locations with Crank Nicolson scheme.

Table 2 .
Estimates results from Richardson-Extrapolation with exact solution.Error can be seen from last column.