Advances in Mathematical Physics
Volume 2016 (2016), Article ID 1836978, 12 pages
http://dx.doi.org/10.1155/2016/1836978

Research Article

Variational Multiscale Element Free Galerkin Method Coupled with Low-Pass Filter for Burgers’ Equation with Small Diffusion

Ping Zhang, Xiaohua Zhang, and Laizhong Song

College of Science, China Three Gorges University, Yichang 443002, China

Received 22 October 2015; Revised 27 January 2016; Accepted 2 February 2016

Academic Editor: Stephen C. Anco

Copyright © 2016 Ping Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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 be an open bounded region with piecewise smooth boundary . The number of space dimensions, , is equal to 1 or 2. Burgers’ equation is written as follows:

where 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 represent the weight function for velocity ; then the standard weak form of (1) is given as follows:

where , , that is, the 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 and its weight function can be decomposed into coarse scale and fine scale, respectively; namely,

where and are the coarse scale and fine scale for velocity and and are the coarse scale and fine scale for weight function , respectively.
Similar to literatures [19, 27], we further assume that and , although nonzero within background integral cell , vanish identically over the boundary when the factor approaches to 1:

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

Next, employing the linearity between coarse scale and fine scale , (7) can be split into coarse scale problem and fine scale problem as follows.
The Coarse Scale Problem . Consider

The Fine Scale Problem . Consider

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 in (8) and (9) by the last converged solution from the fixed point iteration.
The Linearized Coarse Scale Problem . Consider

The Linearized Fine Scale Problem . Consider

Step 2 (solution of the fine scale problem). According to the literatures [19, 27], we next deduce the fine scale solution by using the bubble functions. Let us now consider the fine scale part of the weak form . 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

where . 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 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 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 and are represented via bubble functions over the of the EFG shape function when factor approaches to 1:

where and represent the bubble functions for the trial solution and weight function over the , respectively, and and are their coefficients, respectively, ( 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]:

where is the location of the internal virtual node for the piecewise linear bubble and in the paper .
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]

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:

where and represent the location of the internal virtual node in the background cell and in the paper .
Substituting (13) and (14) into (12), and taking the constant coefficient vector and out of the integral expression, we have

Since is arbitrary, consequently we get

where and are defined as follows:

where is a identity matrix and both and are vectors of gradient of the bubble functions. Substituting (19) into (13), then we can obtain the fine scale over the as follows:

Figure 1: Schematic diagram showing the domains for the bubble functions for a reference triangle cell [21].

Step 3 (solution of the coarse scale problem ). Once we have obtained the fine scale solution , now we turn our attention to the coarse scale problem . Exploiting the linearity of the solution slot between and , (10) can be written as

Using the integration by parts for the terms containing in (22) as well as the assumptions used in (5), (22) has the following forms:

Substituting (18) into (23), subsequently assuming that the coarse scale residual is a constant over the , after rearrangement we get the finally computed formulation:

where

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]

where denotes the given value at point and denotes the value after filtering. With proper choices of coefficients, (22) provides a th-order formulation on a point stencil. The coefficients can be obtained in terms of by Taylor expansion. In the paper, we use the following coefficients of a family of 10th-order filter derived in [24, 28]

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:

with Dirichlet boundary conditions

The exact solution of the problem is given by [10, 11]

where

