The plethora of explicit solutions of the fractional KS equation through liquid–gas bubbles mix under the thermodynamic conditions via Atangana–Baleanu derivative operator

Novel explicit wave solutions are constructed for the Kudryashov–Sinelshchikov (KS) equation through liquid–gas bubbles mix under the thermodynamic conditions. A new fractional definition (Atangana–Baleanu derivative operator) is employed through the modified Khater method to get new wave solutions in distinct types of this model that is used to describe the phenomena of pressure waves through liquid–gas bubbles mix under the thermodynamic conditions. The stability property of the obtained solutions is tested to show the ability of our obtained solutions through the physical experiments. The novelty and advantage of the proposed method are illustrated by applying to this model. Some sketches are plotted to show more about the dynamical behavior of this model.


Introduction
Nowadays, many natural phenomena have been derived in nonlinear partial differential equations with an integer order. These models are included in various and distinct branches of science such as chemistry, physics, biology, engineering, economy, etc. However, using the integer order of these models is not sufficient where the nonlocal property does not appear in these formulas so that many models have been formulated in fractional nonlinear partial differential equations specially to discover that kind of property. Studying these models gives more novel properties of them specially by using the computational and numerical schemes. For using most of these schemes, one needs fractional operators to convert the fractional formulas to nonlinear ordinary differential equations with integer order such as Caputo, Caputo-Fabrizio definition, fractional Riemann-Liouville derivatives, conformable fractional derivative, and so on [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15]. These fractional opera-tors have been employed to investigate the exact and numerical solutions of many phenomena. These solutions have been obtained in explicit formulas by using different analytical schemes such as [16][17][18][19][20][21][22][23][24][25].
This method depends on a new auxiliary equation, which is equal to the auxiliary equation of the main future of the modified mathematical technique [36]. The auxiliary equation of the mK method is given by where δ, , χ , Q are arbitrary constants; whereas the auxiliary equation of the extended exponential-expansion function method is given by where β i (i = 1, 2, 3) are arbitrary constants. So Eqs. (1) and (2) . Using this technique leads to the equality of the mK auxiliary equation with many other analytical methods, but the mK method can obtain more solutions than most of them. This equivalence shows superiority, power, and productivity of the mK method. In this context, the mK method is employed to construct new formulas of solutions for the fractional nonlinear KS equation, which is given by [37][38][39][40][41][42][43]: where [S = S(x, t)] is the function that is used to describe the dynamical behavior of the nonlinear wave processes in a liquid containing gas bubbles. Additionally

Fractional ABR operator
The ABR fractional operator is given by [44][45][46][47][48] where G ϑ is the Mittag-Leffler function, defined by the following formula: and B(ϑ) is a normalization function. Thus For further properties of this fractional operator, you can see [44,49,50]. This leads to where c is an arbitrary constant. This wave transformation converts Eq. (3) to ODE. Integration of the obtained ODEs once with zero constant of the integration gives Calculating the homogeneous balance value in Eq. (6) yields N = 2. Thus, the general formula of solution according to the mK method is given by where a 0 , a 1 , a 2 , b 1 , b 2 are arbitrary constants. The rest of this article is arranged in the following order. In Sect. 2 we apply the mK method to the nonlinear fractional nonlinear (2 + 1)-BLMP equation. Moreover, some sketches are given to show more physical properties of both models. Section 4 discusses the stability property of the obtained solutions. Section 5 gives the conclusion of the whole research.

Abundant wave solutions of the fractional KS equation
Applying the mK method with its auxiliary equation and the suggested general solutions for the fractional KS equation leads to a system of algebraic equations. Using Mathematica 11.2 to find the values of the parameters in this system leads to the following: where (δ = 0, μ = 0, a 2 = 0) .

Stability
This section of our research paper investigates one of the basic properties of any model. It examines the stability property for the fractional nonlinear KS equation by using a Hamiltonian system. The momentum in the Hamiltonian system is given by the following formula: where is an arbitrary constant. Thus, the condition for stability is given in the following condition: where c, are arbitrary constants.
This means that this solution is unstable and, by applying the same steps to other obtained solutions, the stability property of each one of them can be determined.

Conclusion
This research has successfully applied the modified Khater method with a new fractional operator to the fractional nonlinear KS equation that is arising in the nonlinear wave processes in a liquid containing gas bubbles. This new operator is used to avoid the disadvantage of the other fractional operator. Distinct, solitary wave solutions have been obtained for this equation. For more illustrations of the dynamical behavior of this kind of fluid, some solutions have been sketched (Figs. 1, 2, 3) in three different formulas of each figure (two, three-dimensional, and contour plots).