arxiv: v1 [math.ap] 4 Aug 2017

Transcription

1 New Two Step Laplace Adam-Basfort Metod for Integer an Non integer Order Partial Differential Equations arxiv: v1 [mat.ap] 4 Aug 17 Abstract Rodrigue Gnitcogna*, Abdon Atangana** *Department of Matematics, Faculty of Natural and Agricultural Science University of te Free State Bloemfontein. *Department of Matematics, University of Namibia. **Institute of Groundwater Studies, University of te Free State Bloemfontein Sout Africa. Tis paper presents a novel metod tat allows to generalise te use of te Adam-Basfort to Partial Differential Equations wit local and non local operator. Te Metod derives a two step Adam-Basfort numerical sceme in Laplace space and te solution is taken back into te real space via inverse Laplace transform. Te metod yields a powerful numerical algoritm for fractional order derivative were te usually very difficult to manage summation in te numerical sceme disappears. Error Analysis of te metod is also presented. Applications of te metod and numerical simulations are presented on a wave-equation like, and on a fractional order diffusion equation. Keywords: Laplace Transform, Adam-Basfort, PDE, Fractional PDE, Caputo Derivative. 1. Introduction Adam-Basfort metod as been recognised as a powerful numerical tool to solve Partial Differential equations P.D.E [1]. It is a numerical sceme tat is used in many field of applied science. epidemiology, engineering in dynamical systems, in caotic problems []-[6]. Te numerical sceme Preprint submitted to Elsevier October 19, 17

2 is good for bot differential equations wit classical derivatives, and differential equations wit non integer order derivatives. However Te metod is not fully applied to P.D.E wit local and non local operator as it was designed only for Ordinary Differential Equations O.D.E [1]. Various oter metods are used instead for P.D.E wit integer order differentiation, and tose wit real order derivatives[7]-[15]. Noneteless, Due to te accuracy and efficiency of Adam-Basfort tecniques, tere is a need to extend te metods to P.D.E. By eliminating one variable and transforming a P.D.E to an O.D.E via Laplace transform, te newly obtained O.D.E can be analysed in Laplace space, and a furter application of te inverse transform will return te solution to te real space. Tis paper derives a new numerical sceme tat combines te Laplace transform and Adam-Basfort metod to andle P.D.E wit integer order and non integer order derivatives.. Numerical Metod for P.D.E wit Integer order Consider te following general P.D.E ux, t t = Lux, t + Nux, t.1 Were L is a linear operator and N a non linear operator. We start by applying Laplace transform on bot sides of te equation, wit respect to te variable x to obtain ux, t L = L Lux, t + Nux, t = t d up, t = L Lux, t + Nux, t dt d ut = F u, t. dt Were ut = up, t and F u, t = L Lux, t + Nux, t. Next we apply te fundamental teorem of calculus on equation. to obtain ut = ut + t F u, τdτ

3 tis is also Wen t = t n+1 we ave ut = u + t F u, τdτ u n+1 = ut n+1 = u + Wen t = t n u n = ut n = u + it follows tat u n+1 u n = +1 u n+1 u n = +1 F u, τdτ +1 t n F u, τdτ F u, τdτ F u, τdτ F u, τdτ If we approximate F u, t wit te Lagrange polynomial P t F u, t = t t n 1 F u, t n + t t n F u, t n 1 t n t n 1 t n 1 t n we can terefore write u n+1 u n = u n+1 u n = By letting u n+1 u n = F n t n t n 1 P t = t t n 1 F n + t t n F n 1 t n t n 1 t n 1 t n F n t n t n 1 u n+1 u n = +1 t n t n t tn 1 t n t n 1 F n + F u, τdτ F n 1 t tn 1 dt + t n t n 1 t n [ ] tn+1 1 t tt n 1 = t n t n 1 3 t t n F n 1 dt t n 1 t n + F n 1 t n t n 1 t n +1 t n tn 1 t n dt [ 1 t tt n ] tn+1 t n

