Numerical Solution of the Coupled Viscous Burgers’ Equation Using Differential Quadrature Method Based on Fourier Expansion Basis

The differential quadrature method based on Fourier expansion basis is applied in this work to solve coupled viscous Burgers’ equation with appropriate initial and boundary conditions. In the first step for the given problem we have discretized the interval and replaced the differential equation by the Differential quadrature method based on Fourier expansion basis to obtain a system of ordinary differential equation (ODE) then we implement the numerical scheme by computer programing and perform numerical solution. Finally the validation of the present scheme is demonstrated by numerical example and compared with some existing numerical methods in literature. The method is analyzed for stability and convergence. It is found that the proposed numerical scheme produces a good result as compared to other researcher’s result and even generates a value at the nodes or mesh points that the results have not seen yet.


Introduction 1.Background of the Study
In the field of computational mathematics, numerical methods are most widely utilized to solve equations arising in the field of physics, engineering and other sciences.Numerical analysis plays a significant role and helps us to find an ap-Applied Mathematics proximate solution for problems which are difficult to solve analytically.The design and computation of the numerical algorithm is one of the mathematical challenges faced nowadays, but scientists in the field of computational mathematics are trying to develop numerical methods by using of computers for further application.One of those numerical methods is quadrature method.
A differential quadrature method is a numerical method for evaluating derivatives of sufficiently smooth function as proposed by [1].The basic idea of differential quadrature came from Gauss Quadrature [2], which is a useful numerical integration method.Gauss Quadrature is characterized by approximating a definite integral with weighting sum of integrand value of a group at Gauss point.
Extending it to find the derivatives of various orders of sufficiently smooth function gives rise to DQ [1].In other words the derivatives of smooth function are approximated with weighted sum of the function values at a group of so called nodes.Differential quadrature can be formulated through either approximation theory or solving a system of linear equations.However, the rapid development over the recent years on problems involving nonlinear, discontinuity, multiple scales, singularity and irregularity are challenges in the field of computational science and engineering of the various numerical solutions.DQ method has distinguished themselves because of their high accuracy, straight forward implementation and generality in variety of problems [3].

Coupled Viscous Burgers' Equation
Coupled Viscous Burgers' equation is the nonlinear partial differential equation, wildly applicable in the theory of shock waves, mathematical modeling of turbulent fluid and in continuous stochastic processes.
And this equation has a wide application in the various areas of applied mathematics, such as fluid mechanics, nonlinear acoustic gas dynamics and traffic flow.The coupled viscous Burgers' equation was studied for the first time by [4] to model polydisperse sedimentation or evaluate scaled volume concentration of two kinds of particles in fluid suspensions or colloid under the effect of gravity.
The viscous Burgers' equation was presented in 1940 and in 1950 Hopf and in 1951 Cole independently introduced the method that has come to be known as the Cole-Hopf transformation to solve the Viscous Burgers' Equation [3].
The Coupled Viscous Burgers' Equation which is a nonlinear partial differential equation of the form , 0 , , 0 , 0 The boundary condition Differential quadrature method first introduced by [1] is one of the most efficient numerical method to solve partial differential equations.The key procedure in differential quadrature approximation is to determine the weighting coefficient.
[5] and [6] obtain explicit formulations to compute the weighting coefficient of the first and second order derivative.[3] presented a simple algebraic formulation to compute the weighting coefficient of the first order derivative without any restriction on the choice of the grid points and a recurrence relationship to compute the weighting coefficient of the DQ method and its applications were rapidly developed after the late 1980s, thanks to the innovative work in the computation of the weighting coefficients by other researchers and the author.
As a result, the DQ method has emerged as a powerful numerical discretization tool in the past decade.As compared to the conventional low order finite difference and finite element methods, the DQ method can obtain very accurate numerical results using a considerably smaller number of grid points and hence requiring relatively little computational effort.So far, the DQ method has been efficiently employed in a variety of problems in engineering and physical sciences.
Based on the analysis of function approximation and the analysis of linear vector space, the DQ method can be classified as polynomial-based DQ (PDQ) and Fourier series expansion-based DQ (FDQ) methods.PDQ is usually applied to non-periodic problems while FDQ can be applied to both periodic and non-periodic problems.But its performance for periodic problems is much better.The details of PDQ and FDQ methods and their application in engineering can be referred to the book of [3] which is the first book in the area to systematically describe the DQ method and its application in engineering.

