Taylor’s Meshless Petrov-Galerkin Method for the Numerical Solution of Burger’s Equation by Radial Basis Functions

During the last two decades, there has been a considerable interest in developing efficient radial basis functions RBFs algorithms for solving partial differential equations PDEs . In this paper, we introduce the Petrov-Galerkin method for the numerical solution of the one-dimensional nonlinear Burger equation. In this method, the trial space is generated by the multiquadric MQ RBF and the test space is generated by the compactly supported RBF. In the time discretization of the equation, the Taylor series expansion is used. This method is applied on some test experiments, and the numerical results have been comparedwith the exact solutions. The L2, L∞, and root-meansquare RMS errors in the solutions show the efficiency and the accuracy of the method.


Introduction
Burger's equation serves as a useful model for many interesting problems in applied mathematics and related sciences.It is in the models of effectively certain problems of a fluid flow nature, wherein either shocks or viscous dissipation is a significant factor.
The first steady-state solution of Burger's equation were given by Bateman 1 in 1915.However, the equation gets its name from the extensive research of Burgers 2 beginning in 1939.This equation has a large variety of applications in modeling of water in unsaturated soil, dynamics of soil water, statistics of flow problems, mixing and turbulent diffusion, cosmology, and seismology 3-5 .
In one-dimension case, Burger's equation is given by where μ > 0 is a viscosity constant and the subscripts x and t denote the space and time differentiations, respectively.The initial condition and the boundary conditions are given by u x, t 0 u 0 x , 1.2 u a, t g 1 t , u b, t g 2 t .1.3 Various numerical methods have been introduced to solve Burger's equation, such as cubic splines and finite differences 6 , compact differencing method 7, 8 , finite element method FEM 9 , the Tau method 10 , and the method of lines 11 .
Many of these methods need the mesh generation and provide the solution of the problem on mesh points.But in recent years, other methods, namely, meshless methods, have been provided, which do not need to generate the mesh.In a meshless method, a set of scattered nodes is used instead of meshing the domain of the problem.A powerful tool for the scattered data interpolation problem is the RBFs.In the last two decades, the development of the RBFs as a truly meshless method for approximating the solutions of PDEs has drown the attention of many researchers in science and engineering.The initial development was due to the pioneering work of Kansa 12,13 who directly collocated the RBFs for the approximated solution of the equations.
To date, useful meshless methods have been provided for the nonlinear PDEs, wherein the RBFs are used, such as the collocation method for solving Burger's equation 14 , Korteweg-de Vries equation 15, 16 , and sine-Gordon equation 17 .
In this paper, we present the other meshless method, which is called the meshless local Petrov-Galerkin MLPG method.This method was first introduced by Atluri and Zhu 18 .The MLPG method is very general, can be based on the symmetric or unsymmetric local weak forms of the PDEs, and uses a variety of interpolation method for test and trial functions.Hence, we use the local unsymmetric weak form of the problem, and the MQ RBF and the compactly supported RBF are chosen as the trial and the test functions, respectively.
In classical methods, using low order time integration for transient problems may cause loss in the accuracy of the numerical schemes.On the other hand, increasing the order of time discretization can yield better results.The main idea behind the present time integration is to use more time derivatives in the Taylor series expansion.Recently, Da ĝ et al. 19 have developed the Taylor-Galerkin and the Taylor-collocation methods for the numerical solution of Burger's equation by B-splines.The objective of the present paper is to obtain the numerical solution of Burger's equation by combining the MLPG method with the aforesaid time discretization of the Taylor series expansion that according to the best of our knowledge, this is the first demonstration of the application of it.
The paper is organized as follows.A brief knowledge of the RBFs interpolation is given in Section 2. In Section 3, the time discretization of Burger's equation is introduced, and in Section 4, the MLPG method is presented.In Section 5, we give several numerical experiments to demonstrate the accuracy and efficiency of our meshless numerical scheme, and Section 6 contains some conclusion.

