Exact Solutions of the Nonlinear Modiﬁed Benjamin-Bona-Mahony Equation by an Analytical Method

: The current manuscript investigates the exact solutions of the modiﬁed Benjamin-Bona-Mahony (BBM) equation. Due to its efﬁciency and simplicity, the modiﬁed auxiliary equation method is adopted to solve the problem under consideration. As a result, a variety of the exact wave solutions of the modiﬁed BBM equation are obtained. Furthermore, the ﬁndings of the current study remain strong since Jacobi function solutions generate hyperbolic function solutions and trigonometric function solutions, as liming cases of interest. Some of the obtained solutions are illustrated graphically using appropriate values for the parameters.


Introduction
Nonlinear partial differential equations are widely used in science and engineering to model various nonlinear scenarios in real-life applications. For instance, these equations can be used in modeling fluid dynamical problems, wave propagation in corrugated media, the examination of earthquakes and seismic waves, and modeling optical fibers, etc. Despite being an old field of great attention, the field of fluid dynamics is still relevant. Various nonlinear dynamical models have been proposed for modeling the dynamics of surface water waves, such as the Benjamin-Bona-Mahony equation [1], Boussinesq equation [2,3], Korteweg-de-Vries equation and its modifications [4,5], Benney-Luke equation [6] and so on. More specifically, the model of interest in this manuscript is the modified Benjamin-Bona-Mahony (BBM) equation [7][8][9][10] and the references cited therein. In essence, the BBM equation is a nonlinear evolution equation that is used to model the dynamics of unidirectional water waves with small amplitude, in addition to being disturbed by dispersion and nonlinear effects [11]. Moreover, one could find various generalizations and modifications of this model, including its reduction to the Korteweg-de-Vries equation under a special assumption. However, the Korteweg-de-Vries equation is integrable via the inverse scattering approach, as against the current model of interest.
As highlighted above with regard to the unidirectional water waves with small amplitude, such propagation of waves is called a solitary wave [12][13][14]. Moreover, there exist several analytical and computational methods for finding the exact and approximate solitary wave solutions associated with diverse nonlinear evolution equations. In light of this, we mention some of the analytical methods in the literature, such as the Jacobi elliptic functions method [15], the exponential expansion method [16], the tanh-based expansion methods [17], the Kudryashov method [18], the G'/G-expansion method [19], the sub-equation method [20] and the modified exponential rational method [21]. For more on analytical methods, see [22][23][24][25][26][27][28]. Furthermore, we mention some of the well-known computational methods for approximate soliton solutions here, including the famous Adomian decomposition approach [29], the Laplace-homotopy perturbation method [30] and the homotopy analysis approach [31], among others.
The present paper is devoted to examining the modified BBM equation analytically. It is relevant to state here that four types of analytical solutions, including exponential [7], periodic and hyperbolic [8][9][10] and rational [9,10] solutions have been obtained for the model in the literature. This study objectively intends to go beyond what is known in the open literature [7][8][9][10] by constructing more general exact solutions for the model via Jacobi functions. Thus, the modified auxiliary equation method [28] is adopted for the treatment in order to obtain these kinds of solutions. Moreover, the exact solutions are comprehensively analyzed further by determining their constraint conditions (if any). In addition, we aim to examine the obtained solutions graphically, as well as classifying them based on the features of the resultant profiles.
The current paper is organized as follows. Section 2 presents the outline of the methodology of concern. Section 3 applies the method presented in Section 2 on the governing modified BBM model. Section 4 is concerned with the discussion of the obtained results, while Section 5 gives some closing remarks as a concluding note.

