Scholarly article on topic 'Formulation and solution to time-fractional generalized Korteweg-de Vries equation via variational methods'

Formulation and solution to time-fractional generalized Korteweg-de Vries equation via variational methods Academic research paper on "Mathematics"

Share paper
Academic journal
Adv Diff Equ
OECD Field of science

Academic research paper on topic "Formulation and solution to time-fractional generalized Korteweg-de Vries equation via variational methods"

0 Advances in Difference Equations

a SpringerOpen Journal


Open Access

Formulation and solution to time-fractional generalized Korteweg-de Vries equation via variational methods

Youwei Zhang*

Correspondence: Department of Mathematics, Hexi University, Beihuan Road 87, Zhangye, Gansu 734000, China


This paper presents the formulation of the time-fractional generalized Korteweg-de Vries equation (KdV) using the Euler-Lagrange variational technique in the Riemann-Liouville derivative sense and derives an approximate solitary wave solution. Our results show that He's variational-iteration method is a very efficient technique for finding the solution of the proposed equation and extend the existing results.

Keywords: Riemann-Liouville fractional differential operator; Euler-Lagrange equation; Riesz fractional derivative; generalized KdV equations; He's variational-iteration method; solitary wave

£ Springer

1 Introduction

The KdV equation has been used to describe a wide range of physics phenomena of the evolution and interaction of nonlinear waves. It was derived from the propagation of dispersive shallow water waves and is widely used in fluid dynamics, aerodynamics, continuum mechanics, as a model for shock wave formation, solitons, turbulence, boundary layer behavior, mass transport, and the solution representing the water's free surface over a flat bed [1-3]. Camassa and Holm [4] put forward the derivation of solution as a model for dispersive shallow water waves and discovered that it is a formally integrable dimensional Hamiltonian system, and that its solitary waves are solitons. Most of the classical mechanics techniques have been used in studies of conservative systems, but most of the processes observed in the physical real world are nonconservative.

During the past three decades or so, fractional calculus has obtained considerable popularity and importance as generalization of integer-order evolution equations, and one has used it to model some meaningful things, such as fractional calculus in model price volatility in finance [5, 6], in hydrology to model fast spreading of pollutants [7]; the most common hydrologic application of fractional calculus is the generation of fractional Brownian motion as a representation of aquifer material with a long-range correlation structure [8, 9]. The fractional differential equation is used to model the particle motions in a heterogeneous environment and long particle jumps in anomalous diffusion in physics [10,11]. For other exact descriptions of the applications of engineering, mechanics, mathematics, etc., we can refer to [12-16]. If the Lagrangian of conservative system is constructed using fractional derivatives, the resulting equations of motion can be non-conservative. Therefore, in many cases, real physical processes could be modeled in a reli-

©2014 Zhang; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the originalworkis properly cited.

able manner using fractional-order differential equations rather than integer-order equations [17]. Based on the stochastic embedding theory, Cresson [18] defined the fractional embedding of differential operators and provided a fractional Euler-Lagrange equation for Lagrangian systems, then investigated a fractional Noether theorem and a fractional Hamiltonian formulation of fractional Lagrangian systems. Herzallah and Baleanu [19] presented the necessary and sufficient optimality conditions for the Euler-Lagrange fractional equations of fractional variational problems and made an effort in determining in which spaces the functional must exist. Malinowska [20] proposed the Euler-Lagrange equations for fractional variational problems with multiple integrals and proved the fractional Noether-type theorem for conservative and nonconservative generalized physical systems. Riewe [21, 22] formulated a version of the Euler-Lagrange equation for problems of calculus of variation with fractional derivatives. Wu and Baleanu [23] developed some new variational-iteration formulas to find approximate solutions of fractional differential equations and determined the Lagrange multiplier in a more accurate way. For generalized fractional Euler-Lagrange equations and a fractional-order Van der Pol-like oscillator, we can refer to the works by Odzijewicz [24, 25] and Attari etal. [26], respectively. For other the known results refer to Baleanu etal. [27] andlnokuti et al. [28]. In view of the fact that most physical phenomena may be considered as nonconservative, they can be described using fractional-order differential equations. Recently, several methods have been used to solve nonlinear fractional evolution equations using techniques of nonlinear analysis, such as the Adomian decomposition method [29], the homotopy analysis method [30,31], the homotopy perturbation method [32,33] andtheLaguerre spectral method [34-36]. It was mentioned that the variational-iteration method has been used successfully to solve different types of integer and fractional nonlinear evolution equations.

