On the numerical treatment ofthenonlinear partial differentialequation of fractional order

Transcription

1 IOSR Journal of Mahemacs (IOSR-JM) e-iss: , p-iss: X. Volume 2, Issue 6 Ver. I (ov. - Dec.26), PP On he numercal reamen ofhenonlnear paral dfferenalequaon of fraconal order M. M. El-Kojok Deparmen of Mahemacs, College of Scences and Ars, Qassm Unversy, K.S.A Absrac:The exsence and unqueness soluon of a nonlnear paral dfferenal equaon (PDE)of fraconal order are dscussed and proved n a Banach space H due o Pcard s mehod dependng on he properes we expec a soluon o possess. Moreover, some properes concernng he sably of soluons are obaned.a sysemof nonlnear Volerra negral equaons of he second knd (SVIEs) s obaned. The modfed Toeplz marx mehod (MTMM) s used, as a numercal mehod, o oban a nonlnear sysemof algebrac equaons (AS). Also, many mporan heorems relaed o he exsence and unqueness soluon of he algebrac sysem are derved. Fnally,numercal resuls are obaned and he error s calculaed. Keywords: nonlnear paral dfferenal equaon of fraconal order,semgroup, sysem of nonlnear Volerra negral equaons of he second knd,modfed Toeplz marx mehod, nonlnear algebrac sysem. I. Inroducon and Basc formulaons The semgroups mehods play a specal role for paral dfferenal equaons and n applcaons, for example hey descrbe how denses of nal saes evolve n me.moreover, here are equaons whch generae semgroups. These equaons appear n such dverse areas as asrophyscs-flucuaons n he brghness of he Mlky Way [], populaon dynamcs [2,3], and n he heory of jump processes.in [4], El-Bora suded he Cauchy problem n a Banach space E for a lnear fraconal evoluon equaon. In hs paper, he exsence and unqueness of he soluon of he Cauchy problem were dscussed and proved. Also, he soluon was obaned n erms of some probably denses. In [5], El-Bora dscussed he exsence and unqueness soluon of he nonlnear Cauchy problem. Also, some properes concernng he sably of soluons were obaned. In [6] Abdou e al., mproved he work of El-Bora n[4]and obanednumercally he soluon of he Cauchy problem. More nformaon for he PDEsand s soluon can be found n [7-9]. Consder he followng PDE of fraconal order: α u(x, ) α = A u x, + F x,, B u x,,. wh he nal condon: u x, = u x, (.2) n a Banach space H, where u x, s an H-valued funcon on H, T, T <, A s a lnear closed and bounded operaor defned on a dense ses, B,, T s a famly of lnear closed and bounded operaors defned on a dense se S2 S n H no H, F s a gven absrac funcon defned on H, T wh values n H, u x H and < α <. I s assumed ha A generaes an analyc semgroup Q().Ths condon mples: Q() k for, and AQ k for >, (.3) where. s he norm n H and k s a posve consan (see Zadman [7] ). Le us pose ha B g s unformly Hölder connuous n, T, for every g S ; ha s B 2 g B g k ( 2 ) β, (.4) for all 2 >,, 2, T, where k and β are posve consans, β. We pose also ha here exss a number γ (,), such ha B 2 Q k 2 DOI:.979/ Page γ, (.5) where >, 2, T, Hand k 2 s a posve consan (see[4,8,9]). (oce ha Q S for each H and each > ). Also, s assumed ha, he funcon F sasfes he followng condons: () F s unformly Hölder connuous n, T ; ha s F x, 2, W F x,, W l 2 β, (.6) for all 2 >,, 2, T and all x, W H, where land β are posve consans; β, W = B u(x, ) and. s he norm n H. () F sasfes Lpschz condon F x,, W F(x,, W ) x, W W ; x, l, (.7)

