#### Abstract

Nonlinear evolution equations widely describe phenomena in various fields of science, such as plasma, nuclear physics, chemical reactions, optics, shallow water waves, fluid dynamics, signal processing, and image processing. In the present work, the derivation and analysis of Lie symmetries are presented for the time-fractional Benjamin–Bona–Mahony equation (FBBM) with the Riemann–Liouville derivatives. The time FBBM equation is reduced to a nonlinear fractional ordinary differential equation (NLFODE) using its Lie symmetries. These symmetries are derivations using the prolongation theorem. Applying the subequation method, we then use the integrating factor property to solve the NLFODE to obtain a few travelling wave solutions to the time FBBM.

#### 1. Introduction

Partial differential equations running into the thinking of most of the researchers as it represented the importance in several topics of scientific fields as mechanics, optical fibers, medical sciences (as breast cancer), biological science, turbulent bursts, and oceans waves [1–11].

The differential model has a broad application in many phenomena as in [12–14]. Recently, nonlinear fractional differential equations (NLFDEs) show significantly in engineering and applications of other sciences, for example, electrochemistry, physics, electromagnetics, and signal data processing [15–22].

Getting exact solutions for these forms of equations became an important issue; then, most researchers try to achieve this target. The most effective method for obtaining exact solutions for NLFPDEs is the Lie symmetry reduction method. There are many papers for using Lie’s method to obtain explicit solutions for NLFPDEs [23–26].

In our paper, we drive the symmetry vectors for the time FBBM equation and present new closed-form solutions for it. The FBBM equation has many forms [27–30], and we choose to work on the following form:where is the dissipative term.

The manuscript is prearranged as follows. In Section 2, Lie’s group method for FPDEs is exposed. In Section 3, we apply the Lie group reduction method to obtain Lie point symmetry for the time FBBM equation (1). At the end of Section 4, we use these similarity variables to get the reduced equation. In Section 4, we use two methods for solving the resulting ordinary differential equation, the first method is the subequation method and the second method is the integrating factors method to get new solutions that have the properties and form the travelling wave form for the FBBME. In the end, conclusions are written in Section 5.

#### 2. Notations and Introductory

##### 2.1. Fractional Riemann–Liouville Derivative

In this section, we show some definitions for RL fractional derivative [31], which can be considered as follows:where is the total differentiation of integer number of orders , the Gamma function is , and is the (RL) fractional integral of an order of .

*Definition 1. *The partial derivative of order for Riemann–Liouville definition is presented by

##### 2.2. Notations for Lie Symmetry Reduction Method for the Time FPDEs

In this section, we show in detail the main notations and definitions that will be used for obtaining the symmetries of NLFPDEs.

Here, we will consider timing NLFPDEs of the form [31–33]

Assume, equation (2) has a Lie vector X in the formwhere , , and can be called as the infinitesimals of the transformations the independent and the dependent variables , respectively. Let a one-parameter Lie algebra of infinitesimal transformations be of the following form:where « 1 can be defined as a group parameter, in most cases we take it equal one. The explicit expressions of , and , which can be called the prolongation of the infinitesimals and are given bywhere is the total differentiation operator [34] with respect to the independent variables (*i* = 1, 2, then ):

Theorem 1. *Equation (1) concedes a one-parameter group of infinitesimal transformations in equation (2) with the Lie Vector X if and just if the accompanying infinitesimal conditions hold:where and is the 3rd prolongation of the infinitesimal generator X.*

*Definition 2. *Prolonged vector is given by [31]where *q* is the numbers of dependent variables, is the numbers of independent variables, , and PDE involve derivatives up to order *n*. Also, the invariance condition [35] givesThe extended infinitesimal, which deals with fractional derivatives, has the following form [36–38]:whereRemember thatDue to linearization of the infinitesimal in *u* and the presence of , will vanish, where in equation (14).

Lemma 1. *The function can be defined as an invariant solution of (2) if and only if* *(i) is an invariant surface of equation (2)* *(ii)*

#### 3. Lie Symmetry and Reduction of FBBM Equation

In this partition, the Lie symmetry reduction method was applied to find the similarity variables for a one-dimensional time (FBBM) equation. Suppose that (1) is an invariant under (2); we have that

Thus, satisfies equation (1). Applying the third prolongation to (1), we have the accompanying deciding condition, which is given as

Substituting (7) and (8) into (16) and equaling coefficients in derivatives for and power of to zero, the system of equations is obtained:

By solving the obtained equations in (18), we get the following infinitesimal:where are constants. By the previous infinitesimal, equation (1) has two vector fields in the form

*Case 1. *For the infinitesimal generator in (20a), we have a characteristic equation in the following form:By solving the previous equation, we get the variables and Putting into (1), we obtain the following fractional ODE:By solving the above equation, we obtainwhere is constant of integration.

