MPEMathematical Problems in Engineering1563-51471024-123XHindawi Publishing Corporation73046510.1155/2009/730465730465Research ArticleOn the Numerical Solution of Fractional Hyperbolic Partial Differential EquationsAshyralyevAllaberen1, 2DalFadime3PinarZehra3GonçalvesPaulo Batista1Department of MathematicsFatih UniversityBuyukcekmece34500 IstanbulTurkeyfatih.edu.tr2Department of MathematicsITTU744012 AshgabatTurkmenistanittu.edu.tm3Department of MathematicsEge University35100 IzmirTurkeyege.edu.tr200908092009200926022009170720092009Copyright © 2009This 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 stable difference scheme for the numerical solution of the mixed problem for the multidimensional fractional hyperbolic equation is presented. Stability estimates for the solution of this difference scheme and for the first and second orders difference derivatives are obtained. A procedure of modified Gauss elimination method is used for solving this difference scheme in the case of one-dimensional fractional hyperbolic partial differential equations.

1. Introduction

It is known that various problems in fluid mechanics (dynamics, elasticity) and other areas of physics lead to fractional partial differential equations. Methods of solutions of problems for fractional differential equations have been studied extensively by many researchers (see, e.g.,  and the references given therein).

The role played by stability inequalities (well posedness) in the study of boundary-value problems for hyperbolic partial differential equations is well known (see, e.g., ). In the present paper, the mixed boundary value problem for the multidimensional fractional hyperbolic equation

2u(t,x)t2-r=1m(ar(x)uxr)xr+Dt1/2u(t,x)=f(t,x),x=(x1,,xm)Ω,0<t<1,u(0,x)=0,ut(0,x)=0,xΩ¯,u(t,x)=0,xS is considered. Here Dt1/2=D0+1/2 is the standard Riemann-Lioville's derivative of order 1/2 and Ω is the unit open cube in the m-dimensional Euclidean space m:{Ω=x=(x1,,xm):0<xj<1,1jm} with boundary S,Ω¯=ΩS,  ar(x),(xΩ) and f(t,x)(t(0,1),xΩ) are given smooth functions and ar(x)a>0.

The first order of accuracy in t and the second order of accuracy in space variables for the approximate solution of problem (1.1) are presented. The stability estimates for the solution of this difference scheme and its first and second ordes difference derivatives are established. A procedure of modified Gauss elimination method is used for solving this difference scheme in the case of one-dimensional fractional hyperbolic partial differential equations.

2. The Difference Scheme and Stability Estimates

The discretization of problem (1.1) is carried out in two steps. In the first step, let us define the grid space

Ω̃h={x=xr=(h1r1,,hmrm),r=(r1,,rm),0rjNj,hjNj=1,j=1,,m},Ωh=Ω̃hΩ,Sh=Ω̃hS. We introduce the Banach space L2h=L2(Ω̃h) of the grid functions φh(x)={φ(h1r1,,hmrm)} defined on Ω̃h, equipped with the norm

φhL2(Ω̃h)=(xΩh¯|φh(x)|2h1hm)1/2. To the differential operator Ax generated by problem (1.1), we assign the difference operator Ahx by the formula

Ahxuxh=-r=1m(ar(x)ux̅rh)xr,jr acting in the space of grid functions uh(x), satisfying the conditions uh(x)=0 for all xSh. It is known that Ahx is a self-adjoint positive definite operator in L2(Ω̃h). With the help of Ahx we arrive at the initial boundary value problem

d2vh(t,x)dt2+Ahxvh(t,x)+Dt1/2vh(t,x)=fh(t,x),0t1,xΩh,vh(0,x)=0,dvh(0,x)dt=0,xΩ̃ for an finite system of ordinary fractional differential equations.

In the second step, applying the first order of approximation formula (1/π)m=1k(Γ(k-m+(1/2))/(k-m)!)((u(tk)-u(tk-1))/τ1/2)   for Dt1/2u(t) (see ) and using the first order of accuracy stable difference scheme for hyperbolic equations (see ), one can present the first order of acuraccy difference scheme

uk+1h(x)-2ukh(x)+uk-1h(x)τ2+Ahxuk+1h+1πm=1kΓ(k-m+1/2)(k-m)!(umh-um-1h)τ1/2=fkh(x),xΩ̃h,fkh(x)=f(tk,xn),tk=kτ,1kN-1,Nτ=1,u1h(x)-u0h(x)τ=0,u0h(x)=0,xΩ̃h for the approximate solution of problem (2.4). Here Γ(k-m+1/2)=0tk-m-1/2e-tdt.