2 On he numercal reamen of he nonlnear paral dfferenal equaon of fraconal order for allx, W, W H and all, T, where l s a posve consan. () F x,, W l W, l s a consan. (.8) Followng Gelfand and Shlov (see [],[]), we can defne he negral of order α > by I α f = Γ α If < α <, we can defnehe dervave of orderα by d α f() d α = Γ α f (θ) θ α dθ, wherefs an absrac funcon wh values n H. ow, s suable o rewre he PDE((.), (.2))n he form u x, = u x + Γ α θ α A u x, θ dθ + Γ α θ α f θ dθ. f θ = DOI:.979/ Page df θ dθ θ α F x, θ, B θ u x, θ dθ. (.9) Le C H (H, T )be he se of all connuous funcons u(x, ) H, and defne on C H (H, T )a norm by u(x, ) CH (H,T ) = max u(x, ) H, for all, T, x H. x, By a soluon of he PDE((.),(.2)), we mean an absrac funconu(x, ) such ha he followng condons are sasfed: (a)u x, C H H, T and u x, S for all, T, x H. (b) α u (x,) exss and connuous on H, T, where < α <. α (c)u x, sasfes Eq.(.) wh he nal condon (.2) on H, T. Lemma (whou proof):if λ > and < δ <, hen e λη η δ dη λ δ, + δ eλ. (.) II. The exsence and unqueness soluon of PDE of fraconal order Here, he exsence and unqueness soluon of Eq. (.9) and consequenly s equvalen PDE ((.),(.2)),wll be dscussed and proved n a Banach space H by vrue of Pcard s mehod. The formula (.9) s equvalen o he followng negral equaon(see[4]) u x, = ξ α θ Q α θ u x dθ + α θ η α ξ α θ Q η α θ w x, η dθ dη, (2.) whereξ α θ s a probably densy funcon defned on, ), w x, = F x,, B u x, = F x,, W x,, (2.2) and W x, = ξ α θ B Q α θ u x dθ +α θ η α ξ α θ B Q η α θ w x, η dθdη. (2.3) Theorem: The PDE((.),(.2)) has a unque soluon n he Banach spacec H H, T. The proof of hs heorem comes as a resul of he followng lemmas. Lemma 2: Under he condons (.5) and (.7), he negral equaon (2.) has a soluon nc H H, T. Proof: Usng he mehod of successve approxmaons, he formulas (2.2) and (2.3),lead o w n+ x, = F(x,, ξ α θ B Q α θ u x dθ Hence, n vew of he condon (.7), we ge +α θ η α ξ α θ B Q η α θ w n x, η dθ dη ).

3 On he numercal reamen of he nonlnear paral dfferenal equaon of fraconal order w n+ x, w n x, α l θ η α ξ α θ B()Q η α θ ) w n x, η w n x, η dθdη. Usng he condon (.5), we oban w n+ x, w n x, Inroducng (.) n (2.4),we have M max x, e λ +x w n x, w n x, e λ η +x η v dη, v = α γ, M = αl max e λ +x w n+ x, w n x, M x, λ We can choose suffcenly large such ha M + λ v Hence, he above nequaly can be adaped n he form max e λ +x w n+ x, w n x, μ max x, x, By a successve applcaon of he above nequaly, we ge max x, e λ +x w n+ x, w n x, μ max x, μ 2 max x, v v θ γ ξ α θ dθ DOI:.979/ Page + v max x, (2.4) and λ >. (2.5) e λ +x w n x, w n x,. = μ <. (2.6) e λ +x w n x, w n x,. λ +x e w e λ +x w n x, w n 2 x, n x, w n x, μ n max e λ +x w x, w x,, x, wherew x, s he zero approxmaon whch can be aken he zero elemen n hespace H. Thus, he seres n= w n+ x, w n x, n converges unformly n H, T. Snce w n+ x, = = w + x, w x,, follows ha he sequence w n x, converges unformly n he space C H (H, T ) o a connuous funcon F x,, W x, whch sasfes Eq.(2.) for all(x, ) H, T. Consequenly, u(x, ) C H (H, T ), where u x, = ξ θ Q α θ u x dθ + α θ η α ξ α θ Q η α θ F x, η, W(x, η) dθdη. Lemma 3: Under he condons (.5) and (.7), he negral equaon (2.3) has a unquesoluon n C H (H, T ). Proof: Le u x, and u 2 x, be wo soluons of Eq.(2.), hen from he formulas (2.2) and (2.3) wh he ad of condon (.7), we have w 2 x, w x, α l θ η α ξ α θ B Q η α θ w 2 x, η w x, η dθdη. Usng he same argumen of lemma (2), we ge e λ +x w 2 x, w x, H μρ, max x, where ρ = max x, e λ +x w 2 x, w x, H. Thus, from (2.6) we have ρ = max e λ +x w 2 x, w x, H =. x, Ths complees he proof of he lemma. Lemma 4(whou proof): Under he condons (.4), (.5) and (.6), he soluon u x, of Eq.(2.) sasfes a unform Hölder condon.(see [4]) Proof of Theorem :By vrue oflemmas (2), (3)and (4), we deduce ha he soluon u x, of Eq.(2.) represens he unque soluon of he PDE((.),(.2)) n he Banach space C H (H, T ), and u x, S. ow, we wll prove he sably of he soluons of he PDE((.),(.2)).In oher words, we wll show ha he PDE((.), (.2))s correcly formulaed. Theorem 2 : Le u n x, be a sequence of funcons, each of whch s a soluon of Eq.(.) wh he nal condon u n x, = g n x, where g n x S n =,2,. If he sequence g n x converges o an elemen

