Numerical Approximation to Nonlinear One Dimensional Coupled Reaction Diffusion System

This research paper represents a numerical approximation to non-linear coupled one dimension reaction diffusion system, which includes the existence and uniqueness of the time dependent solution with upper and lower bounds of the solution. Also numerical approximation is obtained by finite difference schemes to reach at reasonable level of accuracy, which is magnified by 2 L , L∞ and relative error norms. The accuracy of the approximations is shown by randomly selected grid points along time level and comparison with analytical results. The test example demonstrates the accuracy, efficiency and versatility of the proposed schemes. Moreover, the schemes can be easily applied to a wide class of higher dimension non-linear reaction diffusion equations with a little modifications.


Introduction
Fisher and Kolmogorov Petrovsky Piscounov founded quasilinear partial differential equation which represents reaction diffusion phenomena [1] [2], which play important rule in population dynamics as propagation of beneficial genes; mathematical description to this phenomena is as follows [1] [2] [3], ( ), here, x declares spatial coordinate position with time t, also ( ) , x t ρ is the term used to identify the population density, i is the flux associated to population, and ( ) Q ρ is a limited to originator term that illuminate the net rate of production or growth in the populated area or density associated to population [3].The behaviour of the diffusion is accounted in the flux i, inclines by Fick's law, where the assumption on the diffusion coefficient D, is to be constant [3].It is frequently to imitate the law, which is known as the Pearl-Verhulst logistic law, ( ) related to stranded scale to growth in the population [3].This law explains the that initially population would grow with homogeneous fashion with time of 1 γ − until that growth reach at some stranded level with some saturation in population [3].Fisher and KPP were very much keen about the speed that exists in advantageous genes to propagate and they found a travelling wave solution to Equations ( 1), ( 2), (3) [3], such solution related to velocity ( ) Kametaka (1976) and Uchiyama (1977)] remarks about velocity of ( ) 2 Dγ developed initial conditions (ICs), categorised as [3], , 0 , for 0 , 0 , for , 0 0, for

Governing Equation
In the year of 1948 El' dovich, raised the major features of the diffusion reaction semi-linear equation which make useful contribution to wave of advantageous genes which are now totally classify as in the dynamics of the gas and flame with chemical kinetics [3].Mostly, in process of combustion, the concept of heat which is channelized by the conduction processed by heat related to molecules which moves from the hot part of flame to coolant side [3], by a particular locality, the rate of the chemical reaction such as exothermic chemical reaction (due to presence of heat), totally aware of temperature of the related gas mixture [3].The coefficients associated with molecular diffusion and thermal are at same level of scaled, the temperature and concentration are situated almost isobaric [3].The quantity ρ , then can represent either the temperature or concentration, without confusion, and is often referred to as the "combustion completeness by such coefficient" by Vulis, 1961 [3].Equations ( 1)-( 4) are then hired to modelled this attractive phenomena as a result of travelling wave nature of the flame exists in combustion process or zone, i.e. the transition layer between the region of burnt gas and the gas where combustion has not yet begun [3].Above literature review, leads to the following reaction-diffusion system in two components ( )( ) ( )( ) where in above Equation ( 5), 1 α and 2 α are the linear growth rate whereas

Exact Solution
Exact solution to system in Equation (5) found by E. S. Fahmy [5], he used factorization method, given below:

Smoothness and Uniqueness of the Reaction Diffusion System
It is very important to enforce some assumptions on some physical parameters or related reaction functions, to analysed the smoothness and uniqueness of a positive oriented natural solution by considering lower and upper bound of the solution of the system which mentioned in Equation ( 5) [6].We assume that 1 2 = β β which represents the diffusion coefficients are always positive in the bounded or finite domain such as Ω , whereas initially means at 0 t = the unknowns such as ( ) 0 0 , u v are non-negative [6].The smoothness hypothesis is used to guarantee the existence of the problem of the corresponding linear system [6].Let us look at the complete obtained system in one dimension, imposed the following boundary conditions on Equation (8), which is, imposed the following initial conditions on Equation (8), which is, combined above three equations, to get the model system.

Assumption or Hypothesis (H)
Let us suppose that v f exists and is also bounded which is subsets of domain Ω and there exists a function with . This definition implies that the function f is monotone non-decreasing in v and is uniform bounded for 0 v ≥ [6].Clearly this condition is satisfied by function f [6].Thus this property leads to ( ) ( ) , , , , where above Equations ( 11), ( 12) represents 1 2 , G G are quasi monotone non- increasing and quasi monotone non-decreasing functions in Ω respectively [6].
According to classification of the reaction functions, 1 2 , G G are typed III functions [6].This leads to the following definition of the solutions.

