# New Analytical Solutions of Conformable Time Fractional Bad and Good Modified Boussinesq Equations

Hülya Durur 1 , Orkun Tasbozan 2  and Ali Kurt 3
• 1 Department of Computer Engineering, Faculty of Engineering, Ardahan, Turkey
• 2 Department of Mathematics, Faculty of Science and Art, Hatay, Turkey
• 3 Department of Mathematics, Faculty of Science and Art, Denizli, Turkey
Hülya Durur
• Corresponding author
• Department of Computer Engineering, Faculty of Engineering, Ardahan University, Ardahan, Turkey
• Email
• Search for other articles:
, Orkun Tasbozan
• Department of Mathematics, Faculty of Science and Art, Mustafa Kemal University, Hatay, Turkey
• Search for other articles:
and Ali Kurt
• Department of Mathematics, Faculty of Science and Art, Pamukkale University, Denizli, Turkey
• Search for other articles:

## Abstract

The main purpose of this article is to obtain the new solutions of fractional bad and good modified Boussinesq equations with the aid of auxiliary equation method, which can be considered as a model of shallow water waves. By using the conformable wave transform and chain rule, nonlinear fractional partial differential equations are converted into nonlinear ordinary differential equations. This is an important impact because both Caputo definition and Riemann–Liouville definition do not satisfy the chain rule. By using conformable fractional derivatives, reliable solutions can be achieved for conformable fractional partial differential equations.

## 1 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α(f)(t)=limɛ→0f(t+ɛt1−α)−f(t)ɛ,$
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αa(f)(s)=∫asf(t)t1−α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

$dαdxiαf(x1,...,xn)=limɛ→0f(x1,...,xi−1,xi+ɛxi1−α,...,xn)−f(x1,...,xn)ɛ.$

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

1. Tα(mf + ng) = mTα(f) + nTα(g) for all m,n ∈.
2. Tα(tp) = ptp−α for all p.
3. Tα(f.g) = fTα(g) + gTα(f).
4. $Tα(fg)=gTα(f)−fTα(g)g2.$.
5. Tα(c) = 0 for all constant functions f(t) = c.
6. If f is differentiable, then$Tα(f)(t)=t1−αdf(t)dt$.

## 2 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:

$Dt(2α)u−Dx2u−Dx4u−3Dx2(u2)+3Dx(u2Dxu)=0,$
Good Modified Boussinesq Equation:
$Dt(2α)u−Dx2u+Dx4u−3Dx2(u2)+3Dx(u2Dxu)=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.

## 3 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

$(dzdξ)2=az2(ξ)+bz3(ξ)+cz4(ξ),$
First mentioned by Sirendaoreji [7]. Clearly, we can express the solution procedure of the method step by step as follows.
1. Step 1.The general form of a nonlinear conformable fractional differential equation can be regarded as
$P(u,Dtαu,Dxu,Dt(2α)u,Dx2u,...)=0,$
where the arguments and subscripts of polynomial P show partial derivatives, and $Dt(2α)$ means two times conformable derivative of the function u(x, t).
2. Step 2. Using the wave transformation with chain rule [1]
$u(x,t)=U(ξ),ξ=x−wtαα,$
where w denotes the velocity of the wave. With the aid of this transformation, fractional derivatives of Eq. (5) can be rewritten as:
$∂α(.)∂tα=d(.)dξ,∂(.)∂x=−wd(.)dξ,….$
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 ξ.
3. Step 3.Now, consider U(ξ) is a sum of serial such as
$U(ξ)=∑i=0naizi(ξ),$
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]
4. Step 4.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.
Table 1

Solutions of Eq. (4) with Δ = b2 − 4ac and ε = ±1

