This work is licensed under the Creative Commons Attribution 4.0 International License.
Introduction
In 1695, since L’ Hospital asked the question, what might be a derivative order 1/2. Many researchers tried to find a definition of fractional derivative after this question. Most of the works focused on an integral form of the fractional derivative. Although most famous approaches are the Caputo definition and Riemann–Liouville definition, these two definitions have some drawbacks. For example,
Riemann–Liouville definition does not satisfy Dα1 = 0 when α is not a natural number.
Caputo definition assumes that the function is differentiable.
Both definitions do not satisfy the derivative of the product of two functions.
Both definitions do not satisfy the derivative of the quotient of two functions.
Both of them do not satisfy the chain rule.
Both of them do not satisfy the index rule.
Recently, Khalil et al. [5] introduced a new definition of the integral and conformable fractional derivatives. So using the conformable fractional derivative, we can overcome the aforementioned drawbacks of the existing definitions. Now let us give the definition and some properties of conformable fractional derivative and integral.
Definition 1.1
Let f : [0, ∞) → is a function of α-th order “conformable functional derivate”, which is defined by
{T_\alpha }(f)(t) = \mathop {\lim }\limits_{\varepsilon \to 0} \frac{{f(t + \varepsilon {t^{1 - \alpha }}) - f(t)}}{\varepsilon },
for all t > 0, α ∈ (0, 1).
Definition 1.2
Starting from a ≥ 0, the conformable integral of a function f is defined in [5] as
I_\alpha ^a(f)(s) = \mathop {\mathop {\int }\limits_ }\limits_a \frac{{f(t)}}{{{t^{1 - \alpha }}}}dt.
Similarly, the definition of conformable fractional partial derivative is given by Atangana et al. [2]
.
Definition 1.3
Let f be a function with n variables such as x1,...,xn, the conformable partial derivative of f order α ∈ (0, 1] in xi is defined as follows
\frac{{{d^\alpha }}}{{dx_i^\alpha }}f({x_1},...,{x_n}) = \mathop {\lim }\limits_{\varepsilon \to 0} \frac{{f({x_1},...,{x_{i - 1}},{x_i} + \varepsilon x_i^{1 - \alpha },...,{x_n}) - f({x_1},...,{x_n})}}{\varepsilon }.
Recently, many studies have been performed in various fields such as applied mathematics, physics and engineering related to fractional calculations [12,13, 14, 15, 16]. Phenomena related to nonlinear partial differential equations (NLPDEs) have emerged in many areas such as physics, mechanics and chemistry to investigate the exact solutions for NLPDEs. In recent years, there are a lot of workings with NLPDEs. For example, Whitham [9] studied variational methods and applications on water waves. Sirendaoreji et al. [7] used the auxiliary equation method for solving NLPDEs. Zhang and Xia [11] studied a generalised new auxiliary equation method and its applications to NLPDEs. Tasbozan et al. [8]
studied the Sine-Gordon expansion method to obtain the analytical results for Drinfeld-Sokolov-Wilson system. Yomba [10] discussed the exact results for the nonlinear Klein-Gordon equation and generalised nonlinear Camassa–Holm equation using a generalised auxiliary equation method. Eslami and Mirzazadeh [4]
used the first integral method to obtain the exact solutions of the nonlinear Schrödinger equation.
The properties of this new definition [5] are given below.
Theorem 1.4
Let f, g functions and α ∈ (0, 1] are α-differentiable at a point t > 0, then
If f is differentiable, then{T_\alpha }(f)(t) = {t^{1 - \alpha }}\frac{{df(t)}}{{dt}}.
Fractional Bad–Good Modified Boussinesq Equations
In the 1870s, the Boussinesq equation [3]
was first introduced by Joseph Boussinesq. This equation corresponds to the shallow water wave model that arises narrow channels and coastlines. This equation has been used to describe the motions of long waves in shallow water under gravity forces. Subsequently, it was applied to many other areas of mathematical physics dealing with wave phenomena. In this paper, we consider bad and good modified Boussinesq equations where the fractional terms are in a conformable sense.
Bad Modified Boussinesq Equation:D_t^{(2\alpha )}u - D_x^2u - D_x^4u - 3D_x^2\left( {{u^2}} \right) + 3{D_x}\left( {{u^2}{D_x}u} \right) = 0,Good Modified Boussinesq Equation:D_t^{(2\alpha )}u - D_x^2u + D_x^4u - 3D_x^2\left( {{u^2}} \right) + 3{D_x}\left( {{u^2}{D_x}u} \right) = 0.
In this study, we use the auxiliary equation method to obtain a solution set for the above-mentioned equations using a conformable fractional derivative.
Description of Auxiliary Equation Method
Auxiliary equation method [7] has been used to obtain exact solutions for NLPDEs. Auxiliary equation method can be applied in NLPDEs if the considered equation consists only even-order or only odd-order partial derivative terms. This method depends on the following differential equation
{\left( {\frac{{dz}}{{d\xi }}} \right)^2} = a{z^2}(\xi ) + b{z^3}(\xi ) + c{z^4}(\xi ),
First mentioned by Sirendaoreji [7]. Clearly, we can express the solution procedure of the method step by step as follows.
The general form of a nonlinear conformable fractional differential equation can be regarded as
P\left( {u,D_t^\alpha u,{D_x}u,D_t^{(2\alpha )}u,D_x^2u,...} \right) = 0,
where the arguments and subscripts of polynomial P show partial derivatives, and D_t^{(2\alpha )} means two times conformable derivative of the function u(x, t).
Using the wave transformation with chain rule [1]
u(x,t) = U(\xi ),\xi = x - w\frac{{{t^\alpha }}}{\alpha },
where w denotes the velocity of the wave. With the aid of this transformation, fractional derivatives of Eq. (5) can be rewritten as:
\frac{{{\partial ^\alpha }(.)}}{{\partial {t^\alpha }}} = \frac{{d(.)}}{{d\xi }},\frac{{\partial (.)}}{{\partial x}} = - w\frac{{d(.)}}{{d\xi }}, \ldots .
Using the transformation given in Eq. (6) inside Eq. (5), we obtain the following ordinary differential equation
G(U,{U^'},{U^{''}},{U^{'''}},...) = 0,
where the derivatives are respected to ξ.
Now, consider U(ξ) is a sum of serial such as
U(\xi ) = \sum\limits_{i = 0}^n {a_i}{z^i}(\xi ),
where z(ξ) is the solution of the nonlinear differential equation (4), a, b, c, w, ai are the real constants and n is a positive integer to be determined by the balancing procedure [6]
Balancing the linear and nonlinear terms of the highest order in the ordinary differential equation, Eq. (8) gives the result for the parameter n. Then, we place Eq. (9) into the ordinary differential equation Eq. (8). After this procedure, we get an equation consisting of the powers of z(ξ). All coefficients of z(ξ) are equated to 0 in the final equation. This procedure arouses the system of algebraic equations including a, b, c, w, ai. Solving this system, with respect to these parameters and using the exact solutions of Eq. (4) in Table 1, gives the analytical solutions. Some solutions of Eq. (4) are given in Table 1.
Using the wave transform (6) and integrating both equations twice, the equations turn into the ordinary differential equation as follows.
- U'' + \left( {{w^2} - 1} \right)U - 3{U^2} + {U^3} = 0,U'' + \left( {{w^2} - 1} \right)U - 3{U^2} + {U^3} = 0,
where prime denotes the derivative of the functions with respect to ξ.
Now using the balancing procedure in Eqs. (12) and (13) yields n = 1. Thus, the unknown function U(ξ) can be considered as
U(\xi ) = {a_0} + {a_1}z(\xi ).
Placing Eq. (14) into Eqs (12) and (13) and using Eq. (4) led to an algebraic equation with respect to z(ξ). Equating all the coefficients of same powers of to 0 arouses an algebraic equation system. Solving this system gives the following solution sets.
In this study, an auxiliary equation method is used to obtain the new exact solutions of Bad and Good Boussinesq equations. The study indicates that the auxiliary equation method is direct, effective and understandable and can be used for solving other NLPDEs in mathematical physics. Moreover, using the conformable fractional derivative, one can obtain analytical solutions of the NLPDEs which cannot be solved by Caputo and Riemann–Liouville definitions. In addition, many transactions were made with the program code of the Mathematica program.