4 On he numercal reamen of he nonlnear paral dfferenal equaon of fraconal order u x S, he sequence Ag n x converges and he sequence B()g n x converges unformly on H, T. Then, he sequence of soluons u n x, converges unformly on H, T o a lm funcon u x,, whch s he soluon of he PDE((.),(.2)). Proof: Consder he sequences f n x, and u n x,, where α u n x, α Au n x, = f n x,, u n x, = u n x, g n x, u n x, = g n x, u n x, = α θ η α ξ α θ Q η α θ f n x, η dθdη, and f n x, = F x,, B u n x, + B g n x + Ag n x. Thus, we ge f n x, f m x, Ag n x Ag m x + F x,, B u n x, + B g n x F x,, B u m x, + B g m x. Usng he condon (.7), we oban f n x, f m x, l B u n x, u m x, +l B g n x B g m x + Ag n x Ag m x. Consequenly, f n x, f m x, αl θ η α ξ α θ B Q η α θ f n x, η f m x, η dθdη + l B g n x B g m x + Ag n x Ag m x. In vew of he condons (.5) and (2.5), he above nequaly becomes f n x, f m x, M η v f n x, η f m x, η dη + l B g n x B g m x + A g n x A g m x. Gven ε >, we can fnd a posve neger = (ε) such ha f n x, f m x, M η v f n x, η f m x, η dη + μ 2 ε, for all n, m and (x, ) H, T. Usng (.), he above nequaly akes he form μ e λ(+x) f n x, f m x, μ e λ +x ε. Thus, for suffcenly large λ, we ge max e λ(+x) f n x, f m x, x, ε. Snce H s a complee space, follows ha he sequence f n x, converges unformly on H, T o a connuous funcon f x,, so he sequence u n x, converges unformly on H, T o a connuous funconu x,. I can be proved ha f x, sasfes a unform Hölder condon on, T, hus u x, S. Corollary : The negral equaon (2.) has a unque soluon n he Banach space C R (R, T ). III. The numercal soluon of he PDEof fraconal order In hs secon, we wll usehe (MTMM)o oban numercally, he soluon of he PDE((.),(.2))n he Banach space C R (R, T ), where u x, CR R,T = max u x,,, T, < x <. x, For hs, we wre Eq.(2.2) n he form f x, = u x, = f x, + α p, η Q, η F x, η, B η u x, η dη, ξ α θ Q α θ u x dθ, Q, η = θξ α θ Q η α θ dθ, and he bad kernel p, η = η α, < α <, η T ; T <. (3.) DOI:.979/ Page