Noz(ξ)
1$−absech2(a2ξ)b2−ac(1+ɛtanh(a2ξ))2$a > 0
2$abcsch2(a2ξ)b2−ac(1+ɛcoth(a2ξ))2$a > 0
3$2asech(aξ)ɛΔ−bsech(aξ)$a > 0, Δ > 0
4$2asec(−aξ)ɛΔ−bsec(−aξ)$a < 0, Δ > 0
5$2acsch(aξ)ɛ−Δ−bcsch(aξ)$a > 0, Δ < 0
6$2acsc(−aξ)ɛΔ−bcsc(−aξ)$a < 0, Δ > 0
7$−asech2(a2ξ)b+2ɛactanh(a2ξ)$a > 0, c > 0
8$−asec2(−a2ξ)b+2ɛ−actan(−a2ξ)$a < 0, c > 0
9$acsch2(a2ξ)b+2ɛaccoth(a2ξ)$a > 0, c > 0
10$−acsc2(−a2ξ)b+2ɛ−accot(−a2ξ)$a < 0, c > 0
11$−ab[1+ɛtanh(a2ξ)]$a > 0, Δ = 0
12$−ab[1+ɛcoth(a2ξ)]$a > 0, Δ = 0
13$4aeɛaξ(eɛaξ−b)2−4ac$a > 0
14$±4aeɛaξ1−4ace2ɛaξ$

## 4 Implementation of the Method

We consider the time-fractional Bad–Good modified Boussinesq equations

$Dt(2α)u−Dx2u−Dx4u−3Dx(u2Dxu)=0,$
$Dt(2α)u−Dx2u+Dx4u−3Dx(u2Dxu)=0,$

Using the wave transform (6) and integrating both equations twice, the equations turn into the ordinary differential equation as follows.

$−U′′+(w2−1)U−3U2+U3=0,$
$U′′+(w2−1)U−3U2+U3=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(ξ)=a0+a1z(ξ).$

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.

• Solutions for Bad modified Boussinesq equation:SET:Using obtained constants from Set, with the aid of Table 1, the new wave solutions of Bad modified Boussinesq equation (10) can be given as follows
$u1,2(x,t)=−4(1−w2)sech[12−1+w2(x−tαwα)]28+(1−w2)(1±tanh[12−1+w2(x−tαwα)])2,u3,4(x,t)=−4(−1+w2)csch[12−1+w2(x−tαwα)]28−(−1+w2)(1±coth[12−1+w2(x−tαwα)])2u5,6(x,t)=2(−1+w2)sech[−1+w2(x−tαwα)]±4−2(−1+w2)+2sech[−1+w2(x−tαwα)],u7,8(x,t)=2(−1+w2)sec[1−w2(x−tαwα)]±4−2(−1+w2)+2sec[1−w2(x−tαwα)],u9,10(x,t)=2(−1+w2)csch[−1+w2(x−tαwα)]±−4+2(−1+w2)+2csch[−1+w2(x−tαwα)],u11,12(x,t)=2(−1+w2)csc[1−w2(x−tαwα)]±4−2(−1+w2)+2csc[1−w2(x−tαwα)],u13,14(x,t)=(1−w2)sech[12−1+w2(x−tαwα)]2−2±2(−1+w2)tanh[12−1+w2(x−tαwα)],u15,16(x,t)=(1−w2)sec[121−w2(x−tαwα)]2−2±2(1−w2)tan[121−w2(x−tαwα)],u17,18(x,t)=(−1+w2)csch[12−1+w2(x−tαwα)]2−2±2(−1+w2)coth[12−1+w2(x−tαwα)],u19,20(x,t)=(1−w2)csc[121−w2(x−tαwα)]2−2±2(1−w2)cot[121−w2(x−tαwα)],u21,22(x,t)=4a1e±−1+w2(x−tαwα)(−1+w2)(2a1+e±−1+w2(x−tαwα))2−2a12(−1+w2).$
• Solutions for Good modified Boussinesq equation:SET:
$a0=0, a1=1−w2, b=2a1, c=−a122.$
The new wave solutions of Good modified Boussinesq equation (11) with the aid of Table 1 can be given as follows
$u1,2(x,t)=4(−1+w2)sech[121−w2(x−tαwα)]28+(1−w2)(1±tanh[121−w2(x−tαwα)])2,u3,4(x,t)=4(1−w2)csch[121−w2(x−tαwα)]28+(1−w2)(1±coth[121−w2(x−tαwα)])2,u5,6(x,t)=2(1−w2)sech[1−w2(x−tαwα)]±4+2(1−w2)−2sech[1−w2(x−tαwα)],u7,8(x,t)=2(1−w2)sec[−1+w2(x−tαwα)]±4+2(1−w2)−2sec[−1+w2(x−tαwα)],u9,10(x,t)=2(1−w2)csch[1−w2(x−tαwα)]±−4−2(1−w2)−2csch[1−w2(x−tαwα)],u11,12(x,t)=2(1−w2)csc[−1+w2(x−tαwα)]±4+2(1−w2)−2csc[−1+w2(x−tαwα)],u13,14(x,t)=(−1+w2)sech[121−w2(x−tαwα)]22±2w2−1tanh[121−w2(x−tαwα)],u15,16(x,t)=(−1+w2)sec[12−1+w2(x−tαwα)]22±21−w2tan[12−1+w2(x−tαwα)],u17,18(x,t)=(1−w2)csch[121−w2(x−tαwα)]22±2w2−1coth[121−w2(x−tαwα)],u19,20(x,t)=(−1+w2)csc[12−1+w2(x−tαwα)]22±21−w2cot[12−1+w2(x−tαwα)],u21,22(x,t)=4a1e±1−w2(x−tαwα)(1−w2)(−2a1+e±1−w2(x−tαwα))2+2a12(1−w2).$