4 we ave u n+1 u n = F n 1 t n+1 t n+1 t n 1 1 t n + t n t n 1 F n 1 1 t n+1 t n t n+1 1 t n + t n u n+1 u n = F n 1 t n+1 t n t n+1 +t n t n 1 t n+1 t n F n 1 u n+1 u n = F n 1 t n+1 + t n t n 1 F n 1 1 t n+1 + t n t n 1 1 u n+1 u n = F n t n+1 + t n t n 1 F n 1 t n+1 + t n t n 1 t n+1 t n t n+1 +t n t n t n+1 t n u n+1 u n = F n 1 1 n n n 1 F n 1 n n n u n+1 u n = F n n + 1 n + F n 1 n + 1 n 3 u n+1 = u n + F n 1 F n 1.3 Applying te inverse Laplace transform to return into te real space, we ave: 3 L 1 u n+1 = L [u 1 n + F n 1 ] F n 1 3 ux, t = L [u 1 n + F n 1 ] F n 1 3. New Numerical Metod for P.D.E wit non integer order To illustrate te metod we consider te general fractional P.D.E ux, t t = Lux, t + Nux, t 3.1 4

5 Were L is a linear operator and N a non linear operator. Applying Laplace transform on bot sides of te equation 3.1, we ave ux, t L = L Lux, t + Nux, t t For te Caputo type fractional partial derivative tis will be C ad t up, t = L Lux, t + Nux, t C ad t up, t = F u, t tis is C ad t ut = F u, t 3. were ut = up, t and F u, t = L Lux, t + Nux, t. Te next step is to apply te Caputo fractional Integral operator on equation 3.. Doing so we obtain Wen t = t n+1 Wen t = t n ut ut = 1 Γ u n+1 = ut n+1 = u + 1 Γ u n = ut n = u + 1 Γ t +1 t τ 1 F u, τdτ t n+1 τ 1 F u, τdτ t n τ 1 F u, τdτ u n+1 u n = 1 [ t n+1 ] t n+1 τ 1 F u, τdτ t n τ 1 F u, τdτ Γ +1 n tj+1 t n+1 τ 1 F u, τdτ = t n+1 τ 1 F u, τdτ We approximate F u, t wit te following Lagrange polynomial j= t j 3.3 P t F u, t = t t n 1 F u, t n + t t n F u, t n 1 t n t n 1 t n 1 t n 5

