The Space – Time Kernel-Based Numerical Method for Burgers ’ Equations

Abstract: It is well known that major error occur in the time integration instead of the spatial approximation. In this work, anisotropic kernels are used for temporal as well as spatial approximation to construct a numerical scheme for solving nonlinear Burgers’ equations. The time-dependent PDEs are collocated in both space and time first, contrary to spatial discretization, and time stepping procedures for time integration are then applied. Physically one cannot in general expect that the spatial and temporal features of the solution behaves on the same order. Hence, one should have to incorporate anisotropic kernels. The nonlinear Burgers’ equations are converted by nonlinear transformation to linear equations. The spatial discretizations are carried out to construct differentiation matrices. Comparisons with most available numerical methods are made to solve the Burgers’ equations.


Introduction
Nonlinear partial differential equations can be approximated by various numerical methods such as finite volume methods, finite element methods, meshless methods, boundary element methods, wavelets methods, methods of fundamental solutions, and spectral methods (see [1][2][3][4][5][6][7][8][9] and references therein).In all these approaches, the time derivative is discretized using Runge-Kutta methods, Lie splitting methods, and various explicit and implicit methods.Some work is available in the literature to approximate space derivatives as well as time derivatives.For example, a space-time formulation was developed by Netuzhylov (see [10]); the author used the inverse moving least squares (IMLS) method to solve the problem of moving boundaries.Li and Mao [11] used RBFs on the space-time domain and developed a space-time global RBF method to solve inverse problems.Their resultant system was not square, so they used the least square method to overcome the ill-posedness of the problem.Li's and Mao's techniques were further used in the work [12] to solve other inverse problems.Furthermore, in [13], the authors developed a method to estimate river pollution using global space-time RBFs.Various types of space-time formulations can be found in [14][15][16][17][18]. Recently, the authors in the work [19,20] developed efficient space-time methods for solving time-dependent PDEs.In this work, a space-time numerical scheme is constructed using anisotropic kernels to solve nonlinear Burgers' equations.The time-dependent PDEs are collocated using space-time kernels in both space and time dimensions.We tested our approach to solve 1-D and 2-D Burgers' equations.