Definition
A smooth pair of two vector functions ( ) called upper and lower solutions respectively [6], if they satisfy the following inequalities In the above definitions the smoothness of ( ) , u v is in the sense that these functions are continuously differentiable to the order appeared in Equations (8) and (11), (12) respectively [6].Hypothesis and above definition leads to the following theorem.

Theorem
Let f holds hypothesis (H).If there exist upper and lower solutions ( ) u v converges monotonically from above and below, respectively [6], to a unique solution (u, v) of system (8) [6].Moreover, , , The quality of being of practical use about above theorem, need to construct of lower and upper solutions, with existence problem to be ensured [6], but the asymptotic behaviour of the time-dependent solution and related stability can also be established from the behaviour of the lower and upper solutions [6], according to the definition, so stability of a steady-state solution make according to selection of Lyapunov function [6], be that as it may for the global existence problem or the stability problem with homogeneous boundary conditions [6].

Numerical Methods
Let us apply numerical methods technique, to solve such system which mentioned in Equation ( 5) in finite domain Ω .We partitioned the interval [ ] , a b into n equal parts of width h [7] [8].Place a grid on the rectangle region R by drawing vertical and horizontal lines through the points.Numerical solution is denoted by ( ) , n n m m u v , whereas exact one make in the form of ( )

Forward in Time and Centre 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 system in Equation where Finite difference schemes, must pass certain tests of accuracy, consistency, stability and convergence [9] [10] [11], we discuss these concepts in the following way.

Accuracy of FTCS
Accuracy of the FTCS scheme for system in Equation ( 16), we apply Taylor's series on each term, which is as follows: simplifies above equations, we get the following, dividing above system, by k and simplifies some terms: This leads to principle part of the truncation error (PPTE), which is as follows: which shows that this scheme has 2nd order accuracy in space and first order accuracy in time such as ( )

Consistency of FTCS
According to the definition of the consistency, if the difference between finite difference equation (FDE) and related partial differential equation (PDE), i.e truncation error vanishes as the sizes of the grid spacing go to zero independently, i.e Equation ( 20) is consistent, because the truncation error divided by k tends to zero as h and k tends to zero [11] [12] [13] [14].

Convergence of FTCS
A finite difference method is convergent if the solution of the finite difference equation approaches to exact solution of the partial differential equation as the sizes of the grid spacing h and k tends to zero [13] [14] [15].
0 as , 0 where n m z is called the discretization error.
0 as , 0 where n m Z is called the discretization error.
0, as , 0 0, as , 0 this completes the proof of the convergence.

Stability of FTCS
Another important feature of a finite difference method of solving partial differential equation is the stability of the associated finite difference equation which must be investigated.Let us look again Equation ( 16), Journal of Applied Mathematics and Physics linear form of above Equation (26), is as follows: ( ) ( ) where ( )( ) , from Equations ( 26) and ( 27), we can write in the following way, where ( ) where where matrix G is called the amplification matrix.According to Von-Neumann stability analysis, the condition for stability for system is above Inequality (29) leads to the following two special cases, 1) Special Case 1 • The right hand side of the Inequity (29) gives, (30) which is condition of stability to FTCS scheme in coupled non-linear PDE system.
2) Special Case 2 • The left hand side of the Inequity (29) gives, which is condition of stability to FTCS scheme in coupled non-linear PDE system.Special case (1) and special case (2) lead to the following very interesting results, The Von-Neumann stability analysis is the most common used method of de-

Crank Nicolson Implicit Scheme
Let us apply Crank Nicolson implicit finite difference scheme to Equation ( 5).

Accuracy of CN
Accuracy of the CN scheme to Equation (36), we apply Taylor's series on each term.After some simplification, resultant is as follows, we divide above equation by time step k with using in system (5), to get the accuracy, in the following form, ( )( ) ( )( ) Now principle part of the truncation error (PPTE) is as follows: which shows that this scheme is 2nd order accurate in both time and space, such as , O k h [19] [20].

Consistency of CN
From accuracy, we find principle part of the truncation error along with Equation (42).Which shows that Crank Nicolson scheme is consistent because tends to zero as h and k tends to zero, i.e. 0 n m z → as , 0 h k → .

Stability of CN
Stability of the associated finite difference Equation (36), which is in linear form, ( ) ( ) where , R Q and A are mentioned in Equation ( 16) and Equation ( 27) respectively.From Equation (43), we can write, where ( ) According to Von-Neumann stability analysis, we have where where matrix G is called the amplification matrix.According to Von-Neumann stability analysis (29).We have Above Equation ( 46), satisfies the Von-Neumann stability criterion, which shows that CN scheme for Equation ( 5) is unconditionally stable [21] [22] [23] [24] [25].

Fourth Order Accurate Implicit Scheme
Let us apply another implicit scheme with improved accuracy in space to Equation (5).

