AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 587179 10.1155/2013/587179 587179 Research Article Fractional Subequation Method for Cahn-Hilliard and Klein-Gordon Equations Jafari Hossein 1,2 Tajadodi Haleh 1 Kadkhoda Nematollah 1 0000-0002-0286-7244 Baleanu Dumitru 3, 4, 5 Ahmad Bashir 1 Department of Mathematics University of Mazandaran, P.O. Box 47416-93797 Babolsar Iran umz.ac.ir 2 Department of Mathematical Sciences International Institute for Symmetry Analysis and Mathematical Modelling North-West University, Mafikeng Campus, Private Bag X2046 Mmabatho 2735 South Africa 3 Department of Mathematics and Computer Sciences Faculty of Art and Sciences Çankaya University Ankara Turkey cankaya.edu.tr 4 Department of Chemical and Materials Engineering Faculty of Engineering King Abdulaziz University P.O. Box 80204, Jeddah 21589 Saudi Arabia kau.edu.sa 5 Institute for Space Sciences Măgurele, P.O. Box R 76900, Bucharest Romania spacescience.ro 2013 3 02 2013 2013 10 12 2012 06 01 2013 2013 Copyright © 2013 Hossein Jafari et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The fractional subequation method is applied to solve Cahn-Hilliard and Klein-Gordon equations of fractional order. The accuracy and efficiency of the scheme are discussed for these illustrative examples.

1. Introduction

Fractional calculus deals with fractional integrals and derivatives of any order . Numbers of very interesting and novel applications of fractional partial differential equations (FPDEs) in physics, chemistry, engineering, finance, biology, hydrology, signal processing, viscoelastic materials, fractional variational principles, and so forth, developed mainly in the last few decades , have led recently to an intensive effort to find accurate and stable numerical methods that are also straightforward to be implemented.

Also, the exact solutions of most of the FPDEs cannot be found easily; thus analytical and numerical methods must be used. Some of the numerical methods for solving fractional differential equations (FDE) and FPDEs were discussed in (see [7, 1623] and the references therein).

By taking into account the results from , a new direct method titled fractional subequation method to search for explicit solutions of FPDEs was proposed . We notice that the method relies on the homogeneous balance principle , Jumarie’s modified Riemann-Liouville derivative [27, 28], and the symbolic computation. With the help of this method, some exact solutions of nonlinear time fractional biological population model as well as the (4+1)-dimensional space-time fractional Fokas equation were reported . Recently, the improved fractional subequation method was proposed, and it was used to solve the following two FPDEs in fluid mechanics .

In this paper, we suggest the fractional subequation method and utilize this method to solve the following two FPDEs.

The space-time fractional Cahn-Hilliard equation in the form (1)Dtαu-γDxαu-6u(Dxαu)2-(3u2-1)Dx2αu+Dx4αu=0,

where 0<α1 and u are the functions of (x,t). For the case corresponding to α=1, this equation is related with a number of interesting physical phenomena like the spinodal decomposition, phase separation, and phase ordering dynamics. On the other hand it becomes important in material sciences [30, 31]. However we notice that this equation is very difficult to be solved and several articles investigated it (see, e.g.,  and the references therein).

The nonlinear fractional Klein-Gordon equation  with quadratic nonlinearity reads as (2)Dtt2αu-Dxx2αu+γu-βu2=0,γ,β0.

We notice that the nonlinear fractional Klein-Gordon equation describes many types of nonlinearities. On the other hand the Klein-Gordon equation plays a significant role in several real world applications, for example, the solid state physics, nonlinear optics, and quantum field theory.