In the present paper, He's variational-iteration method [37-40] is applied to solve the time-fractional generalized KdV equation

qDI u(x, t) + aup(x, t)ux(x, t) + buxxx(x, t) = 0,

where a, b are constants, u(x, t) is a field variable, the subscripts denote the partial differentiation of the function u(x, t) with respect to the parameter x and t. x e Q (Q c R)isaspace coordinate in the propagation direction of the field and t e T (= [0, t0] (t0 > 0)) is the time, which occurs in different contexts in mathematical physics. a, b are constant coefficients not equal to zero. The dissipative uxxx term provides damping at small scales, and the nonlinear term upux (p >0) (which has the same form as that in the KdV or 1-dimensional Navier-Stokes equations) stabilizes by transferring energy between large and small scales. For p = 1, we can refer to the known results of the time-fractional KdV equation: formulation and solution using variational methods [41,42]. Forp > 0, p =1, there is little material on the formulation and solution to time-fractional KdV equation. Thus the present paper considers the formulation and solution to a time-fractional KdV equation as the index of the nonlinear term satisfies p > 0, p =1. RDa denotes the Riesz fractional derivative. Making use of the variational-iteration method, this work is devoted to a formulation of a time-fractional generalized KdV equation and derives an approximate solitary wave solution.

This paper is organized as follows: Section 2 states some background material from fractional calculus. Section 3 presents the principle of He's variational-iteration method. Section 4 is devoted to a description of the formulation of the time-fractional generalized

KdV equation using the Euler-Lagrange variational technique and to derive an approximate solitary wave solution. Section 5 presents an analysis for the obtained graphs and figures and discusses the present work.

2 Preliminaries

We recall the necessary definitions for the fractional calculus (see [43-45]) which are used throughout the remaining sections of this paper.

Definition 1 A real multivariable functionf (x, t), t >0 is said to be in the space CY, y e R with respect to t if there exists a real number r (> y), such that f (x, t) = tf1(x, t), where f1(x, t) e C(Q x T). Obviously, Cy c Q if 5 < y.

Definition 2 The left-hand side Riemann-Liouville fractional integral of a functionf e CY (y > -1) is defined by

0Ifx, t) = ^ " >0t e T,

0l0f (x, t)=f (x, t).

Definition 3 The Riemann-Liouville fractional derivatives of the order n -1 < a < n of a functionf e Cy (y > -1) are defined as

1 dn ft f (x, t )

1 dn f ^ M= r^ of- L

r(n - a) dtn J0 (t - t)a+1-n 1 dn ft0 f (x, t )