5 On he numercal reamen of he nonlnear paral dfferenal equaon of fraconal order Here, he unknown funcon u(x, ) C R (R, T ), whle f x,, Q, η and p, η are known funcons and sasfy he followng condons: ( ) f x, s a connuous funcon n (R, T ). (2 ) Q, η wh s paral dervaves are connuous funcon n, T, hence here exss a consan c such ha: Q, η c. (3 ) p, η s a badly behaved funcon of s argumens such ha: (a) for each connuous funcon u x, and 2, he negrals 2 p, η Q, η F x, η, B η u x, η dη, and p, η Q, η F x, η, B η u x, η dη are connuous funcons n (R, T ). (b) p, η s absoluely negrable wh respec o ηfor all T, hus here exss a consan c 2,such ha: p, η dη c 2. (4 )The gven funconf x,, B u x, s connuous nr, T, and sasfes Lpschz condon F x,, B u x, F x,, B v x, x, B u x, v x,, max x, for all u x,, v x, C R (R, T ), wherel s a posve consan. Pung x = x, x =, = x + x, and usng he followng noaons u, x = u, f, x = f, F x, η, B η u x, η = F η, B η u η, he negral equaon (3.) can be ransformed o he followng (SVIEs) u = f + α p, η F η, B η u η d η p, η = p, η Q, η. Remark : le Ξ be he se of all connuous funcons U = u, u 2,, u,, where u ε, T, for all, and defne on Ξ he norm: x, L,, (3.2) Then Ξ s a Banach space. U Ξ = max u = T u,t,. 3.. The exsence of a unque soluon of asvies: In order o guaranee he exsence of a unque soluon of he SVIEs (3.2) n he Banach space Ξ, we wre hs sysem n he negral operaor form L u = f + Lu, (3.3) L u = α p, η F η, B η u η d η. (3.4) Then assume he followng condons: ) max f = f Ξ D, D s consan. T 2 ) The known connuous funcons F, B u for all, sasfy for he consans A, A 2 and A max A, A 2, he followng condons: a max F, B u A U Ξ, T a 2 F, B u () F, B u (2) 2 A 2 u () u (). 3 ) The kernelp, η s a dsconnuous funcon whch sasfes: b foreachconnuous funcon F, B u and 2, he negrals: 2 p, η F η, B η u η d η, and p, η F η, B η u η d η are connuous n, T. (b 2 )p, η s absoluely negrable wh respec o η for all T, hus here exss a consan c, such ha: p, η d η c. Theorem 3 (whou proof): The formula (3.2) has a unque soluon n he space Ξ under he followng condon:δ = α A c <. 3.2.The modfed Toeplz marx mehod (MTMM): DOI:.979/ Page