Theorem 2.1.

Let τ and |h| be sufficiently small numbers. Then, the solutions of difference scheme (2.5) satisfy the following stability estimates: max1kNukhL2h+max1kNukh-uk-1hτL2hC1max1kN-1fkhL2h,max1kN-1τ-2(uk+1h-2ukh+uk-1h)L2h+max1kNr=1m(ukh)x̅rxr,jrL2hC2[f1hL2h+max2kN-1τ-1(fkh-fk-1h)L2h]. Here C1 and C2 do not depend on τ,h, and fkh,1k<N-1.

The proof of Theorem 2.1 is based on the self-adjointness and positive definitness of operator Ahx in L2h and on the following theorem on the coercivity inequality for the solution of the elliptic difference problem in L2h.

Theorem 2.2.

For the solutions of the elliptic difference problem Ahxuh(x)=ωh(x),xΩh,uh(x)=0,xSh, the following coercivity inequality holds : r=1muxrx̅r,jrhL2hCωhL2h.

Remark 2.3.

The stability estimates of Theorem 2.1 are satisfied in the case of operator Au=-k=1nak(x)2uxk2+k=1nbk(x)uxk+c(x)u with Dirichlet condition u=0 in S and Dtα=D0+α is the standard Riemann-Lioville's derivative of order α,0α<1. In this case, A is not self-adjoint operator in H. Nevertheless, Au=A0u+Bu and A0 is a self-adjoint positive definite operator in H and BA0-1 is bounded in H. The proof of this statement is based on the abstract results of  and difference analogy of integral inequality.

Remark 2.4.

The stability estimates of Theorem 2.1 permit us to obtain the estimate of convergence of difference scheme of the first order of accuracy for approximate solutions of the initial-boundary value problem 2u(t,x)t2-r=1nar(x)uxrxr+r=1nbr(x)uxr+Dtαu(t,x)=f(t,x;u(t,x),ut(t,x),ux1(t,x),,uxn(t,x)),x=(x1,,xn)Ω,0<t<1,u(0,x)=0,u(0,x)t=0,xΩ¯,u(t,x)=0,xS for semilinear fractional hyperbolic partial differential equations.

Note that, one has not been able to obtain a sharp estimate for the constants figuring in the stability estimates of Theorem 2.1. Therefore, our interest in the present paper is studying the difference scheme (2.5) by numerical experiments. Applying this difference scheme, the numerical methods are proposed in the following section for solving the one-dimensional fractional hyperbolic partial differential equation. The method is illustrated by numerical experiments.

3. Numerical Results

For the numerical result, the mixed problem

Dt2u(t,x)-uxx(t,x)+Dt1/2u(t,x)=f(t,x),f(t,x)=(2-8t3/23π+(πt)2)sinπx,0<t,x<1,u(0,x)=0,ut(0,x)=0,0x1,u(t,0)=u(t,1)=0,0t1 for the one-dimensional fractional hyperbolic partial differential equation is considered. Applying difference scheme (2.5), we obtain

unk+1-2unk+unk-1τ2-un+1k+1-2unk+1+un-1k+1h2+1Πm=1kΓ(k-m+1/2)(k-m)!(unm-unm-1τ1/2)=φnk,φnk=f(tk,xn),1kN-1,1nM-1,un0=0,τ-1(un1-un0)=0,0nM,u0k=uMk=0,0kN. We get the system of equations in the matrix form

AUn+1+BUn+CUn-1=Dφn,1nM-1,U0=0̃,UM=0,̃ where

0̃=(N+1)×(1),A=[00000000000000a000000a000000a0000000](N+1)×(N+1),B=[b1,100000b2,1b2,20000b3,1b3,2b3,3000b4,1b4,2b4,3b4,400bN,1bN,2bN,3bN,4bN,N0bN+1,1bN+1,2bN+1,3bN+1,4bN+1,NbN+1,N+1](N+1)×(N+1),C=A,D=(N+1)×(N+1),Us=[Us0Us1Us2Us3UsN-1UsN](N+1)×(1),s=n-1,n,n+1,a=-1h2,b1,1=1,b2,1=-1τ,b2,2=1τ,b3,1=1τ2+1τ1/2,b3,2=-2τ2-1τ1/2,b3,2=1τ2+2h1/2,bk+2,1=1πΓ(k-1+1/2)Γ(k)τ1/2,2kN-1,bk+2,k+1=-2τ2-1τ1/2,1kN-1,bk+2,k=1τ2+1π(-Γ(3/2)Γ(2)+Γ(1/2)Γ(1))1τ1/2,2kN-1,bk+2,k+2=1τ2+2h2,1kN-1,bk+2,i+1=1π(-Γ(k-i+1/2)Γ(k-(i-1))+Γ(k-(i+1)+1/2)Γ(k-(i-1)-1))1τ1/2,3kN-1,1ik-2,φnk=(2-8(kτ)3/23π+(πkτ)2)sinπ(nh)),φn=[φn0φn1φn2φnN](N+1)×1. So, we have the second-order difference equation with respect to n matrix coefficients. This type system was developed by Samarskii and Nikolaev . To solve this difference equation we have applied a procedure of modified Gauss elimination method for difference equation with respect to k matrix coefficients. Hence, we seek a solution of the matrix equation in the following form:

