The fractional variational iteration method is used to solve the fractional Cauchy problem. Some examples are given to elucidate the solution procedure and reliability of the obtained results. The variational iteration algorithm leads to exact solutions in the present study.
1. Introduction
Fractional problems have attracted many scholars’ attention as the actual problems can be much better simulated by using the fractional derivatives than by using traditional integral derivatives [1–6].
In [7], the author gave a very lucid as well as elementary discussion of the variational iteration method. The variational iteration algorithm [7–14] is very simple, with results of high precision (sometimes exact solutions can be obtained), and is easy to understand and has been widely applied in various nonlinear problems. In addition, many authors have made a great effort to give sophisticated theoretical verification of the variational iteration method; for example, Odibat [15], Salkuyeh [16], and Tatari and Dehghan [17] proved that the variational iteration algorithm leads to convergent results. Due to its flexibility and ability to solve nonlinear equations accurately and conveniently, the method has been modified or improved to solve nonlinear problems more efficiently [18–21]; furthermore the method has been extended to handle fractional nonlinear models and the fractional variational iteration method (FVIM) has been presented [22].
Many scholars have applied FVIM to solve linear or nonlinear fractional order differential equations [23–28]. For a relatively comprehensive survey on the concepts, theory, and applications of the method, readers are referred to review articles [10, 29].
In [22] the improved VIM is called fractional variational iteration method and the new Lagrange multiplier is determined by using the Laplace transformation. In this paper, we use the fractional variational iteration method to discuss the fractional-order partial differential equation in the form
(1)ut(α)(x,t)+a(x,t)ux(x,t)=ϕ(x),x∈R,t>0(2)u(x,0)=ψ(x),x∈R.
When a(x,t)=a is a constant and ϕ(x)=0, (1) is a linear equation called the fractional transport equation which can describe many interesting phenomena such as the spread of AIDS and the moving of wind. When a(x,t)=u(x,t), the equation is called the nefractional inviscid Burgers’ equation arising in a one-dimensional stream of particles or fluid having zero viscosity. When α=1, (1) is the equation for traditional Cauchy problem [30, 31].
2. The Fractional Variational Iteration Method
To discuss the fractional problems, two definitions are introduced.
Definition 1.
The Caputo derivative is given as
(3)D0Ctαf(t)=1Γ(m-α)∫at(t-τ)m-α-1f(m)(τ)dτ,t>a,m-1<α<m∈Z+.
Definition 2.
The Riemann-Liouville (R-L) integration of f(t) is defined as
(4)Ia,tαf(t)=1Γ(α)∫at(t-τ)α-1f(τ)dτ,α>0.
To illustrate the basic concepts of the fractional variational iteration method, one considers the following general nonlinear fractional system:
(5)0CDtαu+R[u]+N[u]=f(τ),
where R[u] is a linear term and N[u] is a nonlinear one and f(t) is a known analytic function.
According to the fractional variational iteration method, a correction functional can be constructed as follows:
(6)un+1=un+∫0tλ(t,τ)[D0Ctαun+R[un]+N[un]-f(τ)]dτ,
where λ(t,τ) is a general Lagrange multiplier [7–12], which can be identified optimally via the variational theory, the subscript n denotes the nth approximation, and R[un] and N[un] are considered as a restricted variation [7]; that is δR[un]=0, δN[un]=0.
Taking Laplace transform on the correction functional equation established via the R-L integration, Wu and Baleanu [22] present a new way to identify the Lagrange multiplier. The Lagrange multiplier can be identified as
(7)λ(t,τ)=(-1)α(τ-t)α-1Γ(α).
The variational iteration formula (6) can be improved as
(8)un+1=un+∫0t(-1)α(τ-t)α-1Γ(α)[D0Ctαun+R[un]un+∫0t(-1)α(τ-t)α-1Γ(α)+N[un]-f(τ)D0Ctα]dτ.
The initial guess u0 can be freely chosen with possible unknown constants; it can also be solved from its corresponding linear homogeneous equation.
The above iteration formula (8) is also valid for differential equations when α is an arbitrary positive integer.
The fractional variational iteration method can solve effectively, easily, and accurately a large class of nonlinear fractional problems with approximations converging rapidly to the accurate solution.
3. Applications
Since our focus is on the ideas and basic principles, we will consider only the simplest possible equations to clearly illustrate the solution procedure. In particular, we will focus on pure Cauchy problems. These problems are initial value problems.
According to (8), we can construct a correction functional to (1) which reads
(9)un+1(x,t)=un(x,t)+∫0t(-1)α(τ-t)α-1Γ(α)[D0Ctαun(x,τ)+∫0t(-1)α(τ-t)α-1Γ(α)+a(x,τ)ux(x,τ)+∫0t(-1)α(τ-t)α-1Γ(α)-ϕ(x)0CDtα]dτ.
Example 3.
Consider the fractional transport equation:
(10)ut(α)(x,t)+aux(x,t)=0,x∈R,t>0u(x,0)=x2,x∈R.
According to the iteration formula (8), we can obtain
(11)un+1=un-∫0t(τ-t)αΓ(α)[unτ(α)(x,τ)+aunx(x,t)]dτ,unτ(α)(x,τ)+aunx(x,t)(x,t)iiiiiiix∈R,t>0.
The initial iterative value is selected as u0(x,0)=u(x,0)=x2 from the given initial condition. Using (10), we obtain the following successive approximations:
(12)u1(x,t)=x2-2axtαΓ(α+1)u2(x,t)=x2-2axtαΓ(α+1)+2a2t2αΓ(2α+1)u3(x,t)=x2-2axtαΓ(α+1)+2a2t2αΓ(2α+1)⋮un(x,t)=x2-2axtαΓ(α+1)+2a2t2αΓ(2α+1).
So the exact analytical solution is yielded as
(13)u(x,t)=x2-2axtαΓ(α+1)+2a2t2αΓ(2α+1).
Example 4.
Consider the fractional nonlinear Cauchy problem
(14)ut(α)(x,t)+xux(x,t)=0,x∈R,t>0u0(x,0)=u(x,0)=x2,x∈R.
According to the iteration formula (8), we can obtain
(15)un+1=un-∫0t(τ-t)αΓ(α)[ut(α)(x,t)+xux(x,t)]dξ,x∈R,t>0.
The initial iterative value is selected as u0(x,0)=u(x,0)=x2 from the given initial condition. Using (15), we obtain the following successive approximations:
(16)u0(x,t)=x2u1(x,t)=x2-2x2tαΓ(α+1)u2(x,t)=x2-2x2tαΓ(α+1)+x2(2tα)2Γ(2α+1)u3(x,t)=x2-2x2tαΓ(α+1)+x2(2tα)2Γ(2α+1)-x2(2tα)3Γ(3α+1)⋮un(x,t)=x2[1-2tαΓ(α+1)+(2tα)2Γ(2α+1)-(2tα)3Γ(3α+1)+(2tα)4Γ(4α+1)-(2tα)5Γ(5α+1)+⋯].
The VIM admits the use of u=limn→∞un, which gives the exact solution
(17)u(x,t)=x2Eα(-2tα),
where Eα is Mittag-Leffler function
(18)Eα(z)=∑n=0∞znΓ(αn+1).
Example 5.
Consider the following nonhomogeneous fractional Cauchy problem:
(19)ut(α)(x,t)+xux(x,t)=x,x∈R,t>0u(x,0)=ex,x∈R.
According to the iteration formula (9), we can obtain
(20)un+1=un-∫0t(τ-t)αΓ(α)[ut(α)(x,t)+xux(x,t)-x]dξ.
The following successive approximations are obtained by using (20) with the selected initial value u0=ex:
(21)u1(x,t)=ex+(x-ex)tαΓ(α+1)u2(x,t)=x2+(x-ex)tαΓ(α+1)+(ex-1)t2αΓ(2α+1)u3(x,t)=x2+(x-ex)tαΓ(α+1)+(ex-1)t2αΓ(2α+1)-ext3αΓ(3α+1)⋮un(x,t)=[x-Γ(α+1)tαΓ(2α+1)]tαΓ(α+1)+ex[1-tαΓ(α+1)+t2αΓ(2α+1)-t3αΓ(3α+1)+t4αΓ(4α+1)-t5αΓ(5α+1)+⋯].
By using u=limn→∞un, the exact solution of the equation is obtained as
(22)u(x,t)=[x-Γ(α+1)tαΓ(2α+1)]tαΓ(α+1)+exEα(-tα).
Particularly when α=1, the result in (22) is u(x,t)=(x-(t/2))t+ex-t, which is the same as that in [30].
Example 6.
Consider the fractional inviscid Burgers’ equation
(23)ut(α)(x,t)+u(x,t)ux(x,t)=0,x∈R,t>0u(x,0)=x,x∈R.
According to the iteration formula (8), we can obtain
(24)un+1=un-∫0t(τ-t)αΓ(α)[ut(α)(x,τ)+u(x,τ)ux(x,τ)]dτ.
Starting with initial approximation u0(x,0)=u(x,0)=x and proceeding in a similar way illustrated above, we obtain the following successive approximations:
(25)u0(x,t)=1-tαΓ(α+1)u1(x,t)=1-tαΓ(α+1)+2t2αΓ(2α+1)-Γ(2α+1)t3αΓ(3α+1)Γ2(α+1)u2(x,t)=1-tαΓ(α+1)+2t2αΓ(2α+1)+a3t3αΓ(3α+1)+a4t4αΓ(4α+1)+a5t5αΓ(5α+1)+a6t6αΓ(6α+1)+a7t7αΓ(7α+1)⋮,
where
(26)a3=-4Γ(1+3α)-Γ(1+2α)Γ2(1+α)Γ(1+3α)a4=2αΓ(4α)+Γ(2α)α2Γ2(α)Γ(4α)+3Γ(3α)2α2Γ(α)Γ(2α)Γ(4α)a5=-4Γ(1+4α)Γ2(1+2α)-2Γ(1+2α)Γ(1+4α)Γ3(1+α)Γ(1+3α)a6=4Γ(1+5α)Γ2(1+α)Γ(1+3α)(27)a7=-Γ2(1+2α)Γ(1+6α)Γ4(1+α)Γ2(1+3α).
Particularly for the case α=1, the result for u3 in (23) will be
(28)u3=x(1-t+t2-t3+23t4-23t5+19t6-163t7),
which is the same as that in [30].
In order to show the convergence of the iteration solutions, we introduce a function vi(x,t):
(29)vi(x,t)=ui(x,t)x.
It can be seen that vi is a function of t. When α=1, the accurate solution of (23) is u=x/(1+t); then v=1/(1+t).
Figures 1, 2, and 3 show the curves of vi (where i means the iteration times) changing with time t when α=0.5, 1, and 1.5. When α=1, the fifth iteration result v5 is almost the same as the exact solution at t∈[0,1]. When α=1.5, the iteration results converge fast and the fourth iteration result v4 is almost the same as v5 at t∈[0,3].
vi(t)’s curve changing with time t when α=0.5.
vi(t)’s curve changing with time t when α=1.
vi(t)’s curve changing with time t when α=1.5.
4. Conclusions
The variational iteration method has been proved by many authors to be a powerful mathematical tool for various kinds of nonlinear problems.
In this paper, we extend the use of fractional variational iteration method to fractional Cauchy problems and give the numerical examples. Compared with the classical VIM, the modified version method is powerful for solving differential equations with fractional derivatives. The higher order approximate solutions of the Cauchy equation illustrate the method’s efficiency and high accuracy.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The author would like to thank the reviewers for their careful reading and helpful comments. The work is supported by the National Natural Science Foundation of China (Grant no. 11202146).
ChenW.SunH. G.LiX. C.20101stBeijing, ChinaScience Press(Chinese)MillerK. S.RossB.1993New York, NY, USAWileyMR1219954OldhamK. B.SpanierJ.1974New York, NY, USAAcademic PressMR0361633DebnathL.Fractional integrals and fractional differential equations in fluid mechanics20036119155MR2035411ZBL1076.35095PodlubnyI.1999New York, NY, USAAcademic PressMR1658022KilbasA. A.SrivastavaH. M.TrujilloJ. J.2006Amsterdam, The NetherlandsElsevierHeJ. H.Approximate analytical solution for seepage flow with fractional derivatives in porous media19981671-257682-s2.0-003230766110.1016/S0045-7825(98)00108-XMR1665221ZBL0942.76077HeJ. H.Variational iteration method—a kind of non-linear analytical technique: some examples19993446997082-s2.0-0000092673HeJ. H.Variational iteration method—some recent results and new interpretations200720731710.1016/j.cam.2006.07.009MR2332941ZBL1119.65049HeJ. H.WuX. H.Variational iteration method: new development and applications2007547-88818942-s2.0-3474887067710.1016/j.camwa.2006.12.083MR2395625ZBL1141.65372HeJ. H.WuG.-C.AustinF.The variational iteration method which should be followed20101130Jia-QiM.Wei-JiangZ.Xian-FengC.Variational iteration method for solving a class of strongly nonlinear evolution equations20095811739774012-s2.0-71249150404MR2649158ZBL1212.35065NoorM. A.Mohyud-DinS. T.Variational iteration method for solving higher-order nonlinear boundary value problems using He's polynomials2008921411562-s2.0-4374910793610.1016/j.amc.2006.12.071MR2332146ZBL1122.65374WazwazA. M.A study on linear and nonlinear Schrodinger equations by the variational iteration method2008374113611422-s2.0-4064910075110.1016/j.chaos.2006.10.009MR2411539ZBL1148.35353OdibatZ. M.A study on the convergence of variational iteration method2010519-10118111922-s2.0-7764918365010.1016/j.mcm.2009.12.034MR2608904ZBL1198.65147SalkuyehD. K.Convergence of the variational iteration method for solving linear systems of ODEs with constant coefficients2008568202720332-s2.0-5104912429210.1016/j.camwa.2008.03.030MR2466705ZBL1165.65376TatariM.DehghanM.On the convergence of He's variational iteration method200720711211282-s2.0-3425069638310.1016/j.cam.2006.07.017MR2332953ZBL1120.65112GengF.A modified variational iteration method for solving Riccati differential equations2010607186818722-s2.0-7795791121410.1016/j.camwa.2010.07.017MR2719705ZBL1205.65229GhorbaniA.MomaniS.An effective variational iteration algorithm for solving Riccati differential equations20102389229272-s2.0-8485970607910.1016/j.aml.2010.04.012MR2651475ZBL1192.65095TatariM.DehghanM.Improvement of He's variational iteration method for solving systems of differential equations20095811-12216021662-s2.0-7035057045810.1016/j.camwa.2009.03.081MR2557343ZBL1189.65178HesameddiniE.LatifizadehH.Reconstruction of variational iteration algorithms using the Laplace transform20091011-12137713822-s2.0-77950427865WuG.-C.BaleanuD.Variational iteration method for the Burgers' flow with fractional derivatives—new Lagrange multipliers20133796183619010.1016/j.apm.2012.12.018MR3039000KhanY.FarazN.YildirimA.WuQ.Fractional variational iteration method for fractional initial-boundary value problems arising in the application of nonlinear science2011625227322782-s2.0-8005226378310.1016/j.camwa.2011.07.014MR2831688ZBL1231.35288WuG.-C.LeeE. W. M.Fractional variational iteration method and its application201037425250625092-s2.0-7795347899110.1016/j.physleta.2010.04.034MR2640023ZBL1237.34007GuoS.MeiL.The fractional variational iteration method using He's polynomials201137533093132-s2.0-7865088184110.1016/j.physleta.2010.11.047MR2748835ZBL1241.35216WuG.-C.A fractional variational iteration method for solving fractional nonlinear differential equations2011618218621902-s2.0-7995373581310.1016/j.camwa.2010.09.010MR2785581ZBL1219.65085JiJ.ZhangJ.DongY.The fractional variational iteration method improved with the Adomian series201225122223222610.1016/j.aml.2012.06.007MR2967820ZBL1252.65115WuG.-C.Challenge in the variational iteration method—a new approach to identification of the Lagrange multipliers201325217517810.1016/j.jksus.2012.12.002HeJ. H.Asymptotic methods for solitary solutions and compactons20122012130916793MR2999914ZBL1257.3515810.1155/2012/916793ZhouX.-W.YaoL.The variational iteration method for Cauchy problems20106037567602-s2.0-7795542046510.1016/j.camwa.2010.05.022McowenR. C.1996Prentice Hall