The unidirectional propagation of long waves in certain nonlinear dispersive system is explained by the Benjamin-Bona-Mahony-Burger (BBM-Burger) equation. The purpose of this study is to investigate the BBM-Burger equation numerically using Caputo derivative and B-spline basis functions. The fractional derivative is considered in Caputo form, and formula is used for discretization of temporal derivative. The interpolation of space derivative is done with the help of B-spline functions. The effect of and time on solution profile of travelling wave for different domain of is discussed in this paper. The numerical results have been presented to show that the cubic B-spline method is effective and efficient in solving the time fractional Benjamin-Bona-Mahony-Burger (BBM-Burger) equation. Moreover, the convergence and stability of the proposed scheme are analyzed. The error norms are also calculated to check the accuracy of the proposed scheme. The numerical results reflect that the proposed scheme can be used for linear and highly nonlinear models.

1. Introduction

The fractional calculus has a long history, beginning on 30 September, 1695 [1]. Fractional derivatives are an excellent method for explaining the memory and genetic characteristics of various materials and processes. Today, fractional differential equation (FDE) plays an important role in number of fields such as mathematics and dynamic systems [2]. Noninteger-order derivatives and integrals have proved to be more useful in formulating certain electrochemical problems than conventional models [3]. Fractional partial differential equations (FPDEs) provide a powerful tool for modeling a variety of potentially emerging phenomena with a wide range of applications across applied sciences, engineering, and physics. It is the classical generalization of the calculation that involves operators of integration of noninteger orders and differentiation. They are considered to be valuable and efficient tools in both theoretical and technological fields to tackle the complexity and nonlinearity of specific problems, including chemical science, fluid flow, fiber optics, signaling system, polymers, device recognition, and elastic materials [4, 5]. A great deal of attention has been paid to managing fractional differential equations (FDEs). Specifically, fractional-order partial differential equations (PDEs) are increasingly described in viscoelasticity, mathematical biology, finance [68], airplane designing, traffic, population, and particle chemistry [911]. Elasticity, plasma, fluid mechanics, optical fiber, and other disciplines of mathematics and physics involve wave phenomena. The Benjamin-Bona-Mahony-Burger (BBM-Burger) equation, which addresses dissipation as well as nonlinearity and dispersion in the simulation of unidirectional plane waves, is defined as follows [12]: subject to IC: and BCs as follows: where and are positive constants and domain of .

Nonlinearity, dispersive, and dissipative effects are all included in Equation (1). The terms and indicate the effects of dissipation and dispersion, respectively. Cracked rock, thermodynamics, acoustic waves in aharmonic crystals [13], acoustic-gravity waves in fluids, and other applications of the Benjamin-Bona-Mahony-Burger (BBM-Burger) equation and their special cases can be found in engineering and science. Benjamin-Bona-Mahony-Burger (BBM-Burger) equation is used in the analysis of the surface waves of long wavelength in liquids, hydromagnetic waves in cold plasma, and acoustic-gravity waves in compressible fluids.

The water wave model is notoriously difficult to solve analytically. So, several investigations have been done in recent years to determine the numerical solution of water wave model (see, for example, [14, 15] and the references therein). Researchers introduced many different methods to develop an approximate analytical solution for the fractional differential equation and systems, such as the cubic B-spline method (CBS) by [16, 17] and homotopy analysis method (HAM) by [18], unified method (UM) [19], Runge-Kutta method [20], Hermite wavelet technique together with Newton Raphson iteration method [21], and Lie symmetry method [22]. The time fractional Benjamin-Bona-Mahony-Burger (BBM-Burger) equation was proposed to discuss the dynamic physical system. The time fractional operator form of BBM-Burger equation used in [23] is

where is the fractional order of the equation. The initial and boundary conditions are respectively. The exact solution of Equation (2) for is