Radial Basis Function Interpolation
In order to explain multivariate scattered data interpolation by RBFs, suppose a data vector f| X f x 1 , . . ., f x N T ∈ R d of function values, sampled from an unknown function f : , is given.Scattered data interpolation requires computing a suitable interpolant S : R d → R satisfying S| X f| X , that is, S x j f x j , j 1, . . ., N.

2.1
To this end, the RBF interpolation scheme works with a fixed radial function ϕ : 0, ∞ → R, and the interpolant S in 2.1 is assumed to have the form Since enforcing the interpolation conditions in 2.1 leads to a system of N linear equations in the N M unknowns λ j and c i , one usually adds the M additional conditions N j 1 λ j p i x j 0, i 1, . . ., M.

2.3
In general, solving the interpolation problem based on the extended expansion 2.2 , now amounts of solving a system to linear equations of the form: where the pieces are given by A ij ϕ x i − x j , i, j 1, . . ., N, P ik p k x i , i 1, . . ., N, k 1, . . ., M, Λ λ 1 , . . ., λ N T , C c 1 , . . ., c M T , and 0 is a zero vector of length M.
From 20, 21 , we know that we have a unique interpolant S x of f if ϕ r is a conditional positive definite RBF of order m.The commonly used RBFs are MQ, Gaussian, and compactly supported RBFs, see Table 1.The MQ has been found to provide the most accurate approximation in most of the applications of the RBFs 22 .For the MQ RBF, ϕ r √ r 2 c 2 , we have m 1.In this case and for the one-dimensional space, the linear basis function is given by p x 1.Therefore, in one-dimensional and for the MQ approximation scheme, 2.2 yields the following expansion: and 2.3 yields From this equation, we obtain Substituting 2.7 into 2.5 yields or in the compact form, we obtain where the basis functions g j x are given by

2.10
For more details about RBFs, see 23 .
The variable ε in the MQ RBF, ϕ r 1 εr 2 , wherein c is replaced by 1/ε, is known as the shape parameter controls the shape of the functions.The shape parameter ε influences the solution profoundly.Since the exact value of ε cannot be calculated, empirical studies have been made by researchers in the past to obtain a suitable value for the shape parameter ε.There are many methods for choosing the optimal value of ε, such as the brute force method 24 and the leave-one-out cross-validation LOOCV algorithm 25 .Some other methods for finding the optimal shape parameter were produced by Kansa and Hon 26 , Hardy 27 , and Franke 28 .In this paper, we use the value of the shape parameter that is suggested by Hardy 27 .Hardy suggested the value where d i is the distance from the ith center to the nearest neighbor and N is the number of centers.

Time Discretization
For the time discretization of 1.1 , we use the Taylor series expansion.In this approach, the term u n t u t x, t n , t n 1 t n Δt, is arranged with the help of the Taylor series expansion as Differentiating 1.1 with respect to time, u n tt may be written by For the time derivative u n t in 3.2 , using forward difference formula, u n tt can be rewritten as Substituting 3.3 into 3.1 and using the resulting expression in 1.1 yield the following time-discretized form of Burger's equation: 3.4

The Meshless Local Petrov-Galerkin Method
The MLPG approach was first proposed by Atluri and Zhu 18 for solving linear potential problems, by using domain discretization technique.The MLPG approach uses either a local symmetric weak form, or an unsymmetric weak form of the governing equation over the local subdomain.
In the application of a local Petrov-Galerkin scheme to the Burger's equation, 3.4 is multiplied with the weight function w and the resulting equation is integrated over the local subdomain such as Ω s , which is a small region taken for each node in the global domain Ω.Therefore, the local weak form, is obtained which requires u to be twice differentiable and u x to be continuous.If not, u xx will be infinity where u x is discontinuous and 4.1 makes no sense.On the other hand, there is no continuity requirement on the weight function w in 4.1 .Thus, the requirements on u and w are unsymmetric, in order to be admissible in 4.1 .Hence, we denote 4.1 as a local unsymmetric weak form LUWF .The weight function w is often referred to as the test function.In 4.1 , Ω i s is a local subdomain associated with the point i.