Space-Time Kernel-Based Approximation
The time-dependent PDE in the spatial domain Ω d ∈ R d , where (d ≥ 1) and for time interval (0, T f ) satisfies the equation along with some well defined initial and boundary conditions.Here, L is some spatial operator.The time-dependent problem is seen as the problem in d + 1-dimensions.The summarized form of Equation ( 1) along with the boundary and the initial condition is given by the following: These give rise to a well-posed problem (since the resulting system is square) that is used to approximate the solution to u( x) of our problem.The space-time formulation combines space and time variables so that the functional space on Ω = [0, T] × Ω d ⊂ R d+1 contains the unknown function u( x).
The following functional space is incorporated in our work.
where φ is radial kernel, and . is any norm.We have used the Euclidean norm, and points in {x c 1 , x c 2 , ...x c N } ⊂ Ω are centers.The solution u( x) at the point x ∈ Ω is then defined by The space-time kernel-based method requires the approximated unknown to satisfy equations as well as initial and boundary conditions at all nodes in domain Ω.For simplicity, we take N I internal nodes, N B boundary nodes, and N T initial time (t = 0) nodes.The unknown u( x) = u(t, x) is used to satisfy the equation by putting Equations ( 2) and (3).We then have N ∑ j=1 a j Bφ( Equations ( 7)- (9) give rise to a linear algebraic system of equations with N × N equations and N unknowns.In a more compact form, we obtain or The system of linear Equation ( 11) can be solved for unknown coefficient a using an equation solver.The solution to Equations ( 2)-( 4) can be obtained from Equation (6).

Stability of the Numerical Scheme
To discuss the stability of our numerical scheme expressed in Equation (11), which is given by where M is N × N differentiation matrix, the stability constant of Equation ( 12) is defined by The value c s is bounded for every discrete norms .on R N .Hence, we have Similarly, in the case of pseudoinverse M † of M, we obtain and we thus write Hence, Equations ( 14) and ( 16) provide the bounds for stability constant c s .

Numerical Experiments
To validate the accuracy and robustness of our numerical scheme, we solve a nonlinear Burgers' equation .The Hopf-Cole transformation [21] is used to transform a nonlinear Burgers' equation into a heat equation.The solution accuracy greatly depends on the scale factor in the kernel function (for example, see [22,23] and references therein) The space-time method is applied with anisotropic radial kernels, namely anisotropic C 2 Matern.To achieve a desirable accuracy, we used two scale factors: one for the time variable and one for the space variable, as involved in the anisotropic kernels.
The anisotropic kernel C 2 Matern and its derivative, which has enough smoothness, are defined by where ε x is the space shape parameter, ε t is the temporal shape parameter, and eps is a Matlab function.The accuracy of the space-time method is tested in terms of the maximum absolute error (MAE) defined by where u and u are the exact and the approximate solutions, respectively.

Example 1
We consider a Burgers' equation defined by with the following boundary conditions: and the initial condition Using the Hopf-Cole transformation [21,24], The nonlinear Equation ( 18) is converted into a linear heat equation with the boundary and initial conditions given by Hence, the solution u of the heat equation expressed in Equation ( 22) with conditions expressed in Equations ( 23) and ( 24) can be obtained using the space-time method discussed above.The solution to the Burgers' equation can then be recovered from the transformation expressed in Equation (21).The exact solution to the linearized problem is defined by Equations ( 22)-( 24) [25]: where the coefficients are defined by and Hence, using the Hopf-Cole transformation given by Equation ( 21), the (exact) Fourier solution to the problem given by Equations ( 18)-( 20) is obtained as where α 0 and α j (j = 1, 2, 3, ....) are defined by Equations ( 26) and (27), respectively.All results are obtained by the space-time meshless method using the anisotropic kernel with two scale factors: one in a time variable and one in a space variable.The results are displayed in Figures 1 and 2 and Table 1, respectively.The results obtained by the space-time kernel-based method are comparable to the exact solution and the results in [25].The main advantage of the current method is the avoidance of a time stepping procedure that requires a small time step for higher accuracy and stability.On the contrary, here the dimensions of our problem is increased by 1, yet it faces no difficulty, as these kernel-based methods are designed for multi-dimensional domains.The second advantage is that we can apply the method for irregular domains as well.Table 1.The space-time method results and the results in [25] where N denotes the number of collocation points, corresponding to Example 1.

Space-Time Solution Exact Solution [25]
x N =

Example 2
As a second test problem, we considered the Burgers' Equation in Equation ( 18) with the exact solution The initial and boundary conditions can be extracted from the exact solution.The Hopf-Cole transformation converts the Burgers' equation to the heat equation The solution to the heat equation for the specified initial condition is given by where b > |a| to ensure that u(t, x) > 0 for all time.The solution to the Burgers equation can now be easily obtained using the Hopf-Cole transformation The space-time meshless method is used to solve the test problem given in Equation ( 2) incorporating anisotropic kernel with two scale factors: one in the time variable and one in the space variable, as shown in Figures 3 and 4 and Table 2, respectively.It is observed that the space-time kernel-based method results are comparable to exact solutions.

Example 3
In the last test example, we consider a 2D Burgers' equation that can be transformed into the 2D heat equation using Hopf-Cole transformation [26].We consider the transformed 2D heat equation to validate the space-time kernel-based method for anisotropic kernels: with initial and boundary conditions defined by u(0, x, y) = sin(πx) sin(πy), 0 < with an exact solution given by u(t, x, y) = sin(πx) sin(πy)e −2π 2 t , 0 < x, y < 1, 0 < t < 1.
We applied the space-time method to solve a 2D heat equation in the space-time domain [0, 1] 2 × [0, 1].The results are displayed in Figures 5 and 6 and Table 3. Once again we obtained very accurate results with ease and stability.The present method can work in multi-dimensions in irregular domains as well.

Concluding Remarks
In the present work, we constructed a numerical scheme for Burgers equations.The scheme is based on space-time anisotropic kernels.The time variable is also considered as a space variable, and this increases the space dimension by 1.This space-time method avoids time-integration methods such as θ-weighted schemes and implicit and explicit numerical schemes.Such methods require a very small time step for stability and higher accuracy.The present numerical scheme reduces the cost of computation by avoiding the need to recompute the matrix for each time level, as contrary to methods such as the RK4 method for evolutionary PDEs.

Figure 1 .Figure 2 .
Figure 1.Distribution of centers in the space-time domain: Initial time centers (red) at t = 0; boundaries centers at x = 0 (blue) and at x = 1 (green).Others are internal centers.

Figure 3 .Figure 4 .
Figure 3.The first two plots show the solution to heat equations, while the second two plots show the solution to Burgers equation in the space-time domain [0, 1] × [0, 1] corresponding to Example 2.

Table 2 .
Numerical results when t = 1, a = 1, and b = 1 in the space-time domain[0, 1] × [0, 1],where N x points in the x− direction and N t in the t− direction denotes the number of collocation points in the space-time domain, corresponding to Example 2.