Research Article  Open Access
A Meshfree QuasiInterpolation Method for Solving Burgers’ Equation
Abstract
The main aim of this work is to consider a meshfree algorithm for solving Burgers’ equation with the quartic Bspline quasiinterpolation. Quasiinterpolation is very useful in the study of approximation theory and its applications, since it can yield solutions directly without the need to solve any linear system of equations and overcome the illconditioning problem resulting from using the Bspline as a global interpolant. The numerical scheme is presented, by using the derivative of the quasiinterpolation to approximate the spatial derivative of the dependent variable and a low order forward difference to approximate the time derivative of the dependent variable. Compared to other numerical methods, the main advantages of our scheme are higher accuracy and lower computational complexity. Meanwhile, the algorithm is very simple and easy to implement and the numerical experiments show that it is feasible and valid.
1. Introduction
Burgers’ equation plays a significant role in various fields, such as turbulence problems, heat conduction, shock waves, continuous stochastic processes, number theory, gas dynamics, and propagation of elastic waves [1–5]. The onedimensional Burgers’ equation first suggested by Bateman [6] and later treated by Burgers [1] has the form where is the coefficient of kinematic viscosity and the subscripts and denote space and time derivatives. Initial and boundary conditions are where , , and will be chosen in a later section.
Burgers’ equation is a quasilinear parabolic partial differential equation, whose analytic solutions can be constructed from a linear partial differential equation by using HopfCole transformation [1, 2, 7]. But some analytic solutions consist of infinite series, converging very slowly for small viscosity coefficient . Thus, many researchers have spent a great deal of effort to compute the solution of Burgers’ equation using various numerical methods. Finite difference methods were presented to solve the numerical solution of Burgers’ equation in [8–11]. Finite element methods for the solution of Burgers’ equation were introduced in [12–15]. Recently, various powerful mathematical methods such as Galerkin finite element method [16, 17], spectral collocation method [18, 19], sinc differential quadrature method [20], factorized diagonal padé approximation [21], Bspline collocation method [22], and reproducing kernel function method [23] have also been used in attempting to solve the equation.
In 1968 Hardy proposed the multiquadric (MQ) which is a kind of radial basis function (RBF). In Franke’s review paper, the MQ was rated as one of the best methods among 29 scattered data interpolation and ease of implementation. Since Kansa successfully applied MQ for solving partial differential equation, more and more reasearchers have been attracted by this meshfree, scattered data approximation scheme [24]. The meshfree method uses a set of scattered nodes, instead of meshing the domain of the problem. It has been successfully applied to solve many physical and engineering problems with only a minimum of meshing or no meshing at all [25–30]. In recent years, many meshfree metheods have been developed, such as the elementfree Galerkin method [31], the smooth particle hydrodynamics method [32], the elementfree kpRitz method [33–36], the meshless local PetrovGalerkin method [37], and the reproducing kernel particle method [38].
With the use of univariate multiquadric (MQ) quasiinterpolation, solution of Burgers’ equations was obtained by Chen and Wu [39]. Moreover, Hon and Mao [40] developed an efficient numerical scheme for Burgers’ equation applying the MQ as a spatial approximate scheme and a low order explicit finite difference approximation to the time derivation. Zhu and Wang [41] presented the numerical scheme for solving the Burgers’ equation, by using the derivative of the cubic Bspline quasiinterpolation to approximate the time derivative of the dependent variable and a low order forward difference to approximate the time derivative of the dependent variable. In this paper, we provide a numerical scheme to solve Bugers’ equation using the quartic Bspline quasiinterpolation. Then we do not require to solve any linear system of equation so that we do not meet the question of the illcondition of the matrix.Therefore, we can solve the computational time and decrease the numerical error.
This paper is arranged as follows. In Section 2, the definition of quartic Bspline has been described and univariate quartic Bspline quasiinterpolants have been presented. In Section 3, we mainly propose the numerical techniques using quartic Bspline interpolation to solve Burgers’ equation. In Section 4, numerical examples of Burgers’ equation are presented and compared with those obtained with some previous results. At last, we conclude the paper in Section 5.
2. Univariate Quartic BSpline QuasiInterpolant
For an interval , we introduce a set of equallyspaced knots of partition . We assume that , , , and . Let be the space of continuouslydifferentiable, piecewise, quarticdegree polynomials on . A detailed description of Bspline functions generated by subdivision regarding the Bsplines basis in can be found in [45].
The zero degree Bspline is defined as and, for positive constant , it is defined in the following recursive form:
We apply this recursion to get the quartic Bspline , which is defined in as follows:As usual, we add multiple knots at the endpoints: and .
In [24], univariate quartic Bspline quasiinterpolants can be defined as operators of the form The coefficients are listed as follows: and , , . For , we have the error estimate
We use to denote the space of polynomials of the total degree at most . In general, we impose that is exact on the space ; that is, for all . As a consequence of this property, the approximation order of is on smooth functions. In this paper, the coefficient is a linear combination of discrete values of at some points. The main advantage of is that they have a direct construction without solving any system of linear equations. Moreover, they are local in the sense that values of depend only on values of in a neighborhood of . Finally, they have a rather small infinity norm and, therefore, are nearly optimal approximant.
Differentiating interpolation polynomials leads to the classic finite difference for the approximate computation of derivatives. Therefore, we can draw a conclusion of approximating derivatives of by derivatives of . The general theory will be developed elsewhere. We can evaluate the value of at by and . and can be computed by the formula of Bspline’s derivatives as follows: where
By some trivial computations, we can obtain the value of at the knots, which are illustrated in Table 1. Then, we get the differential formulas for quartic Bspline QIs as