Fourier Expansion Basis
The polynomial approximation is suitable for most of the engineering problems, but for some problems, especially for those with periodic behaviors, Fourier series expansion could be a better choice for the true solution instead of polynomial approximate the interval [0, 2π] [7] and [4], the Fourier series expansion can be given by ( ) ( ) where the coefficient 0 , k a c and k d are expressed as ( ) For practical application, the truncated Fourier series expansion is usually used.Thus Equation (1.4) can be described as ( ) ( ) ( ) Moreover, the convergence of the above Equation (1.8) of ( ) f x as N → ∞ is guaranteed by Weirstrass's second theorem.i.e.
The function ( ) f x be continuous on the interval [ ] 0, 2π .Then for any 0 ε > there exists an integer n and a trigonometric S n such that the inequality is satisfied for all values of x. (1.9) where ( )

Statement of the Problem
Consider the Coupled Viscous burgers' Equation in Equations (1.1)-(1.3) is a nonlinear partial differential equation, by using Fourier basis based on weighted average differential quadrature method one can expect to find its' solution numerically by using some approximation method.
In this regard [4] presented numerical simulations for the Coupled Viscous Burgers' Equation and compared the results with experimental data recently; [9] proposed a Fourier Pseudospectral method for solving Coupled Viscous Burgers' Equation; [10] and [11] applied differential quadrature method to solve the viscous Burgers' Equation; [7] use finite difference and cubic spline finite element methods to solve Burgers' equation; [10] and [11] used polynomial Differential quadrature method for numerical solution of coupled viscous Burger' equations; [12] transform the Burgers' equation to linear heat equation using Hopf-Cole transformation and then use explicit finite difference and exact explicit finite difference method to solve the transformed linear heat equation with Neumann boundary condition; [13] solve one dimensional Burgers' Equation by using differential quadrature method based on Fourier Expansion basis; [10] used cubic B-spline collocation scheme based on Crank-Nicolson formulation for time integration and cubic B-spline functions for space integration by linearizing the nonlinear terms to solve coupled viscous Burgers' Equation.
In this study, we look for the solution of coupled viscous Burger' equation by applying the differential quadrature method based on Fourier Expansion basis.
As a result, this study attempted to answer the following basic research questions.
1) How do we describe the differential quadrature method based on Fourier Expansion basis for coupled viscous Burgers' equation?
2) How the present method is applied to solve coupled viscous Burgers' equation?
3) To what extent the proposed method is approximate the exact solution?

Objectives of the Study 1.3.1. General Objective
The general objective of this study is to find the numerical solution of coupled viscous burger equation using differential quadrature method based on Fourier expansion basis.

Specific Objectives
The specific objectives of the present study are: 1) To describe the differential quadrature method based on fourier expansion basis in solving coupled viscous Burgers' equation numerically.
2) To solve the coupled viscous Burgers' equation using differential quadrature method based on a Fourier expansion.
3) To determine the accuracy that the proposed method approximate the exist solutions.

Significance of the Study
The outcomes of this study have the following importance to: 1) Find an alternative numerical solution of coupled viscous Burgers' Equation.
2) Apply differential quadrature method based on Fourier expansion basis to find the solution for some practical problems like coupled viscous Burgers' equation.

Delimitation of the Study
This study was delimited to the numerical solution of the coupled viscous Burgers' equation.
The study also delimited to the practically most important (linear) differential quadrature based on Fourier expansion that are of the same type as differential quadrature and Fourier expansion formulas, i.e., linear combinations of weighted differential function evaluations and their error analysis for quadrature.

Study Area and Period
The study was conducted in Jimma University under the department of mathematics from September 2016 G.C. to September 2017 G.C. Conceptually the study focus on the solution of coupled viscous Burgers' equation using differential quadrature method based on Fourier expansion basis.

Study Design
This study was employed mixed-design (documentary review design and experimental design) on the coupled viscous Burgers' equation type.

Source of Information
The relevant sources of information for this study is books, published articles DOI: 10.4236/am.2018.97057826 Applied Mathematics and related studies from internet and the experimental results were obtained by writing MATLAB code for the present methods.