n«f( a 1 9 n f f (x, t ) . , ^

tDa f (x, t) = —-- — ---—■— dx, t e T.

' r(n - a) dtn Jt (t - t)a+1-n

Lemma 1 The integration formula of Riemann-Liouville fractional derivativefor the order 0 < a <1

if (x, t)0Dg(x, t) dt = / g(x, t)tD0f (x, t) dt

is valid under the assumptions thatf, g e C(Q x T) and that for arbitrary x e Q, tD^0f, oD^g exist at every point t e T and are continuous in t.

Definition 4 The Riesz fractional integral of the order n -1 < a < n of a function f e CY (y > -1) is defined as

Riaf (x, t) =1 (0iff (x, t)+tiaf (x, t)) = —^ f0 it - t ry(x, t) dt,

2s 0 2T(a) J0

where 0If and tI^0 are, respectively, the left- and right-hand side Riemann-Liouville fractional integral operators.

Definition 5 The Riesz fractional derivative of the order n -1 <a < n of a functionf e CY (y > -1) is defined by

RDf (x, t) = 1 (0Df (x, t) + (-1)^/ (x, t))

1 dn t0

2T(n -a) dtn

f° it - T in-"-1f (x, T) dt, 0

where oD'à and are, respectively, the left- and right-hand side Riemann-Liouville fractional differential operators.

Lemma 2 Let a > 0 and ¡3 > 0 be such that n -1 < a < n, m -1 < P < m and a + ¡3 < n, and letf e L1(^ x T) and 0im-af e ACm(Q x T). Then we have the following index rule:

m t-a-i

Ri^a CRrfift., A\ Rr\a+Pf(„A X^Rr^P -if

RdhRdj (x, t)) = R0Dat+Pf(x, t) -j2R0Df-f (x, t)|t=o

i=1 r(1- a - i)

Remark 1 One can express the Riesz fractional differential operator RDa-1 of the order 0 < a < 1 as the Riesz fractional integral operator RlT~a, i.e.

R0Da-1f (x, t)=U-af (x, t), t e T.

3 Variational-iteration method

The variational-iteration method provides an effective procedure for explicit and solitary wave solutions of a wide and general class of differential systems representing real physical problems. Moreover, the variational-iteration method can overcome the foregoing restrictions and limitations of approximate techniques so that it provides us with a possibility to analyze strongly nonlinear evolution equations. Therefore, we extend this method to solve the time-fractional KdV equation. The basic features of the variational-iteration method are outlined as follows.

We consider a nonlinear evolution equation that consists of a linear part Lu(x, t), a nonlinear part Nu(x, t), and a free term g(x, t) represented as

Lu(x, t) + Nu(x, t) = g(x, t). (1)

According to the variational-iteration method, the (n + 1)th approximation solution of (1) can be read using an iteration correction functional as

un+1(x,t) = un(x,t)+ / k(t)(Lu(x,t)+Nu(x,t)-g(x,t))dt, (2)

where X(t) is a Lagrangian multiplier and u(x, t) is considered as a restricted variation function, i.e., su(x, t) = 0. Extremizing the variation of the correction functional (2) leads to the Lagrangian multiplier X(t). The initial iteration u0(x, t) can be used as the initial value u(x, 0). As n tends to infinity, the iteration leads to the solitary wave solution of (1), i.e.

u(x,t) = lim un(x,t).

4 Time-fractional generalized KdV equation

The generalized KdV equation in (1 +1) dimensions is given by

ut (x, t) + aup(x, t)ux(x, t) + buxxx(x, t) = 0, (3)

wherep > 0, a, b are constants, u(x, t) is a field variable, x e ^ is a space coordinate in the propagation direction of the field and t e T is the time. Employing a potential function

v(x, t) on the field variable, setting u(x, t) = vx(x, t), yields the potential equation of the generalized KdV equation (3) in the form

vxt (x, t) + avp(x, t)vxx(x, t) + bvxxxx(x, t) = 0. (4)

The Lagrangian of this generalized KdV equation (3) can be defined using the semi-inverse method [46, 47] as follows. The functional of the potential equation (4) can be represented as

J(v) = dx (v(x, t) (C1 vxt(x, t) + C2avp(x, t)vxx(x, t) + C3bvxxxx(x, t))) dt, (5)

with ci (i = 1,2,3) being unknown constants to be determined later. Integrating (5) by parts and taking vt |q = vx|q = vx|t = 0 yield

J(v)= / dx [-civx(x, t)vt(x, t)- -2a-vp+2(x, t)

IQ JT\ p + 1

- c3bvx(x, t)vxxx(x, t) dt. (6)

The constants ci (i = 1,2,..., 6) can be determined taking the variation of the functional (6) to make it optimal. By applying the variation of the functional, integrating each term by parts, and making use of the variation optimum condition of the functional J(v), it yields the following expression:

-2c1vxt(x, t) - c2a(p + 2)vp(x, t)vxx(x, t) - 2c3bv*xxx(x, t) = 0. (7)

We notice that the obtained result (7) is equivalent to (4), so one sees that the constants ci (i = 1,2,..., 6) are, respectively,

c1 = --, c2 = 7, c3 = -. 2 p + 2 2

In addition, the functional expression given by (6) obtains directly the Lagrangian form of the generalized KdV equation,

L(vt, vx, vxxx) = 1 vx(x, t)vt (x, t) + ---—yx+2(x, t)

2 (p + 1)(p + 2)

+ 2 vx(x, t)vxxx(x, t).

Similarly, the Lagrangian of the time-fractional version of the generalized KdV equation could be read as

F (0Datv, vx, vxxx) = 1 vx(x, t)0Datv(x, t) + ----vp+2(x, t)

2 (p + 1)(p + 2)

+ 2vx(x, t)vxxx(x, t), a e ]0,1]. (8)

