An O ( k 2 + kh 2 + h 4 ) Accurate Two-level Implicit Cubic Spline Method for One Space Dimensional Quasi-linear Parabolic Equations

In this piece of work, using three spatial grid points, we discuss a new two-level implicit cubic spline method of O(k + kh + h) for the solution of quasi-linear parabolic equation   , , , , xx x t u f x t u u u  , 0< x <1, t > 0 subject to appropriate initial and Dirichlet boundary conditions, where h > 0, k > 0 are grid sizes in space and time-directions, respectively. The cubic spline approximation produces at each time level a spline function which may be used to obtain the solution at any point in the range of the space variable. The proposed cubic spline method is applicable to parabolic equations having singularity. The stability analysis for diffusionconvection equation shows the unconditionally stable character of the cubic spline method. The numerical tests are performed and comparative results are provided to illustrate the usefulness of the proposed method.


Introduction
The use of cubic splines for the numerical solution of linear two point boundary value problems has been discussed by Bickley [1], Fyfe [2], Albasiny and Hoskins [3] and Rubin and Khosla [4].Later, Chawla et al [5,6] have developed fourth order accurate cubic spline methods for singular two point boundary value problems.In 1983, Jain and Aziz [7] have derived fourth order cubic spline method for the solution of general two point non-linear boundary value problems.Khan and Aziz [8] have used parametric cubic spline approach for the solution of a system of second order boundary value problems.Recently, Monoj Kumar et al [9][10][11], and Rashidinia et al [12][13] have discussed higher order cubic spline finite difference method for singular two point boundary value problems.A cubic spline technique for the solution of 1D heat equation was discussed by Papamichael and Whiteman [14].This technique was extended by Fleck Jr [15] and Raggett and Wilson [16] for solving one-dimensional wave equation.Archer [17] has developed a fourth order collocation method for quasi-linear parabolic equation.Jain and Lohar [18] have solved non-linear parabolic equations by using second order cubic spline method.Recently, Rashidinia et al [19] have studied nonpolynomial cubic spline methods for the solution of parabolic equations.High order finite difference methods for the solution of non-linear parabolic equations have been discussed by Jain et al [20] and Mohanty et al [21][22].In the present paper, using three spatial grid points (see Figure 1), a new cubic spline technique similar to that of Jain and Aziz [7] is developed for the solution of one dimensional general quasi-linear parabolic equation.We use the cubic spline approximation in the space direction together with a finite difference approximation in the time direction.This approximation produces at each time level a cubic spline function, which may be used to obtain the solution at any point in the range of the space variable.In next section, we discuss the cubic spline method.In Section 3, we give the complete derivation of the method.Stability of a linear difference scheme, which is consistent with diffusion-convection equation, is discussed in Section 4 and it is shown that the linear scheme is unconditionally stable.Difficulties were experienced in the past for the high order cubic spline solution of parabolic equations with jth -level (j + 1)th -level singular coefficients.The solution usually deteriorates in the vicinity of the singularity.In this paper, we refine our technique in such a way that the solution retains its order and accuracy everywhere in the solution region.In Section 5, we compared the proposed cubic spline method with the corresponding finite difference methods discussed in [20][21][22].Finally, concluding remarks are given in Section 6.

The Two Level Implicit Cubic Spline Method
Consider the quasi-linear parabolic equation in which the function u(x,t) satisfies , , , , , , where the solution space is defined by and the boundary conditions are given by where The region  is covered in the usual manner by a  , where (N + 1)h = 1, and h > 0 and k > 0 are grid sizes in x-and t-directions, respectively.The mesh ratio parameter is given by U be the exact solution value of u(x,t) and j l u denotes a discrete approximation to u(x,t) at the grid point   , l j x t .We let   j S x denote the cubic spline interpolating the value j l u at the jth time level, and is given by which satisfies at jth-level the following properties: 1)   j S x coincides with a polynomial of degree three  , and where , , , , , We consider the following approximations:     , , , , , , , , , Further, we define   ˆ, , , , Then at each grid point   , l j x t , the cubic spline method with accuracy of O(k 2 + kh 2 + h 4 ) for the solution of differential Equation (1) may be written as where   and Note that, the initial and Dirichlet boundary conditions are given by ( 2) and (3), respectively.Incorporating the initial and boundary conditions, we can write the cubic spline method (12) in a tri-diagonal matrix form.If the differential Equation ( 1) is linear, we can solve the linear system using a tri-diagonal solver; in the non-linear case, we can use generalized Newton-Raphson method to solve the non-linear system (see Kelly [23], and Hageman and Young [24]).

