بهبود جواب معادله غیرخطی تعمیم یافته بنجامین بااستفاده از یک روش بدون شبکه توابع پایه ای شعاعی
Subject Areas : Numerical Analysis
مهران
نعمتی
1
(گروه ریاضی ، واحد رودبار، دانشگاه آزاد اسلامی، رودبار، ایران)
سیده فائزه
تیموری
2
(گروه کامپیوتر، واحد رودبار، دانشگاه آزاد اسلامی، رودبار، ایران)
Keywords: تفاضل متناهی, توابع پایه ای شعاعی, روش کرانک نیکلسون, معادله غیر خطی تعمیم یافته بنجامین,
Abstract :
در این مقاله براساس روش توابع پایه ای شعاعی ، معادله دیفرانسیل غیرخطی تعمیم یافته بنجامین حل شده است.برای گسسته سازی قسمت زمانی از تفاضل متناهی و کرانک نیکلسون استفاده کردیم. و قسمت فضایی با استفاده از درونیابی توابع پایه ای شعاعی تقریب زده شده است. در نتیجه یک دستگاه معادلات جبری خطی حاصل می شود که باحل این دستگاه جواب های تقریبی بدست می ایند.در ادامه باحل مثال عددی نشان داده می شود که روش پیشنهادی کارا هست واینکه خطا بهبود یافته است.
Improved solution to nonlinear generalized Benjamin–Bona–Mahony–Burgers (GBBMB) equation by a meshless RBFs method
Abstract
In this paper, based on the RBF collocation method and finite differences, a numerical method is proposed to solve nonlinear generalized Benjamin–Bona–Mahony–Burgers (GBBMB) equation. First order finite differences and Crank-Nicolson method are applied to discretize the temporal parts. The spatial parts are approximated by MQ-RBF interpolation which results in a linear system of algebraic equations. Approximate solutions are determined by solving such a system. The proposed scheme is verified by solving some test problems and computing error norms and
. Results show the efficiency of the suggested method and the error has been improved.
Key words: Radial basis functions (RBFs), Finite differences, Crank-Nicolson scheme, Nonlinear Generalized Benjamin–Bona–Mahony–Burgers (GBBMB) equation.
1. Introduction
Most phenomena in applied sciences and engineering are modeled by nonlinear partial differential equations (PDEs). Since their exact solutions cannot be found simply, studying those types of equations is often hard. Hence, numerical methods for approximating nonlinear PDEs have been regarded by many researchers and they have been successfully applied to a lot of real-world problems (e.g., [1-5]). One well-known class of numerical methods is meshless methods. These schemes have attracted great attention in recent decades and they have been applied as powerful tools especially for problems in computational mechanics. The advantage of meshless methods over the traditional numerical techniques including finite difference method (FDM), finite element method (FEM), and finite volume method (FVM), is that they do not require mesh generation and domain or surface discretization.
There exist several types of meshless methods including moving least square meshless method [6], meshless local Petrov- Galerkin method (MLPG) [7], smooth-particle hydrodynamics [8], the reproduced kernel particle method (RKPM) [9], the finite point method [10], mesh-free weak-strong form (MWS) [11], the diffuse element method (DEM) [12], and radial basis functions (RBFs) method. Each of these approaches has specific advantages for certain classes of problems. Among them, RBFs method is the simplest and the most efficient one.
RBFs method was introduced by Ronald Hardy, an Iowa State geodesist, for the first time in a paper appearing in 1971 [13]. He proposed this method for efficient interpolation of the scattered data on a topographic surface. A set of N distinct points called centers is used in RBFs interpolation. There are no constraints on the geometry of domains of the problem or on the position of the centers [14]. There are different types of RBFs such as Multiquadric (MQ), Thin Plate Spline (TPS), Gaussian, Linear, Inverse quadric (IQ), Inverse Multiquadric (IMQ). Hardy used the MQ function for his interpolation scheme. Later, Duchon proposed TPS for data interpolation in 1975 [15]. The use of radial basis functions was limited to scattered data interpolation until 1990, a physicist named Edward Kansa used RBFs to solve PDEs for the first time [16,17]. His method leads to an ill-conditioned matrix for a large number of nodes because of having asymmetrical nature of the interpolation matrix. Fasshauer modified the Kansa's method in 1996 [18] by proposing a Hermite based approach in which the collocation matrices are symmetric in nature, and their condition number is smaller. In the continuation, many efforts have done by numerous researchers to improve the method and also propose new kinds of it (e.g., [19-24]).
During recent years, RBFs method has been considered as an efficient tool for solving different kinds of problems including PDEs [25-28], integral equations [29,30], and fractional equations [31,32]. The current work is devoted to the numerical study of the following two-dimensional nonlinear PDE by the RBFs method.
(1)
with following initial and boundary conditions.
(2)
where ,
is the vector valued function,
,
,
are gradient and Laplacian operators respectively. This equation is known as the nonlinear generalized Benjamin–Bona–Mahony–Burgers (GBBMB) equation.
The GBBMB equation is used in different scientific fields such as the analysis of the surface waves of long wavelength in liquids, hydromagnetic waves in a cold plasma, acoustic gravity waves incompressible fluids, and acoustic waves in a harmonic crystal. The GBBMB equation and its modified or generalized form have been investigated by many authors both numerically and analytically. For further information refer to [33-44]. The nonlinear GBBMB has been solved using different analytical and numerical methods that in the following, some of them are reported. Tari and Ganji [45] approximated solutions of nonlinear GBBMB equation by He’s methods analytically. In [46], an exponential function method was presented to solve especial type of nonlinear GBBMB equation. Dehghan et al. studied nonlinear GBBMB by an RBF meshless method in [47]. Haq et al. proposed a numerical method based on Haar wavelets and finite differences to solve nonlinear GBBMB in [54]. Hajiketabi et al. introduced a new numerical method for solving the nonlinear GBBMB equation in [55]. Recently, Ali Ebrahimijahan and Mehdi Dehghan [48] have proposed a numerical technique based on the integrated radial basis functions (IRBFs) for solving the nonlinear GBBMB and regularized long-wave equations. The present study aims to apply a numerical method based on the RBF collocation method and finite differences to solve the nonlinear GBBMB equation. Finite differences and Crank-Nicolson scheme are employed for temporal part discretization, while two-dimensional RBF interpolation is implemented for approximating the spatial parts. The MQ-RBF is the focus of this paper because of its popularity in applications and its good approximation characteristic.
The manuscript is organized as following. In section 2, the RBF method is explained briefly by presenting basic concepts and definitions. In Section 3, RBF collocation method is implemented for the time discretized nonlinear GBBMB equation. Then, the proposed method is applied for some test problems, and consequently the results are reported in Section 4. Finally, a conclusion is presented in Section 5.
2. A brief review of RBFs method
In this section, some basic concepts and definitions are expressed for the radial basis functions interpolation.
Definition 1. Let be d-dimensional Euclidean space and
. A radial basis function is a function which is both continuous and multivariable like
that its value at any point
, is dependent on the distance from a certain point
. This function could be written as
where
and
is the Euclidean norm on
. The function
is an univariable function in r and
is a center of RBF
.
Definition 2. Given the data , with
,
, and
, the scattered data interpolation problem is defined as finding a smooth function s such that
, for
Function s is called an interpolant.
A radial basis function interpolant u at centers assumes the following form.
(3)
where ,
is a radial basis function, coefficients
are constants to be determined such that the following interpolation condition at the set of N centers,
is hold.
(4)
Imposing the interpolation condition (4) to (3) leads to a linear system as follows.
(5)
where ,
, and B is a
matrix called the interpolation matrix or the system matrix defined as follows.
where
According to Definition 1, RBF is independent of the spatial dimension. This property makes it possible to transform a multivariable problem into a one-variable problem easily. This is the main advantage of the RBF interpolation scheme over the other classical methods.
Generally, RBFs are sorted into two major categories: Infinitely smooth and piecewise smooth. Some well-known RBFs are listed in Table 1. Infinitely smooth RBFs include a free parameter which is called shape parameter, denoted by . Although this parameter can be chosen arbitrarily, a proper choice of its value is necessary. Because in an infinitely smooth RBF interpolation, the value of a shape parameter influences the accuracy of the scheme [49]. Piecewise smooth RBFs have algebraic convergence rates while infinitely smooth RBFs achieve spectral or exponential convergence rates [50, 51].
Table 1. Some well-known RBFs
Category | Name of the function | Definition |
| Multiquadric (MQ) |
|
| Inverse Multiquadric (IMQ) |
|
Infinitely smooth RBFs | Inverse Quadric (IQ) |
|
| Gaussian |
|
Piecewise smooth RBFs | Linear | r |
Cubic |
| |
Thin Plate Spline (TPS) |
|
3. Time discretization of nonlinear GBBM equation
In this section, Let be an arbitrary interval in
. by applying forward finite differences and also Crank-Nicolson scheme, the time variable is discretized for the first-order time derivative as follows.
(6)
Now, choose . Then,
Approximated non-linear terms by the following formulas
(7)
yields to
This equation can be simplified as follows.
(8)
Discretizing Eq. (8) in space by RBF expansion (3) results in
Considering N collocation points in
leads to the following linear system.
Where ,
,
,
and
.
and
are matrices of second derivative of the system matrix, B, respectively in
, and
.
and
are matrices of first derivative of the system matrix, B, respectively in
, and
.
Let
and
By the assumption is non-singular,
is obtained as follows.
Recalling that , the approximate PDE solution at
is obtained as follows.
(9)
Where and
.
4. Numerical experiments
In this section, the following some test problems are numerically solved for the purpose of verifying the ability of the proposed method with regards to the nonlinear GBBMB equations. Among all of the RBFs, MQ, the most popular RBF, is used in computations due to the rapid convergent rate. Here, following MQ radial basis function is used.
Where is the shape parameter.
Meshless methods which are based on radial basis functions (RBFs) contain a free shape parameter that plays an important role for the accuracy and condition number of the coefficient matrix of the method. Most authors use the trial and error method for obtaining a good shape parameter that results in best accuracy. Here, the shape parameter is chosen by trial and error method.
The domain is chosen as the unit region, i.e.
. In order to test the accuracy, two error norms,
and
defined as follows are computed.
where and
denote the approximate and exact solutions, respectively.
4.1. Test Problem 1
Consider the nonlinear GBBMB equation as follows.
with conditions
and the source term
The exact solution of the problem is .
Numerical solutions are calculated at various values of time variable with
and time step
and shape parameter= 1.4. Consequently, the results are provided in Table 2, showing that the proposed method is accurate sufficiently. The numerical results at
,
and different time steps
and the numerical results presented in [54], are all provided in Table 3. Values of shape parameter are derived by trial and error method. Comparison the results show that better approximations are obtained by the proposed scheme.
Approximate and exact solutions and also absolute error are illustrated in Fig.1,. According to this figure, it can be seen that approximate solutions are very close to the exact ones.
Table 2. Error norms of test problem 1 at and
T |
|
|
0.1 | 6.7949e-06 | 3.4022e-05 |
0.2 | 1.2265e-05 | 6.1129e-05 |
0.5 | 2.1897e-05 | 1.0869e-04 |
1 | 2.9094e-05 | 1.2306e-04 |
Fig 1. Graphs of approximate and exact solutions and absolute error of test problem 1 at and
Table 3. Error norms of test problem 1 at
|
|
|
|
| CPU time | Shape parameter |
1/10 | 1.9000e-3 | 3.4000e-3 | 4.1369e - 03 | 1.0243e - 02 | 0.1550 | 1.5 |
1/20 | 9.9446e-04 | 1.8000e-3 | 2.4115e - 03 | 5.8698e - 03 | 0.1718 | 1.1 |
1/40 | 5.7353e-04 | 9.8716e-04 | 1.5207e - 03 | 3.6346e - 03 | 0.1747 | 0.8 |
1/80 | 3.1700e-04 | 5.5014e-04 | 1.0682e - 03 | 2.5188 - 03 | 0.1958 | 0.6 |
1/160 | 1.2022e-04 | 1.7277e-04 | 8.4022e - 04 | 1.9708 - 03 | 0.2301 | 0.6 |
1/320 | 5.0058e-05 | 8.2924e-05 | 7.2574e - 04 | 1.7034 - 03 | 0.3083 | 0.5 |
1/640 | 2.0853e-05 | 3.9993e-05 | 6.7980e - 04 | 1.5724 - 03 | 0.47338 | 0.4 |
1/1280 | 3.2936e-05 | 5.0085e-05 | 6.6576e - 04 | 1.5079 - 03 | 0.7473 | 0.3 |
1/2560 | 2.4328e-05 | 4.8605e-05 | 6.5874e - 04 | 1.4759 - 03 | 1.3567 | 0.3 |
4.2. Test Problem 2
Consider the nonlinear GBBMB equation as follows.
With initial and boundary conditions
and the source term
with the exact solution .
In table 4, Numerical solutions are calculated at various values of time variable with
and time step
and shape parameter= 1.4. The results are provided in Table 2 show that the proposed method is accurate sufficiently.
Table 4. Error norms of test problem 2 at and
T |
|
|
0.1 | 3.7096e-06 | 1.5701e-05 |
0.2 | 7.4830e-05 | 2.8961e-05 |
0.5 | 1.9129e-05 | 6.5406e-05 |
1 | 3.9390e-05 | 1.5772e-04 |
5. Conclusion
In this study, an PDE called nonlinear generalized Benjamin–Bona–Mahony–Burgers (GBBMB) equation was studied numerically. The finite difference formula and Crank Nicolson technique were implemented to discretized the temporal parts. As a result, a time semi-discrete formula was obtained. After that, a fully discrete formula was achieved by approximating the spatial terms using RBF interpolation. Numerical results show that the suggested method has better accuracy and the error has been improved.
References
[1] H.M. Patil, R. Maniyeri, Finite difference method based analysis of bio-heat transfer in human breast cyst, Thermal Science and Engineering Progress (2019).
[2] Jiang, W., Spencer, B.W., Dolbow, J.E., Ceramic Nuclear Fuel Fracture Modeling with the Extended Finite Element Method, Engineering Fracture Mechanics (2019), doi: https://doi.org/10.1016/ j.engfracmech.2019.106713
[3] E. Rossi, M. Paloni, G. Storti, R. Rota, Modeling Dual Reflux-Pressure Swing Adsorption processes: numerical solution based on the Finite Volume Method, Chemical Engineering Science (2019), doi: https://doi.org/10.1016/j.ces.2019.03.055
[4] Longfei Gao, David Keyes, Combining finite element and finite difference methods for isotropic elastic wave simulations in an energy-conserving manner, Journal of Computational Physics 378 (2019) 665–685
[5] Jithin Jose, Sung-Jin Choi, Knut Erik Teigen Giljarhus, Ove Tobias Gudmestad, A comparison of numerical simulations of breaking wave forces on a monopile structure using two different numerical models based on finite difference and finite volume methods, Ocean Engineering 137 (2017) 78–88.
[6] E.Dabboura, H.Sadat, C.Prax, A moving least squares meshless method for solving the generalized Kuramoto-Sivashinsky equation, Alexandria Engineering Journal, Volume 55, Issue 3, September 2016, Pages 2783-2787.
[7] S. N. Atluri and T. -L. Zhu, The meshless local Petrov-Galerkin (MLPG) approach for solving problems in elasto-statics, Computational Mechanics volume 25, pages169–179(2000).
[8] Wang Zhi-Bin, ,Chen Rong, WangHong, Liao Qiang,,ZhuXun,, LiShu-Zhe, An overview of smoothed particle hydrodynamics for simulating multiphase flow, Applied Mathematical Modelling, Volume 40, Issues 23–24, December 2016, Pages 9625-9655.
[9] W. K. Liu, S. Jun, and Y. F. Zhang, “Reproducing kernel particle methods,” International Journal for Numerical Methods in Fluids, vol. 20, no. 8-9, pp. 1081–1106, 1995.
[10] E. Onate, S. Idelsohn, O. C. Zienkiewicz, and R. L. Taylor, “A finite point method in computational mechanics. Applications to convective transport and fluid flow,” International Journal for Numerical Methods in Engineering, vol. 39, no. 22, pp. 3839–3866, 1996.
[11] G. R. Liu and Y. T. Gu, “A meshfree method: meshfree weak-strong (MWS) form method, for 2-D solids,” Computational Mechanics, vol. 33, no. 1, pp. 2–14, 2003.
[12] B. Nayroles, G. Touzot, and P. Villon, “Generalizing the finite element method: diffuse approximation and diffuse elements,” Computational Mechanics, vol. 10, no. 5, pp. 307–318, 1992.
[13] R. L. Hardy, Multiquadric equations of topography and other irregular surfaces, Journal of Geophysical Research 76 (1971), no. 8, 1905–1915.
[14] Sarra, S A 2017 The Matlab Radial Basis Function Toolbox. Journal of Open Research Software, 5: 8, DOI: https://doi.org/10.5334/jors.131
[15] J. Duchon, Splines minimizing rotation-invariant semi-norms in sobolev spaces, in: W. Schempp, K. Zeller (Eds.), Constructive Theory of Functions of Several Variables, Springer, Berlin Heidelberg, 1977, pp. 85–100.
[16] E. J. Kansa, Multiquadrics – A scattered data approximation scheme with applications to computational fluid dynamics I: Surface approximations and partial derivative estimates, Computers and Mathematics with Applications 19 (1990), no. 8/9, 127–145.
[17] E.J. Kansa, Multiquadrics – A scattered data approximation scheme with applications to computational fluid-dynamics – II solutions to parabolic, hyperbolic and elliptic partial differential equations, Comput. Math. Appl. 19 (8) (1990) 147–161.
[18] Fasshauer GE. Solving partial differential equations by collocation with radial basis functions, In: Mehaute A, Rabut C, Schumaker L L, editors. Surface Fitting and Multiresolution Methods, Vanderbilt University Press, Nashville. 1997; 131–138.
[19] Rosales AH, Rocca A La, Power H. Radial basis function Hermite collocation approach for the numerical simulation of the effect of precipitation inhibitor on the crystallization process of an over-saturated solution. Numerical Methods for Partial Differential Equations. 2006; 22(2):361–380.
[20] Li J C, Chen CS. Some observations on unsymmetric radial basis function collocation methods for convectiondiffusion problems. International Journal for Numerical Methods in Engineering. 2003; 57(8):1085–1094.
[21] Wendland H. Fast evaluation of radial basis functions methods based on partition of unity. In: Approximation Theory X (St. Louis, MO, 2001), Vanderbilt University Press, Nashville, TN, 2002; 473–483.
[22] Fornberg B, Piret C. A stable algorithm for flat radial basis functions on a sphere. SIAM Journal on Scientific Computing. 2007; 30:60–80.
[23] Šarler B, Vertnik R. Meshfree explicit local radial basis function collocation method for diffusion problems. Computers and Mathematics with Applications. 2006; 51(8):1269–1282.
[24] Ling L, Kansa EJ. A least-squares preconditioner for radial basis functions collocation methods. Advances in Computational Mathematics. 2005; 23:31–54.
[25] S.Kazem, J.A.Rad, K.Parand, Radial basis functions methods for solving Fokker–Planck equation, Engineering Analysis with Boundary Elements, Volume 36, Issue 2, February 2012, Pages 181-189.
[26] Siraj-ul-Islam, R.Vertnik, B.Šarler, Local radial basis function collocation method along with explicit time stepping for hyperbolic partial differential equations, Applied Numerical Mathematics, Volume 67, May 2013, Pages 136-151
[27] Mohan K. Kadalbajoo, Alpesh Kumar, Lok Pati Tripathi, A radial basis functions based finite differences method for wave equation with an integral condition, Applied Mathematics and Computation, Volume 253, 15 February 2015, Pages 8-16.
[28] Pedro González Casanova, Christian Gout, Jorge Zavaleta, Radial basis function methods for optimal control of the convection–diffusion equation: A numerical study, Engineering Analysis with Boundary Elements, Volume 108, November 2019, Pages 201-209
[29] H.Laeli Dastjerdi, M. Nili Ahmadabadi, The numerical solution of nonlinear two-dimensional Volterra–Fredholm integral equations of the second kind based on the radial basis functions approximation with error analysis, Applied Mathematics and Computation, Volume 293, 15 January 2017, Pages 545-554
[30] Pouria Assari, Mehdi Dehghan, The approximate solution of nonlinear Volterra integral equations of the second kind using radial basis functions, Applied Numerical Mathematics, Volume 131, September 2018, Pages 140-157.
[31] Chandhini G., Prashanthi K.S., V. AntonyVijesh, A radial basis function method for fractional Darboux problems Engineering , Analysis with Boundary Elements, Volume 86, January 2018, Pages 1-18
[32] CécilePiret, EmmanuelHanert, A radial basis functions method for fractional diffusion equations, Journal of Computational PhysicsVolume 238, 1 April 2013, Pages 71-81.
[33] Gomez CA, Salas AH, Frias BA (2010) New periodic and soliton solutions for the generalized BBM and BBM-Burgers equations. Appl Math Comput 217:1430–1434
[34] Guo C, Fang S (2012) Optimal decay rates of solutions for a multidimensional generalized Benjamin–Bona–Mahony equation. Nonlinear Anal Theor Appl 75:3385–3392
[35] Kadri T, Khiari N, Abidi F, Omrani K (2008) Methods for the numerical solution of the Benjamin–Bona–Mahony–Burgers equation. Numer Methods Part Differ Equ 24(6):1501–1516
[36] Noor MA, Noor KI, Waheed A, Al-Said EA (2011) Some new solitonary solutions of the modified Benjamin–Bona–Mahony equation. Comput Math Appl 62:2126–2131
[37] Abbasbandy S, Shirzadi A (2010) The first integral method for modified Benjamin–Bona–Mahony equation. Commun Nonlinear Sci Numer Simul 15:1759–1764
[38] Abdollahzadeh M, Hosseini M, Ghanbarpour M, Kashani S (2011) Exact travelling solutions for Benjamin–Bona–Mahony–Burgers equations by G G′—expansion method. Int J Appl Math Comput3:70–76
[39] Omrani K, Ayadi M (2008) Finite difference discretization of the Benjamin–Bona–Mahony–Burgers equation. Numer Methods Part Differ Equ 24(1):239–248
[40] Omrani K (2006) The convergence of fully discrete Galerkin approximations for the Benjamin–Bona–Mahony (BBM) equation. Appl Math Comput 180:614–621
[41] Qinghua X, Zheng CZ (2012) Degenerate boundary layer solutions to the generalized Benjamin–Bona–Mahony–Burgers equation. Acta Math Sci 32:1743–1758
[42] Xiao Q, Zhao H (2013) Nonlinear stability of generalized Benjamin–Bona–Mahony–Burgers shock profiles in several dimensions. J Math Anal Appl 406:165–187.
[43] Yin H, Hu J (2010) Exponential decay rate of solutions toward traveling waves for the Cauchy problem of generalized Benjamin–Bona–Mahony–Burgers equations. Nonlinear Anal Theory Methods Appl 73:1729–1738.
[44] M. Dehghan, M. Abbaszadeh, A. Mohebbi, The use of interpolating element-free Galerkin technique for solving 2D generalized Benjamin-BonaMahony-Burgers and regularized long-wave equations on non-rectangular domains with error estimate, Comput. Math. Appl. 286 (2015) 211–231.
[45] H. Tari, D.D. Ganji, Approximate explicit solutions of nonlinear BBMB equations by He’s methods and comparison with the exact solution, Phys. Lett. A 367 (2007) 95–101.
[46] Z.Z. Ganji, D.D. Ganji, H. Bararnia, Approximate general and explicit solutions of nonlinear BBMB equations by exp-function method, Appl. Math. Model. 33 (2009) 1836–1841.
[47] Mehdi Dehghan, Mostafa Abbaszadeh, Akbar Mohebbi, The numerical solution of nonlinear high dimensional generalized Benjamin–Bona–Mahony–Burgers equation via the meshless method of radial basis functions, Computers and Mathematics with Applications 68 (2014) 212–237.
[48] Ali Ebrahimijahan and Mehdi Dehghan, The numerical solution of nonlinear generalized Benjamin–Bona–Mahony–Burgers and regularized long‑wave equations via the meshless method of integrated radial basis functions, Engineering with Computers (2019). https://doi.org/10.1007/s00366-019-00811-3
[49] S. Rippa, An algorithm for selecting a good value for the parameter c in radial basis function interpolation, Adv. Comput. Math. 11 (1999) 193–210.
[50] I. R. H. Jackson, Convergence properties of radial basis functions, Constructive Approximation 4 (1988), no. 1, 243–264.
[51] R. B. Platte, How fast do radial basis function interpolants of analytic functions converge?, IMA Journal of Numerical Analysis 31 (2011), no. 4, 1578–1597.
[52] Q. Liu, Y. Gu, P. Zhuang, F. Liu, Y. Nie, An implicit RBF meshless approach for time fractional diffusion equations, Comput. Mech. 48 (2011) 1–12.
[53] A. Quarteroni, A. Valli, Numerical Approximation of Partial Differential Equations, Springer-Verlag, New York, 1997.
[54] Sirajul Haq, Abdul Ghafoor, Manzoor Hussain, Shamsul Arifeen, Numerical solutions of two dimensional Sobolev and generalized Benjamin–Bona–Mahony–Burgers equations via Haar wavelets, Computers & Mathematics with Applications Volume 77, Issue 2, 15 January 2019, Pages 565-575.
[55] M. Hajiketabi, S. Abbasbandy, F. Casas, The Lie-group method based on radial basis functions for solving nonlinear high dimensional generalized Benjamin–Bona–Mahony–Burgers equation in arbitrary domains, Applied Mathematics and Computation 321 (2018) 223–243.