Uj=αj+1Uj+1+βj+1,n=M-1,,2,1, where αj(j=1,,M) are (N+1)×(N+1) square matrices and βj(j=1,,M) are (N+1)×1 column matrices defined by

αn+1=-(B+Cαn)-1A,βn+1=(B+Cαn)-1(Dφn-Cβn),n=1,2,,M-1, where

α1=(N+1)×(N+1),β1=(N+1)×1. Now, we will give the results of the numerical analysis. First, as we noted above one can not obtain a sharp estimate for the constants C1 and C2 figuring in the stability estimates of Theorem 2.1. We have

C1=maxf,u(Ct1),C2=maxf,u(Ct2), where

Ct1=max1kNukhL2h+max1kNτ-1(ukh-uk-1h)L2h(max1kN-1fkhL2h)-1,Ct2=[max1kN-1τ-2(uk+1h-2ukh+uk-1h)L2h+max1kN-1r=1n(uk+1h)x¯r,xr,jrL2h]×(max2kN-1τ-1(fkh-fk-1h)L2h+f1hL2h)-1. The constants Ct1 and Ct2 in the case of numerical solution of initial-boundary value problem (3.1) are computed.

The numerical solutions are recorded for different values of N and M,unk represents the numerical solutions of this difference scheme at (tk,xn). The constants Ct1 and Ct2 are given in Table 1 for N=20,40,80, and M=80, respectively.

The difference scheme.

 M=80 M=80 M=80 N=20 N=40 N=80 The values of Ct1 1.0379 1.0667 1.0800 The values of Ct2 0.6265 0.6186 0.6140

Recall that we have not been able to obtain a sharp estimate for the constants C1 and C2 figuring in the stability estimates. The numerical results in the Tables 1 and 2 give Ct11.00 and Ct20.62, respectively. That means the constants C1 and C2 figuring in the stability estimates in the case of numerical solution of initial-boundary value problem (3.1) of this difference scheme is stable with no large constants.

The difference scheme.

 Method M=80 M=80 M=80 N=20 N=40 N=80 Comparison of errors (E0) for approximate solutions 0.0452 0.0235 0.0125 Comparison of errors (E1) for approximate solutions 0.2702 0.1452 0.0779

Second, for the accurate comparison of the difference scheme considered, the errors computed by

E0=max1kN-1(n=1M-1|u(tk,xn)-unk|2h)1/2,E1=max1kN-1(n=1M-1|utt(tk,xn)-unk+1-2unk+unk-1τ2|2h)1/2 of the numerical solution are recorded for higher values of N and M, where u(tk,xn) represents the exact solution and unk represents the numerical solution at (tk,xn). The errors E0 and E1 results are shown in Table 2 for N=20,40,80 and M=80, respectively.

The exact solution.

Difference Scheme.