Accuracy of Douglas Scheme
Accuracy of the Douglas scheme to Equation (36), we apply Taylor's series on each term.After some simplification, resultant is as follows, Dividing above Equations ( 54) and (55) by k and take into account Equation ( 5), so resultants are as follows: [ ] Now principle part of the truncation error (PPTE) is as follows: which shows that this scheme is 4th order accurate in space, such as , O k h [26] [27] [28] [29] [30].

Consistency of Douglas Scheme
From accuracy, we find principle part of the truncation error along with Equation (58).Which shows that Douglas scheme is consistent because n m z PPTE = tends to zero as h and k tends to zero, i.e.

Stability of Douglas Scheme
Stability of the associated finite difference Equation (51), which is in linear form, is Consider Equations (44) (60), in the following few important steps.

Richardson Extrapolation Technique
Richardson extrapolation method lead to considerable improvement of numerical results which solving the partial differential equation system by finite difference method.Richardson's extrapolation formulae are [38] [39], ( ) ( ) above formula leads to get fourth order accuracy [39] [40], 3

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 2 L and L ∞ , which are outlined below:

Results
Numerical computations have been performed using the uniform grid [14] [21] [25] [29].For the test problem (5), we analysed results with various finite difference schemes.In Tables 1-3, we fixed domain as 8 8 x − ≤ ≤ with k = 0.0001 and find results at 0.5 t = . In these three tables, we obtain insight error.6.All these six tables show results for three implicit scheme to analyse error and performance of the schemes.According to results mentioned in tables, these scheme play very important role for convergence although non-linear system is difficult to stabilized.In Figure 1, FTCS conditionally stable finite difference scheme is used at different grid, also we did in Figure 2 with Crank Nicolson implicit scheme.Fixed some parameters in CN scheme to get common region at different time level, as we mentioned in Figure 3.In

Conclusions
In this chapter, the solution to one dimensional coupled Fisher KPP system is successfully approximated by a various numerical finite difference schemes.
Explicit FTCS is conditionally stable, and we give more attention to parameter R 1 and R 2 , which can be used to stabilized the results as we can see from Figure 1.Crank Nicolson and Douglas schemes are implicit schemes with unconditionally stable nature, this guarantee is given by Von-Neumann stability analysis and ( ) , v x t are the concentrations of the chemicals in the reaction diffusion process.

Figure 4 and
Figure 4 and Figure 5, we did comparison of results with two different schemes at various time levels.In Figure 6 and Figure 7 attractive common region at different time level is obtained.In Figures 8-10, comparison of results with Richardson, Douglas and exact in ( ) , u x t component and ( ) ( ) , , u x t v x t

Figure 1 .
Figure 1.Shows results using FTCS scheme at two different grids, keeping the conditionally stable term.Area under the curves is the common region for two different grids.

Figure 2 .
Figure 2. Shows results using Crank Nicolson scheme at two different grids.

Figure 3 .
Figure 3. Shows results using Crank Nicolson scheme at two different time level.Common region is the wave front that we study in literature [35] [39] [41].

Figure 4 .
Figure 4. Shows results using Crank Nicolson scheme at various time levels.We did comparison with exact solution for ( ) , u x t and ( ) , v x t components.

Figure 5 .
Figure 5. Shows results using fourth order implicit scheme at various time levels.We did comparison with exact solution for ( ) , u x t and ( ) , v x t components.

Figure 6 .
Figure 6.Shows fourth order implicit scheme at three different time level.Common region is the wave front that we study in literature [17] [23] [29] [41].

Figure 7 .
Figure 7. Shows results using fourth order implicit scheme at various time levels, with common region.

Figure 8 .
Figure 8. Shows results using Richardson Extrapolation and fourth order implicit scheme, with comparison to exact solution for u component.

Figure 9 .
Figure 9. Shows results using Richardson Extrapolation and fourth order implicit scheme, with comparison to exact solution for ( ) ( ) , , u x t v x t form.

Figure 10 .
Figure 10.Shows results using Richardson Extrapolation at different time levels, along fourth order implicit scheme.

Figure 11 .
Figure 11.Shows a comparison of two schemes in u(x,t)/v(x,t) form, with exact solution.
stability criterion as it is generally the easiest to apply.It can only be )

Table 4 and
Table 5 show results for 2 , L L ∞ and

Table 2 .
Estimates of results using Fourth Order Implicit Scheme.

Table 3 .
Estimates of results using Sixth Order Extrapolation Method.

Table 4 .
Estimates of results, with comparison of errors as mentioned in formulae.Fixed some parameters such as

Table 5 .
Estimates of results, with two different Implicit FD schemes.Relative error values give encouragement to our results.Fixed some parameters as schemes CN and Douglas at different grids and time levels.Rate of convergence is defined in Table

Table 6 .
Estimates the rate of convergence for two implicit FD schemes for different grids.