The Modified Auxiliary Equation Method
In this section, we give the outline of the methodology of concern, the modified auxiliary equation method [28]. Thus, considering the partial differential equation below, we present the main steps of the method as follows: (1) Step 1. First, we use the transformation where c is a constant to be determined. Then, upon employing the transformation given in Equation (2) in Equation (1), the following nonlinear ordinary differential equation is obtained: Step 2. We assume that the solution of Equation (3) is given by where n is a positive integer, and η i are arbitrary constants to be determined. Moreover, where ν 0 , ν 1 , ν 2 are arbitrary constants. Moreover, Equation (5) admits the following solutions: (5) has a solution ψ(ξ) = sn(ξ, k), where sn(ξ, k) defines the Jacobi function and k denotes the elliptic modulus such that 0 < k < 1. (5) has a solution ψ(ξ) = cn(ξ, k), and cn(ξ, k) defines the Jacobi function and k is as defined earlier.
Step 3. The value of n in Equation (4) can be determined via the application of the homogeneous balancing principle [28].
Step 4. After substitution of Equation (4) together with Equation (5) into Equation (3) and setting all terms with the same power of ψ(ξ) equal to zero, we find a set of overdetermined equations for η i . Consequently, a solution of Equation (1) is obtained.

Analytical Solutions of the Modified BBM Equation
In this section, the modified auxiliary equation method will be employed to analytically treat the governing modified BBM equation. Therefore, we consider the modified BBM equation, which is given by [7][8][9][10] where α is a nonzero constant. Hence, we shall use the transformation where c is a constant. Next, the governing model given in Equation (6) is written via Equation (7) in the following way: Using the homogeneous balancing principle on Equation (8), we find n = 1. Thus, the solution of Equation (8) is written as Then, after substituting Equation (9) with the use of Equation (5) into Equation (8) and thereafter setting the coefficients of ψ(ξ) equal to zero, we obtain a system of algebraic equations. Solving this system for η 0 , η 1 , η −1 and c, the following sets of solutions are obtained: Set 1.
On substituting these values into Equation (9), various exact solutions can be constructed as the following cases: Case 1. If ν 0 = 1, ν 1 = −(1 + k 2 ), ν 2 = k 2 , then the modified BBM equation in Equation (6) has a solution in the form The solution in Equation (11) further leads to when k → 1. Figure 1 represents the positive solutions given in Equations (11) and (12) when α = 1. From this figure, it is can be observed that Figure 1a shows a singular partial periodicity, while Figure 1b shows a singular behavior.

Set 3.
Then, substituting these results into Equation (9), we obtain the following solution cases. Figure 4 represents the solution determined in Equation (27) when α = 1 and k = 0.4. The 3D plot in Figure 4 shows a singular solution. Further, the same solution expressed in Equation (27) reduces to the following expression
By substituting these results into Equation (9), we find the following solution cases.
In the same manner, the following solution cases are obtained.

Discussion of the Results
The current manuscript takes into consideration one of the most famous evolution equations that plays an important role in the dynamics of fluid flows: the modified BBM equation [7][8][9][10]. Moreover, the model is realized in various nonlinear physical processes. A reliable analytical integration scheme obtained by the named modified auxiliary equation method [28] has been used for the analytical examination. Importantly, this reliable method leads to six (6) different sets of exact solutions, in addition to the various solution cases associated with each solution set. Furthermore, diverse Jacobi functions have been obtained by the present reliable method, in addition to the provision of various periodic and hyperbolic function solutions, as limiting cases of interest. Finally, it is highly recommended to make use of the deployed method to study other forms of evolution equations, such as higher-order evolution equations involving different forms of nonlinearities.

Conclusions
In conclusion, the current study examines the modified BBM equation by constructing a variety of exact solutions. We have used the modified auxiliary equation method due to its documented applications in tackling diverse evolution equations. Comparing our results with the results in the literature, our solutions are found for the first time since they are given in terms of Jacobi functions. Moreover, certain solution cases of interest have been graphed, showing different forms of behavior, including periodicity, singular periodicity and kick-type, among others. Lastly, the present study may give more insights into the dynamicity of the modified BBM equation with regard to its response to spatial and temporal evolution, as shown graphically. As a future perspective, the current study can be extended to other forms of evolution equations featuring different forms of nonlinearities.