Then the functional of the time-fractional generalized KdV equation will take the form of the expression

J(v)= / dx F(0Dat v, vx, VxxXjdt, jq jt

where the time-fractional Lagrangian F(0Da v, vx, vxxx) is given by (8). Following Agrawal's method [48-50], the variation of the functional (9) with respect to v(x, t) leads to

С Ci д F дF

SJ(v)= dx ( -—tS(oDat v(x, t)) + — 5vx(x, t) Jq Jt\д oDa v dvx

+ --s vxxx (x, t) ) dt. (0)

d Vxxx /

By Lemma 1, upon integrating the right-hand side of (10), one has

s/(v)=L dxL(M mv) - s( £) - U 3)svdt,

noting that Sv|t = Sv|s = Svx\Q = 0.

Obviously, optimizing the variation of the functional/(v), i.e., S/(v) = 0, yields the Euler-Lagrange equation for time-fractional generalized KdV equation as in the following expression:

dt (JL- ) - A ( dL ) - 4 (JL_) = 0. (11)

V d 0-a Vj 8x\d Vx) 3x3\ d Vxxx)

Substituting the Lagrangian of the time-fractional generalized KdV equation (8) into the Euler-Lagrange formula (11) one obtains

tDT^2Vx(x,t)^ -0D^2Vx(x,t)^ -avp(x,t)vxx(x,t)-bvxxxx(x,t) = 0.

Once again, substituting the potential function vx(x, t) for u(x, t) yields the time-fractional generalized KdV equation for the state function u(x, t),

1 (0Dau(x, t) - tDTu(x, t)) + aup(x, t)ux(x, t) + buxxx(x, t) = 0. (12)

According to the Riesz fractional derivative RDau(x, t), the time-fractional generalized KdV equation represented in (12) can be written as

RDa u(x, t) + aup(x, t)ux(x, t) + buxxx(x, t) = 0. (13)

Acting from the left-hand side by the Riesz fractional operator RD1-a on (13) leads to

d ta-2

— u(x, t)-QDa-1u(x, t) |t=0 —:-r

d t 0 t 0 r(a -1)

