# A Novel Analytical Technique to Obtain Kink Solutions for Higher Order Nonlinear Fractional Evolution Equations

Authors: Qazi Mahmood Ul Hassan; Jamshad Ahmad; Muhammad Shakeel
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

We use the fractional derivatives in Caputo’s sense to construct exact solutions to fractional fifth order nonlinear evolution equations. A generalized fractional complex transform is appropriately used to convert this equation to ordinary differential equation which subsequently resulted in a number of exact solutions.


---

## Body

## 1. Introduction

The concept of differentiation and integration to noninteger order is not new in any case. The notion of fractional calculus emerged when the ideas of classical calculus were proposed by Leibniz, who mentioned it in a letter to L’Hospital in 1695. The foundation of the earliest, more or less, systematic studies can be traced back to the beginning and middle of the 19th century by Liouville in 1832, Riemann in 1853, and Holmgren in 1864, although Euler in 1730, Lagrange in 1772, and others also made contributions. Recently, it has turned out those differential equations involving derivatives of noninteger [1, 2]. For example, the nonlinear oscillation of earthquakes can be modeled with fractional derivatives [3]. There has been some attempt to solve linear problems with multiple fractional derivatives [3, 4]. Not much work has been done on nonlinear problems and only a few numerical schemes have been proposed for solving nonlinear fractional differential equations [5, 6]. More recently, applications have included classes of nonlinear equation with multiorder fractional derivatives. The generalized fractional complex transform was applied in [7–13] to convert fractional order differential equation to ordinary differential equation. Finally, by using Exp-function method [14–25] we obtain generalized solitary solutions and periodic solutions. Recently the theory of local fractional integrals and derivatives [26–28] is one of useful tools to handle the fractal and continuously nondifferentiable functions. It is to be tinted that that c=d and p=q are the only relations that can be obtained by applying Exp-function method [29] to any nonlinear ordinary differential equation. Most scientific problems and phenomena in different fields of sciences and engineering occur nonlinearly. Except in a limited number of these problems are linear, this method has been effectively and accurately shown to solve a large class of nonlinear problems. The solution procedure of this method, with the aid of Maple, is of utter simplicity and this method can easily extend to other kinds of nonlinear evolution equations. In engineering and science, scientific phenomena give a variety of solutions that are characterized by distinct features. Traveling waves appear in many distinct physical structures in solitary wave theory [30] such as solitons, kinks, peakons, cuspons, compactons, and many others. Solitons are localized traveling waves which are asymptotically zero at large distances. In other words, solitons are localized wave packets with exponential wings or tails. Solitons are generated from robust balance between nonlinearity and dispersion. Solitons exhibit properties typically associated with particles. Kink waves [30, 31] are solitons that rise or descend from one asymptotic state to another and hence another type of traveling waves as in the case of the Burgers hierarchy. Peakons, that are peaked solitary wave solutions, are another type of travelling waves as in the case of Camassa-Holm equation. For peakons, the traveling wave solutions are smooth except for a peak at a corner of its crest. Peakons are the points at which spatial derivative changes sign so that peakons have a finite jump in 1st derivative of the solution. Cuspons are other forms of solitons where solution exhibits cusps at their crests. Unlike peakons where the derivatives at the peak differ only by a sign, the derivatives at the jump of a cuspon diverge. The compactons are solitons with compact spatial support such that each compacton is a soliton confined to a finite core or a soliton without exponential tails or wings. Other types of travelling waves arise in science such as negatons, positons, and complexitons. In this research, we use the Exp-function method along with generalized fractional complex transform to obtain new Kink waves’ solutions for [30–32].


## 2. Preliminaries and Notation [1, 2]

In this section, we give some basic definitions and properties of the fractional calculus theory [1, 2] which will be used further in this work. For the finite derivative in [a,b] we define the following fractional integral and derivatives.Definition 1.
A real functionf(x), x>0, is said to be in the space Cμ, μ∈R if there exists a real number (p>μ) such that f(x)=xpf1(x), where f1(x)=C(0,∞), and it is said to be in the space Cμmμ if fm∈Cμ, m∈N.Definition 2.
The Riemann-Liouville fractional integral operator of orderα≥0 of a function f∈Cμ, μ≥-1, is defined as
(1)Jα(x)=1Γ(α)∫0x(x-t)α-1f(t)dt,α>0,x>0,J0(x)=f(x).
Properties of the operator Ja can be found in [1]; we mention only the following.
Forf∈Cμ, μ≥-1, α,β≥0, and γ≥-1(2)JαJβf(x)=Jα+Bf(x),JαJBf(x)=JBJαf(x),Jαxγ=Γ(γ+1)Γ(α+γ+1)xα+γ.
The Riemann-Liouville derivative has certain disadvantages when trying to model real-world phenomena with fractional differential equations. Therefore, we will introduce a modified fractional differential operator proposed by M. Caputo in his work on the theory of viscoelasticity [2].Definition 3.
Form to be the smallest integer that exceeds α, the Caputo time fractional derivative operator of order α>0 defined as
(3)DaCtαf(t)=Jm-αDmf(t)=1Γ(m-α)∫0t(t-τ)m-α-1fm(t)dt,
for m-1<α≤1m, m∈N, t>0, f∈C-1m.