*Case 2. *For in equation (20b), the similarity variables for the infinitesimal generator can be obtained from the equation:The previous equation is called the characteristic equation; by solving it, we have the similarity variable as a result in the form:The group invariant solutionwhere is a new arbitrary function of and By using equation (26), equation (1) is transformed into FODE.

Theorem 3. *The transformation in (25) and (26), which is obtained from the similarity group method, reduces equation (1) to NLFODE as below:*

Using the operator EK fractional differential operator [32, 34],where

*Proof. *Let Depending on the Riemann–Liouville (RL) derivatives, definitions, and similarity variables in (25) and (26), we obtainLet and . Thus, (31) becomesSubstitute the EK fractional operator in (30) into (32), we haveFor simplicity, let and ; we acquire.

Hence, equation (33) will be rewritten asRepeating times, we haveUsing the definition of EK fractional differential operator in (28) to rewrite (35), we obtain

*Remark 1. *The fractional derivative must achieve the linearization property [37, 39]:Using the invariant group solution in (26), (37), and (38), we obtainHence,Thus, (1) can be reduced toand the theorem is totally proofed.

#### 4. Explicit Solutions for FBBM Equation

##### 4.1. Clarifications for the Subequation Method

The subequation method [39] is presented in this section. Consider the NLFPDE in the formwhere is a dependent variable, is a series of and its fractional derivatives, and are the Riemann–Liouville (RL) derivatives of w.r.t . Here, we present the principles for the subequation technique. By using the d’Alembert transformation,where is constant that will be determined later, and we can rewrite (41) as NLFODE:

According to the subequation procedure, assume that the wave solution will be written in the following form:where are constants, which will be determined later, belongs to integers numbers, which are determined by equaling the highest order derivatives and nonlinear terms in (44) together, and the function achieves the Riccati equation of fractional orderwhere is a constant. Some trigonometric solutions of the fractional Riccati equation (46) are

By substituting forms (45) into (44) and setting the coefficients of to be zero, we obtain an algebraic system in . By solving the determinate system, we obtain the constants . Substituting these constants and the solutions of (47) into (45), we obtain the closed form solutions of (42).

##### 4.2. Applying the Subequation Method to the Time FBBM Equation

We now implement a subequation method to (1). We will use the transformationwhere *c* is a constant, and this will transform (1) into an NLFODE:

We now assume that (49) has the solution in the formwhere (*i* = 1, ..., *n*) are constants, which will be determined, and achieves equation (46).

Balancing the highest order derivative terms with nonlinear terms in equation (49), we obtain *n* = 2. Hence,

We then substitute (51) along with (46) into (49), then collect the coefficients of , and set them to equal zero. A set of algebraic equations are obtained in knowns *c*, . Solving these algebraic equations with the help of the software program (Maple), we get the following values.

Thus, from (47), we obtain five forms of explicit travelling wave solutions of (1), namely,

Thus, from (47), we obtain five forms of explicit travelling wave solutions of (1), namely,where is arbitrary constant. We plot the result in equation (57) in the three dimensions, contour plot, and density plot, as shown in Figures 1–3, respectively.

##### 4.3. Applying Simple Transformation

We solve the conformable FBBM equation using simple transformation to change the fraction order in partial derivative to nonsolvable ODE. For the reduction of (1) to ODE, we use the following transformation:where and are arbitrary constants; we can rewrite (1) as NLODE:

This equation has no implicit solution but possesses two integrating factors. We apply the integrating factor technique to obtain an analytical solution for (59).

Equation (59) has two integrating factors (IF) as follows:

Using these integrating factors by the same steps in [39] and neglecting the constants of integration, equation (59) will be reduced to

By solving this equation, we obtain travelling wave solution for (1):

Replacing ,

In other manner, equation (59) have two Lie vectors. The first one of them reduces it to

Equation (64) has closed form solution, but, in the back substitution step, we are unable to get even if we neglect the values of constants. So, from here, we can say the integrating factor method for reducing and solve ODEs, occasionally, more effectiveness than the Lie reduction method. Result obtained in (63) is plotted in Figure 4 at different values of . We observe that, by decreasing the value of , the top of the wave has a parabolic shape.

Comparing our result in (63) with results in [5], specially equation (17), we find that the two solutions are travelling wave solutions, but the amplitude and direction of flow are different.

#### 5. Conclusions

In this paper, we show the importance and the effective of the Lie symmetry reduction method on the FBBM equations. We obtain time FBBM equation’s Lie symmetry generators and then reduce the equation to FODE using these symmetry vectors. The projected analysis is extremely effective and dependable for getting similarity solutions for fractional differential equations. New travelling solutions were derived for the FBBM equation using the subequation method.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.