3. Numerical scheme Using the Meshfree QuasiInterpolation
In this section, we present the numerical scheme for solving Burgers’ equation based on the quartic Bspline quasiinterpolation.
Discretizing the Burgers’ equation in time with meshlength , we get
We can get where is the approximation of the value of at the point . Then, we can use the derivatives of the quartic Bspline quasiinterpolant to approximate and . To dump the dispersion of the scheme, we define a switch function , whose values are and at the discrete points , as follows: where . Thus, the resulting numerical scheme is
Starting from the initial condition, we can compute the numerical solution of Burgers’ equation step by step using the Bspline quasiinterpolation scheme (16) and formulas (11).
4. Numerical Results
To investigate the applicability of the quasiinterpolation method to Burgers’ equation, four selected example problems are studied. To show the efficiency of the present method for our problem in comparison with the exact solution, we use the following norms to assess the performance of our scheme:
Example 1. Burgers’ equation is solved over the region and the initial and boundary conditions are given in Asaithambi [42]: and the exact solution of this problem has the following nice compact closedform, as given by Wood [46]:
In this computational study, we set , , . The comparison of the numerical solutions obtained by the present method, at the different coefficient of kinematic viscosity , are presented with the solutions obtained by Asaithambi [42] and the exact solution in Table 2.

Example 2. In this example, we consider the exact solution of Burgers’ equation [47]: where , , , and are constants. The boundary conditions are and initial condition is used for the exact solution at .
We solve the problem with , , and by our method. In Table 3, and errors at the time level are compared with the error obtained by Chen and Wu [39], Zhu and Wang [41], Da et al. [17], and Saka and Da [43]. For comparison, the parameters are adopted as time step , space step , and viscosity coefficient . From Table 3, we can find that our method provides better accuracy than most methods through the and error norms. The profiles of initial wave and its propagation are depicted at some times in Figure 1.
Example 3. Consider Burgers’ equation with the initial condition and the boundary conditions
The analytical solution of this problem was given by Cole [2] in the term of an infinite series as with the Fourier coefficients
In Table 4, we have computed the numerical solutions of this example at differential time levels with parameter values , , and . The comparison of our results with the exact solutions as well as the solutions obtained in [11, 15, 44] is reported in Table 4. From Table 4, we can find that the presented scheme provides better accuracy. Moreover, in Tables 5, 6 and 7, we compare our method with Hon and Mao’s scheme, Chen and Wu’s MQQI method, and Zhu’s BSQI method at with , for , respectively. For the MQQI method, the shape parameter , , for Table 5, respectively, as [39]. Solutions found with the present method are in good agreement with the result and better than other methods. These show that the method works well.