## 5 Conclusion

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.

## References

• [1]

Abdeljawad, T., (2015), On conformable fractional calculus, Journal of Computational and Applied Mathematics, 279, 57–66.

• [2]

Atangana, A., Baleanu, D., Alsaedi, A., (2015), New properties of conformable derivative. Open Mathematics, 13, 1–10.

• [3]

Boussinesq, J. (1872), Thorie des ondes et des remous qui se propagent le long d’un canal rectangulaire horizontal, en communiquant au liquide contenu dans ce canal des vitesses sensiblement pareilles de la surface au fond, Journal de Mathmatiques Pures et Appliques, 55–108.

• [4]

Eslami, M., Mirzazadeh, M., (2013), Topological 1-soliton solution of nonlinear Schrödinger equation with dual-power law nonlinearity in nonlinear optical fibers, The European Physical Journal Plus, 128(11), 140.

• [5]

Khalil, R., Al Horani, M., Yousef, A., Sababheh, M., (2014), A new definition of fractional derivative, Journal of Computational and Applied Mathematics, 264, 65–70.

• [6]

Malfliet, W., (1992), Solitary wave solutions of nonlinear wave equations, American Journal of Physics, 60(7), 650–654.

• [7]

Sirendaoreji, Jiong, S., (2003), Auxiliary equation method for solving nonlinear partial differential equations, Physics Letters A, 309, 387–396.

• [8]

Tasbozan, O., Şenol, M., Kurt, A., Özkan, O., (2018), New solutions of fractional Drinfeld-Sokolov-Wilson system in shallow water waves, Ocean Engineering, 161, 62–68.

• [9]

Whitham, G. B., (1967), Variational Methods and Applications to Water Waves, In Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences (Vol. 299, pp. 6–25).

• [10]

Yomba, E., (2008), A generalized auxiliary equation method and its application to nonlinear Klein-Gordon and generalized nonlinear Camassa-Holm equations, Physics Letters A, 372(7), 1048–1060.

• [11]

Zhang, S., Xia, T., (2007), A generalized new auxiliary equation method and its applications to nonlinear partial differential equations, Physics Letters A, 363, 356–360.

• [12]

Yokuş, A. (2018), Comparison of Caputo and conformable derivatives for time-fractional Korteweg–de Vries equation via the finite difference method. International Journal of Modern Physics B, 32(29), 1850365

• [13]