The paper suggests a fractional subequation method to find the exact analytical solutions of nonlinear fractional partial differential equations with the Jumarie’s modified Riemann-Liouville derivative of order α which is defined as  (3)Dxαf(x)={1Γ(1-α)0x(x-ξ)-α-1[f(ξ)-f(0)],α<0,1Γ(1-α)ddx0x(x-ξ)-α[f(ξ)-f(0)],0<α<1,[f(α-n)(x)](n),nα<n+1,n1. As pointed out by Kolwankar and Gangal , even though the variable t is taking all real positive values the actual evolution takes place only for values of t in the fractal set C. We take χ(t)=1 which is a flag function. We conclude that, from the viewpoint of the Kolwankar-Gangal’s local fractional derivative, the parameter α is the fractal dimension of time. Thus, the approximate solution is generated by some distribute function defined over the fractal sets in some closed interval [0,1]. They are continuous but not differentiable functions with respect to t.

The organization of the manuscript is as follows. In Section 2, we briefly explain the fractional subequation method for solving fractional partial differential equations. In Section 3, we extend the application of the proposed method to two nonlinear equations. Finally, Section 4 is devoted to our conclusions.

2. The Method

The fundamental ingredients of the fractional subequation method for solving fractional partial differential equations are described in . The starting point is to consider a given nonlinear fractional partial differential equation in u(x,t)(4)p(u,ux,ut,Dtαu,Dxαu,)=0,0<α<1, where Dtαu and Dxαu are Jumarie’s modified Riemann-Liouville derivatives of u, u=u(x,t) is an unknown function, and P is a polynomial in u and its various partial derivatives, in which the highest order derivatives and nonlinear terms are involved.

To specify u explicitly, we use in this paper the proposal in four basic steps proposed in [29, 35]; namely, we reduce, by using the traveling wave transformation, the given nonlinear FPDE to a nonlinear fractional differential equation (FDE). After that we assume that the reduced equation obtained previously admits the following solution (5)u(ξ)=i=0naiφi, where ai (i=0,1,,n-1,n) are constants to be found, n denotes a positive integer determined by balancing the highest order derivatives with the highest nonlinear terms in (4) or the modified one (see  for more details), and the new variable φ=φ(ξ) fulfilling the fractional Riccati equation: (6)Dξαφ=σ+φ2,0<α1. The next step is to substitute (5) along with (6) into the modified version of the equation and to use the properties of Jumarie’s modified Riemann-Liouville derivative, in order to get a polynomial in φ(ξ). Requesting all coefficients of φk (k=0,1,2,) to be zero, we end up to a set of overdetermined nonlinear algebraic equations for c,  k,  ai  (i=0,1,,n-1,n).

Finally, assuming that c,  k,  ai  (i=0,1,,n-1,n) are obtained by solving the algebraic equations in the previous step, and substituting these constants and the solutions of (6) into (5), we get the explicit solutions of (4).

3. Main Results

In this section, we apply the method presented in Section 2 for solving the FPDEs (1) and (2), respectively.

Example 1.

We consider the space-time fractional Cahn-Hilliard equation as (7)Dtαu-γDxαu-6u(Dxαu)2-(3u2-1)Dx2αu+Dx4αu=0.

Making use of the travelling wave transformation (8)u=u(ξ),ξ=kx+ct. Equation (7) is reduced into a nonlinear FDE easy to solve, namely, (9)cαDξαu-γkαDξαu-6u(kαDξαu)2-(3u2-1)k2αDξ2αu+k4αDξ4αu=0. Next we suppose that (9) has a solution in the form given below (10)u=i=0naiφi, where φ obeys the subequation (6).

By balancing the highest order derivative terms and nonlinear terms in (9), gives the value of  n=1, we substitute (10), along with (6), into (9), and then setting the coefficients of φj  (j=0,1,,5) to zero, we finally end up with a system of algebraic equations, namely, (11)a1σcα-6a0a12σ2k2α-a1γσkα=0,-6a13σ2k2α-16a1σ2k4α-6a02a1σk2α+2a1σk2α=0,a1cα-24a0a12σk2α-a1γkα=0,-18a13σk2α-40a1σk4α-6a02a1k2α+2a1k2α=0,-18a0a12k2α=0,-12a13k2α-24a1k4α=0. Solving the set of algebraic equations yields (12)a0=0,a1=±i2k2α, where k2α=1/2σ,   cα=kα and σ denotes an arbitrary constant.

By using (8)–(12) after some tedious calculations, the exact solutions of (7), namely, generalized hyperbolic function solutions (see  for their definitions) and generalized trigonometric function solutions are obtained as (13)u={±i2k2α(-σtanhα(-σ(kx+ct))),σ<0,±i2k2α(-σcothα(-σ(kx+ct))),σ<0,±2k2α(σtanα(σ(kx+ct))),σ>0,±i2k2α(σcotα(σ(kx+ct))),σ>0. We stress on the fact that when α1 these obtained exact solutions give the ones of the standard form equation of the space-time fractional Cahn-Hilliard equation (7).

Example 2.

The next step is to investigate the fractional nonlinear Klein-Gordon equation in the following form: (14)Dtt2αu-Dxx2αu+γu-βu2=0.

To solve (14), we perform the traveling wave transformation (15)u=u(ξ),ξ=kx+ct; therefore (14) is reduced to the following nonlinear fractional ODE, namely, (16)c2αDξ2αu-k2αD2αu+γu-βu2=0. Next, we assume that (16) admits a solution in the form (17)u=i=0naiφi. At this stage we apply the same technique as in the case of the previous example. Namely, by balancing the highest order derivative terms and nonlinear terms in (16), then substituting (17), with  n=2, with (6) into (16), we finally obtain the corresponding system of algebraic equations as (18)-a02β+a0γ+2a2σ2c2α-2a2σ2k2α=0,a1γ-2a0a1β+2a1σc2α-2a1σk2α=0,-a12β+a2γ-2a0a2β+8a2σc2α-8a2σk2α=0,-2a1a2β+2a1c2α-2a1k2α=0,-a22β+6a2c2α-6a2k2α=0. After using the Mathematica to solve (18) the following solutions are reported: (19)a0=γ+8σc2α-8σk2α2β,a1=0,a2=-6(k2α-c2α)β, where σ denotes an arbitrary constant. Finally, from (15)–(19) we obtain the following generalized hyperbolic function solutions, generalized trigonometric function solutions, and the rational solution of (14) as (20)u={γ+8σc2α-8σk2α2β+6σ(k2α-c2α)β×(tanhα2(-σξ)),σ<0,γ+8σc2α-8σk2α2β+6σ(k2α-c2α)β×(cothα2(-σξ)),σ<0,γ+8σc2α-8σk2α2β-6σ(k2α-c2α)β×(tanα(σξ)),σ>0,γ+8σc2α-8σk2α2β-6σ(k2α-c2α)β×(cotα(σξ)),σ>0,γ2β-6(k2α-c2α)Γ2(1+α)β(ξα+ω)2,σ=0, where ξ=kx+ct.

As α1 (20) the results obtained above become the ones of (14).

4. Conclusions

In this paper, a fractional subequation method is used to construct the exact analytical solutions of the space-time fractional Cahn-Hilliard (1) and the fractional nonlinear Klein-Gordon equation (2). These solutions include the generalized hyperbolic function solutions, generalized trigonometric function solutions, and rational function solutions, which may be very useful to further understand the mechanisms of the complicated nonlinear physical phenomena and FPDEs. Also, this method help us to find all exact solutions of the Fan subequations involving all possible parameters, it is concise and efficient. Mathematica has been used for computations and programming in this work.