Example 4. We consider particular solution of Burgers’ equation: where . Initial condition is obtained from when is used. Boundary conditions are . Analytical solution represents shocklike solution of the onedimensional Burgers’ equation. Parameters and are selected for comparison over the domain . Accuracy of our method is shown by calculating the error norms. These together with some previous results are given in Table 8. Table 8 shows that our method provides better accuracy than MQQI method and BSQI method. Although the accuracy is not higher than that of QBCM method, we know that, at each time step, the complexity of our method is lower than theirs. The numerical solutions are depicted with , , and for in Figure 2.

5. Conclusion
Following the recent development of the quasiinterpolation method for scattered data interpolation and the meshfree method for solving partial differential equations, this paper combines these ideas and proposes a new meshfree quasiinterpolation method for Burgers’ equation. The method does not require solving a large size matrix equation and, hence, the illconditioning problem from using Bspline functions as global interpolants can be avoided. We have made comparison studies between the present results and the exact solutions. The agreement of our numerical results with those exact solutions is excellent. For the highdimensional Burgers’ equations, we believe our scheme can also be applicable. In this case, we would use multivariate spline quasiinterpolation instead of univariate spline quasiinterpolation. We will consider these problems in our future work.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by the Disciplinary Construction Guide Foundation of Harbin Institute of Technology at Weihai (no. WH20140206) and the Scientific Research Foundation of Harbin Institute of Technology at Weihai (no. HIT(WH)201319).
References
 J. M. Burgers, “A mathematical model illustrating the theory of turbulence,” Advances in Applied Mechanics, vol. 1, pp. 171–199, 1948. View at: Google Scholar
 J. D. Cole, “On a quasilinear parabolic equation occurring in aerodynamics,” Quarterly of Applied Mathematics, vol. 9, no. 3, pp. 225–236, 1951. View at: Google Scholar  MathSciNet
 M. J. Lighthill, “Viscosity effects in sound waves of finite amplitude,” in Surveys in Mechanics, pp. 250–351, 1956. View at: Google Scholar  MathSciNet
 L. A. Pospelov, “Propagation of finite amplitude elastic waves(Longitudinal elastic wave of finite amplitude propagation in isotropic solid),” Soviet PhysicsAcoustics, vol. 11, pp. 302–304, 1966. View at: Google Scholar
 B. van der Pol, “On a nonlinear partial differential equation satisfied by the logarithm of the Jacobian thetafunctions, with arithmetical applications. I, II,” Indagationes Mathematicae, vol. 13, pp. 261–284, 1951. View at: Google Scholar  MathSciNet
 H. Bateman, “Some recent researches on the motion of fluids,” Monthly Weather Review, vol. 43, no. 4, pp. 163–170, 1915. View at: Google Scholar
 E. Hopf, “The partial differential equation ${U}_{x}+{UU}_{t}={U}_{xx}$,” Communications on Pure and Applied Mathematics, vol. 3, no. 3, pp. 201–230, 1950. View at: Publisher Site  Google Scholar  MathSciNet
 I. A. Hassanien, A. A. Salama, and H. A. Hosham, “Fourthorder finite difference method for solving Burgers' equation,” Applied Mathematics and Computation, vol. 170, no. 2, pp. 781–800, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Ciment, S. H. Leventhal, and B. Weinberg, “The operator compact implicit method for parabolic equations,” Journal of Computational Physics, vol. 28, no. 2, pp. 135–166, 1978. View at: Publisher Site  Google Scholar  MathSciNet
 R. S. Hirsh, “Higher order accurate difference solutions of fluid mechanics problems by a compact differencing technique,” Journal of Computational Physics, vol. 19, no. 1, pp. 90–109, 1975. View at: Publisher Site  Google Scholar
 S. Kutluay, A. R. Bahadir, and A. Özdeş, “Numerical solution of onedimensional Burgers equation: explicit and exactexplicit finite difference methods,” Journal of Computational and Applied Mathematics, vol. 103, no. 2, pp. 251–261, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 P. Arminjon and C. Beauchamp, “A finite element method for Burgers’ equation in hydrodynamics,” International Journal for Numerical Methods in Engineering, vol. 12, no. 3, pp. 415–428, 1978. View at: Publisher Site  Google Scholar
 L. Iskandar and A. Mohsen, “Some numerical experiments on the splitting of Burgers' equation,” Numerical Methods for Partial Differential Equations, vol. 8, no. 3, pp. 267–276, 1992. View at: Publisher Site  Google Scholar
 P. C. Jain and M. Raja, “Splittingup technique for Burgers’ equations,” Indian Journal of Pure and Applied Mathematics, vol. 10, pp. 1543–1551, 1979. View at: Google Scholar
 T. Öziş, E. N. Aksan, and A. Özdeş, “A finite element approach for solution of Burgers' equation,” Applied Mathematics and Computation, vol. 139, no. 23, pp. 417–428, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 A. Dogan, “A Galerkin finite element approach to Burgers' equation,” Applied Mathematics and Computation, vol. 157, no. 2, pp. 331–346, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 İ. Dağ, B. Saka, and A. Boz, “Bspline Galerkin methods for numerical solutions of the Burgers' equation,” Applied Mathematics and Computation, vol. 166, no. 3, pp. 506–522, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 A. H. Khater, R. S. Temsah, and M. M. Hassan, “A Chebyshev spectral collocation method for solving Burgers'type equations,” Journal of Computational and Applied Mathematics, vol. 222, no. 2, pp. 333–350, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. K. Khalifa, K. I. Noor, and M. A. Noor, “Some numerical methods for solving Burgers equation,” International Journal of Physical Sciences, vol. 6, no. 7, pp. 1702–1710, 2011. View at: Google Scholar
 A. Korkmaz and I. Daǧ, “Shock wave simulations using sinc differential quadrature method,” Engineering Computations, vol. 28, no. 6, pp. 654–674, 2011. View at: Publisher Site  Google Scholar
 K. Altiparmak and T. Özis, “Numerical solution of Burgers' equation with factorized diagonal Padé approximation,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 21, no. 34, pp. 310–319, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 İ. Dağ, D. Irk, and A. Şahin, “Bspline collocation methods for numerical solutions of the Burgers' equation,” Mathematical Problems in Engineering, vol. 2005, no. 5, pp. 521–538, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 S. S. Xie, S. Heo, S. Kim, G. Woo, and S. Yi, “Numerical solution of onedimensional Burgers' equation using reproducing kernel function,” Journal of Computational and Applied Mathematics, vol. 214, no. 2, pp. 417–434, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 P. Sablonniere, “Univariate spline quasiinterpolants and applications to numerical analysis,” Rendiconti del Seminario Matematico, vol. 63, no. 3, pp. 211–222, 2005. View at: Google Scholar
 Y. Ren and X. Li, “A meshfree method for Signorini problems using boundary integral equations,” Mathematical Problems in Engineering, vol. 2014, Article ID 490127, 12 pages, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 R. J. Cheng, L. W. Zhang, and K. M. Liew, “Modeling of biological population problems using the elementfree kpRitz method,” Applied Mathematics and Computation, vol. 227, pp. 274–290, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 L. W. Zhang, Y. J. Deng, and K. M. Liew, “An improved elementfree Galerkin method for numerical modeling of the biological population problems,” Engineering Analysis with Boundary Elements, vol. 40, pp. 181–188, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 L. W. Zhang, P. Zhu, and K. M. Liew, “Thermal buckling of functionally graded plates using a local Kriging meshless method,” Composite Structures, vol. 108, pp. 472–492, 2014. View at: Publisher Site  Google Scholar
 Z. X. Lei, L. W. Zhang, K. M. Liew, and J. L. Yu, “Dynamic stability analysis of carbon n anotubereinforced functionally graded cylindrical panels using the elementfree KpRitz method,” Composite Structures, vol. 113, pp. 328–338, 2014. View at: Publisher Site  Google Scholar
 K. M. Liew, Z. X. Lei, J. L. Yu, and L. W. Zhang, “Postbuckling of carbon nanotubereinforced functionally graded cylindrical panels under axial compression using a meshless approach,” Computer Methods in Applied Mechanics and Engineering, vol. 268, pp. 1–17, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 T. Belytschko, Y. Y. Lu, and L. Gu, “Elementfree Galerkin methods,” International Journal for Numerical Methods in Engineering, vol. 37, no. 2, pp. 229–256, 1994. View at: Publisher Site  Google Scholar  MathSciNet
 J. J. Monaghan, “An introduction to SPH,” Computer Physics Communications, vol. 48, no. 1, pp. 89–96, 1988. View at: Publisher Site  Google Scholar
 K. M. Liew, X. Zhao, and T. Y. Ng, “The elementfree KpRitz method for vibration of laminated rotating cylindrical panels,” International Journal of Structural Stability and Dynamics, vol. 2, no. 4, pp. 523–558, 2002. View at: Publisher Site  Google Scholar
 K. M. Liew, H. Y. Wu, and T. Y. Ng, “Meshless method for modeling of human proximal femur: treatment of nonconvex boundaries and stress analysis,” Computational Mechanics, vol. 28, no. 5, pp. 390–400, 2002. View at: Publisher Site  Google Scholar
 K. M. Liew, Y. C. Wu, G. P. Zou, and T. Y. Ng, “Elastoplasticity revisited: Numerical analysis via reproducing kernel particle method and parametric quadratic programming,” International Journal for Numerical Methods in Engineering, vol. 55, no. 6, pp. 669–683, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. M. Liew, X. L. Chen, and J. N. Reddy, “Meshfree radial basis function method for buckling analysis of nonuniformly loaded arbitrarily shaped shear deformable plates,” Computer Methods in Applied Mechanics and Engineering, vol. 193, no. 3–5, pp. 205–224, 2004. View at: Publisher Site  Google Scholar
 S. N. Atluri and T. Zhu, “A new Meshless Local PetrovGalerkin (MLPG) approach in computational mechanics,” Computational Mechanics, vol. 22, no. 2, pp. 117–127, 1998. View at: Publisher Site  Google Scholar  MathSciNet
 W. K. Liu, S. Jun, and Y. F. Zhang, “Reproducing kernel particle methods,” International Journal for Numerical Methods in Fluids, vol. 20, no. 89, pp. 1081–1106, 1995. View at: Publisher Site  Google Scholar  MathSciNet
 R. Chen and Z. Wu, “Applying multiquadratic quasiinterpolation to solve Burgers' equation,” Applied Mathematics and Computation, vol. 172, no. 1, pp. 472–484, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 C. Zhu and R. Wang, “Numerical solution of Burgers' equation by cubic BSpline quasiinterpolation,” Applied Mathematics and Computation, vol. 208, no. 1, pp. 260–272, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. Asaithambi, “Numerical solution of the Burgers' equation by automatic differentiation,” Applied Mathematics and Computation, vol. 216, no. 9, pp. 2700–2708, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 B. Saka and I. Dağ, “Quartic Bspline collocation method to the numerical solutions of the Burgers' equation,” Chaos, Solitons and Fractals, vol. 32, no. 3, pp. 1125–1137, 2007. View at: Publisher Site  Google Scholar
 I. A. Hassanien, A. A. Salama, and H. A. Hosham, “Fourthorder finite difference method for solving Burgers' equation,” Applied Mathematics and Computation, vol. 170, no. 2, pp. 781–800, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 C. de Boor, A Practical Guide to Splines, Springer, New York, NY, USA, 1978. View at: MathSciNet
 W. L. Wood, “An exact solution for Burgers' equation,” Communications in Numerical Methods in Engineering, vol. 22, no. 7, pp. 797–798, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 I. Christie, D. F. Griffiths, and A. R. Mitchell, “Product approximation for nonlinear problems in the finite element method,” IMA Journal of Numerical Analysis, vol. 1, no. 3, pp. 253–266, 1981. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2014 Mingzhu Li 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.