On abundant new solutions of two fractional complex models

We use an analytical scheme to construct distinct novel solutions of two well-known fractional complex models (the fractional Korteweg–de Vries equation (KdV) equation and the fractional Zakharov–Kuznetsov–Benjamin–Bona–Mahony (ZKBBM) equation). A new fractional definition is used to covert the fractional formula of these equations into integer-order ordinary differential equations. We obtain solitons, rational functions, the trigonometric functions, the hyperbolic functions, and many other explicit wave solutions. We illustrate physical explanations of these solutions by different types of sketches.


Introduction
Fractional nonlinear evolution equation is one of the noticeable branches of science in recent years. Fractional calculus has a great profound physical background able to formulate many various phenomena in distinct fields such as physics, mechanical engineering, economics, chemistry, signal processing, food supplement, applied mathematics, quasichaotic dynamical systems, hydrodynamics, system identification, statistics, finance, fluid mechanics, solid-state biology, dynamical systems with chaotic dynamical behavior, optical fibers, electric control theory, economics, and diffusion problems. Mathematical modeling of these phenomena contains fractional derivatives, which provide a great explanation of the nonlocal property of these models since they depend on both historical and current states of the problem in contrast to the classical calculus depending on the current state only. Based on the importance of this kind of calculus, many definitions were derived such as conformable fractional, fractional Riemann-Liouville, Caputo, and Caputo-Fabrizio derivatives [7,8,23,24,41,43,50]. These definitions are employed to convert fractional nonlinear partial differential equations to nonlinear integer-order ordinary differential equations, and then computational and numerical schemes can be applied to get various types of solutions for these models and examples of these schemes [3, 9, 11-19, 21, 22, 25, 32, 35, 36, 39, 40, 42, 44, 45, 51-53, 57].
This method depends on a new auxiliary Riccati equation [47]. The auxiliary equation of the mK method is given by where [δ, , χ, Q] are arbitrary constants such that [Q = 0, Q = 1], whereas the Riccati equation is given by where [E 0 , E 1 , E 2 ] are arbitrary constants. So Eqs. (1) and (2) . Using this technique leads to the mK auxiliary equation, which includes many other analytical methods, but the mK method can obtain more solutions than most of them. This shows the superiority, power, and productivity of the mK method.
In this context, we employ the mK method to construct new formulas of solutions for the fractional KdV and ZKBBM equations, which are given, respectively, by [20,26,37,46,49,54,55] where [λ, ν, μ] are arbitrary constants. The KdV model is one of the essential models in studying the shallow-water waves, and it has a strong physical impact in describing the interaction of two long waves with various dispersion relations. It is used only for the instant of time (local property); that is why the solitary wave in the soliton solutions of it may behave not very well, whereas the fractional KdV is used to estimate the effect of higher-order dispersion of the regular KdV equation to increasing the amplitude of the soliton. On the other hand, the fractional ZKBBM equation is used to investigate the gravity water waves in the long-wave regime.
In this research, we use new fractional derivative operator defined as follows.

B(α) being a normalization function. Thus
Applying this definition of ABR fractional operator to Eqs. (3) and (4), respectively, with the wave transformation [K( , where k, ω are arbitrary constants, leads to conversion of Eq. (3) and (4) into the corresponding ODEs. Integration of the obtained ODEs with zero constant of integration gives Calculating the homogeneous balance value in Eqs. (8) and (9) yields N = 1. Thus both equations have the same general formula of solution given according to the mK method by The rest of the paper is organized as follows. In Sect. 2, we apply the mK method to the nonlinear fractional Kdv and ZKBBM equations. Moreover, we give some sketches to show more physical properties of both models. In Sect. 4, we discuss the obtained computational results and compare them with those obtained in previous works. Moreover, we compare the obtained numerical results. In Sect. 5, we give the conclusion of the whole research.

Abundant wave solutions of the fractional KdV and ZKBBM equations
In this section, we apply an analytical scheme to the nonlinear fractional KdV and ZKBBM equations and show physical properties of the two models.

The fractional KdV equation
Applying the mK method with its auxiliary equation and the suggested general solutions of the fractional KdV equation leads to a system of algebraic equations. Using Mathematica 11.2, we find the values of the parameters in this system, which lead to two families of solutions.
Family I Consequently, the closed-form solutions for the fractional KdV models are given as follows. When When [χ = 0 & = δ], When [χ 2 -4δ = 0], Family II Consequently, the closed-form solutions for the fractional KdV models are given as follows.

The fractional ZKBBM equation
Applying the mK method with its auxiliary equation and the suggested general solutions for the fractional ZKBBM equation leads to a system of algebraic equations. Using Mathematica 11.2, we find the values of the parameters in this system, which lead to the following families of solutions.
Family I Consequently, the closed-form solutions for the fractional ZKBBM models are given as follows.

Interpretation of figures
In this section, we give a physical interpretation of the shown figures. All our obtained solutions are considered as traveling wave solutions. We further give a physical interpretation of the shown figures: 1. Fig. 1 shows the bright cone wave solution (13) in the three-dimensional plot (a) to explain the perspective view of the solution, the two-dimensional plot (b) to explain the wave propagation pattern of the wave along the x-axis, and the contour plot (c)  Fig. 2 shows the dark cone wave solution (14) in the three-dimensional plot (a) to explain the perspective view of the solution, the two-dimensional plot (b) to explain the wave propagation pattern of the wave along the x-axis, and the contour plot (c) to explain the overhead view of the solution when [α = 1 2 , δ = 6, λ = 3, m = 1, n = 1, χ = 5, = 1]. 3. Fig. 3 shows the periodic bright cone-wave solution (40) in the three-dimensional plot (a) to explain the perspective view of the solution, the two-dimensional plot (b) to explain the wave propagation pattern of the wave along the x-axis, and the

Results and discussion
This section is divided into two main parts. In the first part, we studyg the obtained computational solutions for the two fractional suggested models. whereas in the second part, we compare them with the other obtained results in previous works. • In [56], two analytical methods are applied to three different models involving two our models. However, they use two schemes, but a very few special solutions are obtained. • Two analytical schemes in [56] are just particular cases of the mK method when

Conclusion
In our paper, we solved the flaws and disadvantages of the ( G G )-expansion method that is used by Ali Akbar et al. [56] because, as shown in the previous section, it is just a particular case of our method. Moreover, we use a new definition of fractional derivative, which successfully converts the fractional-order differential equations from our models to integer-order ordinary differential equations. Abundant new solutions for both models were obtained, and to further clarify the physical meaning of these solutions, some plots are sketched in three-and two-dimensional and contour plots (Figs. 1, 2, 3, 4).