For this problem, extensive numerical results indicated that the solutions are smooth when [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 norm error is used:

where is the space step, is the number of node, is the numerical solution, and is the exact solution. Table 1 shows the 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 on a log-log scale. It can be noted that the order of convergence of the VMEFG method with filter is linear.

Table 1: The errors obtained by VMEFG method with filter at , 0.3, and 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

(a)

(b)

Figure 2: 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.

(a)

(b)

(c)

Figure 3: Comparison of VMEFG method with filter and exact solution at , 0.3, and 0.5.

Figure 4: The errors from Table 1 plotted against on a log-log scale.

Example 2. In this example, we consider the 1D Burgers’ equation

with the initial condition

and the boundary conditions

The exact solution of this problem is [10, 30]

where the Fourier coefficients are

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 . Figure 5 presents the numerical solutions obtained by EFG method and VMEFG method for Re = 10³. 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 = 10⁴. It can be observed that the profiles of develop extremely sharp gradients near the right-hand boundary for , and the numerical results of VMEFG method without filter show the nonphysical oscillations near the right-hand boundary for , 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 , especially for , are well distinguished.
Table 2 shows the results obtained by VMEFG method with low-pass filter at for Re = 10⁴ and those available in the literatures. It can be found that our results are in agreement with the results presented in [30] for Re = 10⁴. 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 () is very difficult due to the presence of shock [21]. A direct application of the VMEFG method for Re = 10⁵ leads to highly oscillatory results as shown in Figure 7. It can be seen that the oscillation starts at 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 = 10⁹ further. Figure 8 shows the results obtained without and with low-pass filter for Re = 10⁹. 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.

Table 2: Comparison between the numerical solutions at for Re = 10⁴.

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

(a)

(b)

Figure 5: The numerical solutions of Re = 10³ via EFG and VMEFG, respectively. (a) is obtained by EFG and (b) is obtained by VMEFG.

(a)

(b)

Figure 6: Effects of filtering for solving one-dimensional Burgers’ equation for Re = 10⁴ by the VMEFG method. (a) is obtained without filter and (b) is with low-pass filter.

(a)

(b)

Figure 7: Effects of filtering for solving one-dimensional Burgers’ equation for Re = 10⁵ by the VMEFG method. (a) is obtained without filter and (b) is with low-pass filter.

(a)

(b)

Figure 8: Effects of filtering for solving one-dimensional Burgers’ equation for Re = 10⁹ 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 and defined as the following:

with the initial conditions

and the boundary conditions

The above problem has various symmetries. As Zhang et al. [14] pointed out, there exist and over the global domain and and along the transverse section of the domain. Thus, in the following, we only present the results of .
This example was solved by many researchers and in their studies, the maximum value of Re is the order of 10³. Meanwhile, from those results available in the literature [14], we can see that spurious oscillation still appeared for Re = 10³.
In this study, 51 × 51 nodes are distributed uniformly in the domain and . 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 along the domain diagonal () 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 = 10⁴ 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 () and (). 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 () for Re = 10⁴. It can be noted obviously that the solutions are also smooth and the sharp shock profiles are distinguished.

(a)

(b)

Figure 9: 2D Burgers’ equation for Re = 10³. (a) are computed without filter; (b) are with the low-pass filter.

(a)

(b)

Figure 10: Solutions of the VMEFG wihout filter (a) and with low-pass filter (b) at different times for Re = 10³.

(a)

(b)

Figure 11: 2D Burgers’ equation for Re = 10⁴. (a) are computed without filter; (b) are with the low-pass filter.

Figure 12: Solutions of the VMEFG with low-pass filter at different times for Re = 10⁴.

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).

References

P.-G. Zhang and J.-P. Wang, “A predictor-corrector compact finite difference scheme for Burgers' equation,” Applied Mathematics and Computation, vol. 219, no. 3, pp. 892–898, 2012.

M. Yousuf, “On the class of high order time stepping schemes based on Padé approximations for the numerical solution of Burgers' equation,” Applied Mathematics and Computation, vol. 205, no. 1, pp. 442–453, 2008.

Y. C. Hon and X. Z. Mao, “An efficient numerical scheme for Burgers' equation,” Applied Mathematics and Computation, vol. 95, no. 1, pp. 37–50, 1998.

T. Özis, E. N. Aksan, and A. Özdes, “A finite element approach for solution of Burgers' equation,” Applied Mathematics and Computation, vol. 139, pp. 417–428, 2003.

A. Dogan, “A Galerkin finite element approach to Burgers' equation,” Applied Mathematics and Computation, vol. 157, no. 2, pp. 331–346, 2004.

A. Huerta, B. Roig, and J. Donea, “Time-accurate solution of stabilized convection–diffusion–reaction equations: II—accuracy analysis and examples,” Communications in Numerical Methods in Engineering, vol. 18, no. 8, pp. 575–584, 2002.

J. Donea, B. Roig, and A. Huerta, “High-order accurate time-stepping schemes for convection-diffusion problems,” Computer Methods in Applied Mechanics and Engineering, vol. 182, no. 3-4, pp. 249–275, 2000.

Z. Guo-Zhong, Y. Xi-Jun, and W. Di, “Numerical solution of the Burgers' equation by local discontinuous Galerkin method,” Applied Mathematics and Computation, vol. 216, no. 12, pp. 3671–3679, 2010.

G. Arora and B. K. Singh, “Numerical solution of Burgers' equation with modified cubic B-Spline differential quadrature method,” Applied Mathematics and Computation, vol. 224, pp. 166–177, 2013.

S. Dhawan, S. Kapoor, S. Kumar, and S. Rawat, “Contemporary review of techniques for the solution of nonlinear Burgers equation,” Journal of Computational Science, vol. 3, no. 5, pp. 405–419, 2012.

