# A Novel Method for Analytical Solutions of Fractional Partial Differential Equations

Authors: Mehmet Ali Akinlar; Muhammet Kurulay
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

A new solution technique for analytical solutions of fractional partial differential equations (FPDEs) is presented. The solutions are expressed as a finite sum of a vector type functional. By employing MAPLE software, it is shown that the solutions might be extended to an arbitrary degree which makes the present method not only different from the others in the literature but also quite efficient. The method is applied to special Bagley-Torvik and Diethelm fractional differential equations as well as a more general fractional differential equation.


---

## Body

## 1. Introduction

Fractional calculus is a significantly important and useful branch of mathematics having a broad range of applications at almost any branch of science. Techniques of fractional calculus have been employed at the modeling of many different phenomena in engineering, physics, and mathematics. Problems in fractional calculus are not only important but also quite challenging which usually involves hard mathematical solution techniques (see, e.g., [1]). Unfortunately, a general solution theory for almost each problem in this area has yet to be established. Each application venue has developed its own approaches and implementations. As a consequence, a single standard method for the problems in fractional calculus has not emerged. Therefore, finding reliable and efficient solution techniques along with fast implementation methods is a significantly important and active research area.In the literature, a number of methods have been developed for the numerical or analytical solutions for FPDEs. We can list some of these methods as follows: Adomian decomposition method [2], the collocation method [3], the fractional differential transform method [4], homotopy analysis method [5], homotopy perturbation method [6], and some other methods [7, 8] listed on the references of these papers. In this paper, we present a new method for the analytical solutions of FPDEs. The solutions are expressed as a finite sum of a vector type functional. By employing MAPLE software, it is shown that the solutions might be extended to an arbitrary degree which makes the present method not only different from the others in the literature but also quite efficient. The method is applied to special Bagley-Torvik and Diethelm fractional differential equations as well as a more general fractional differential equation.Now let us briefly review some significant concepts in fractional calculus. The fractional calculus is a name for the theory of integrals and derivatives of arbitrary order, which unifies and generalizes the notions of integer-order differentiation andn-fold integration. We have well-known definitions of a fractional derivative of order α>0 such as Riemann-Liouville, Grunwald-Letnikow, Caputo, and generalized functions approach [9]. The most commonly used definitions of fractional derivative belongs to Riemann-Liouville and Caputo. For the purpose of this paper, Caputo’s definition of fractional differentiation will be used, taking the advantage of Caputo’s approach that the initial conditions for fractional differential equations with Caputo’s derivatives take on the traditional form as for integer-order differential equations. We give some basic definitions and properties of the fractional calculus theory which were used through paper.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)J0vf(x)=1Γ(v)∫0x‍(x-t)v-1f(t)dt,v>0,J0f(x)=f(x).
It has the following properties.
Forf∈Cμ, μ≥-1, α,β≥0, and γ>1,(1)
JαJβf(x)=Jα+βf(x),(2)
JαJβf(x)=JβJαf(x),(3)
Jαxγ=(Γ(γ+1)/Γ(α+γ+1))xα+γ.The Riemann-Liouville fractional derivative is mostly used by mathematicians but this approach is not suitable for the physical problems of the real world since it requires the definition of fractional order initial conditions, which have no physically meaningful explanation yet. Caputo introduced an alternative definition, which has the advantage of defining integer order initial conditions for fractional order differential equations.Definition 3.
The fractional derivative off(x) in the Caputo sense is defined as
(2)D*vf(x)=Jam-vDmf(x)=1Γ(m-v)∫0x‍(x-t)m-v-1f(m)(t)dt,
for m-1<v<m, m∈N, x>0, f∈C-1m.Lemma 4.
Ifm-1<α<m, m∈N, and f∈Cμm, μ≥-1, then
(3)D*αJαf(x)=f(x),JαD*vf(x)=f(x)-∑k=0m-1fk(0+)xkk!,x>0.
The Caputo fractional derivative is considered here because it allows traditional initial and boundary conditions to be included in the formulation of the problem.Definition 5.
Form to be the smallest integer that exceeds α, the Caputo time-fractional derivative operator of order α>0 is defined as
(4)D*tαu(x,t)=∂αu(x,t)∂tα={1Γ(m-α)∫0t‍(t-ξ)m-α-1∂mu(x,ξ)∂ξmdξ,form-1<α<m,∂mu(x,t)∂tm,forα=m∈N,
and the space-fractional derivative operator of order β>0 is defined as
(5)D*xαu(x,t)=∂βu(x,t)∂xβ={1Γ(m-β)∫0x‍(x-θ)m-β-1∂mu(θ,t)∂θmdθ,form-1<β<m,∂mu(x,t)∂xm,forβ=m∈N.Another concept which plays a very significant role in the fractional calculus is the Gamma function. Next we briefly overview the definition and some important properties of Gamma function.Definition 6.
For0<n<∞, the integral (sometimes known as Euler integral) defined as
(6)Γ(n)=∫0∞‍(-log(t))(n-1)dt
is said to be Gamma function. From (6), it is clear that for n>0,
(7)Γ(n)=∫0∞‍t(n-1)e-tdt.
Again for n>0, the function defined as
(8)Γ(n)=2∫0∞‍t(2n-1)e-t2dt
is known as the Gamma function. Detailed information about the significant features of the Gamma function might be obtained from any well-written differential equation book.
Organization of the paper is in the following way. Firstly, we overview basic concepts of fractional derivative. Because we employ Caputo sense derivative, we describe it in detail. Secondly, we introduce a new method for analytical solutions of FPDEs. In the third section, we illustrate three computational examples as the application of the present method and complete the paper with a discussion section.


## 2. A Novel Method for Analytical Solutions of FPDEs

Let us consider the FPDE given as(9)Dnx(t)=f(t,x(t),Dα1x(t),Dα2x(t),…,Dαmx(t)),t>0,
where αi is a real number for every i and n is a natural number. Initial conditions for this equation are given by
(10)Djx(0)=aj,j=0,1,…,n-1,0<α1<α2<⋯<αm<n.
Writing
(11)yi(t)=Dαix(t),
we can express (9) as
(12)Dnx(t)=f(t,y0(t),y1(t),…,ym(t)),t>0.
Now let us assume that the solution of (12) is given by
(13)y=y0+y1xα1+y2xα2+⋯+exαn,
where y is an m dimensional vector function and e is also a vector type function having the same dimension with y. If the solution (13) is plugged into (12) and the higher order derivatives are ignored, a linear equation expressed in terms of e is obtained. As a consequence, e is plugged in (13) and solution is obtained. If the same process is repeated for the higher order derivatives, series solution of (9) (or (13)) from an arbitrary degree is obtained.Now in order to express the solution in a new power series form, let us explain the application of this method to the power series. When (13) is plugged into (9), we can write
(14)f(x)=f0+f1xα1+f2xα2+⋯+(fn+p1e1+p2e2+⋯+pmem)xαn,
where p1,p2,…,pm are constant real numbers, e1,e2,…,em are the bases of e, and m is the dimension of the vector e. It is clear that using (13) it is possible to define that
(15)yi=yi,0+yi,1xα1+yi,2xα2+⋯+eixαn.
If y and its derivatives are written at (14), one obtains
(16)fi(x)=fi,n+pi,1e1+⋯+(pi,mem)xαn-j+O(xαn-j+1),
where fi is an element of f(y,y′,x) appearing at (9) and j is an index which is only 1 less than the degree of the highest derivative given in the initial problem. By solving (16), one obtains ei, for each i=1,2,…,m. By writing these ei values at (15), for each i, the polynomial yi having the degree of n is obtained. This gives us power series solution of (9) from an arbitrary order. Finally let us suppose that step size of x is h. If this value is plugged in y and y′, y and y′ are obtained at the point of x=x0+h. Repeating the same process, one obtains solution of (9).In the next section, we illustrate the application of this new and novel method to the analytical solutions of some FPDEs.


## 3. Computational Applications

Example 7.
In the first example, we consider a special case of Bagley-Torvik equation(17)d2ydx2+d3/2ydx3/2+y=1+x
with the initial conditions
(18)y(0)=1,y′(0)=1.
Now bearing in mind the aforementioned solution procedure, let us assume that the solution of (17) is given by
(19)y=y0+y1x1/α+y2x2/α+y3x3/α+y4x4/α.
Forα=2, from (13) we can write that
(20)y=1+x+ex2.
Taking the derivatives of y, we get
(21)Γ(3)Γ(1.5)ex0.5+1+x+ex2=1+x
which implies that e=0. In a similar way, one obtains
(22)y=1+x+ex5/α.
Again by taking the derivatives of y, one obtains that e=0. In the next step, we have
(23)y=1+x+y6x6/α
which implies that e=0. Repeating this procedure, one obtains the solution of (17) as
(24)y=1+x.Example 8.
In this example, we consider the initial value problem studied by Diethelm and given as(25)D0.5y(x)=y(x)+x2+2Γ(2.5)x1.5
with initial conditions
(26)y(0)=0.
Let us suppose that the solution of (25) is given by y=ex. When we write this solution in the original equation (25), we get e=0. In the same way, if we consider the solution
(27)y=ex2,
we obtain e=1. Again by considering the solution
(28)y=x2+ex3,
we obtain e=0. Repeating the procedure in this way, we obtain the solution of the initial value problem (25) as
(29)y=x2
which is also exact solution of the problem.
In the last example, we consider a more general example to illustrate the application of the novel method.Example 9.
Let us consider the equation(30)a(t)D2x(t)+b(t)Dα2x(t)+c(t)Dx(t)+e(t)Dα1x(t)+k(t)x(t)=f(t),0<α1<α2<1,
where
(31)f(t)=-a-b(t)Γ(3-α2)t2-α2-c(t)t-e(t)Γ(3-α1)t2-α1+k(t)(2-12t2),
with the initial conditions
(32)x(0)=2,x′(0)=0.
Following the steps of aforementioned solution algorithm, one can obtain the solution of (30) as
(33)x(t)=2-12t2,
which is also the exact solution of this problem.


## 4. Conclusion

A new technique for the analytical solutions of FPDEs has been successfully developed in this paper. By employing MAPLE software, it is shown that the solutions might be extended to an arbitrary degree which makes the present method not only different from the others in the literature but also quite efficient. The method is applied to special Bagley-Torvik and Diethelm fractional partial differential equations as well as a more general fractional differential equation. Experimental results prove that the present method is a useful and highly efficient technique.



---

*Source: 195708-2013-07-16.xml*