+ 0 D1fa(aup (x, t)—u(x, t) + b u(x, t)^ = 0, (14)

\ dx dx3 /

from Lemma 2. In view of the variational-iteration method, combining with (14), the (n + l)th approximate solution of (13) can be read using the iteration correction functional as

un+i(x, t) = un(x, t)+ I X(x) Jo

— Un(x, T)~lDar-1Un(x, T)|t=0 -7T

dT I (a -1)

+ oD1T-a ( aUpn(x, t)—Un(x, t) + bUn(x, t )) dT, n > 0, (15)

dx dx3 /

where the function Un(x, t) is considered as a restricted variation function, i.e., &Un(x, t) = 0. The extremum of the variation of (15) subject to the restricted variation function straightforwardly yields

Sun+1(x, t) = Sun(x, t)+ I X(x)S—un(x, t) dx

= Sun(x, t) +X(x)Sun(x, x)|t=t - I —X(x)Sun(x, x) dx

This expression reduces to the following stationary conditions: d

—X(x ) = 0, 1 + X(x ) = 0, dx

which are converted to the Lagrangian multiplier at X(x) = -1. Therefore, the correction functional (15) takes the following form:

un+1(x, t) = un(x, t) -

x, t)- /

— un(x, T)-oIt1-"un(x, T)|t=0"^--T

dx I (a - 1)

d d3 + °D1T-a I aupn(x, t) — un(x, t) + b—run(x, x) dx dx3

dx, n > 0, (16)

and since a -1 < 0, the fractional derivative operator RDa-1 reduces to the fractional integral operator Rl1-a by Remark 1.

The right-hand side Riemann-Liouville fractional derivative is interpreted as a future state of the process in physics. For this reason, the right-derivative is usually neglected in applications, when the present state of the process does not depend on the results of the future development, and so the right-derivative is used equal to zero in the following calculations. The zeroth order solitary wave solution can be taken as the initial value of the state variable, which is taken in this case as

2 ( P \

U0(x, t) = u(x, 0) = ksechp I —— (x + no) I, \2v/b /

where k = ((p+1)|p+2))p, n0 is a constant.

Substituting this zeroth order approximation solitary wave solution into (16) and using Definition 5 lead to the first order approximate solitary wave solution