Derivation of the Cubic Spline Method
For the derivation of the cubic spline method (12) for the solution of the quasi-linear parabolic Equation (1), we simply follow the approaches given by Jain and Aziz [7]. , Further at the grid point   , l j x t , we may write the differential Equation (1) as

Cubic Spline Scheme for Parabolic Equation with Singular Coefficients
Consider the linear parabolic equation of the form subject to appropriate initial and Dirichlet boundary conditions are prescribed by ( 2) and (3), respectively.Assume that the functions E x x   the above parabolic equation becomes a cylindrical and spherical problem for 1   and 2, respectively.
Applying the formula (12) to the differential Equation (26), we get where , , , At the grid point   Note that, the cubic spline sheme ( 27) is of O k kh h   accuracy for the approximate solution of Equation ( 26).However, the scheme (27) fails to compute at l = 1 when the coefficients In order to get a meaningful cubic spline scheme of O k kh h   accuracy in compact operator form, we need the following approximations: Now, by the help of the approximations ( 5)-( 10c) and (29a)-(29c), neglecting high order terms, we can re-write the scheme (27) as a two-level implicit cubic spline scheme in operator compact form where are averaging and central difference operators with respect to x-direction etc.The cubic spline scheme (30) has a local truncation error of   and is free from the terms   1 1 l  , and hence, it can be solved for l = 1(1)N in the region 0 < x < 1, t > 0. Now consider the linear singular parabolic equation O k kh h   for the solution of linear singular parabolic Equation (31).Now, consider the linear diffusion-convection parabolic equation The constant terms 0   and 0   are called the dif- fusivity and convective terms, respectively.For 0     , the problem is said to be convection dominated problem.Substituting , 0 O k h  and consistent with the differential Equation (32).Further note that, the scheme (33) is identical with the scheme discussed by Jain et al [20].It has been shown that, the scheme (33) is unconditionally stable and using this scheme computational result for the solution of Equation ( 32) is reported in [20].

Numerical Results
In order to demonstrate the application of the proposed cubic spline method, we have solved the following four problems whose analytical solutions are known to us.The right-hand side homogeneous functions, initial and boundary conditions can be obtained using the exact solution as a test procedure.We have also compared the proposed cubic spline method with the corresponding finite difference method of   O k kh h   discussed in [20][21][22].We have solved the system of linear difference equations using a tri-diagonal solver and the system of non-linear difference equations by Newton-Raphson method (see Kelly [23], and Hageman and Young [24]).All computations were performed using double length arithmetic.
Example 1: (Linear parabolic equation with singular coefficients) The exact solution is given by u(x,t) = exp(-εt) sinh x.The root mean square (RMS) errors at t = 1.0 are tabulated in Table 1 The exact solution is given by The exact solution is given by u(x,t)=exp(-t) sin( x  ).
The RMS errors at t = 1.0 are tabulated in

Final Remarks
We presented a new cubic spline discretization strategy of   O k kh h   for the solution of one-space dimensional quasi-linear parabolic partial differential equations.The proposed method with a little modification is applicable to singular parabolic problems.For singular problems, our numerical results show that the new cubic spline method may be advantageous compared with the corresponding finite difference method discussed in [20][21][22].Numerical results for non-linear problems also show better as compared to the method discussed in [20][21][22], and numerical oscillation do not appear for high Reynolds number.

For 1   and 2 ,
the equation above represents the parabolic equation in cylindrical and spherical coordinates respectively.Replacing the variable x by r and substituting