Variational multiscale element free Galerkin (VMEFG) method is applied to Burgers’ equation. It can be found that, for the very small diffusivity coefficients, VMEFG method still suffers from instability in the presence of boundary or interior layers. In order to overcome this problem, the high order low-pass filter is used to smooth the solution. Three test examples with very small diffusion are presented and the solutions obtained are compared with exact solutions and some other numerical methods. The numerical results are found in which the VMEFG coupled with low-pass filter works very well for Burgers’ equation with very small diffusivity coefficients.
1. Introduction
Burgers’ equation, which is a nonlinear partial differential equation, can model several physical phenomena such as traffic, shock waves, gas dynamics, longitudinal elastic waves in an isotropic solid, and turbulence problems. In general, the exact solution of Burgers’ equation can be only obtained in special case. Moreover, in many cases these solutions involve infinite series which may converge very slowly for small diffusivity coefficients. Therefore, finding accurate and efficient numerical methods for solving Burgers’ equation has been an attractive research undertaking. The main challenge of numerical solution of Burgers’ equation lies in nonphysical oscillations for convection-dominated problem. Thus, many numerical schemes are constructed or developed in the last few years to solve Burgers’ equation, such as finite difference method (FDM) [1–3], finite element method (FEM) [4–8], B-spline methods [9–11], and lattice Boltzmann method (LBM) [12]. A comprehensive review of the different methods for solving Burgers’ equation can refer to [10].
In recent years, meshless methods emerge as a new numerical method for solving various mathematical-physical problems. Compared with the traditional computational methods based on mesh, meshless methods have many outstanding advantages. In these methods, the approximation is built only based on nodes and no predefined nodal connectivity is required; meanwhile, the removal and the addition of nodes in the domain are easily performed.
Previous studies [13, 14] show that the numerical solution may be corrupted by nonphysical oscillations when meshless methods are directly used to solving Burgers’ equation with small diffusion. Usually, spatial stabilization technique is required for the meshless methods and many meshless stabilized methods are developed in recent years. Most of these stabilized methods are a generalization of the stabilized FEM [13–17], such as streamline-upwind Petrov-Galerkin (SUPG), Galerkin least-square (GLS), and subgrid scales (SGS). As Fries and Matthies [16] had pointed out, standard stabilization approaches in FEM can be directly applied to the meshless methods; however, the choice of the stabilization parameter required special attention. But in practice, the choice of the stabilization parameter is more or less dependent on the researches’ experience. In order to avoid the selection of the stabilization parameter, Zhang and his coworkers [18] presented variational multiscale element free Galerkin method (VMEFG) which is free of user-defined stability parameters and a consistent stabilized method. Subsequently, they applied the method to solve Stokes problems [18], water wave problems [19], MHD flows problems [20], and convection-diffusion-reaction problem [21] successfully. But, as far as the solution of Burgers’ equation is concerned, for the very small diffusivity coefficients, the VMEFG method also produces the numerical pseudo oscillations as we will see in the paper; therefore, we must find another technique to deal with it.
In recent years, high order low-pass filters have been proposed to treat hyperbolic conservation laws [22, 23] and some experiences show that a low-pass filter is ideal for suppressing the amplitude of undesirable high frequency components and does not affect the remaining components of the solution [24]. To the best of our knowledge, most numerical applications are restricted to problem involving moderate and larger diffusion. Moreover, no published results with meshless methods in conjunction with low-pass filter are used to solve Burgers’ equations with very small diffusivity coefficients. Therefore, the main goal of the paper is to evaluate VMEFG coupled with high order low-pass filter for Burgers’ equation with very small diffusion.
This paper is organized as follows. In Section 2, we list some elementary knowledge about the VMEFG method and low-pass filter. In Section 3, numerical experiments are performed using our proposed method. And the conclusions are presented in the Section 4.
2. Numerical Algorithms
Let Ω⊂Rnd be an open bounded region with piecewise smooth boundary Γ. The number of space dimensions, nd, is equal to 1 or 2. Burgers’ equation is written as follows:(1)∂u∂t+u·∇u-1Re∇·∇u=0inΩ×0,T,where u is the dependent variable resembling the flow velocity and Re is the Reynolds number characterizing the size of viscosity. From (1), it can be obviously seen that Burgers’ equation is a simplified form of the Navier-Stokes equation having viscosity and nonlinear convection term. When the diffusivity coefficient is very small, Burgers’ equation models convection-dominated problem.
Let w represent the weight function for velocity u; then the standard weak form of (1) is given as follows:(2)w,u˙+w,u·∇u+∇w,1Re∇u=0,where u˙=∂u/∂t, (·,·)=∫Ω(·)dΩ, that is, the L2 product of the indicated arguments over domain Ω.
In the paper, the moving least-square (MLS) method is used to construct the shape functions of meshless method and the detailed MLS approximation can refer to [25]. In general, owing to that the MLS approximation does not possess the Kronecker Delta condition property, this makes the imposition of essential boundary conditions more complicated and time-consuming than that of FEM. Here we used a simple technique proposed by Zhang et al. [19] to enforce the essential boundary conditions; that is, the nodal influence domain of meshless method is extended to have arbitrary polygon shape. Their numerical results indicated that the shape functions almost possess the Kronecker delta function property when the dimensionless size of the nodal influence domain α approaches to 1.
2.1. Variational Multiscale Element Free Galerkin Method
The VMEFG method is inheritance of the variational multiscale FEM which was offered by Hughes et al. [26] and the main idea of the method is to divide the unknown physical quantity into fine scale and coarse scale. Then it can use bubble functions to determine fine scale solution analytically and gets stabilization parameter naturally. Last, it substitutes the fine scale solution into coarse scale problem and then obtains the coarse scale solution numerically. In the following, the brief introduction of VMEFG is presented.
Step 1 (multiscale decomposition).
Assume that the velocity u and its weight function w can be decomposed into coarse scale and fine scale, respectively; namely,(3)u=u-+u^,(4)w=w-+w^,where u- and u^ are the coarse scale and fine scale for velocity u and w- and w^ are the coarse scale and fine scale for weight function w, respectively.
Similar to literatures [19, 27], we further assume that u^ and w^, although nonzero within background integral cell ΩK, vanish identically over the ΩK boundary when the factor α approaches to 1:(5)u^=0on∂ΩK,(6)w^=0on∂ΩK.We now substitute the trial solutions equation (3) and the weighting functions equation (4) in the standard variational equation (2), and this becomes the point of departure from the conventional Galerkin formulations(7)w-+w^,u-˙+w-+w^,u-+u^·∇u-+u^+∇w-+w^,1Re∇u-+u^=0.Next, employing the linearity between coarse scale w- and fine scale w^, (7) can be split into coarse scale problem W- and fine scale problem W^ as follows.
The Coarse Scale Problem W-. Consider(8)w-,u-˙+w-,u-+u^·∇u-+u^+∇w-,1Re∇u-+u^=0.
The Fine Scale Problem W^. Consider(9)w^,u-˙+w^,u-+u^·∇u-+u^+∇w^,1Re∇u-+u^=0.Obviously, the coarse and fine scale equations are in fact nonlinear equations because of the convection term. In general, to solve nonlinear equations we need to linearize them. Here, we substitute the nonlinear convective coefficient u-+u^ in (8) and (9) by the last converged solution uc from the fixed point iteration.
The Linearized Coarse Scale Problem W-. Consider(10)w-,u-˙+w-,uc·∇u-+u^+∇w-,1Re∇u-+u^=0.
The Linearized Fine Scale Problem W^. Consider(11)w^,u-˙+w^,uc·∇u-+u^+∇w^,1Re∇u-+u^=0.
Step 2 (solution of the fine scale problem).
According to the literatures [19, 27], we next deduce the fine scale solution u^ by using the bubble functions. Let us now consider the fine scale part of the weak form W^. Exploiting linearity of the solution slot in the first term and employing integration by parts to the second term on the left-hand side in (11), we get(12)w^,uc·∇u^+∇w^,1Re∇u^=w^,-u-˙-uc·∇u-+∇·1Re∇u-=w^,r-,where r-=-u-˙-uc·∇u-+∇·1/Re∇u-. In this process, we have applied the integration by parts as well as the relevant assumptions such as (6). In addition, it is important to note that r- is the residual of Burgers’ equation for the coarse scale. This is a crucial ingredient of the VMEFG method and ensures that the resultant formulation yields a consistent stabilized method.
In order to obtain fine scale solution u^ analytically from (12) and deduce the structure of the stability parameter τ, we use bubble functions similar to FEM context [27]. Without loss of generality, we assume that the fine scale u^ and w^ are represented via bubble functions over the ΩK of the EFG shape function when factor α approaches to 1:(13)u^K=b1Kξβ⟶u^iK=b1Kξβi,onΩK,(14)w^K=b2Kξγ⟶w^iK=b2Kξγi,onΩK,where b1K and b2K represent the bubble functions for the trial solution u^ and weight function w^ over the ΩK, respectively, and β and γ are their coefficients, respectively, i=1,…,nd (nd indicates the dimension of the problem).
In the paper, for the 1D problem, the bubble functions fine scale trial solution and weighting function are given as follows [27]:(15)b1K=1-x2,-1≤x≤1,b2Kx=11+Xbx+1,-1≤x≤Xb,-11-Xbx-Xb+1,-1≤x≤Xb,where Xb is the location of the internal virtual node for the piecewise linear bubble and in the paper Xb=-0.5.
For the 2D problem, when the triangle background cells are used and in a reference background cell, the bubble for the fine scale trial solution is a quadratic bubble defined as [21](16)b1Kx,y=271-x-y.For the fine scale weighting function, the reference background cell is divided into three regions (see Figure 1) and the bubble functions on these regions are given as follows:(17)b2Kx,y=xyXbYbforx,yinregion1,1-x-y1-Xb-Ybforx,yinregion2,xyXbYbforx,yinregion3,where Xb and Yb represent the location of the internal virtual node in the background cell and in the paper (Xb,Yb)=(0.34,0.34).
Substituting (13) and (14) into (12), and taking the constant coefficient vector β and γ out of the integral expression, we have(18)γT∫Kb2Kuc·∇b1KdΩI+∫K1Re∇b1K:∇b2KdΩI+∫K1Re∇b1K⊗∇b2KdΩβ=γTb2K,r-.Since γ is arbitrary, consequently we get(19)β=S-1R,where S and R are defined as follows:(20)S=∫Kb2Kuc·∇b1KdΩI+∫K1Re∇b1K:∇b2KdΩI+∫K1Re∇b1K⊗∇b2KdΩ,R=∫Kb2Kr-dΩ,where I is a nd×nd identity matrix and both ∇b1K and ∇b2K are nd×1 vectors of gradient of the bubble functions. Substituting (19) into (13), then we can obtain the fine scale u^ over the ΩK as follows:(21)u^K=b1Kξβ=b1KξS-1R.
Schematic diagram showing the domains for the bubble functions for a reference triangle cell [21].
Step 3 (solution of the coarse scale problem W-).
Once we have obtained the fine scale solution u^, now we turn our attention to the coarse scale problem W-. Exploiting the linearity of the solution slot between u- and u^, (10) can be written as(22)w-,u-˙+w-,uc·∇u-+w-,uc·∇u^+∇w-,1Re∇u-+∇w-,1Re∇u^=0.Using the integration by parts for the terms containing u^ in (22) as well as the assumptions used in (5), (22) has the following forms:(23)w-,u-˙+w-,uc·∇u--uc·∇w-,u^+∇w-,1Re∇u--∇·1Re∇w-,u^=0.Substituting (18) into (23), subsequently assuming that the coarse scale residual r- is a constant over the ΩK, after rearrangement we get the finally computed formulation:(24)w-,u-˙+w-,uc·∇u-+uc·∇w-+∇·1Re∇w-,τu-˙+uc·∇u--∇·1Re∇u-=0,where(25)τ=b1K∫Kb2KdΩ∫Kb2Kuc·∇b1KdΩI+∫K1Re∇b1K:∇b2KdΩI+∫K1Re∇b1K⊗∇b2KdΩ-1.Compared with the standard variational form (2), it is important to note that the third term on the left appears in (24) owing to the assumption that there exists fine scale in the problem, which in fact models the effect of the fine scale.
2.2. Low-Pass Filter
Although the VMEFG method can play a role in stability for some moderate Re values [21], but for very large Re values, it will still produce spurious oscillations as we will see in the next section. Previous researches [22–24, 28] indicate that a low-pass filter can suppress the spurious oscillations; thus in the paper we use the high order low-pass filter to smooth the solutions for very large Re values at each time loop.
For one-dimensional problems, a high order tridiagonal filter can be formulated as [23, 24](26)αfϕ^i-1+ϕ^i+αfϕ^i+1=∑m=1Man2ϕi+m+ϕi-m,where ϕi denotes the given value at point i and ϕ^i denotes the value after filtering. With proper choices of coefficients, (22) provides a 2Mth-order formulation on a 2M+1 point stencil. The M+1 coefficients a0,a1,…,aM can be obtained in terms of αf by Taylor expansion. In the paper, we use the following coefficients of a family of 10th-order filter derived in [24, 28](27)a0=193+126αf256,a1=105+302αf256,a2=15-1+αf64,a3=451-2αf512,a4=5-1+2αf256,a5=1-2αf512.
For two-dimensional problems, the explicit high order filter can refer to [29].
3. Numerical Examples
To show the effectiveness of the combination VMEFG and low-pass filter, numerical solutions of Burgers’ equation with small diffusivity coefficients for three problems will be considered.
Example 1.
Consider the following 1D Burgers’ problem:(28)∂u∂t+u∂u∂x=1Re∂2u∂x2,0<x<1,t>0,with Dirichlet boundary conditions(29)u0,t=u0,u1,t=u-,t>0.
The exact solution of the problem is given by [10, 11](30)ux,t=0.1e-A+0.5e-B+e-Ce-A+e-B+e-C,where(31)A=0.05Rex-0.5+4.95t,B=0.25Rex-0.5+0.75t,C=0.5Rex-0.375.
For this problem, extensive numerical results indicated that the solutions are smooth when Re≤200 [11]; as the Re number increases, the solutions will appear discontinuous.
Figure 2 shows the VMEFG solutions without filter and with filter at different times and Re = 2000. It can be seen that VMEFG solutions have some spurious oscillation under this larger Re number, while the solutions obtained by VMEFG method with filter are still smooth. This example also illustrates that VMEFG method produces the numerical pseudo oscillations for very small diffusivity coefficient. Thus, other techniques need to overcome the difficulty.
Figure 3 presents the comparison of solutions obtained by VMEFG method with filter and exact solutions at different times. It can be seen that the VMEFG method with filter solutions are in well agreement with the exact solutions. From Figures 2 and 3, it can be seen that the VMEFG method coupled high order low-pass filter is one of the techniques which can work out the difficulty well.
In order to investigate the accuracy of the VMEFG method with filter, the following L2 norm error is used:(32)e2=h∑i=1Nuinum-uiexact21/2,where h is the space step, N is the number of node, unum is the numerical solution, and uexact is the exact solution. Table 1 shows the L2 norm errors by taking different numbers of node at different times for the simulation and observes that the numerical results are satisfactory. Figure 4 illustrates these errors plotted against h on a log-log scale. It can be noted that the order of convergence of the VMEFG method with filter is linear.
The L2 errors obtained by VMEFG method with filter at t=0.1, 0.3, and 0.5.
h
t=0.1
t=0.3
t=0.5
1/200
0.032125
0.046050
0.064297
1/400
0.016440
0.024454
0.034351
1/800
0.009747
0.015153
0.020432
1/1600
0.004423
0.006895
0.009768
The numerical solutions of Re = 2000 obtained by VMEFG method with and without filter, respectively. (a) is obtained without filter and (b) is with low-pass filter.
Comparison of VMEFG method with filter and exact solution at t=0.1, 0.3, and 0.5.
t=0.1
t=0.3
t=0.5
The errors from Table 1 plotted against h on a log-log scale.
Example 2.
In this example, we consider the 1D Burgers’ equation(33)∂u∂t+u∂u∂x=1Re∂2u∂x2,0<x<1,t>0,with the initial condition(34)ux,0=sinπx,0<x<1,and the boundary conditions(35)u0,t=u1,t=0.0,t>0.
The exact solution of this problem is [10, 30](36)ux,t=2π1Re∑n=1∞anexp-n2π2t/Rensinnπxa0+∑n=1∞anexp-n2π2t/Recosnπx,where the Fourier coefficients are(37)a0=∫01exp-Re2π1-cosπxdx,an=∫01exp-Re2π1-cosπxcosnπxdxn=1,2,3,….
Although the example has the exact solution, the convergence is very slow for larger value of Re owing to that the exact solution of the equation involves Fourier series solutions. Therefore, in order to show the computational accuracy of the proposed method for small diffusivity coefficients, that is, large Re values, the numerical solutions are compared with the results available in the literatures [3, 30].
In this example, 200 nodes are distributed uniformly in the domain and Δt=0.001. Figure 5 presents the numerical solutions obtained by EFG method and VMEFG method for Re = 103. It can be seen that there exists some oscillation in the solution of EFG, but, in contrast, the solutions of VMEFG are always smooth, which indicates that VMEFG method has the ability to accommodate the exact solution to the problem. Meanwhile, this demonstrates that the VMEFG has better stability and accuracy than EFG method for moderate large Re number indeed.
Figure 6 illustrates the numerical results obtained by VMEFG without and with low-pass filter at different time for Re = 104. It can be observed that the profiles of u develop extremely sharp gradients near the right-hand boundary for t≥0.4, and the numerical results of VMEFG method without filter show the nonphysical oscillations near the right-hand boundary for t≥0.6, which indicates that the VMEFG method still has spurious oscillations for very large Re numbers, whereas the numerical results obtained by VMEFG method with low-pass filter are still smooth and meanwhile the sharp gradients of the profiles near x=1, especially for t=0.6, are well distinguished.
Table 2 shows the results obtained by VMEFG method with low-pass filter at t=1.0 for Re = 104 and those available in the literatures. It can be found that our results are in agreement with the results presented in [30] for Re = 104. From Figure 6 and Table 2, we can conclude that the VMEFG method with low-pass filter gives remarkable accuracy for large values of Re.
The numerical solution of Burgers’ equation at a high Reynolds number (Re≥105) is very difficult due to the presence of shock [21]. A direct application of the VMEFG method for Re = 105 leads to highly oscillatory results as shown in Figure 7. It can be seen that the oscillation starts at t=0.4 and is accumulated and amplified in the following times. Figure 7 also shows the results obtained by using the low-pass filter under the same condition. It can be found that the oscillations are eliminated and sharp shock profiles are resolved. Meanwhile, we consider Re = 109 further. Figure 8 shows the results obtained without and with low-pass filter for Re = 109. Figure 8(a) is obtained without filter; we see some numerical oscillations obviously. But the low-pass filter eliminates those high frequency oscillations very well.
Comparison between the numerical solutions at t=1.0 for Re = 104.
x
Reference [3]
Reference [30]
VMEFG with filter
0.05
0.0424
0.0379
0.037716792
0.11
0.0843
0.0834
0.082959609
0.16
0.1263
0.1213
0.120632606
0.22
0.1684
0.1667
0.165785694
0.27
0.2103
0.2044
0.203352212
0.33
0.2522
0.2497
0.248337909
0.38
0.2939
0.2872
0.285730246
0.44
0.3355
0.3322
0.330462517
0.50
0.3769
0.3769
0.375015652
0.55
0.4182
0.4141
0.411982983
0.61
0.4592
0.4584
0.456120122
0.66
0.4999
0.4951
0.492687106
0.72
0.5404
0.5388
0.536270966
0.77
0.5805
0.5749
0.572307656
0.83
0.620
0.6179
0.615158458
0.88
0.6600
0.6533
0.650490215
0.94
0.6957
0.6952
0.692359873
0.96
—
0.7090
0.706171615
0.98
—
0.7228
0.719903495
0.99
—
0.7296
0.728178293
The numerical solutions of Re = 103 via EFG and VMEFG, respectively. (a) is obtained by EFG and (b) is obtained by VMEFG.
Effects of filtering for solving one-dimensional Burgers’ equation for Re = 104 by the VMEFG method. (a) is obtained without filter and (b) is with low-pass filter.
Effects of filtering for solving one-dimensional Burgers’ equation for Re = 105 by the VMEFG method. (a) is obtained without filter and (b) is with low-pass filter.
Effects of filtering for solving one-dimensional Burgers’ equation for Re = 109 by the VMEFG method. (a) is obtained without filter and (b) is with low-pass filter.
Example 3.
For the third example, we consider the coupled nonlinear system of Burgers’ equation which has been studied by Zhang et al. [14] and Donea et al. [7]. The coupled nonlinear system of Burgers’ equation is considered over the square domain Ω=[0,1]×[0,1] and defined as the following:(38)∂u∂t+u∂u∂x+v∂u∂y=1Re∂2u∂x2+∂2u∂y2,∂v∂t+u∂v∂x+v∂v∂y=1Re∂2v∂x2+∂2v∂y2,with the initial conditions(39)ux,y,0=sinπxcosπy,x,y∈Ω,vx,y,0=cosπxsinπy,x,y∈Ω,and the boundary conditions(40)u0,y,t=u1,y,t=0.0,t>0,vx,0,t=vx,1,t=0.0t>0,∂u∂nx,0,t=∂u∂nx,1,t=0.0t>0,∂v∂n0,y,t=∂v∂n1,y,t=0.0t>0.
The above problem has various symmetries. As Zhang et al. [14] pointed out, there exist u(x,y,t)=v(y,x,t) and u(x,y,t)=-u(1-x,1-y,t) over the global domain and u(x,x,t)=v(x,x,t) and u(x,x,t)=-u(1-x,1-x,t) along the transverse section of the domain. Thus, in the following, we only present the results of u(x,y,t).
This example was solved by many researchers and in their studies, the maximum value of Re is the order of 103. Meanwhile, from those results available in the literature [14], we can see that spurious oscillation still appeared for Re = 103.
In this study, 51 × 51 nodes are distributed uniformly in the domain and Δt=0.01. Figure 9 illustrates the comparison for effects of low-pass filter. It can be seen that the spurious oscillation also exists in our results obtained by using VMEFG method without filter, whereas the results obtained with low-pass filter are excellent; clearly, all oscillations are effectively removed and the shock front is very sharp. To quantify the computational accuracy and stabilization of VMEFG with low-pass filter, the results of u along the domain diagonal (y=x) are presented in Figure 10. It can be seen obviously that the results obtained by VMEFG have some spurious oscillation, while the results obtained by VMEFG with low-pass filter have no spurious oscillation and the sharp shock profiles are resolved.
To test the VMEFG method with low-pass filter for solving this coupled Burgers’ equataion further, we next consider the Re = 104 and the results are presented in Figure 11. To the best of our knowledge, there have been no results reported on so large Re value. The top four figuers are obtained without filter; we see that some numerical oscillations appear along the diagnonal of the domain passing through the points (0,1) and (1,0). But the low-pass filter eliminates those very well and the shocks along the diagonal of the domain are also very sharp. Figure 12 shows the results obtained by VMEFG with low-pass filter along the domain diagonal (y=x) for Re = 104. It can be noted obviously that the solutions are also smooth and the sharp shock profiles are distinguished.
2D Burgers’ equation for Re = 103. (a) are computed without filter; (b) are with the low-pass filter.
Solutions u(x,x,t) of the VMEFG wihout filter (a) and with low-pass filter (b) at different times for Re = 103.
2D Burgers’ equation for Re = 104. (a) are computed without filter; (b) are with the low-pass filter.
Solutions u(x,x,t) of the VMEFG with low-pass filter at different times for Re = 104.
4. Conclusions
In this paper, VMEFG method in conjunction with low-pass filter has been applied for the numerical solution of Burgers’ equations with small diffusivity coefficients, which lead to convection dominated problem. To the best of our knowledge, this is the first time for meshless methods to be used for solving Burgers’ equation with such small diffusivity coefficients. The accuracy and stabilization of the proposed method are validated by three numerical examples. Compared with numerical results available in the literatures, it can be seen that the VMEFG method coupled with low-pass filter produces better results with ease of implementation for Burgers’ equation with very large Re numbers.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is financially supported by the Natural Sciences Foundation of China (nos. 11102101, 11171181, and 61273183) and the Academic Mainstay Foundation of Hubei Province of China (no. Q20111208).
ZhangP.-G.WangJ.-P.A predictor-corrector compact finite difference scheme for Burgers' equation2012219389289810.1016/j.amc.2012.06.064MR29812812-s2.0-84867337996YousufM.On the class of high order time stepping schemes based on Padé approximations for the numerical solution of Burgers' equation2008205144245310.1016/j.amc.2008.08.025MR24666502-s2.0-54049108803HonY. C.MaoX. Z.An efficient numerical scheme for Burgers' equation1998951375010.1016/S0096-3003(97)10060-1MR16302602-s2.0-0002879831ÖzisT.AksanE. N.ÖzdesA.A finite element approach for solution of Burgers' equation2003139417428DoganA.A Galerkin finite element approach to Burgers' equation2004157233134610.1016/j.amc.2003.08.037MR20882572-s2.0-4444378330HuertaA.RoigB.DoneaJ.Time-accurate solution of stabilized convection–diffusion–reaction equations: II—accuracy analysis and examples200218857558410.1002/cnm.518MR19183192-s2.0-0036684315DoneaJ.RoigB.HuertaA.High-order accurate time-stepping schemes for convection-diffusion problems20001823-424927510.1016/s0045-7825(99)00193-0MR17442492-s2.0-0034681589Guo-ZhongZ.Xi-JunY.DiW.Numerical solution of the Burgers' equation by local discontinuous Galerkin method2010216123671367910.1016/j.amc.2010.05.021MR26617322-s2.0-77953609346AroraG.SinghB. K.Numerical solution of Burgers' equation with modified cubic B-Spline differential quadrature method201322416617710.1016/j.amc.2013.08.071MR31276092-s2.0-84884693508DhawanS.KapoorS.KumarS.RawatS.Contemporary review of techniques for the solution of nonlinear Burgers equation20123540541910.1016/j.jocs.2012.06.0032-s2.0-84865388484KapoorS.DhawanS.A computational technique for the solution of Burgers' equation200868495GaoY.LeL.-H.ShiB.-C.Numerical solution of Burgers' equation by lattice Boltzmann method2013219147685769210.1016/j.amc.2013.01.056MR30326082-s2.0-84875928875HuertaA.MéndezS. F.Time accurate consistently stabilized mesh-free methods for convection dominated problems20035691225124210.1002/nme.6022-s2.0-0037423941ZhangX. H.OuyangJ.ZhangL.Element-free characteristic Galerkin method for Burgers' equation200933335636210.1016/j.enganabound.2008.07.001MR24811252-s2.0-58149503030GuY. T.LiuG. R.Meshless techniques for convection dominated problems200638217118210.1007/s00466-005-0736-8MR22292632-s2.0-33646422497FriesT. P.MatthiesH. G.A stabilized and coupled meshfree/meshbased method for the incompressible Navier-Stokes equations-part I: stabilization200619544–476205622410.1016/j.cma.2005.12.002FriesT. P.MatthiesH. G.A stabilized and coupled meshfree/meshbased method for the incompressible Navier-Stokes equations—part II: coupling200619561916204ZhangL.OuyangJ.ZhangX. H.ZhangW.On a multi-scale element-free Galerkin method for the Stokes problem2008203274575310.1016/j.amc.2008.05.081MR24589902-s2.0-51549088767ZhangL.OuyangJ.JiangT.RuanC. L.Variational multiscale element free Galerkin method for the water wave problems2011230125045506010.1016/j.jcp.2011.03.026MR27959952-s2.0-79955643974ZhangL.OuyangJ.ZhangX. H.The variational multisclae element free Galerkin method for MHD flows at high Hartmann numbers201318411061118ZhangX. H.XiangH.Variational multiscale element free Galerkin method for convection-diffusion-reaction equation with small diffusion201446859210.1016/j.enganabound.2014.05.010MR32295792-s2.0-84902689668WeiG. W.GuY.Conjugate filter approach for solving Burgers' equation2002149243945610.1016/s0377-0427(02)00488-0MR19372942-s2.0-0037114626ZhouY. C.WeiG. W.High resolution conjugate filters for the simulation of flows2003189115917910.1016/S0021-9991(03)00206-7MR1988145ZBL1097.765812-s2.0-0038454600LiJ.VisbalM. R.High-order compact schemes for nonlinear dispersive waves200626112310.1007/s10915-004-4797-1MR2221177ZBL1089.760432-s2.0-32944461101NguyenV. P.RabczukT.BordasS.DuflotM.Meshless methods: a review and computer implementation aspects200879376381310.1016/j.matcom.2008.01.003MR2477562ZBL1152.740552-s2.0-55949119340HughesT. J.FeijooG. R.MazzeiL.QuincyJ.-B.The variational multiscale method—a paradigm for computational mechanics19981661-232410.1016/s0045-7825(98)00079-6MR16601412-s2.0-0032203197MasudA.KhurramR. A.A multiscale finite element method for the incompressible Navier-Stokes equations200619513–161750177710.1016/j.cma.2005.05.048MR22039912-s2.0-30944434045LeleS. K.Compact finite difference schemes with spectral-like resolution19921031164210.1016/0021-9991(92)90324-rMR11880882-s2.0-9144220381SenguptaT. K.BhumkarY. G.New explicit two-dimensional higher order filters201039101848186310.1016/j.compfluid.2010.06.014ZBL1245.761282-s2.0-77956338423HashemianA.ShodjaH. M.A meshless approach for solution of Burgers' equation20082201-222623910.1016/j.cam.2007.08.014MR2444167ZBL1149.650792-s2.0-47849110861