u1(x,t) = ksechp( (x + n0) 11 1 + ^-tanh( —^(x + n0))

\2Vb J\ Vbr(a + 1) \2Vb JJ

Substituting the first order approximate solitary wave solution into (16), using Definition 5 then leads to the second order approximate solitary wave solution in the following form:

u2(x, t) = ksechp I —— (x + n0) \2v/b

r(a + 1) k

akp+1 2+2/ p / \ (

sechp | —— (x + n0) tanh

2Vb J \

. 4b p

(x + n0)

+ 4= sechp . _

4b \24b

(x + no) I tanh

3kp 24b

sech p

x tanhi -P

(x + n0)

(x + -

kp2 24b

(x + n0)

sech p

(x + n0)

r(2a + 1)

akp+1 2+2/ p —-—sech p —— (x + n0) b

x( I -3-p+. tanh^ -4 2 2 \2^b

24b (x + n0)

3kp (p + 2) sechp /(x + n0)) tanh2 8b \2vb

3kp2(p + 2) sechp2+2(—4—

8b V^vb

3kp(p + 2)2

(x + n0)

sechp | —4= (x + n0)) tanh2

(x + n0)

(x + n0)

kp3 2 / p

+ SechV24(x + n0)

k(p + 2)3 2/ p \ .(v

1 -sechp' 1 ^ . ^

3k(p + 2)2 8b

(x + no) ) tanh

(x + n0)

sech p

(x + n0) tanh

(x + n0)

kp2(p + p ,.....t ,2/ p

, ■sechp 4b

r(2a + 1)t3a

r(3a + 1)(r(a + 1))2

' apkp+1 . 2b4b

(x + n0) tanh

(x + n0)

2+2 sech p

x tanh

(x + n0) p -

_P 24b

(x + n0) (x + n0)

Making use of Definition 5 and the Maple package or Mathematics, we obtain u3(x, t), u4(x, t) and so on, and substituting the n - 1 order approximate solitary wave solution into (16) leads to the n order approximate solitary wave solution. As n tends to infinity, the iteration leads to the solitary wave solution of the time-fractional generalized KdV equation

2 / p u(x, t) = k sechp (

(x -1 + n0) .

5 Discussion

The target of the present work is the effect of the fractional-order derivative on the structure and propagation of the resulting solitary waves obtained from time-fractional generalized KdV equations. We derive the Lagrangian of the generalized KdV equations by the semi-inverse method, and then we take a similar form for the Lagrangian to the time-fractional generalized KdV equations. Using the Euler-Lagrange variational technique, we continue our calculations until the third order iteration. Meanwhile, our approximate calculations are carried out on the solution of the time-fractional generalized KdV equations taking into account the values of the coefficients and some meaningful values namely, a = 1, j, 2 and jr. The solitary wave solution of time-fractional generalized KdV equations are obtained. In addition, a 3-dimensional representation of the solution u(x, t) for the time-fractional generalized KdV equations with space x and time t for different values of the order a is presented, respectively. In Figure 1, the solution u is still a single soliton wave solution for all values of the order a. It shows that the balancing scenario between nonlinearity and dispersion is still valid. Figure 2 presents the change of amplitude and width of the soliton due to the variation of the order a, 2- and 3-dimensional graphs depicted the behavior of the solution u(x, t) at time t =1 corresponding to different values of the order a. This behavior indicates that the increasing of the value a is uniform in both the height and the width of the solitary wave solution. That is, the order a can be

(c)«-l/2 (d) a=l>5

Figure 1 The function u as a 3-dimensions graph for different order a when p = 3, a =10, b = 1, no =0.

used to modify the shape of the solitary wave without change of the nonlinearity and the dispersion effects in the medium. Figure 3 is devoted to the study of the expression of the relation between the amplitude of the soliton and the fractional order at different time values. These figures show that at the same time, the increasing of the fractional a increases the amplitude of the solitary wave to some value of a.

Competing interests

The author declares that they have no competing interests. Author's contributions

The author wrote this paper carefully, gave a rigorous derivation process, and read and approved the finalmanuscript.

Received: 7 November 2013 Accepted: 13 January 2014 Published: 14 Feb 2014


1. Camassa, R, Holm, D, Hyman, J: A new integrable shallow water equation. Adv. Appl. Mech. 31,1-33 (1994)

2. Johnson, RS: Camassa-Holm, Korteweg-de Vries and related models for waterwaves. J. Fluid Mech. 455,63-82 (2002)

3. Fokas, A, Fuchssteiner, B: Symplectic structures, their Backlund transformation and hereditary symmetries. Physica D 4,47-66(1981)

4. Camassa, R, Holm, D: An integrable shallow water equation with peaked solutions. Phys. Rev. Lett. 71,1661-1664 (1993)

5. Gorenflo, R, Mainardi, F, Scalas, E, Raberto, M: Fractionalcalculus and continuous-time finance, III. The diffusion limit. In: MathematicalFinance (Konstanz, 2000). Trends Math., pp. 171-180 (2001)

6. Sabatelli, L, Keating, S, Dudley, J, Richmond, P: Waiting time distributions in financialmarkets. Eur. Phys. J. B 27, 273-275 (2002)

7. Schumer, R, Benson, DA, Meerschaert, MM, Wheatcraft, SW: Eulerian derivation of the fractionaladvection-dispersion equation. J. Contam. Hydrol. 48,69-88 (2001)

8. Benson, DA, Meerschaert, MM, Revielle, J: Fractionalcalculus in hydrologic modeling: a numericalperspective. Adv. Water Resour. 51, 479-497 (2013)

9. Roop, JP: Computationalaspects of FEM approximation of fractionaladvection dispersion equations on bounded domains in R2. J. Comput. Appl. Math. 193, 243-268 (2006)

10. Hilfer, R: Applications of FractionalCalculus in Physics. World Scientific, Singapore (2000)

11. Meerschaert, MM, Benson, DA, Baeumer, B: Multidimensionaladvection and fractionaldispersion. Phys. Rev. E 59, 5026-5028(1999)

12. Lundstrom, B, Higgs, M, Spain, W, Fairhall, A: Fractionaldifferentiation by neocorticalpyramidalneurons. Nat. Neurosci. 11, 1335-1342(2008)

13. Metzler, R, Klafter, J: Boundary value problems for fractional diffusion equations. Phys. A, Stat. Mech. Appl. 278, 107-125 (2000)

14. Rossikhin, YA, Shitikova, MV: Application of fractionalderivatives to the analysis of damped vibrations of viscoelastic single mass systems. Acta Mech. 120,109-125 (1997)

15. Schumer, R, Benson, DA, Meerschaert, MM, Baeumer, B: Multiscaling fractionaladvection-dispersion equations and their solutions. Water Resour. Res. 39,1022-1032 (2003)

16. Yuste, SB, Acedo, L, Lindenberg, K: Reaction front in an A + B ^ C reaction subdiffusion process. Phys. Rev. E 69, 036126(2004)

17. Tavazoei, MS, Haeri, M: Describing function based methods for predicting chaos in a class of fractionalorder differentialequations. Nonlinear Dyn. 57, 363-373 (2009)

18. Cresson, J: Fractionalembedding of differentialoperators and Lagrangian systems. J. Math. Phys. 48, 033504 (2007)

19. Herzallah, MAE, Baleanu, D: FractionalEuler-Lagrange equations revisited. Nonlinear Dyn. 69, 977-982 (2012)

20. Malinowska, AB: A formulation of the fractionalNoether-type theorem for multidimensionalLagrangians. Appl. Math. Lett. 25, 1941-1946(2012)

21. Riewe, F: Nonconservative Lagrangian and Hamiltonian mechanics. Phys. Rev. E 53,1890-1899 (1996)

22. Riewe, F: Mechanics with fractionalderivatives. Phys. Rev. E 55, 3581-3592 (1997)

23. Wu, GC, Baleanu, D: Variationaliteration method for the Burgers' flow with fractionalderivatives - new Lagrange multipliers. Appl. Math. Model. 37,6183-6190 (2013)

24. Odzijewicz, T, Malinowska, AB, Torres, DFM: Generalized fractionalcalculus with applications to the calculus of variations. Comput. Math. Appl. 64, 3351-3366 (2012)

25. Odzijewicz, T, Malinowska, AB, Torres, DFM: Fractionalvariationalcalculus with classicaland combined Caputo derivatives. Nonlinear Anal. TMA 75,1507-1515 (2012)

26. Attari, M, Haeri, M, Tavazoei, MS: Analysis of a fractionalorder Van der Pol-like oscillator via describing function method. Nonlinear Dyn. 61, 265-274 (2010)

27. Baleanu, D, Muslih, SI: Lagrangian formulation of classicalfields within Riemann-Liouville fractionalderivatives. Phys. Scr. 72, 119-123(2005)

28. Inokuti, M, Sekine, H, Mura, T: Generaluse of the Lagrange multiplier in non-linear mathematicalphysics. In: Nemat-Nasser, S (ed.) VariationalMethod in the Mechanics of Solids. Pergamon, Oxford (1978)

29. Saha Ray, S, Bera, R: An approximate solution of a nonlinear fractionaldifferentialequation by Adomian decomposition method. Appl. Math. Comput. 167, 561-571 (2005)

30. Cang, J, Tan, Y, Xu, H, Liao, S: Series solutions of nonlinear fractionalRiccati differentialequations. Chaos Solitons Fractals 40,1-9 (2009)

31. Liao, S: The proposed homotopy analysis technique for the solution of nonlinear problems. PhD thesis, Shanghai Jiao Tong University (1992)

32. Sweilam, NH, Khader, MM, Al-Bar, RF: Numericalstudies for a multi-order fractionaldifferentialequation. Phys. Lett. A 371 , 26-33 (2007)

33. Golmankhaneh, AK, Golmankhaneh, AK, Baleanu, D: Homotopy perturbation method for solving a system of Schrodinger-Korteweg-de Vries equation. Rom. Rep. Phys. 63,609-623 (2011)

34. Baleanu, D, Bhrawy, AH, Taha, TM: Two efficient generalized Laguerre spectralalgorithms for fractionalinitialvalue problems. Abstr. Appl. Anal. 2013, Article ID 546502 (2013)

35. Baleanu, D, Bhrawy, AH, Taha, TM: A modified generalized Laguerre spectralmethod for fractionaldifferential equations on the half line. Abstr. Appl. Anal. 2013, Article ID 413529 (2013)

36. Bhrawy, AH, Baleanu, D: A spectralLegendre-Gauss-Lobatto collocation method for a space-fractionaladvection diffusion equations with variable coefficients. Rep. Math. Phys. 72, 219-233 (2013)

37. He, J: A new approach to nonlinear partial differentialequations. Commun. Nonlinear Sci. Numer. Simul. 2, 230-235 (1997)

38. He, J: Approximate analyticalsolution for seepage flow with fractionalderivatives in porous media. Comput. Methods Appl. Mech. Eng. 167,57-68(1998)

39. Molliq, RY, Noorani, MSM, Hashim, I: Variationaliteration method for fractionalheat- and wave-like equations. Nonlinear Anal., RealWorld Appl. 10,1854-1869 (2009)

40. Momani, S, Odibat, Z, Alawnah, A: Variationaliteration method for solving the space- and time-fractionalKdV equation. Numer. Methods PartialDiffer. Equ. 24, 261-271 (2008)

41. El-Wakil, S, Abulwafa, E, Zahran, M, Mahmoud, A: Time-fractionalKdV equation: formulation and solution using variationalmethods. Nonlinear Dyn. 65, 55-63 (2011)

42. Atangana, A, Secer, A: The time-fractionalcoupled-Korteweg-de-Vries equations. Abstr. Appl. Anal. 2013, Article ID 947986(2013)

Page 12 of 12

43. Kilbas, AA, Srivastava, HM, Trujillo, JJ: Theory and Applications of FractionalDifferentialEquations. North-Holland Mathematics Studies, vol. 204. Elsevier, Amsterdam (2006)

44. Podlubny, I: FractionalDifferentialEquations. Academic Press, San Diego (1999)

45. Samko, SG, Kilbas, AA, Marichev, OI: Fractionallntegrals and Derivatives: Theory and Applications. Gordon & Breach, New York (1993)

46. He, J: Semi-inverse method of establishing generalized variationalprinciples for fluid mechanics with emphasis on turbo-machinery aerodynamics. Int. J. Turbo Jet-Engines 14, 23-28 (1997)

47. He, J: Variationaliteration method - a kind of nonlinear analyticaltechnique: some examples. Int. J. Non-Linear Mech. 34,699-708(1999)

48. Agrawal, OP: Formulation of Euler-Lagrange equations for fractionalvariationalproblems. J. Math. Anal. Appl. 272, 368-379 (2002)

49. Agrawal, OP: A generalformulation and solution scheme for fractionaloptimalcontrolproblems. Nonlinear Dyn. 38, 323-337 (2004)

50. Agrawal, OP: Fractionalvariationalcalculus in terms of Riesz fractionalderivatives. J. Phys. A, Math. Theor. 40,62-87 (2007)


Cite this article as: Zhang: Formulation and solution to time-fractional generalized Korteweg-de Vries equation via variational methods. Advances in Difference Equations 2014, 2014:65

Submit your manuscript to a SpringerOpen journal and benefit from:

► Convenient online submission

► Rigorous peer review

► Immediate publication on acceptance

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article

Submit your next manuscript at ►