S. Kapoor and S. Dhawan, “A computational technique for the solution of Burgers' equation,” International Journal of Applied Mathematics and Mechanics, vol. 6, pp. 84–95, 2008.

Y. Gao, L.-H. Le, and B.-C. Shi, “Numerical solution of Burgers' equation by lattice Boltzmann method,” Applied Mathematics and Computation, vol. 219, no. 14, pp. 7685–7692, 2013.

A. Huerta and S. F. Méndez, “Time accurate consistently stabilized mesh-free methods for convection dominated problems,” International Journal for Numerical Methods in Engineering, vol. 56, no. 9, pp. 1225–1242, 2003.

X. H. Zhang, J. Ouyang, and L. Zhang, “Element-free characteristic Galerkin method for Burgers' equation,” Engineering Analysis with Boundary Elements, vol. 33, no. 3, pp. 356–362, 2009.

Y. T. Gu and G. R. Liu, “Meshless techniques for convection dominated problems,” Computational Mechanics, vol. 38, no. 2, pp. 171–182, 2006.

T. P. Fries and H. G. Matthies, “A stabilized and coupled meshfree/meshbased method for the incompressible Navier-Stokes equations-part I: stabilization,” Computer Methods in Applied Mechanics and Engineering, vol. 195, no. 44–47, pp. 6205–6224, 2006.

T. P. Fries and H. G. Matthies, “A stabilized and coupled meshfree/meshbased method for the incompressible Navier-Stokes equations—part II: coupling,” Computer Methods in Applied Mechanics and Engineering, vol. 195, pp. 6191–6204, 2006.

L. Zhang, J. Ouyang, X. H. Zhang, and W. Zhang, “On a multi-scale element-free Galerkin method for the Stokes problem,” Applied Mathematics and Computation, vol. 203, no. 2, pp. 745–753, 2008.

L. Zhang, J. Ouyang, T. Jiang, and C. L. Ruan, “Variational multiscale element free Galerkin method for the water wave problems,” Journal of Computational Physics, vol. 230, no. 12, pp. 5045–5060, 2011.

L. Zhang, J. Ouyang, and X. H. Zhang, “The variational multisclae element free Galerkin method for MHD flows at high Hartmann numbers,” Computer Physics Communications, vol. 184, pp. 1106–1118, 2013.

X. H. Zhang and H. Xiang, “Variational multiscale element free Galerkin method for convection-diffusion-reaction equation with small diffusion,” Engineering Analysis with Boundary Elements, vol. 46, pp. 85–92, 2014.

G. W. Wei and Y. Gu, “Conjugate filter approach for solving Burgers' equation,” Journal of Computational and Applied Mathematics, vol. 149, no. 2, pp. 439–456, 2002.

Y. C. Zhou and G. W. Wei, “High resolution conjugate filters for the simulation of flows,” Journal of Computational Physics, vol. 189, no. 1, pp. 159–179, 2003.

J. Li and M. R. Visbal, “High-order compact schemes for nonlinear dispersive waves,” Journal of Scientific Computing, vol. 26, no. 1, pp. 1–23, 2006.

V. P. Nguyen, T. Rabczuk, S. Bordas, and M. Duflot, “Meshless methods: a review and computer implementation aspects,” Mathematics and Computers in Simulation, vol. 79, no. 3, pp. 763–813, 2008.

T. J. Hughes, G. R. Feijoo, L. Mazzei, and J.-B. Quincy, “The variational multiscale method—a paradigm for computational mechanics,” Computer Methods in Applied Mechanics and Engineering, vol. 166, no. 1-2, pp. 3–24, 1998.

A. Masud and R. A. Khurram, “A multiscale finite element method for the incompressible Navier-Stokes equations,” Computer Methods in Applied Mechanics and Engineering, vol. 195, no. 13–16, pp. 1750–1777, 2006.

S. K. Lele, “Compact finite difference schemes with spectral-like resolution,” Journal of Computational Physics, vol. 103, no. 1, pp. 16–42, 1992.

T. K. Sengupta and Y. G. Bhumkar, “New explicit two-dimensional higher order filters,” Computers and Fluids, vol. 39, no. 10, pp. 1848–1863, 2010.

A. Hashemian and H. M. Shodja, “A meshless approach for solution of Burgers' equation,” Journal of Computational and Applied Mathematics, vol. 220, no. 1-2, pp. 226–239, 2008.