6 On he numercal reamen of he nonlnear paral dfferenal equaon of fraconal order Here, we presen he MTMM o oban he numercal soluon of a VIE of he second knd wh sngular kernel. So, we assume he VIE: u = f + α p, η F η, B η u η dη. (3.5) Followng he same way of Abdou eal.,(see [2],[3]), we can apply he MTMM for Volerra erm o oban he followng equaon u α n= D n F n, B n u n = f. (3.6) Pung = m, = T, n 3.6 and usng he followng noaons u m = u m, D n m =, f m = f m, F m, B m u m = F m B m u m, (3.7) we ge he followng AS and u m α = n= F n B n u n = f m, n m, (3.8) A m, n = A n m + B n m, < n < B m, n =, DOI:.979/ Page (3.9) A n = F n +, n + B n + I F(n +, B n + J(), (3.) n + B n = F n, B n J() F n, n B n I(), (3.) = F n, B n F (n +, n + B(n + )) F (n, n B(n )) F(n +, B(n + ), (3.2) n +. (3.3) I = P, η F η, B η dη, J = p, η F η, ηb η dη n n The marx can be wren n he Toeplz marx form: = G mn E mn. Here, he marx G mn = A n m + B n m, n m, (3.4) s called he Toeplz marx of order + and m, n = E mn = B < n < (3.5) A m, n =, represens a marx of order ( + ) whose elemens are zeros excep he frs and he las rows (columns). Defnon : The Toeplz marx mehod s sad o be convergen of order r n he nerval [, T], f and only f for suffcenly large, here exss a consan d > ndependen on such ha u u d r. (3.6) Defnon 2: The esmae local error R akes he form R = P, η F (η, B (η) u η dη F n (B n u n ). (3.7) n= Lemma 5: If he kernelp, η of Eq. (3.5) sasfes condon (3 ) and he followng condon lm p, η p, η dη = ;,, T, (3.8) hen n= exss, and lm m m n= D m n =.

7 On he numercal reamen of he nonlnear paral dfferenal equaon of fraconal order Proof: From he formulas (3.) and (3.3), we ge A n () F n +, n + B n + n + n + n p (, η) F(η, B(η) dη + F(n +, B(n + ) p (, η) F(η, η B(η) dη. n Summng fromn = o n =, hen akng n accoun he connuy of he funcon F, B u nerval, T and fnally usng he condon (3 ),here exss a small consan M, Snce, each erm of n= A n () M,. M = 2L c ; F(, B u() L n= A n s bounded above, hence for = m, we deduce A n (m ) n= Smlarly, from he formulas (3.) and (3.3), we have n= B n m n he M. (3.9) M. (3.2) In he lgh of (3.9), and wh he help of (3.9) and (3.2), here exss a small consan M 2, such ha Hence, n= n= exss. n= A n m + n= By vrue of he formulas (3.) and (3.3), we have for, ϵ, T A n A n F n +, n + B n + n + B n m p, η p, η F η, η B η dη n + n + M 2, M 2 = 2M. F n +, B n + p, η p, η F(η, η B η ) dη. n Summng from n = o n =, and akng n accoun he connuy of he funcon F, he above nequaly can be adaped n he form n= A n A n 2 L 2 Pung = m, = m, hen usng he condon (3.9), we ge lm m m Smlarly, n vew of he formulas (3.) and (3.3), we can prove ha lm m m n= n= Fnally, wh he ad of (3.9), (3.2) and (3.22), we have lm m m P, η P, η dη. A n m A n m =. (3.2) n= B n m B n m =. (3.22) D m n =. IV. The exsence of a unque soluon of AS The SVIEs(3.2) afer usng MTMM akes he form DOI:.979/ Page

8 On he numercal reamen of he nonlnear paral dfferenal equaon of fraconal order, u,m = f,m + α F,n B,n u,n. n= Lemma 6 (whou proof):if he kernel p, η of Eq. (3.2) sasfes he condons (3 ) and (3.8), hen we have n= exss, and lm m m D n= m D n mn =. Accordng o lemma (5) and n order o guaranee he exsence of a unque soluon of he AS (4.) n he Banach space l, we wre hs sysem n he followng operaor form Vu,m = f,m + Vu,m, (, m ) (4.2) Then we assume he followng condons: I),m Vu,m = α n= f,m H, H s a consan. II), n= (4.) F,n B,n u,n. (4.3) e, e s a consan. III) The known connuous funcons F,n B,n u,n sasfy (, n)and for he consans q, q 2 and q max q, q 2 he followng condons: (a ) F,n B,n u,n q u,n,n l, 2 2 (b ) F,n B,n u,n F,n B,n u,n q 2 u,n u,n. Theorem 4 (whou proof): The formula (4.2) has a unque soluon n he Banach space l under he followng condon:σ = α q e <. Defnon 3: The esmae oal error R j s deermned by he relaon when j = max,, hen he sum R j = p, η F x, η, B η u x, η dη F,n n= F,n B,n u,n n= B,n u,n, ends o P, η F x, η, B η u x, η dη. Theorem 5: If he sequence of connuous funcons f n x, converges unformly o he funcon f x,, and he funcons Q, η, p, η and F x,, B u x, sasfy, respecvely, he condons (2 ), (3 b) and (4 ). Then, he sequence of approxmae soluons u n x, converges unformly o he exac soluon of Eq.(3.) n he Banach space C R (R, T ). Proof: The formula (3.) wh s approxmae soluon gve u x, u n x, max f x, f n x, max x, +α x, p, η Q, η max x,η F, η, B η u x, η F, η, B η u n x, η dη, η T, < x <. Snce B s a bounded operaor, here exss a posve consan M, such ha B u x, M u x,. (4.4) In vew of he condons (4.4), (2 ), (3 b), and (4 ), he above nequaly can be adaped n he form u x, u n x, CR R,T D f x, f n x, CR R,T ; D = αc c 2 ML. Snce f x, f n x, CR R,T as n, so ha u x, u n x, CR R,T. Theorem 6 The oal error R j sasfes lm R j =. j Proof: From he defnon of R j, we have. DOI:.979/ Page

9 On he numercal reamen of he nonlnear paral dfferenal equaon of fraconal order R j,j u m u m j +, n=,j Usng he condons (I) and (III - b ), we ge F n, B n u n F n, B n u n R j + eq u m u m j l, j. Snce each erm of R j s bounded above, hence for = m, we deduce R j + eq max u u j j = + eq u u j Ξ. The above nequaly can be adaped n he form R j l + eq u x, u x, j. C R (R,T ) j. Snce u x, u x, j C R (R,T ) as j, hen R j l, and consequenly R j. V. Applcaon In Eq. (3.) le: < α <, Q, η =, F x,, B u x, = x , hus, we ge a VIE of he second knd wh Abel kernel u x, = f x, + α η α x 2 + η 2 3 dη, (5.) where he exac soluon: u x, = x Usng Maple 2, he resuls are obaned numercally for x =,.8, wh α =.98,.8,.6,.4,. and.2.thenerval,.8 s dvded no = 4 unes. The followng ables and dagrams are seleced from a large amoun of daa o compare beweenheexacsoluonof Eq. (5.) (Exac.sol.) and s numercal soluon (Approx.sol.). x = u E α =.98 α =.8 α =.4 u T E T u T E T u T E T E E E E E E E E Table () Dg. () x = u E α =.6 α =. α =.2 u T E T u T E T u T E T E E E E E E E E E E E E DOI:.979/ Page

10 On he numercal reamen of he nonlnear paral dfferenal equaon of fraconal order Table (2) Dg. (2) VI. Conclusons From hs paper, we can conclude he followng pons: ) The PDE ((.),(.2)) of fraconal order s equvalen o he followng VIE: u x, = u x + Γ α θ α A u x, θ dθ + Γ α θ α F x, θ, B θ u x, θ dθ. From he above negral equaon, we can dscuss many cases for he nonlnear negro- dfferenal equaon f for example F x,, B u x, = f ( u( x, ), u ( x, ), u ( x, )). 2) The MTMM, as a bes mehod o solve he sngular negral equaons, s used o oban a AS, and many heorems are derved o prove he exsence and unqueness of heas. 3) From he numercal resuls, we esablsh he followng: a) For fxed values of α, he error values are ncreasng wh he ncrease values of x and. b) For fxed values of x and,he error values are slowly ncreasng wh he decrease values of α.so, he change n he values of α s lghly effecve n he numercal calculaons. c) The error s a x = = for all cases we have suded. d) The maxmum value of he error a α =.98 s , whle he maxmum value of he error a α =.2 s , for x = =.8. References []. S.Chandrasekhar and G.Münch, The heory of flucuaon n brghness of he Mlky-Way, Asrophy.J.25, [2]. M.C.Mackey and Rudnck, Global sably n a delayed paral dfferenal equaon descrbng cellular replcaon, J.Mah.Bol,33(994), [3]. J.A.J.Mez and O.Dekmann, The Dynamc of Physologcally Srucured Populaons Sprnger lecure oes n Bomahemacs, Vol.68, Sprnger-Verlag, ew York,986. [4]. M.M.EL-Bora, Some probably denses and fundamenal soluons of fraconal evoluon equaon, Chaos, Solons & Fracals 4(22), [5]. M.M.El-Bora, Semgroups and some nonlnear fraconal dfferenal equaons, J.Appl.Mah. Compu. 49 (24), [6]. M.A.Abdou, M.M.El- Kojok &S.A.Raad, Analyc and umerc Soluon of Lnear Paral Dfferenal Equaon of Fraconal Order,Global.J.Inc.Vol.3(23), [7]. S.D. Zadman, Absrac Dfferenal Equaons, Research oes n Mahemacs,Pman Advanced Publshng Program, San Francsco, London Melbourne, 979. [8]. [8] M.M.El-Bora, Semgroups and ncorrec problems for evoluon equaons, 6 h IMACS WORLD COGRESS LAUSAE, Augus 2.Froner Research" -F Vo. 3 Issue 3 Verson. Aprl 23. [9]. M.M.El-Bora, O.L. Mousafa and F.H.Mcheal, On he correc formulaon of a nonlnear dfferenal equaons n Banach space, In.J.Mah. 22 () (999). []. I.M. Gelfand and G.E.Shlov, Generalzed Funcons, Vol., auka, Moscow,959. []. W.R.Schneder and W.Wayes, Fraconal dffuson and wave equaon,j.mah. Phys. 3 (989). [2]. M.A. Abdou, M.M. El-Bora, M.M. El-Kojok, Toeplz marx mehod and nonlnear negral equaon of Hammersen ype, j.com. Appl. Mah. 223 (29) [3]. M.A. Abdou, A.A. Badr, M.M. El-Kojok, On he soluon of a mxed nonlnear negral equaon, j.com. Appl. Mah. 27 (2) DOI:.979/ Page