PodlubnyI.Fractional Differential Equations1999198San Diego, Calif, USAAcademic Pressxxiv+340Mathematics in Science and EngineeringMR1658022ZBL1056.93542SamkoS. G.KilbasA. A.MarichevO. I.Fractional Integrals and Derivatives1993Yverdon, SwitzerlandGordon and Breachxxxvi+976MR1347689ZBL0924.44003LavoieJ.-L.OslerT. J.TremblayR.Fractional derivatives and special functionsSIAM Review1976182240268MR040991410.1137/1018042ZBL0324.44002TarasovV. E.Fractional derivative as fractional power of derivativeInternational Journal of Mathematics2007183281299MR231461310.1142/S0129167X07004102ZBL1119.26011El-MesiryA. E. M.El-SayedA. M. A.El-SakaH. A. A.Numerical methods for multi-term fractional (arbitrary) orders differential equationsApplied Mathematics and Computation20051603683699MR2113113ZBL1062.65073El-SayedA. M. A.GaafarF. M.Fractional-order differential equations with memory and fractional-order relaxation-oscillation modelPure Mathematics and Applications2001123296310MR1942895El-SayedA. M. A.El-MesiryA. E. M.El-SakaH. A. A.Numerical solution for multi-term fractional (arbitrary) orders differential equationsComputational & Applied Mathematics20042313354MR2111630GorenfloR.MainardiF.CarpinteriA.MainardiF.Fractional calculus: integral and differential equations of fractional orderFractals and Fractional Calculus in Continuum Mechanics1997378Vienna, AustriaSpringer223276CISM Courses and LecturesMR1611585MatignonD.Stability results for fractional differential equations with applications to control processingComputational Engineering in System Application19962Lille, FranceAshyralyevA.A note on fractional derivatives and fractional powers of operatorsJournal of Mathematical Analysis and Applications20093571232236MR252682310.1016/j.jmaa.2009.04.012PodlubnyI.El-SayedA. M. A.On Two Definitions of Fractional Calculus1996Slovak Academy of Science-Institute of Experimental PhysicsKreinS. G.Linear Differential Equations in a Banach Space1967Moscow, RussiaNaukaMR0247239SobolevskiiP. E.ChebotaryevaL. M.Approximate solution by method of lines of the Cauchy problem for an abstract hyperbolic equationsIzvestiya Vysshikh Uchebnykh Zavedenii. Matematika19775103116AshyralyevA.MartinezM.PasterJ.PiskarevS.Weak maximal regularity for abstract hyperbolic problems in function spaces, further progress in analysisProceedings of the 6th International ISAAC CongressAugust 2007Ankara, TurkeyWorld Scientific679689AshyralyevA.AggezN.A note on the difference schemes of the nonlocal boundary value problems for hyperbolic equationsNumerical Functional Analysis and Optimization2004255-6439462MR210626910.1081/NFA-200041711ZBL1065.35021AshyralyevA.MuradovI.On difference schemes a second order of accuracy for hyperbolic equationsModelling Processes of Explotation of Gas Places and Applied Problems of Theoretical Gasohydrodynamics1998Ashgabat, TurkmenistanIlim127138AshyralyevA.SobolevskiiP. E.New Difference Schemes for Partial Differential Equations2004148Basel, SwitzerlandBirkhäuserx+443Operator Theory: Advances and ApplicationsMR2075306ZBL1069.65097AshyralyevA.OzdemirY.On nonlocal boundary value problems for hyperbolic-parabolic equationsTaiwanese Journal of Mathematics200711410751089MR2348553ZBL1140.65039AshyralyevA.YildirimO.On multipoint nonlocal boundary value problems for hyperbolic differential and difference equationsTaiwanese Journal of Mathematics20091322SamarskiiA. A.GavrilyukI. P.MakarovV. L.Stability and regularization of three-level difference schemes with unbounded operator coefficients in Banach spacesSIAM Journal on Numerical Analysis2001392708723MR186025610.1137/S0036142999357221ZBL1002.65101AshyralyevA.SobolevskiiP. E.Two new approaches for construction of the high order of accuracy difference schemes for hyperbolic differential equationsDiscrete Dynamics in Nature and Society200520052183213MR220345610.1155/DDNS.2005.183ZBL1094.65077AshyralyevA.KoksalM. E.On the second order of accuracy difference scheme for hyperbolic equations in a Hilbert spaceNumerical Functional Analysis and Optimization2005267-8739772MR219202110.1080/01630560500431068ZBL1098.65055AshyralyevA.KoksalM. E.On the stability of the second order of accuracy difference scheme for hyperbolic equations in a Hilbert spaceDiscrete Dynamics in Nature and Society200720072657491AshyraliyevM.A note on the stability of the integral-differential equation of the hyperbolic type in a Hilbert spaceNumerical Functional Analysis and Optimization2008297-8750769MR245018910.1080/01630560802292069ZBL1146.45001AshyralyevA.SobolevskiiP. E.A note on the difference schemes for hyperbolic equationsAbstract and Applied Analysis2001626370MR185875010.1155/S1085337501000501ZBL1007.65064SobolevskiiP. E.Difference Methods for the Approximate Solution of Differential Equations1975Voronezh, RussiaIzdatelstvo Voronezhskogo Gosud UniversitetaSamarskiiA. A.NikolaevE. S.Numerical Methods for Grid Equations19892Basel, SwitzerlandBirkhäuserxvi+502Iterative MethodsMR1004469