Mathematical Procedures
The study is an experimental as it evolves entirely laboratory work with the help of computer and MATLAB software.Farther, important materials for the study were collected by the researcher using the documentary analysis.The required numerical data was collected by coding and running using MATLAB software to get the numerical results and table of some examples that have exact solution, to show the validity and efficiency of the method.Hence, in order to achieve the stated objectives, the study was followed the procedures 1) Problem preparation or formulation.
3) Replacing the partial differential equation (or the Coupled Viscous Burgers' Equation) by using differential quadrature method based on Fourier expansion basis to obtain a system of first order ordinary differential equations.
4) The obtained systems of ordinary differential equations was solved by classical fourth order Runge-Kutta method.
5) Validating the schemes using numerical examples.
6) Writing MATLAB code for the method to solve the systems obtained.

Differential Quadrature Method
In seeking an efficient discretization technique to obtain accurate numerical solution using a considerably small number of grid points, [1] and [14] introduced the method of differential quadrature (DQ).where a partial derivative of a function with respect to a coordinate direction is expressed as a linear weighted sum of all the functional values at all mesh points along that direction.The key to DQ is to determine the weighting coefficient for the discretization of a derivative of any order.[1] suggested two methods to determine the weighting coefficients of first order derivative.The first method is based on an ill-conditioned algebraic equation system and the second method uses a simple algebraic formulation, but the coordinates of the grid points are fixed by the roots of the shifted Legendre polynomial.In earlier applications of the DQ method, Bellman's first method was usually used because it allows the use of arbitrary grid points distribution.
However, since the algebraic equation system of this method is ill-conditioned, the number of the grid points usually used is less than 13 [15].After that, [5] obtain explicit formulations to compute the weighting coefficient for the first and second order derivative.[16] presented a simple algebraic formulation to compute the weighting coefficient of the first order derivative without any restriction on the choice of the grid points and a recurrence relationship to compute the weighting coefficient of the DQ method and [3] and [17]  developed some simple algebraic formulation to compute the weighting coefficient.Recently, the most frequently used DQ procedures to solve one and two dimensional differential equations are Lagrange interpolation polynomials based differential quadrature method (PDQM).The DQ method approximates the derivative of a smooth function at a grid point by a linear weighted summation of all the functional value in the whole computational domain.
The coupled viscous Burgers' equation is a nonlinear partial differential equation of the form , 0 , , 0 , 0 The boundary condition where η, α, ε and β are arbitrary constants depending on the system parameters.

