The main aim of this work is to consider a meshfree algorithm for solving Burgers’ equation with the quartic B-spline quasi-interpolation. Quasi-interpolation 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 ill-conditioning problem resulting from using the B-spline as a global interpolant. The numerical scheme is presented, by using the derivative of the quasi-interpolation 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 one-dimensional Burgers’ equation first suggested by Bateman [6] and later treated by Burgers [1] has the form
(1)Ut+UUx-λUxx=0,
where λ>0 is the coefficient of kinematic viscosity and the subscripts x and t denote space and time derivatives. Initial and boundary conditions are
(2)U(x,0)=f(x),a≤x≤b,U(a,t)=β1,u(b,t)=β2,t≥0,
where β1, β2, and f(x) will be chosen in a later section.
Burgers’ equation is a quasi-linear parabolic partial differential equation, whose analytic solutions can be constructed from a linear partial differential equation by using Hopf-Cole 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], B-spline 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 element-free Galerkin method [31], the smooth particle hydrodynamics method [32], the element-free kp-Ritz method [33–36], the meshless local Petrov-Galerkin method [37], and the reproducing kernel particle method [38].
With the use of univariate multiquadric (MQ) quasi-interpolation, 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 B-spline quasi-interpolation 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 B-spline quasi-interpolation. Then we do not require to solve any linear system of equation so that we do not meet the question of the ill-condition 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 B-spline has been described and univariate quartic B-spline quasi-interpolants have been presented. In Section 3, we mainly propose the numerical techniques using quartic B-spline 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 B-Spline Quasi-Interpolant
For an interval I=[a,b], we introduce a set of equally-spaced knots of partition Ω={x0,x1,…,xn}. We assume that n≥5, xi=a+ih(i=0,1,…,n), x0=a, and xn=b. Let S4[π] be the space of continuously-differentiable, piecewise, quartic-degree polynomials on π. A detailed description of B-spline functions generated by subdivision regarding the B-splines basis in S4[π] can be found in [45].
The zero degree B-spline is defined as
(3)Ni,0(x)={1,x∈[xi,xi+1],0,otherwise,
and, for positive constant p, it is defined in the following recursive form:
(4)Ni,p=x-xixi+p-xiNi,p-1(x)+xi+p+1-xxi+P+1-xi+1Ni+1,p-1,p≥1.
We apply this recursion to get the quartic B-spline Ni,4(x), which is defined in S4(π) as follows:(5)Ni,4(x)=124h4{(x-xi-2)4,x∈[xi-2,xi-1],(x-xi-2)4-5(x-xi-1)4,x∈[xi-1,xi],(x-xi-2)4-5(x-xi-1)4+10(x-xi)4,x∈[xi,xi+1],(x-xi+3)4-5(x-xi+2)4,x∈[xi+1,xi+2],(x-xi+3)4,x∈[xi+2,xi+3],0,otherwise.As usual, we add multiple knots at the endpoints: a=x-4=x-3=⋯=x0 and b=xn=xn+1=⋯=xn+4.
In [24], univariate quartic B-spline quasi-interpolants (abbr.QIs) can be defined as operators of the form
(6)Q4(f)=∑j=1n+4μjNj,4.
The coefficients are listed as follows:
(7)μ1(f)=f1,μ2(f)=17105f1+3532f2-3596f3+21160f4-5244f5,μ3(f)=-1945f1+377288f2+61288f3-59480f4+7288f5,μ4(f)=47315f1-77144f2+251144f3-97240f4+471008f5,μj(f)=471152(fj-4+fj+1)-107288(fj-3+fj-1)μj(f)=+319192fj-2,j=5,…,n,μn+1(f)=47315fn+2-77144fn+1+251144fn-97240fn-1μn+1(f)=+471008fn-2,μn+2(f)=-1945fn+2+377288fn+1+61288fn-59480fn-1μn+2(f)=+7288fn-2,μn+3(f)=17105fn+2+3532fn+1-3596fn+21160fn-1μn+3(f)=-5244fn-2,μn+4(f)=fn+2,
and fi=f(ti), ti=(1/2)(xi-2+xi-1), i=1,…,n+2. For f∈C5(I), we have the error estimate
(8)∥f-Q4(f)∥∞=O(h5).
We use ∏4 to denote the space of polynomials of the total degree at most 4. In general, we impose that Q4 is exact on the space ∏4; that is, Q4(p)=p for all p∈∏4. As a consequence of this property, the approximation order of Q4 is O(h5) on smooth functions. In this paper, the coefficient μj is a linear combination of discrete values of f at some points. The main advantage of QIs is that they have a direct construction without solving any system of linear equations. Moreover, they are local in the sense that values of Q4f(x) depend only on values of f in a neighborhood of x. 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 f by derivatives of Q4f. The general theory will be developed elsewhere. We can evaluate the value of f at xi by (Q4f)′=∑j=1n+4μj(f)Nj,4′ and (Q4f)′′=∑j=1n+4μj(f)Nj,4′′. Nj,4′ and Nj,4′′ can be computed by the formula of B-spline’s derivatives as follows:
(9)Ni,4(k)=4!(4-k)!∑j=1nαk,jNi+j,4-k,
where
(10)α0,0=1,αk,0=αk-1,0xi+3-k-xi,αk,k=-αk-1,k-1xi+5-xi+k,αk,j=αk-j,j-αk-1,j-1xi+j+5-k-xi+j.
By some trivial computations, we can obtain the value of Ni,4(k)(k=0,1,2,3) at the knots, which are illustrated in Table 1. Then, we get the differential formulas for quartic B-spline QIs as
(11)f′=∑j=1n+4μj(f)Nj,4′,f′′=∑j=1n+4μj(f)Nj,4′′.
The values of Ni,4(k)(x) at the knots.
xi-1
xi
xi+1
xi+2
Otherwise
Ni,4(x)
124
1124
1124
124
0
Ni,4′(x)
16h
36h
-36h
-16h
0
Ni,4′′(x)
12h2
-12h2
-12h2
-12h2
0
Ni,4′′′(x)
1h3
-3h3
3h3
-1h3
0
3. Numerical scheme Using the Meshfree Quasi-Interpolation
In this section, we present the numerical scheme for solving Burgers’ equation based on the quartic B-spline quasi-interpolation.
Discretizing the Burgers’ equation
(12)Ut+UUx-λUxx=0,
in time with meshlength τ, we get
(13)Ujk+1-Ujkτ+Ujk(Ux)jk-λ(Uxx)jk=0.
We can get
(14)Ujk+1=Ujk+τUjk(Ux)jk-τλ(Uxx)jk,
where Ujk is the approximation of the value of U(x,t) at the point (xj,tk). Then, we can use the derivatives of the quartic B-spline quasi-interpolant Q4U(xj,tk) to approximate (Ux)jk and (Uxx)jk. To dump the dispersion of the scheme, we define a switch function g(x,t), whose values are 0 and 1 at the discrete points (xj,tk), as follows:
(15)g(xj,tk)=max{0,1+min{0,sign((Ux)jk·(Ux)lk)}},
where l=j-sign(Ujk). Thus, the resulting numerical scheme is
(16)Ujk+1=Ujk+τUjk(Ux)jkg(xj,tk)-τλ(Uxx)jk.
Starting from the initial condition, we can compute the numerical solution of Burgers’ equation step by step using the B-spline quasi-interpolation scheme (16) and formulas (11).
4. Numerical Results
To investigate the applicability of the quasi-interpolation 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:
(17)L∞=maxj|Ujexact-Ujnum|,L2=h∑j=1n(Ujexact-Ujnum)2.
Example 1.
Burgers’ equation is solved over the region [0,1] and the initial and boundary conditions are given in Asaithambi [42]:
(18)U(x,0)=2λπsinπxα+cosπx(α>1),U(0,t)=0,U(1,t)=0,t>0,
and the exact solution of this problem has the following nice compact closed-form, as given by Wood [46]:
(19)U(x,t)=2λπe-π2λtsinπxα+e-π2λtcosπx(α>1).
In this computational study, we set α=2, h=0.025, Δτ=0.0001. 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.
Comparison of exact and numerical solution at t=0.001.
x
λ=1
λ=0.5
Our method
Asai [42]
Exact
Our method
Asai [42]
Exact
0.1
0.653563
0.653589
0.653544
0.327870
0.327874
0.327870
0.2
1.305519
1.305611
1.305534
0.655028
0.655078
0.655069
0.3
1.949321
1.949485
1.949364
0.978449
0.978427
0.978413
0.4
2.565977
2.566103
2.565925
1.288417
1.288485
1.288463
0.5
3.110769
3.110992
3.110739
1.563014
1.563096
1.563064
0.6
3.492902
3.493222
3.492866
1.756653
1.756691
1.756642
0.7
3.549538
3.550079
3.549595
1.787184
1.787281
1.787206
0.8
3.050089
3.050702
3.050134
1.537658
1.537794
1.537694
0.9
1.816492
1.817077
1.816660
0.916795
0.916941
0.916860
Example 2.
In this example, we consider the exact solution of Burgers’ equation [47]:
(20)U(x,t)=α+μ+(μ-α)expη1+expη,0≤x≤1,t≥0,
where η=(α(x-μt-γ))/λ, α, μ, and γ are constants. The boundary conditions are
(21)U(0,t)=1,U(1,t)=0.2,t≥0,
and initial condition is used for the exact solution at t=0.
We solve the problem with α=0.4, μ=0.6, and γ=0.125 by our method. In Table 3, L2 and L∞ errors at the time level t=0.5 are compared with the error obtained by Chen and Wu [39], Zhu and Wang [41], Dag˘ et al. [17], and Saka and Dag˘ [43]. For comparison, the parameters are adopted as time step τ=0.01, space step h=1/36, and viscosity coefficient λ=0.01. From Table 3, we can find that our method provides better accuracy than most methods through the L2 and L∞ error norms. The profiles of initial wave and its propagation are depicted at some times in Figure 1.
The computational results at t=0.5 for λ=0.01 with h=1/36 and τ=0.01.
BSQI [41]
MQQI [39]
QBCM I [17]
QBGM I [43]
Our method
L2×103
3.43253
5.77786
0.77033
1.92558
0.85269
L∞×103
9.26698
20.8467
3.03817
6.35489
3.79716
The numerical solutions with h=1/36, τ=0.01 for λ=0.01.
Example 3.
Consider Burgers’ equation with the initial condition
(22)U(x,0)=sin(πx),0≤x≤1
and the boundary conditions
(23)U(0,t)=U(1,t)=0.
The analytical solution of this problem was given by Cole [2] in the term of an infinite series as
(24)U(x,t)=2πλ∑k=1∞kAksin(kπx)exp(-k2π2λt)A0+∑k=1∞Akcos(kπx)exp(-k2π2λt)
with the Fourier coefficients
(25)A0=∫01exp{-(2πλ)-1(1-cos(πx))}dx,Ak=2∫01exp{-(2πλ)-1(1-cos(πx))}cos(kπx)dx,k≥1.
In Table 4, we have computed the numerical solutions of this example at differential time levels with parameter values λ=0.1, h=0.025, and τ=0.0001. 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 t=1 with τ=0.001, h=0.01 for λ=0.1,0.01,0.0001, respectively. For the MQQI method, the shape parameter c=7.2×10-3, 2.9×10-3, 1.43×10-4 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.
Comparison of results at different time for λ=0.1 with h=0.025 and τ=0.0001.
x
t
Hassanien [44]
Kutluay [11]
Ozis [15]
Our method
Exact
0.25
0.4
0.3089
0.3083
0.3143
0.3089
0.3089
0.6
0.2407
0.2404
0.2437
0.2407
0.2407
0.8
0.1957
0.1954
0.1976
0.1957
0.1957
1.0
0.1626
0.1624
0.1639
0.1626
0.1626
3.0
0.0272
0.0272
0.0274
0.0272
0.0272
0.50
0.4
0.5696
0.5691
0.5764
0.5696
0.5696
0.6
0.4472
0.4468
0.4517
0.4472
0.4472
0.8
0.3592
0.3589
0.3625
0.3592
0.3592
1.0
0.2919
0.2916
0.2944
0.2919
0.2919
3.0
0.0402
0.0402
0.0406
0.0402
0.0402
0.75
0.4
0.6254
0.6256
0.6259
0.6257
0.6254
0.6
0.4872
0.4570
0.4903
0.4872
0.4872
0.8
0.3739
0.3737
0.3771
0.3739
0.3739
1.0
0.2875
0.2872
0.2902
0.2875
0.2875
3.0
0.0298
0.0297
0.0133
0.0298
0.0298
Comparison of results at t=1 for λ=0.1.
x
Hon and Mao [40]
MQQI [39]
BSQI [41]
Our method
Exact
0.1
0.0664
0.07124
0.06628
0.06630
0.06632
0.2
0.1313
0.13431
0.13115
0.13119
0.13121
0.3
0.1928
0.19339
0.19269
0.19271
0.19279
0.4
0.2481
0.24538
0.24792
0.24797
0.24803
0.5
0.2919
0.28517
0.29175
0.29185
0.29191
0.6
0.3159
0.30473
0.31580
0.31598
0.31607
0.7
0.3079
0.29288
0.30791
0.30800
0.30810
0.8
0.2534
0.23784
0.25337
0.25344
0.25372
0.9
0.1459
0.13542
0.14583
0.14587
0.14606
Comparison of results at t=1 for λ=0.01.
x
Hon and Mao [40]
MQQI [39]
BSQI [41]
Our method
Exact
0.1
0.0755
0.07868
0.07530
0.07538
0.0754
0.2
0.1507
0.15202
0.15049
0.15066
0.1506
0.3
0.2257
0.22554
0.22544
0.22573
0.2257
0.4
0.3003
0.29904
0.30002
0.30028
0.3003
0.5
0.3744
0.37226
0.37407
0.37437
0.3744
0.6
0.4478
0.44484
0.44742
0.44778
0.4478
0.7
0.5202
0.51643
0.51985
0.52038
0.5203
0.8
0.5913
0.58622
0.59106
0.59151
0.5915
0.9
0.6607
0.62956
0.65964
0.66007
0.6600
Comparison of results at t=1 for λ=0.0001.
x
Hon and Mao [40]
MQQI [39]
BSQI [41]
Our method
0.05
0.0422
0.0422
0.0422
0.0422
0.16
0.1263
0.1263
0.1262
0.1262
0.27
0.2103
0.2103
0.2096
0.2103
0.38
0.2939
0.2939
0.2928
0.2939
0.50
0.3769
0.3769
0.3754
0.3769
0.61
0.4592
0.4592
0.4573
0.4592
0.72
0.5404
0.5404
0.5381
0.5404
0.83
0.6203
0.6201
0.6174
0.6203
0.94
0.6983
0.6957
0.6947
0.6983
Example 4.
We consider particular solution of Burgers’ equation:
(26)U(x,t)=x/t1+t/t0exp(x2/4λt),t≥1,0≤x≤1,
where t0=exp(1/8λ). Initial condition is obtained from when t=1 is used. Boundary conditions are U(0,t)=U(1.2,t)=0. Analytical solution represents shock-like solution of the one-dimensional Burgers’ equation. Parameters h=0.02,0.005 and λ=0.005,0.01 are selected for comparison over the domain [0,1]. 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 h=0.02, τ=0.001, and λ=0.005 for t≤4 in Figure 2.
Comparison of results at different times for τ=0.01.
L2×103
L∞×103
L2×103
L∞×103
L2×103
L∞×103
h=0.02, λ=0.005
t=1.8
t=1.8
t=2.4
t=2.4
t=3.2
t=3.2
BSQI [41]
1.66464
5.12020
2.06695
6.31491
2.36889
6.85425
QBCM I [17]
0.19127
0.54058
0.14246
0.39241
0.93617
5.54899
QBCM II [17]
0.49130
1.16930
0.41864
0.93664
1.28863
7.49147
MQQI [39]
6.88480
25.6767
7.89738
27.2424
8.56856
2.68122
Our method
0.6642
0.91725
0.7573
1.1465
0.8592
1.2103
h=0.02, λ=0.01
t=1.8
t=1.8
t=2.4
t=2.4
t=3.2
t=3.2
BSQI [41]
0.82751
2.59444
0.98595
2.35031
1.58264
5.73827
QBCM I [17]
0.17014
0.40431
0.20476
0.86363
1.29951
6.69425
QBCM II [17]
0.24003
0.48800
0.30849
1.14760
1.57548
8.06798
MQQI [39]
5.89555
14.7550
6.64358
15.9892
6.90385
16.3403
Our method
0.82751
0.50367
0.46281
1.05625
0.88261
4.73827
The numerical solutions for t≤4.
5. Conclusion
Following the recent development of the quasi-interpolation method for scattered data interpolation and the meshfree method for solving partial differential equations, this paper combines these ideas and proposes a new meshfree quasi-interpolation method for Burgers’ equation. The method does not require solving a large size matrix equation and, hence, the ill-conditioning problem from using B-spline 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 high-dimensional Burgers’ equations, we believe our scheme can also be applicable. In this case, we would use multivariate spline quasi-interpolation instead of univariate spline quasi-interpolation. 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).
BurgersJ. M.A mathematical model illustrating the theory of turbulence19481171199ColeJ. D.On a quasi-linear parabolic equation occurring in aerodynamics195193225236MR0042889LighthillM. J.Viscosity effects in sound waves of finite amplitude1956250351MR0077346PospelovL. A.Propagation of finite amplitude elastic waves(Longitudinal elastic wave of finite amplitude propagation in isotropic solid)1966113023042-s2.0-11644278702van der PolB.On a non-linear partial differential equation satisfied by the logarithm of the Jacobian theta-functions, with arithmetical applications. I, II195113261284MR0042599BatemanH.Some recent researches on the motion of fluids1915434163170HopfE.The partial differential equation Ux+UUt=Uxx19503320123010.1002/cpa.3160030302MR0047234HassanienI. A.SalamaA. A.HoshamH. A.Fourth-order finite difference method for solving Burgers' equation2005170278180010.1016/j.amc.2004.12.052MR2175245ZBL1084.650782-s2.0-27144455801CimentM.LeventhalS. H.WeinbergB.The operator compact implicit method for parabolic equations197828213516610.1016/0021-9991(78)90031-1MR5055882-s2.0-0008466688HirshR. S.Higher order accurate difference solutions of fluid mechanics problems by a compact differencing technique19751919010910.1016/0021-9991(75)90118-7KutluayS.BahadirA. R.ÖzdeşA.Numerical solution of one-dimensional Burgers equation: explicit and exact-explicit finite difference methods1999103225126110.1016/S0377-0427(98)00261-1MR16776532-s2.0-0039728591ArminjonP.BeauchampC.A finite element method for Burgers’ equation in hydrodynamics197812341542810.1002/nme.16201203042-s2.0-0017919256IskandarL.MohsenA.Some numerical experiments on the splitting of Burgers' equation19928326727610.1002/num.16900803032-s2.0-0026866334JainP. C.RajaM.Splitting-up technique for Burgers’ equations19791015431551Zbl0432.35071ÖzişT.AksanE. N.ÖzdeşA.A finite element approach for solution of Burgers' equation20031392-341742810.1016/S0096-3003(02)00204-7MR19486512-s2.0-0037448173DoganA.A Galerkin finite element approach to Burgers' equation2004157233134610.1016/j.amc.2003.08.037MR2088257ZBL1054.651032-s2.0-4444378330Dağİ.SakaB.BozA.B-spline Galerkin methods for numerical solutions of the Burgers' equation2005166350652210.1016/j.amc.2004.06.078MR21504862-s2.0-20444477122KhaterA. H.TemsahR. S.HassanM. M.A Chebyshev spectral collocation method for solving Burgers'-type equations2008222233335010.1016/j.cam.2007.11.007MR2474624ZBL1153.651022-s2.0-53449090789KhalifaA. K.NoorK. I.NoorM. A.Some numerical methods for solving Burgers equation201167170217102-s2.0-79960752482KorkmazA.DaǧI.Shock wave simulations using sinc differential quadrature method201128665467410.1108/026444011111546192-s2.0-80052213097AltiparmakK.ÖzisT.Numerical solution of Burgers' equation with factorized diagonal Padé approximation2011213-431031910.1108/09615531111108486MR28478192-s2.0-79955703532Dağİ.IrkD.ŞahinA.B-spline collocation methods for numerical solutions of the Burgers' equation20052005552153810.1155/MPE.2005.521MR21946512-s2.0-31044432795XieS. S.HeoS.KimS.WooG.YiS.Numerical solution of one-dimensional Burgers' equation using reproducing kernel function2008214241743410.1016/j.cam.2007.03.010MR23983432-s2.0-39049173059SablonniereP.Univariate spline quasi-interpolants and applications to numerical analysis2005633211222RenY.LiX.A meshfree method for Signorini problems using boundary integral equations201420141249012710.1155/2014/490127MR3173374ChengR. J.ZhangL. W.LiewK. M.Modeling of biological population problems using the element-free kp-Ritz method201422727429010.1016/j.amc.2013.11.033MR3146315ZhangL. W.DengY. J.LiewK. M.An improved element-free Galerkin method for numerical modeling of the biological population problems20144018118810.1016/j.enganabound.2013.12.008MR3161279ZhangL. W.ZhuP.LiewK. M.Thermal buckling of functionally graded plates using a local Kriging meshless method201410847249210.1016/j.compstruct.2013.09.0432-s2.0-80052384467LeiZ. X.ZhangL. W.LiewK. M.YuJ. L.Dynamic stability analysis of carbon n anotube-reinforced functionally graded cylindrical panels using the element-free Kp-Ritz method201411332833810.1016/j.compstruct.2014.03.035LiewK. M.LeiZ. X.YuJ. L.ZhangL. W.Postbuckling of carbon nanotube-reinforced functionally graded cylindrical panels under axial compression using a meshless approach201426811710.1016/j.cma.2013.09.001MR3133485BelytschkoT.LuY. Y.GuL.Element-free Galerkin methods199437222925610.1002/nme.1620370205MR12568182-s2.0-0028259955MonaghanJ. J.An introduction to SPH1988481899610.1016/0010-4655(88)90026-42-s2.0-0023165466LiewK. M.ZhaoX.NgT. Y.The element-free Kp-Ritz method for vibration of laminated rotating cylindrical panels20022452355810.1142/S0219455402000701LiewK. M.WuH. Y.NgT. Y.Meshless method for modeling of human proximal femur: treatment of nonconvex boundaries and stress analysis200228539040010.1007/s00466-002-0303-52-s2.0-0036577031LiewK. M.WuY. C.ZouG. P.NgT. Y.Elasto-plasticity revisited: Numerical analysis via reproducing kernel particle method and parametric quadratic programming200255666968310.1002/nme.523ZBL1033.740502-s2.0-0037202023LiewK. M.ChenX. L.ReddyJ. N.Mesh-free radial basis function method for buckling analysis of non-uniformly loaded arbitrarily shaped shear deformable plates20041933–520522410.1016/j.cma.2003.10.0022-s2.0-0347300455AtluriS. N.ZhuT.A new Meshless Local Petrov-Galerkin (MLPG) approach in computational mechanics199822211712710.1007/s004660050346MR16494202-s2.0-0032136132LiuW. K.JunS.ZhangY. F.Reproducing kernel particle methods1995208-91081110610.1002/fld.1650200824MR13339172-s2.0-0029102512ChenR.WuZ.Applying multiquadratic quasi-interpolation to solve Burgers' equation2006172147248410.1016/j.amc.2005.02.027MR21979162-s2.0-31144472588HonY. C.MaoX. Z.An efficient numerical scheme for Burgers' equation1998951375010.1016/S0096-3003(97)10060-1MR16302602-s2.0-0002879831ZhuC.WangR.Numerical solution of Burgers' equation by cubic B-Spline quasi-interpolation2009208126027210.1016/j.amc.2008.11.045MR2490794ZBL1159.650872-s2.0-58549091278AsaithambiA.Numerical solution of the Burgers' equation by automatic differentiation201021692700270810.1016/j.amc.2010.03.115MR2653083ZBL1193.651542-s2.0-77953133901SakaB.DağI.Quartic B-spline collocation method to the numerical solutions of the Burgers' equation20073231125113710.1016/j.chaos.2005.11.0372-s2.0-33845310038HassanienI. A.SalamaA. A.HoshamH. A.Fourth-order finite difference method for solving Burgers' equation2005170278180010.1016/j.amc.2004.12.052MR21752452-s2.0-27144455801de BoorC.1978New York, NY, USASpringerMR507062WoodW. L.An exact solution for Burgers' equation200622779779810.1002/cnm.850MR22449572-s2.0-33745845521ChristieI.GriffithsD. F.MitchellA. R.Product approximation for nonlinear problems in the finite element method19811325326610.1093/imanum/1.3.253MR6413092-s2.0-0003991061