Many scholars contributed by developing new analytical and numerical schemes to solve fractional model of BBM-Burger. Zhang et al. [23] implemented modified residual power series method (RPSM) for the simulations of BBM-Burger equation. expansion method has been implemented to investigate the exact solution of time fractional model of BBM-Burger equation in [24]. The authors in [25] utilized Lie symmetry approach to develop solution of fractional BBM equation. New soliton solutions for -dimensional Extended Date–Jimbo–Kashiwara–Miwa Equation has been established by means of Hirota method together with quadratic test functions in [26]. In [27], the authors studied numerical solution of -dimensional shallow water wave equation and explicit traveling wave solutions for Zoomeron evolution equation using Lie vector technique. Ali et al. [28, 29] developed new method based on hybrid orthonormal Bernstein and block-pulse function wavelet method (HOBW) for the solution of fractional voltera-Fredholm differential equation and for the solution of 2D Bratu problem. Runge-Kutta (RK4) method and Keller–Box scheme have been implemented for the analysis of the melting heat transport of blood with time-dependent cross nanofluid model by [30]. The authors in [31] developed the Gudermannian neural networks to solve the nonlinear fifth-order induction motor model (FO-IMM). Ali [32, 33] analyzed heat transfer of Ostwald-de Waele fluid using method of lines and penalty-artificial compressibility technique (PACT) for studying free convection of magnetohydrodynamic unsteady incompressible flow in a square cavity, respectively. New homotopy analysis transform method has been established in [34] to constitute analytical solution of BBM-Burger equation. The exact solution of BBM-Burger equation of fractional order has been developed by HAM in [12]. The author in [35] investigated the solution of Benjamin-Bona-Mahony equation using linearized difference scheme. The linear difference scheme with Crank-Nicolson has been implemented in [36] for the treatment of Benjamin-Bona-Mahony equation with a time fractional nonlocal viscous term. Majeed et al. [37] employed cubic B-spline scheme for the approximate solution of time fractional nonhomogeneous BBM-Burger equation. [38] used Pohozhaev nonlinear capacity method to investigate the blowing-up solutions of different nonlinear fractional partial differential equations like KdV, time fractional modified KdV-Burger, BBM-Burger, Rosenau, Ostrovsky, and Degasperis-Procesi. [39] has investigated the analytical solution of time fractional Benjamin-Bona-Mahony equation using method based on Laplace Transformation and Adomian decomposition method. Two problems have been solved with different initial conditions. In [40], the authors have derived the Benjamin-Bona-Mahony and KdV equations from water wave model. Further travelling wave equations of these models have been derived by two different methods. Space-time fractional BBM-Burger models have been solved analytically using Adomian decomposition method by [41]. The energy method is used by [42] to prove the uniqueness and existence of BBM-Burger model. The authors also used the Fourier Spectral method to get numerical solution of this model. In [43], the authors used EMAEM method to get different solutions of 3D fractional WBBM models. The sine-cosine function method is used in [44] to investigate the analytical solution of Benjamin-Bona-Mahony and MRLW models. The modified simple method and expansion method have been used in [45] to find the exact solution of TFBBM equations.

In this study, cubic B-spline method along with Caputo derivative will be used to solve time fractional BBM-Burger equation numerically. formula will be used to tackle the noninteger-order derivative, while B-spline basis functions will be used to discretize the space derivative. The advantages of B-spline collocation method over the existing schemes are that it efficiently delivers a piecewise-continuous, closed form solution and it is simpler and can be used to a wide range of problems involving partial and fractional partial differential equations. The advantage of this method over others is that once the solution has been computed, the information needed for interpolation at different locations in the interval is available. The key benefit of this method is that it can approximate the analytical curve to a given degree of smoothness. As a result, the spline approach, when compared to previous methods, has the flexibility to obtain an estimate at any point in the domain with more accurate results. The behavior of solution profile for different domains of will be studied, and the effect of parameter will be discussed for various values. The objective of this study is to solve nonlinear fractional-order partial differential equation with mixed partial derivative numerically using scheme based on cubic B-spline bases. We proved that the proposed scheme is unconditionally stable.

The organization of the paper is as follows: Section 2 gives brief introduction and preliminaries. Section 3 presents the algorithm and implementation of the suggested scheme to solve the fractional model of BBM-Burger. The stability and convergence of the proposed scheme are analyzed in Sections 4 and 5, respectively. Discussion on results is presented in Section 6, and Section 7 reviews the concluding remarks.

2. Basic Definitions

Definition 1. The -order Caputo derivative [3] is defined as

Definition 2. The Caputo derivative is discretized by formula [46]

Lemma 3. The coefficients in Equation (6) satisfy the following properties [16]: (i)(ii)(iii)

Definition 4. The generalized linearization formula used in [47] is applied in this work. Here, is for positive integer.

3. Proposed Scheme

To implement the numerical schemes, let be the domain for time and is splitted up into uniformly sized subintervals : where The time fractional derivative in Equation (5) for can be evaluated at by the formula [46]. Using Equation (6) together with forward difference scheme and Crank-Nicolson for time derivative, Equation (2) becomes the following: where and . Replacing Equation (7) to Equation (8), we get

To apply the proposed scheme, the domain of spatial variable is partitioned into subdomain by , where and The solution of Equation (9) at time step is given by where are unknowns to be determined and are spline function basis of third degree defined as follows [16]:

Different derivatives like , , and of Equation (10) in terms of at knot points can be expressed as follows:

The variation of over the typical component is given by

On substituting Equation (10) and its necessary derivatives of Equation (12) in Equation (9), we obtain the following set of equations: where

Equation (14) will be solved for first, and then, loop will be executed for , which gives system of linear equations with unknowns . To get consistent solution of this system, we need two more equations, for this derivative equations will be used.

Equations (14) and (16) have the matrix form as follows: where

The approximate solution for time level can be achieved for by using vector in Equation (10).

3.1. Initial Value

To compute the initial value , the enforced conditions of the given problem will help. So, Equation (10) becomes

To determine , the derivatives of the given auxiliary conditions are approximated as follows [48]. (i) for (ii) for (iii) for

This gives matrix equation.

The numerical simulations have been done with the help of the MATLAB 2015b software.

4. Stability Analysis

The most frequently used formula to determine the stability is von Neumann Fourier formula [47]. Solution will be considered in single Fourier mode.

Using Equation (21) in Equation (14), we have