Differential Quadrature Method Based on the Fourier Expansion basis
Let the first and second order derivative of ( ) u x at a point i x be approx- imated by the following equations , for 1, 2, , , for 1, 2, , where ij a and ij b represents the weighting coefficients of first and second de- rivative respectively, N is the number of grid points.And the key procedure in DQ is to determine the weighting coefficient.
For, the solution of any partial differential equation has to be convergent and bounded, so the solution ( ) u x of the viscous Burger equation be approximated by a Fourier expansion of the form ( ) ( ) ( ) u x in this equation constitutes an ( ) 1, sin π , cos π ,sin 2π , cos 2π , ,sin π and cos π are linearly independent vectors.Thus, we can consider (4.4a) as a set of base vectors.
Actually FDQ uses two sets of base functions, one is the base functions (4.4a) and the other is the terms of the Lagrange interpolating trigonometric poly-nomial given by where ( ) ( ) By using the above two sets of base vectors to drive explicit formulations to compute the weighting coefficients of the first and second order derivatives, for the non-diagonal weighting coefficient we need to use the second set of base vectors.For simplicity, we set where is the Kronecker delta operator.
Then Equation (4.4b) can be reduced to Using the same approach as in a polynomial differential quadrature method [18] we let all the base vectors given by Equation (4.10) satisfy Equations ((4.1) and (4.2)) and obtain ( ) ( ) It is observed from Equations ((4.11) and (4.12)), the computation of ij a and ij b is equivalent to the evaluations of ( ) ( ) can be calculated by Equation (4.9).We successively differentiate Equation (4.8) From the above equation, we can obtain the following result Substituting Equations ((4.16) and (4.17)) in to Equation ( Similarly substituting Equations ((4.19) and (4.20)) in to (4.12) and using Equations ((4.17), (4.18)), we obtain ( ) From Equations ((4.21) and (4.23)), ij a and ij b ( i j ≠ ) can be obtained.
However, the calculation of ii a from Equation of the first and second order derivatives can be calculated with FDQ.It should be indicated that these equations can be applied to the periodic problems and the non-periodic problems.For the non-periodic problems, the range in the computational domain is 0 π x ≤ ≤ , while for the periodic problems, the range in the computational domain is 0 2π x ≤ ≤ .
Where N is the number of grid points, ij a and ij b are the weighting coeffi- cients in the x direction, and when FDQ method is used ij a and ij b are com- puted by Equations ((4.21), (4.23)).The resulting system of ODEs (4.27) can now be solved by using the classical fourth order Runge-Kutta method.

Numerical Experiments
In this section, we consider two examples to test the efficiency of the proposed method.For describing the error, we consider maximum error norm for u defined as follows where ( ) To illustrate the efficiency of the proposed numerical scheme, we solve two test examples and throughout the numerical experiment, we consider step length in time space t ∆ .
and 0 , , a α β are arbitrary constants.In Table 3 and Table 4, we present comparison between the numerical solutions obtained by the present method and the method proposed in literature.Table 3 and Table 4 shows the maximum absolute errors for 1, 2 η = = and various values of α and β .
From Table 1 and Table 2 it can be seen that the solution of Maximum absolute norm u and v decrease to zero at a time t as the value of , η ε increase it can also observed that the Differential quadrature method based on Fourier Expansion basis method is capable of finding numerical solution for larger value of , η ε .
Here the parameters α and β determine the size of the boundary layer near the edges of the domain.We use the value of 1, 2 η ε = = in our numerical experiment.
Table 3.Comparison of numerical results of Example 2 with the results obtained by [9] and [11] for the variable u with 0 0.05 a = , N = 21 with 0.01  [9] and [11] for the variable v with 0 0.05 a = , N = 21 with 0.01 In this case we solve the PDE on the random point distribution.The optimal shape parameter α and β is also taken the same as we achieved in the linear PDE case.The error norm in Table 2 shows as differential quadrature method based on Fourier basis.Also works well for this nonlinear problem.However, it is shown that the method can be fluctuated as applied to a wide class of higher-dimension, nonlinear partial differential equations with a little modification depending on the system parameters.In Table 3 and Table 4 we compare the differential quadrature method results with those available ones in the literature.As seen from Table 3 and Table 4 the differential quadrature method based on Fourier expansion is quite comparable with those methods in the literature.

Conclusions
This study demonstrates the application of FDQ method to solve one-dimensional coupled viscous Burgers' equation.This method produces a system of first order ordinary differential equations, which can be solved numerically by the classical fourth order Runge-Kutta method.To verify the method two test problems have been considered and the results obtained are compared with some other methods which are available in the literature.The numerical solutions obtained show that the present method produces reasonably accurate numerical solution of coupled viscous Burgers' equation.
Therefore differential quadrature method based on Fourier expansion basis is a reliable method to obtain the numerical solution of some physically important nonlinear problems governed by partial differential equation; this is because it brings the problem into a smooth and periodic one.It is necessary to postulate a shock across which u changes discontinuously to have a unique solution and so a physically result.
It is shown that the present numerical scheme gives better solutions on most of the grid points.However, it is shown that the method can be fluctuated as applied to a wide class of higher-dimension, nonlinear partial differential equations with a little modification depending on the system parameters.

Future Scope
In the future by applying different differential quadrature method, I will try to solve different partial differential equations numerically.
If the viscous term is dropped from the Burgers' equation the nonlinearity allows discontinuous solutions to develop.A wave is convicting from left to right and solutions for successive times are indicated.Points on the wave with larger values of u convict faster and consequently overtake parts of the wave convicting with smaller values of u .Coupled viscous Burgers' equation is the nonlinear partial differential equation, that encounters in the theory of shock waves, mathematical modeling of turbulent fluid and in continuous stochastic processes.

Table 2 .
Maximum absolute norm of u and v of Example 1 for different values of t ∆ the proposed method, at t = 1.

Table 4 .
Comparison of numerical results of Example 2 with the results obtained by