ISRN Applied Mathematics
In this paper, the compactly supported RBFs are chosen as the test functions.In this case, Ω i s is support of the compactly supported RBFs.The compactly supported RBFs are generally expressed in the form ψ r 1 − r n p r Wu 29 and Wendland 30 , and where p r is a prescribed polynomial.By replacing r with r/δ for δ > 0, the basis function has support on 0, δ .In the numerical results for compactly supported RBF, the Wendland function is used as follow: ψ r 1 − r 8 32r 3 25r 2 8r 1 .

4.3
In LUWF 4.1 , the collocation approach is used to impose both the essential as well as natural boundary conditions.The 10-point Gauss quadrature rule 31 is used for the numerical integration of 4.1 .

Numerical Experiments
Three test experiments are studied to investigate the robustness and the accuracy of the proposed method.The L 2 , L ∞ , and RMS errors, which are defined by the parameters μ 0.005 and Δt 0.01 over the solution domain 0, 1.2 .The computational results are listed in Tables 2, 3, and 4. We also plot the profiles of the solutions at t 1.3, 1.6, 2.1, and 2.4 in Figure 1. Figure 2 represents the percentage absolute error at different time levels.Figure 3 represents the percentage absolute error for various values of N at the time t 1.7 of Experiment 1.
According to Table 3, we observe that first the value of errors decrease, but for N 40 these values increase.When we add center points in order to improve the accuracy, the condition number of the interpolation matrix grows and then the problem becomes illconditioned.The ill-conditioning grows due to the decrease in the distance between points, and not to the increase in the number of center points 32 .the exact solution.We study this experiment with parameters α 0.4, β 0.6, λ 0.125, and μ 1.The computational results are listed in Tables 5, 6, and 7. We also plot profiles of the percentage absolute error at different time levels in Figure 4.The initial function u 0 x in 1.2 is obtained from 5.4 at time t 0. The boundary functions g 1 t and g 2 t in 1.3 can be obtained from the exact solution.The computational results are listed in Tables 8, 9, and 10. Figure 5 represents the L 2 and L ∞ norm errors for t 0.7 with various values of μ in Experiment 3.

Conclusion
In this paper, we present a new meshless method, the MLPG method, for solving Burger's equation.The Wendland compactly supported RBF is used as the test function in   interpolation.In the weak form, integrals are evaluated over the local subdomain instead of the global domain.This method can be extended to solve Burger's equation in the higherorder dimensions because the MQ RBF can easily be extended to higher dimensions.Also, because of the infinite differentiability of the RBFs, we can use the MLPG method for the PDEs with higher-order derivatives with respect to x, for example, the Korteweg-de Vries equation.

Figure 3 :
Figure 3: The percentage absolute error for μ 0.005, Δt 0.01, and various values of N at t 1.7 of Experiment 1.

Experiment 3 .
In this experiment, we study 1.1 with the exact solution 34 , u x, t 0.1e −A 0.5e −B e −C e −A e −B e −C , 0 x 1, t 0, 5.4

Figure 4 :
Figure 4: The percentage absolute error for μ 1, Δt 0.01, and N 11 at different times of Experiment 2.

Figure 5 :
Figure 5: L 2 and L ∞ norm errors for t 0.7 with various values of μ of Experiment 3.

Table 1 :
Different radial basis functions.

Table 3 :
Errors for μ 0.005, Δt 0.01, and various values of N at t 1.7 of Experiment 1.

Table 6 :
Comparison of results for N 11 at t 1, of Experiment 2.

Table 7 :
Errors at different times for μ 1, Δt 0.01, and various values of N of Experiment 2.
the local weak form, and in time discretization, the Taylor series expansion is used.Results of numerical experiments indicate that it is an accurate and efficient numerical scheme.The proposed method is a truly meshless method, which does not require domain elements in the

Table 10 :
Errors at different times for various values of parameter μ with N 20 of Experiment 3.