Burgers’ equation with as coefficient of viscosity can be defined as
with Dirichlet boundary conditions (BCs),
and initial conditions (ICs),
Linearized form of Burgers’ equation (by using Hopf-Cole transformation) is given as
with the Neumann boundary conditions
and the initial conditions
The study of Burger’s equation is popular among the scientific community as it is very simple form of Naviers’ Stokes equation and due to it’s appearence in various field of applied mathematics and physics such as in the context of gas dynamics, in the theory of shock waves, traffic flow, mathematical modeling of turbulent fluid and in continuous stochastic processes. In , it was first introduced by Bateman . Later in , it was introduced by Burger [6, 7] as a class of equation which delineate the mathematical model of turbulence. Recently in , Ryu etc.  propose some nowcasting rainfall models based on Burger’s equation. This equation has been solved analytically for some initial condition and solution is represented in the form of Fourier series expansion which converges slowly for small values of viscosity. Exact solution does not work very well for small values of viscosity and hence it always attracts researchers to test newly devloped numerical method on this nonlinear parabolic PDEs.
Recently, with development in computer speed several numerical schemes based on finite difference method, finite element method, spectral method, differential quadrature method, decomposition method, moving least squares particle method, Haar wavelet quasilinearization approach etc., have been developed to solve the Burger’s equation [8, 37, 2, 35, 31, 19, 47, 26, 1, 15, 3, 23, 24, 25, 33, 32, 17, 21, 22, 4, 18, 39, 40, 38, 16, 30, 45, 46].
Crank-Nicolson (CN) method [14, 41, 42, 28] is a second order method which is based on Trapezoidal formula which is A-stable but not L-stable. In the presence of inconsistencies  CN produces unwanted oscillations. Chawla etc.  produces generalised Trapezoidal formula (GTF()), where which is L-stable and gives a quite stable result. Chawla etc.  proposed a modified Simpson’s rule (ASIMP) which is A-stable and used to give fourth-order time integration formula but it produces unwanted oscillations like CN due to lack of L-stability. To remove this oscillation, Chawla etc.  produced L-stable version of Simpson’s rule and implemented it to derive a third-order time integration formula for the diffusion equation which gives stable and accurate result. Lajja  proposed L-stable derivative free error corrected Trapezoidal rule (LSDFECT). Verma etc.  developed a fifth order time integration formula for the diffusion equation which is weakly L-stable.
Here we derive order time integration formula which is weakly -stable and generalize above mentioned existing results. The issue of slow convergence of series solution for small forces analytical solution of Burgers’ equation to deviate from the exact solution. So, it is not easy to compute the solution for small values of . The newly developed method computes the solution even for small values of . To compute the numerical solution we use Mathematica 11.3 and find out that numerical solutions are in good agreement for small values of . The result are in good agreement with exact solution when inconsistencies are present in the initial and boundary condition.
The paper is organized as follows. In section 2, we give close form solution which we use to compute exact solution. In section 3, we derive higher order integration method in time for . In section 4, we use this technique combined with finite difference to solve Burgers’ equation and demonstrate the stability. In section 5, we illustrate the numerical results with tables and 2D-3D graphs.
2 Close Form Solution
with the Neumann boundary condition (BC)
and the initial condition (IC)
by non-linear transformation
The analytical solution of the linearized heat equation (2.1) is given by
where and are Fourier coefficient and is given by
The analytical solution by Hopf-Cole transformation is
3 Illustration of the proposed method
We consider the initial value problem
The Newton Cotes time integration formula is given by
Now, we use the fifth order Hermite approximation for , , , , which are given by
and sixth order Taylor’s approximation
Now, we define
Now, the time integral formula (3.2) for the interval becomes
3.1 Local trunction error
Using Taylor’s series expansion, we have
then it follows that
Also, we have
From all of the above, we deduce that
Thus the scheme (3.19) is seventh order convergent.
3.2 Stability of the formula (3.19)
Consider the test problem
and assume , then we have
From figure 1 it can be seen that and hence our scheme is not -stable. Since as and hence scheme is weakly -stable.
3.3 Stability Region for the formula (3.19)
We use the boundary locus method [29, p.64, chapter 7] and determine the boundary of the region. It can be easily seen that outside of the region it is unconditionally stable.
4 Solution of the Burgers’ equation
4.1 The final scheme
We discretize the solution space with uniform mesh expressed as For that, we partition the interval in to equal sub intervals with the spatial grid , where is a positive integer and is the spatial step.
Also partition the interval in to equal subintervals with the temporal grid where and is a positive integer.
Now define and consider linearized Burger’s equation (2.1) and compute the solution for a given and for on . Then we use (2.4)-(2.5) to deduce the following formula for computing the which is the solution of the nonlinear Burgers’ equation (1.1),
Here we approximate second order spatial derivative by fourth order finite difference ratio which is given by
and convert the linearized Burgers’ equation into an initial value problem in vector form.
Now, we apply finite difference discretization on (2.1) with the Neumann boundary conditions
we get the following equation
where and is the pentadiagonal matrix given by
Let , then applying the time integration formula on the initial value problem (4.1), we get
Now, use in Equation (4.3), we get
This method is BY using (4.3) we can compute and hence is computed at different ’s for a given time level . Physical properties of the solutions are discussed later in form of figures and tables.
4.2 Unconditional stability of the scheme for the heat equation
Equation (4.3) can be written as
The matrix is similar to a symmetric matrix.
Let us introduce a diagonal matrix
i.e., is similar to a symmetric matrix .
Now, we will show that is similar to symmetric matrix. Let
but matrices and are symmetric and commute and therefore is similar to a symmetric matrix and therefore all the eigen values of the matrix are real. ∎
Eigen values of the matrix are
Let be the eigen vectors of the matrix corresponding to the eigen value Then, we have
We set then we get fourth order difference equation
with the boundary conditions . The characteristc equation of the equation (4.11) is . Assume are the characteristc roots then we have