6 P t = t t n 1 F n + t t n F n 1 t n t n 1 t n 1 t n Te first fractional integral in equation 3.3 can ten be expressed as +1 = = n [ t n+1 τ 1 F u, τdτ = F n t n t n 1 j= n [ tj+1 Fn j= t j tj+1 n j= tj+1 t j t j t n+1 t 1 t t n 1 dt+ F n 1 t n 1 t n t n+1 t 1 t t n 1 dt F n 1 t t n+1 t 1 tn 1 F n + t t n F n 1 dt t n t n 1 t n 1 t n tj+1 We now implement te following cange of variable. We let y = t n+1 t, dt = dy, t = t n+1 y t j tj+1 t j ] t n+1 t 1 t t n dt ] t n+1 t 1 t t n dt tj+1 t j t n+1 t 1 t t n 1 dt = +1 t j+1 t n+1 t j y 1 y + t n+1 t n 1 dy = +1 t j+1 t n+1 t j y y 1 dy = 1 ] [y +1 tn+1 t j+1 ] [y tn+1 t j+1 t n+1 t j t n+1 t j = 1 t n+1 t j+1 +1 t n+1 t j +1 t n+1 t j+1 t n+1 t j. On te oter and tj+1 t j t n+1 t 1 t t n dt = +1 t j+1 t n+1 t j y 1 y + t n+1 t n dy = 1 = +1 t j+1 t n+1 t j y y 1 dy = 1 ] [y +1 tn+1 t j+1 ] [y tn+1 t j+1 t n+1 t j t n+1 t j t n+1 t j+1 t n+1 t j t n+1 t j+1 +1 t n+1 t j +1 6

7 It will follow tat +1 t n+1 τ 1 F u, τdτ = F n { n j= [ 1 t ] n+1 t j+1 +1 t n+1 t j +1 n [ ]} t n+1 t j+1 t n+1 t j j= F n 1 { n j= [ 1 t ] n+1 t j+1 +1 t n+1 t j +1 n [ ]} t n+1 t j+1 t n+1 t j j= = F n 1 t n+1 t +1 t n+1 t F n 1 1 t n+1 t +1 t n+1 t = F n n F n 1 n n n + 1 [ = n + 1 n n + 1 n F n F n 1 ]. 7

8 Te second fractional integral in equation 3.3 will similarly be evaluated as n 1 t n τ 1 F u, τdτ = j= tj+1 t n τ 1 F u, τdτ = F n 1 n t j j= t t n t 1 tn 1 F n + t t n F n 1 dt t n t n 1 t n 1 t n tj+1 t n t 1 t t n 1 dt F n 1 t j We also implement te following cange of variable. We let y = t n t, dt = dy, t = t n y t n τ 1 F u, τdτ = F n 1 n j= t n τ 1 F u, τdτ = F n 1 n j= j= t j+1 n 1 tj+1 j= y 1 t n t n 1 ydy F n 1 t n t j t j+1 y y 1 dy F n 1 t n t j t n τ 1 F u, τdτ = F n 1 [ n y +1 ]tn tj+1 y F n 1 t n t j t j t n t 1 t t n dt n 1 tj+1 j= n 1 tj+1 j= t n t j n 1 [ y +1 j= y dy t n t j ]t n t j+1 t n t j y dy = F n n 1 tn t j+1 +1 j= t n t j+1 t n t j +1 + t n t j F n 1 n 1 tn t j+1 +1 j= t n t j +1 = F n { n 1 tn t j+1 +1 j= t n t j +1 n 1 F n 1 j= t n t j+1 t n t j } n 1 tn t j+1 +1 j= t n t j +1 = F { n t n t +1 } t n t F n 1 tn t +1 8

9 = F n t n t +1 + t n t + F n 1 t n t +1 = F n n n +1 + n+1 +1 F n 1 Terefore we can ten write n t n τ 1 F u, τdτ = F n+1 n + n+1 F n 1 Rewriting equation 3.3 substituting in te later results we can ten ave [ u n+1 u n = n + 1 n F n Γ n + 1 n n ] F n 1 F n+1 n + n+1 F n 1 u n+1 u n = [ n + 1 n Γ n n+1 n n+1 n F n 3.4 F n 1 ] Remark 3.1. For = 1 we recover te classic Adam-Masfort Numerical sceme. In fact u n+1 u n = n n + 1 n + n + Fn n + 1+n n + 1 F n 1 Γ1 u n+1 u n = n n 1 n n 1 Fn n + 1+n n n 1 u n+1 u n = 3 F n 1 F n 1 To find te numerical Sceme in te real space we need to apply te inverse Laplace 3.4. We obtain te following iterative sceme in te real space: { [ u n+1 x, t = L 1 u n + n + 1 n + n+1 n F n Γ n n+1 n Fn 1 F n 1 ] }

10 Te above equation can be discretised in x using any classical metod, including but not limited to forward, backward difference, Crank Nicolson. 4. Error Analysis of te Laplace Adam-Basfort Metod Let C D ux, t = Lux, t + Nux, t 4.1 be a general fractional partial differential equation. As we establised earlier te numerical solutions using Laplace Adam-Basfort metod is given as { [ u n+1 x, t = L 1 u n + n + 1 n + n+1 n F n 4. Γ n n+1 n ] } F n 1 + Rn Were R n < Proof. Following te derivation presented earlier u n+1 u n = 1 [ t n+1 t n+1 τ 1 F u, τdτ Γ Were F u, ζ = t t n 1 F n + t t n F n 1 + F u, ζ t n t n 1 t n 1 t n! u n+1 u n = 1 Γ 1 Γ + 1 Γ +1 ] t n τ 1 F u, τdτ 1 t t i i= t t n+1 t 1 tn 1 F n + t t n F n 1 dt t n t n 1 t n 1 t n t t n t 1 tn 1 F n + t t n F n 1 dt t n t n 1 t n 1 t n +1 F u, ζ! 1 F u, ζ Γ! 1 t t i t n+1 t 1 dt i= 1 t t i t n t 1 dt i= 1

11 Tis is equal to [ u n+1 u n = n + 1 n + n+1 n F n Γ n n+1 n ] F n 1 + Rn Terefore one can easily deduce R n = 1 Γ +1 Rn 1 +1 Γ R n 8Γ + F u, ζ! F u, ζ! F u, ζ! F u, ζ! 1 t t i t n+1 t 1 dt i= 1 i= t t i t n t 1 dt 1 t t i t n+1 t 1 dt i= 1 t t i t n t 1 dt i= +1 max {F u, ζ} ζ,t n+1 + tn+1 t 1 dt t n t 1 dt Rn t max 8Γ {F u, ζ} n+1 + t n ζ,t n+1 Rn max 8Γ {F u, ζ}n n < + ζ,t n+1 5. Applications In tis section we apply te metod to some integer and non integer order P.D.E. We start wit te classical wave equation. 11

12 5.1. Example 1: Integer Order P.D.E Consider te following wave-equation like: ux, t t = c ux, t x Applying Laplace transform on bot sides we get te following 5.1 dux, t dt = c p up, t u, t Silencing te variable p writing up, t = ut te later equation can be rewritten as dut = F t, ut dt Te later equation is.. We proved earlier tat te solution in te Laplace space to te previous equation is given by.3. 3 u n+1 = u n + F n 1 F n 1 Applying te inverse transform to.3 we ave : ux, t n+1 = ux, t n + 3 c ux, t n x Discretising in te space variable we ave ux i, t n+1 = ux i, t n + 3 c [ uxi+1, t n ux i, t n x Leting ux i, t n = u n i, and x = l c ux, t n 1 x 5. ] [ ] c uxi+1, t n 1 ux i, t n 1 x u n+1 i u n+1 i = u n i + 3 cun i+1 u n i l 1 cun 1 i+1 un 1 i l = 1 3c l un i + 3c l un i+1 c l un 1 i+1 + c l un 1 i 5.3 1

13 Stability Analysis Assume we ave a Fourier expansion in space of ux, t = f ût expjf x Equation 5.3 becomes û n+1 e jifl = 1 3c l û ne jifl + 3c l ûne ji+1fl c l ûn 1e ji+1fl + c l tis is û n+1 û n û n+1 = 1 3c l û n + 3c l ûne jfl c l ûn 1e jfl + c l ûn 1 û n+1 = 1 3c l û n+1 = 1 3c û n l + 3c Finding te ration we ave: û n+1 û n û n+1 û n = 1 3c l = 1 3c 3c 1 cosfl +j l = 1 3c l + 3c l ejfl û n + c l c l ejfl û n 1 l ejfl c l + 3c l ejfl + c l sinfl+c l fl sin + c fl sin l û n 1 û n l c l e jfl + c l ejfl ûn 1 û n û n 1 û n ûn 1e jifl 1 cos jflûn 1 j c sin flûn 1 û n l û n ûn 1 Let us prove tat n u n < u. For n = û 1 = 1 3c fl sin l û Tis means tat we ave û 1 3c fl û < 1 1 < 1 sin < 1 l = < 3c fl sin < l 13 û n + j c ûn 1 sinfl3 l û n

14 = < 3c l fl sin < = < 3c fl sin < 1 4l Te previous condition will certainly be acieved if since. < 3c 4l < 1 < 3c fl sin 4l 3c < 1 = 4l 3c 4l. l < 4 3 c 5.4 Now let us assume j n, û n < û and prove tat û n+1 < û Te condition 5.4 ensures tat te coefficients in te following equation obtained earlier above are positive. û n+1 = 1 3c l û n + 3c l ûne jfl c l ûn 1e jfl + c. Te later ten implies tat l ûn 1 û n+1 < 1 3c l û n + 3c l û n e jfl c l û n 1 e jfl + c l û n 1. Using te induction ypotesis tis becomes û n+1 < 1 3c l û + 3c l û c l û + c l û. and terefore we obtain û n+1 < 1 3c + 3c c l l l + c û l. û n+1 < û. We can terefore conclude tat te numerical sceme solution presented above is stable. 14

15 5.. Example: Fractional Order P.D.E Consider te following fractional Order P.D.E, te fractional derivative is given in te Caputo Sense. ux, t t = d ux, t x 5.5 Applying Laplace transform on bot sides we get te following { C D ux, t } t up, t = d L x = d L p up, t p u, t u, t Silencing te variable p writing up, t = ut te later equation can be rewritten as C D t up, t = d F t, ut Te later equation is 3.. We proved earlier tat te solution in te Laplace space to te previous equation is given by 3.4 u n+1 u n = [ n + 1 n Γ n n+1 n Applying te inverse transform to 3.4 we ave : ux, t n+1 = ux, t n + + n+1 n F n F n 1 ] Γ δ nd ux, t n x Γ δ,1 n d ux, t n 1 x were and δ n = n + 1 n δ,1 n = n n+1 n n+1 n

16 Discretising in te space variable we ave ux i, t n+1 = ux i, t n + Γ δ nd ux i+1, t n ux i, t n + ux i 1, t n x Γ δ,1 n d ux i+1, t n 1 ux i, t n 1 + ux i 1, t n 1 x Leting ux i, t n = u n i, and x = l te equation becomes Tis is u n+1 i = u n+1 i = u n i + Γ δ nd un i+1 u n i + u n i 1 l Γ δ,1 n 1 d l Γ δ n u n i + d l Γ δ,1 n u n 1 i+1 + d l Γ δ,1 n d un 1 i+1 un 1 i + u n 1 i 1 l d l Γ δ nu n i+1 + u n 1 d l Γ δ nu n i 1 i d l Γ δ,1 n u n 1 i Stability Analysis of te numerical sceme solution fractional Order P.D.E Assume we ave a Fourier expansion in space of ux, t = f ût expjf x wile leting u n i Equation 5.6 becomes û n+1 e jfil = Tis is 1 d l Γ δ n = û n exp jfi x = û n exp jfil û n e jfil + d l Γ δ nû n e jfi+1l + d l Γ δ nû n e jfi 1l d l Γ δ,1 n û n 1 e jfi+1l + d l Γ δ,1 n û n 1 e jfil d l Γ δ,1 n û n 1 e jfi 1l û n+1 = 1 d l Γ δ n û n + d l Γ δ nû n e jfl + d l Γ δ nû n e jfl d l Γ δ,1 n û n 1 e jfl + d l Γ δ,1 n û n 1 d l Γ δ,1 n û n 1 e jfl 16

17 Wic implies û n+1 = 1 d l Γ δ n û n + d l Γ δ,1 n û n 1 + d l Γ δ nû n e jfl + e jfl d l Γ δ,1 n û n 1 e jfl + e jfl Furter we ave û n+1 = 1 d l Γ δ n û n + d l Γ δ,1 n û n 1 + d l Γ δ nû n cos fl d l Γ δ,1 n û n 1 cos fl and û n+1 = 1 d l Γ δ n + d d l Γ δ n cos fl û n + l Γ δ,1 n d l Γ δ,1 n cos fl û n 1 Tus we can write û n+1 = 1 d d l Γ δ n 1 cos fl û n + l Γ δ,1 n 1 cosfl û n From equation 5.7 we ave û n+1 û n = 1 d d û n 1 l Γ δ n 1 cos fl + l Γ δ,1 n 1 cosfl û n = 1 4 d l Γ δ n sin fl 4 d + l Γ δ,1 n fl ûn 1 sin û n For n = û 1 û û 1 = 1 4 d fl û l Γ δ n sin < 1 1 < 1 4 d fl l Γ δ n sin < 1 = 17

18 as < 4 d fl l Γ δ n sin Tis will certainly be acieved if Terefore < = < d l Γ δ n sin < d l Γ δ n < 1 d fl l Γ δ n sin d l Γ δ n. < d l Γ δ n < 1 = l fl < 1 < Γ. 5.8 dδn We will now prove tat n u n < u, we proved already tat u 1 < u, let us assume tat u j < u, j n and prove tat u n+1 < u û n+1 = 1 4 d l Γ δ n sin fl 4 d fl û n + l Γ δ,1 n sin û n 1 Let A 1 = 1 4 d l Γ δ n sin fl 4 d fl, A = l Γ δ,1 n sin Because of te condition 5.8 we ave bot A 1, A >. û n+1 = 1 4 d l Γ δ n sin fl 4 d fl û n + l Γ δ,1 n sin û n 1 By induction ypotesis we ave = A 1 û n + A û n 1 û n+1 < A 1 û n + A û n 1 û n+1 < A 1 û + A û Tis means û n+1 < A 1 + A û 18

19 But A 1 + A = 1 4 d l Γ δ n sin fl 4 d + l Γ δ,1 n fl sin = 1 Terefore û n+1 < û Tis prove tat n u n < u te numerical sceme solution to te fractional Order P.D.E is stable. 6. Grapical Simulations 19

20 Figure 1: Exact Solution of equation 5.1 wit ux,=expx; u,t=expct.

21 Figure : Approximate Solution of equation 5.1 given by 5.3, wit ux,=expx; u,t=expct. 1

22 Figure 3: Exact Solution of equation 5.1 wit ux,=cosx; u,t=cosct.

23 Figure 4: Approximate Solution of equation 5.1 given by 5.3, wit ux,=cosx; u,t=cosct. 3

24 7. Conclusion In order to extend te well-known Adams-Basfort numerical sceme to partial differential equations wit integer and non-integer order derivatives, we introduced a new reliable and efficient numerical sceme. Te metod is a combination of Laplace transform, Adams-Basfort and forward or backward numerical sceme. We developed te metod for general partial differential equations wit local and non-local differentiation. We presented in detail te error analysis and te convergence of te metod. In te case of fractional partial differential equations, te metod provides a numerical algoritm tat is easier to implement. Unlike te conventional metods, forward, backward, Crank-Nicolson, te cumbersome summation tat always appears in te additional term of teir numerical algoritm for te case of fractional partial differential equations, does not exist wit our metod. Tis leads to an easier proof of stability and convergence. We illustrated te metod by solving two partial differential equations including wave equation for te local case and a diffusion equation for fractional case. We studied te stability of eac example. Te proof sows witout doubt tat our metod is very stable and also converges very quickly to te exact solution. We believe tis metod will turn out to be a very useful numerical sceme tat will elp solving nonlinear and linear partial differential equations wit local and non-local operators. 4

25 Conflict of Interests Disclosure Te autors declare no conflict of Interests. References References [1] J.C. Butcer, Numerical metods for ordinary differential equations in te t century, Journal of Computational and Applied Matematics, Volume 15, Issue 1,, Pages 1-9, ISSN , ttp://dx.doi.org/1.116/s [] Armenio, V. 1997, AN IMPROVED MAC METHOD SIMAC FOR UNSTEADY HIGH-REYNOLDS FREE SURFACE FLOWS. Int. J. Numer. Met. Fluids, 4: doi:1.1/sici : 185::AID-FLD CO-Q. [3] Dietelm K, Freed AD 1999b On te solution of nonlinear fractional differential equations used in te modeling of viscoplasticity. In: Scientific computing in cemical engineering II: computational fluid dynamics, reaction engineering, and molecular properties, pp 17-4 [4] Donald N. Slinn, J.J. Riley, A Model for te Simulation of Turbulent Boundary Layers in an Incompressible Stratified Flow, Journal of Computational Pysics, Volume 144, Issue, 1998, Pages 55-6, ISSN , ttp://dx.doi.org/1.16/jcp [5] J. Peinado, J. Ibez, E. Arias, V. Hernndez, Adams-Basfort and Adams-Moulton metods for solving differential Riccati equations, Computers & Matematics wit Applications, Volume 6, Issue 11, 1, Pages , ISSN , ttp://dx.doi.org/1.116/j.camwa [6] Durran, D. R., 1991: Te tird-order Adams?Batfort metod: An attractive alternative to leapfrog time differencing. Mon. Wea. Rev., 119, 7?7, doi:ttps://doi.org/1.1175/ :TTOABM..CO. 5

26 [7] Y. He, W. Sun Stability and convergence of te Crank-Nicolson/Adams- Basfort sceme for te time-dependent Navier?Stokes equations SIAM J. Numer. Anal., 45 7, pp [8] G.C. Wu, D. Baleanu and Z.G. Deng, Variational iteration metod as a kernel constructive tecnique, Applied Matematical Modelling, 3915, , 15. [9] Atangana, Abdon, and Juan Jose Nieto. Numerical solution for te model of RLC circuit via te fractional derivative witout singular kernel. Advances in Mecanical Engineering : [1] Atangana, Abdon, and Dumitru Baleanu. Caputo-Fabrizio derivative applied to groundwater flow witin confined aquifer. Journal of Engineering Mecanics : D4165. [11] Gmez-Aguilar, J. F., et al. Atangana-Baleanu fractional derivative applied to electromagnetic waves in dielectric media. Journal of Electromagnetic Waves and Applications : [1] Yang, Xiao-Jun, and Dumitru Baleanu. Fractal eat conduction problem solved by local fractional variation iteration metod. Termal Science : [13] Yang, Xiao Jun, Dumitru Baleanu, and Hari Moan Srivastava. Local fractional integral transforms and teir applications. Academic Press, 15. [14] Zang, Yu, Hari M. Srivastava, and Miaela-Cristina Baleanu. Local fractional variational iteration algoritm II for non-omogeneous model associated wit te non-differentiable eat flow. Advances in Mecanical Engineering : [15] Gnitcogna, R.B., Atangana, A. Comparison of two iteration metods for solving nonlinear fractional partial differential equations 15 International Journal of Matematical Models and Metods in Applied Sciences, 9, pp