## 3. Chain Rule for Fractional Calculus and Fractional Complex Transform

In [7], the authors used the following chain rule ∂αu/∂tα=(∂u/∂s)(∂αs/∂tα) to convert a fractional differential equation with Jumarie’s modification of Riemann-Liouville derivative into its classical differential partner. In [10], the authors showed that this chain rule is invalid and showed following relation [8]:
(4)Dtau=σt′dudηDtaη,Dxau=σx′dudηDxaη.
To determine σs, consider a special case as follows:
(5)s=tα,u=sm
and we have
(6)∂αu∂tα=Γ(1+mα)tmα-αΓ(1+mα-α)=σ·∂u∂s=σmtmα-α.
Thus one can calculate σs as
(7)σs=Γ(1+mα)Γ(1+mα-α).
Other fractional indexes (σx′,σy′,σz′) can determine in similar way. Li and He [3, 7–9] proposed the following fractional complex transform for converting fractional differential equations into ordinary differential equations, so that all analytical methods for advanced calculus can be easily applied to fractional calculus:
(8)u(x,t)=u(η),η=kxβΓ(1+β)+ωtαΓ(1+α)+MxγΓ(1+γ),
where k, ω, and M are constants.


## 4. Exp-Function Method [33–36]

Consider the general nonlinear partial differential equation of fractional order:(9)P(u,ut,ux,uxx,uxxx,…,Dtαu,Dxαu,Dxxαu,…)=0,0<α≤1,
where Dtαu, Dxαu, and Dxxαu are the fractional derivative of u with respect to t, x, xx, respectively.Use(10)u(x,t)=u(η),η=kxβΓ(1+β)+ωtαΓ(1+α)+MxγΓ(1+γ),
where k, ω, and M are constants.Then (9) becomes
(11)Q(u,u′,u′′,u′′′,uiv)=0,
where the prime denotes derivative with respect to η. In accordance with Exp-function method, we assume that the wave solution can be expressed in the following form:
(12)u(η)=∑n-cdanexp⁡⁡[nη]∑m-pqbmexp⁡⁡[mη],
where p, q, c, and d are positive integers which are known to be further determined and an and bm are unknown constants. Equation (8) can be rewritten as
(13)u(η)=acexp⁡⁡(cη)+⋯+a-dexp⁡⁡(-dη)bpexp⁡⁡(pη)+⋯+b-qexp⁡⁡(-qη).
This equivalent formulation plays an important and fundamental for finding the analytic solution of problems. c and p can be determined by [29].


## 5. Solution Procedure