Yokus, A., Kaya, D., (2017), Numerical and exact solutions for time fractional Burgers’ equation. Journal of Nonlinear Sciences and Applications, 10(7), 3419–3428.

• [14]

Yokuş, A., Kaya, D., (2015), Conservation laws and a new expansion method for sixth order Boussinesq equation. In AIP Conference Proceedings (Vol. 1676, No. 1, p. 020062). AIP Publishing.

• [15]

Yokus, A., (2018), Numerical solution for space and time fractional order Burger type equation. Alexandria engineering journal, 57(3), 2085–2091.

• [16]

Yokus, A., Bulut, H., (2018), On the numerical investigations to the Cahn-Allen equation by using finite difference method. An International Journal of Optimization and Control: Theories & Applications (IJOCTA), 9(1), 18–23.

If the inline PDF is not rendering correctly, you can download the PDF file here.

• [1]

Abdeljawad, T., (2015), On conformable fractional calculus, Journal of Computational and Applied Mathematics, 279, 57–66.

• [2]

Atangana, A., Baleanu, D., Alsaedi, A., (2015), New properties of conformable derivative. Open Mathematics, 13, 1–10.

• [3]

Boussinesq, J. (1872), Thorie des ondes et des remous qui se propagent le long d’un canal rectangulaire horizontal, en communiquant au liquide contenu dans ce canal des vitesses sensiblement pareilles de la surface au fond, Journal de Mathmatiques Pures et Appliques, 55–108.

• [4]

Eslami, M., Mirzazadeh, M., (2013), Topological 1-soliton solution of nonlinear Schrödinger equation with dual-power law nonlinearity in nonlinear optical fibers, The European Physical Journal Plus, 128(11), 140.

• [5]

Khalil, R., Al Horani, M., Yousef, A., Sababheh, M., (2014), A new definition of fractional derivative, Journal of Computational and Applied Mathematics, 264, 65–70.

• [6]

Malfliet, W., (1992), Solitary wave solutions of nonlinear wave equations, American Journal of Physics, 60(7), 650–654.

• [7]

Sirendaoreji, Jiong, S., (2003), Auxiliary equation method for solving nonlinear partial differential equations, Physics Letters A, 309, 387–396.

• [8]

Tasbozan, O., Şenol, M., Kurt, A., Özkan, O., (2018), New solutions of fractional Drinfeld-Sokolov-Wilson system in shallow water waves, Ocean Engineering, 161, 62–68.

• [9]

Whitham, G. B., (1967), Variational Methods and Applications to Water Waves, In Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences (Vol. 299, pp. 6–25).

• [10]

Yomba, E., (2008), A generalized auxiliary equation method and its application to nonlinear Klein-Gordon and generalized nonlinear Camassa-Holm equations, Physics Letters A, 372(7), 1048–1060.

• [11]

Zhang, S., Xia, T., (2007), A generalized new auxiliary equation method and its applications to nonlinear partial differential equations, Physics Letters A, 363, 356–360.

• [12]

Yokuş, A. (2018), Comparison of Caputo and conformable derivatives for time-fractional Korteweg–de Vries equation via the finite difference method. International Journal of Modern Physics B, 32(29), 1850365

• [13]

Yokus, A., Kaya, D., (2017), Numerical and exact solutions for time fractional Burgers’ equation. Journal of Nonlinear Sciences and Applications, 10(7), 3419–3428.

• [14]

Yokuş, A., Kaya, D., (2015), Conservation laws and a new expansion method for sixth order Boussinesq equation. In AIP Conference Proceedings (Vol. 1676, No. 1, p. 020062). AIP Publishing.

• [15]

Yokus, A., (2018), Numerical solution for space and time fractional order Burger type equation. Alexandria engineering journal, 57(3), 2085–2091.

• [16]

Yokus, A., Bulut, H., (2018), On the numerical investigations to the Cahn-Allen equation by using finite difference method. An International Journal of Optimization and Control: Theories & Applications (IJOCTA), 9(1), 18–23.

OPEN ACCESS