By putting values of and and after simplification, we have where , , and were defined above. If augment factor , then the suggested scheme will be stable. Otherwise, it will be unstable.

The above result indicates that the proposed scheme is stable unconditionally.

5. Convergence Analysis

Let be the exact solution of Equation (2). Assume that which takes different values at various locations as in [49] and is independent of , , , and with the smoothness conditions that

We first linearize the term by taking as a constant “.” The linearized form of the proposed scheme is then where .

After simplifying Equation (26), we get where , , and Using Equation (25), for the scheme Equation (27), we have where is the exact solution at point . where and is a constant.

Theorem 5. If and are the solutions of given equation and Equation (27), respectively, and if satisfies smoothness condition Equation (25), then for sufficiently small and ,

Proof. Subtracting Equation (27) from Equation (28), we obtain where . Multiplying the above equation by and summing over from to , we obtain Using the inequality , we obtain Since , therefore Dividing Equation (33) throughout by , we get where . Note that we have used . This completes the proof.

6. Discussion

The numerical results obtained by the suggested scheme for various parameters are reported in this part in the form of tables and figures. and norms are computed and defined as

The approximate solution (10) can be written in piecewise form as follows:

The numerical solution of travelling wave models is discussed graphically and numerically in our computational work, the values of the parameter , step length , and are chosen, and for to computation, it was carried out over different finite space range of , , and . Figure 1 shows the wave evolution for different time with . Figure 1(a) depicts the behavior of travelling wave at , while Figure 1(b) demonstrates the identical patron of travelling waves for different time , 0.75, 2.5, and 5. The domain of varies from -5 to to observe the behavior of wave, and it is clearly shown that both solutions are in good consensus with each other with identical patron of waves.

The graphical illustration displayed in Figure 2 expresses the behavior of solution profile for different time steps , 2.5, and 5 in broader domain of and keeping , , and fixed. Figure 2 express the behavior of wave evolution at , 2.5, and 5 in wide domain of , while keeping , , and fixed in this case. As the number of iterations is raised, the solution profile produces better results, and the traveling wave gets smoother in a larger domain. The wave evolution is also observed at different values of . It is observed that as the value of approaches to 1, the numerical solution tends toward the solution, and if we use the alpha value near 0, numerical value moves far away from the solution as shown in Figure 3. The 3D form of numerical solution with different time level is shown in Figure 4 effect of is observed. The efficiency of the proposed scheme in tabular form for this model can be observed by Table 1. In this table, numerical and solutions are compared. Table 2 presents the and norms.

To study the effect of solution profile on a larger domain, we extended the space domain from -20 to 20 in Figure 5. The parameters for this simulation were chosen , , and for different times , 0.75, 2.5, and . Figure 4 displays the solution profile for time , while the effect of different time on solution profile of travelling wave is shown in Figure 4. The graphics illustrate the smoothness of the solution profile. Figure 6 shows how different fractional parameters, such as , affect the solution profile.

The wave swings for small values of in the broader domain and then smooths out and approaches the classical solution for . The numerical solution of wave model in 3D form can be seen in Figure 7. Both solutions are extremely similar to each other for different parameters, as seen in the graphics. The and approximate outcomes are compared in Table 3. Both solutions are comparable up to four decimal places, as may be shown from tabular data. The error norms obtained at are shown in Table 4. The behavior of solution profile for broader domain has been studied and illustrated in the figures. These graphics demonstrate that solution curves become smooth and swing as we extend the domain, and both results express the identical behavior. So, it is observed that the proposed scheme gives better results for smaller and larger domains but if the spatial derivative in the model is of third or fourth order, then the proposed scheme has to be modified.

7. Concluding Remarks

In this work, an algorithm based on B-spline functions and Caputo derivative has been successfully implemented on fractional BBM-Burger equation. The temporal derivative is considered in Caputo’s form, and discretization has been done by formula. The space derivative has been interpolated by means of Crank-Nicolson and B-spline function. Applicability of the scheme is examined on test problem, and the following points are extracted. (i)The behavior of travelling wave is studied for different domains of , and it is observed that as we extend the domain of , the solution profile improves and wave becomes more smooth(ii)The effect of the fractional parameter on various spatial ranges has been studied, and it has been discovered that the wave varies for smaller values of before becoming smooth for larger values of (iii)The approximation solution tends to overlap the classical solution, as may be seen in visuals and tabular data(iv)The study’s findings on the stability of the suggested scheme reveal that it is unconditionally stable(v)The error norms for different domains of are presented in the form of table which reflects that the proposed scheme is suitable for any domain of

In general, cubic B-spline scheme is an effective and convenient method in finding approximate solution for the time fractional BBM-Burger equation and other long waves in certain nonlinear dispersive systems.

Data Availability

The experimental data used to support the findings of this study are available within this paper.

Conflicts of Interest

The authors declare that they have no conflicts of interest to report regarding the present study.

Authors’ Contributions

All authors equally contributed to this work. All authors read and approved the final manuscript.


We thank Dr. Muhammad Kashif Iqbal for his assistance in proofreading of the manuscript.