Consider the following new fifth order nonlinear(2+1)-dimensional evolution equations of fractional order:
(14)Dt3αu-(Dtαu)xxxx-(Dtαu)xx-4(ux(Dtαu)x)x=0.
Using (8) in (14) then it can be converted to an ordinary differential equation. Consider
(15)-ω3u⃛+ωk4u(v)+k2ωu⃛+4ωk3u˙u⃛=0,
where the prime denotes the derivative with respect to η. The solution of (15) can be expressed in form (13). To determine the value of c and p, by using [26],
(16)p=c,q=d.Case 1.
We can freely choose the values ofc and d, but we will illustrate that the final solution does not strongly depend upon the choice of values of c and d. For simplicity, we set p=c=1 and q=d=1 (15) reduces to
(17)u(η)=a1exp⁡[η]+a0+a-1exp⁡[-η]b1exp⁡[η]+a0+b-1exp⁡[-η].
Substituting (17) into (15), we have
(18)1A[c4a1exp⁡[4η]=c3exp⁡[3η]+c2exp⁡[2η]+c1exp⁡[η]+c0+c-1exp⁡[-η]+c-2exp⁡[-2η]+c-3exp⁡[-3η]+c-4exp⁡[-4η]]=0,
where A=(b1exp⁡(η)+b0+b-1exp⁡(-η))4 and ci are constants obtained by Maple software 16. Equating the coefficients of exp⁡⁡(nη) to be zero, we obtain
(19)c-4=0,c-3=0,c-2=0,c-1=0,c0=0,c1=0,c2=0,c3=0,c4=0.
For solution of (19) we have five solution sets satisfying the given (15).
1st Solution Set. Consider
(20)ω=-k2+1k,a-1=b-1(3kb0+a0)b0,a0=a0,a1=0,b-1=b-1,b0=b0,b1=0.b-1(3kb0+a0)b0
We, therefore, obtained the following generalized solitary solution u(x,t) of (14) (Figure 1):
(21)u(x,t)=b-1(3kb0+a0)e-kx-(σk2+1ktα/Γ(1+α))/b0+a0b-1e-kx-(σk2+1ktα/Γ(1+α))+b0.
2nd Solution Set. Consider
(22)ω=4k2+1k,a-1=b-1(6kb1+a1)b1,a0=a0,a1=a1,b-1=b-1,b0=0,b1=b1.
We, therefore, obtained the following generalized solitary solutionu(x,t) of (14) (Figure 2):
(23)u(x,t)=(b-1(6kb1+a1)e-kx+(σ4k2+1ktα/Γ(1+α))/b1+a1ekx-(σ4k2+1ktα/Γ(1+α)))×(b-1e-kx+(σ4k2+1ktα/Γ(1+α))+b1ekx-(σ4k2+1ktα/Γ(1+α)))-1.
3rd Solution Set. Consider
(24)ω=k2+1k,a-1=b-1(3kb0+a0)b0,a0=a0,a1=0,b-1=b-1,b0=b0,b1=0.
We, therefore, obtained the following generalized solitary solution u(x,t) of (14) (Figure 3):
(25)u(x,t)=b-1(3kb0+a0)e-kx+(σk2+1ktα/Γ(1+α))/b0+a0b-1e-kx+(σk2+1ktα/Γ(1+α))+b0.
4th Solution Set. Consider
(26)ω=k2+1k,a-1=191k2b14(9k2a0b0b13-9k2a1b02b12-3ka02b13+9ka0a1b0b12-6ka12b02b1-a02a1b12+2a0a12b0b1-a13b02),b1=b1,b0=b0,a0=a0,a1=a1,b-1=193ka0b0b12-3ka1b02b1-a02b12+2a0a1b0b1-a12b02k2b13.
We, therefore, obtained the following generalized solitary solutionu(x,t) of (14) (Figure 4):
(27)u(x,t)=(191k2b14(e-kx+(σk2+1ktα/Γ(1+α))(9k2a0b0b13-9k2a1b02b12-3ka02b13+9ka0a1b0b12-6ka12b02b1-a02a1b12+2a0a12b0b1-a13b02)×e-kx+(σk2+1ktα/Γ(1+α)))+a0+a1ekx-(σk2+1ktα/Γ(1+α))191k2b14)×(b1ekx-(σk2+1ktα/Γ(1+α))19(e-kx+(σk2+1ktα/Γ(1+α))(3ka0b0b12-3ka1b02b1-a02b12+2a0a1b0b1-a12b02)×e-kx+(σk2+1ktα/Γ(1+α)))(k2b13)-1+b0+b1ekx-(σk2+1ktα/Γ(1+α))19)-1.
5th Solution Set. Consider
(28)ω=k2+1k,a-1=b-1(6kb1+a1)b1,a0=0,a1=a1,b-1=b-1,b0=0,b1=b1.
We, therefore, obtained the following generalized solitary solution u(x,t) of (14) (Figure 5):
(29)u(x,t)=b-1e-kx+(σk2+1ktα/Γ(1+α))+a1ekx-(σk2+1ktα/Γ(1+α))b-1e-kx+(σk2+1ktα/Γ(1+α))+b1ekx-(σk2+1ktα/Γ(1+α)).Figure 1
Kink waves’ solutions of (14) for 1st solution set.Figure 2
Kink waves’ solutions of (14) for 2nd solution set.Figure 3
Kink waves’ solutions of (14) for 3rd solution set.Figure 4
Kink waves’ solutions of (14) for 4th solution set.Figure 5
Kink waves’ solutions of (14) for 5th solution set.Case 2.
Ifp=c=2 and q=d=1 then trial solution(14) reduces to
(30)u(η)=a2exp⁡⁡[2η]+a1exp⁡⁡[η]+a0+a-1exp⁡⁡[-η]b2exp⁡⁡[2η]+b1exp⁡⁡[η]+b0+b-1exp⁡⁡[-η].
Proceeding as before, we obtain the following.1st Solution. Consider
(31)a-1=a-1,a0=a-1b0b-1,a1=a-1b1b-1,a2=a-1b2b-1,b-1=b-1,b0=b0,b1=b1,b2=b2.
Hence we get the generalized solitary wave solution u(x,t) of (14) (Figure 6):
(32)u(x,t)=(a-1e-kx+(σωktα/Γ(1+α))+a-1b0b-1+a-1b0b-1ekx-(σωktα/Γ(1+α))+a-1b2b-1e2kx-2(σωktα/Γ(1+α)))×(b-1e-kx+(σωktα/Γ(1+α))+b0+b1ekx-(σωktα/Γ(1+α))+b2e2kx-2(σωktα/Γ(1+α)))-1.
In both cases, for different choices of c, p, d, and q we get the same soliton solutions which clearly illustrate that final solution does not strongly depend on these parameters.Figure 6
Kink waves’ solutions of (14) for 1st solution set of Case 2.


## 6. Conclusions

Exp-function method is applied to construct solitary solutions of the nonlinear new fifth order evolution equations of fractional orders. The reliability of proposed algorithm is fully supported by the computational work, the subsequent results, and graphical representations. It is observed that Exp-function method is very convenient to apply and is very useful for finding solutions of a wide class of nonlinear problems of fractional orders.



---

*Source: 213482-2014-07-01.xml*