Analytic solution of homogeneous time-invariant fractional IVP

Transcription

1 Jaradat et al. Advances in Difference Equations ( :13 R E S E A R C H Open Access Analytic solution of homogeneous time-invariant fractional IVP Imad Jaradat *, Marwan Alquran and Mohammad Al-Dolat * Correspondence: iajaradat@just.edu.jo Department of Mathematics & Statistics, Jordan University of Science and Technology, Irbid, Jordan Abstract In this work, we propose a novel analytical solution approach for solving a general homogeneous time-invariant fractional initial value problem in the normal form [ ] u(x, t = F(u(x, t, 0 t < R, D t u(x,0=f(x, where D t is the Caputo fractional operator with 0 < 1. The solution is given analytically in the form of a convergent multi-fractional power series without using any particular treatments for the nonlinear terms. The new approach is taken to search patterns for compacton solutions of several nonlinear time-fractional dispersive equations, namely K (,, ZK (,, DD (1,,,andK (,,1.Remarkably, the graphical analysis showed that the n-term approximate memory solutions, labeled by the memory parameter 0 < 1, are continuously homotopic as they reflect, in some sense, some memory and heredity properties. MSC: 6A33; 35R11; 35F5; 35C10; 0C15 Keywords: Caputo derivative; Fractional power series; Fractional PDE; Power series method 1 Introduction The importance of fractional differential equations (FDEs is stimulated by the appearance of many scientific models that have a nonlocal dynamical property. It has been observed that such models possess a continuum flow due to the fractional derivative effect that can be characterized by a long-term memory in time (memory index. This memory index can be interpreted as a chaotic or bifurcation behavior of a certain phenomenon for a short time subject to past circumstances. For example, it has been observed that the universal electromagnetic, acoustic, and mechanical responses are influenced by a remnant memory which can be accurately modeled by the fractional diffusion wave equations [1]. Also, it has been shown that the propagation of mechanical diffusive waves in viscoelastic media can be identified by the fractional wave equation []. More applications that relate to Newtonian mechanics, electromagnetism, quantum mechanics, electrochemistry, signal and image processing, and biomedical engineering can be found in [3 10]. In general, it is not an easy task to extract an analytical solution for nonlinear fractional differential equations. Almost all attempts were developed by either finding nu- The Author(s 018. This article is distributed under the terms of the Creative Commons Attribution.0 International License ( which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s and the source, provide a link to the Creative Commons license, and indicate if changes were made.

2 Jaradat et al. Advances in Difference Equations ( :13 Page of 1 merical solutions over a specific range or considering few terms of an iterative computational series solution as an approximate. Such available methods are (He s variational iteration methods [11, 1], the iterative Laplace transform method [13], Adomian s decomposition method [1], the differential transform method [15, 16], homotopy analysis/perturbation methods [17, 18], the homotopy analysis Laplace transform method [19], Chebyshev/Jacobi/Legendre operational matrix methods [0], the fractional Lie groupmethod [1], the generalized Taylor power series method [, 3], and the residual power series method [ 8]. The motivation of the current work is to explore a closed-form solution of a general homogeneous time-invariant fractional initial value problem in the normal form, D [ ] ( t u(x, t = F u(x, t, 0 t < R, u(x,0=f (x, where D t is the Caputo fractional operator with (0, 1], u(x, t is an unknown function, F is an analytic differential operator in the variables x = x 1,...,x m that involves both linear and nonlinear terms, R R,andf (x C (R m. Classes of these equations include, but are not limited to, the Schrödinger equation, Korteweg de Vries equation, Burgers Fisher equation, Cauchy reaction diffusion equation, Boussinesq equation, and Sharma Tasso Olver equation. In particular, we have successfully applied the present approach to finding closed-form solutions for various nonlinear time-fractional dispersive equations, namely K (,, ZK (,, DD (1,,, and K (,,1 equations. It should be pointed out here that the existence of mild solutions of the nonlocal problem of time-fractional evolution equations is extensively studied in [9 38]. The organization of the current paper is as follows: In Sect., we recall some necessary definitions and theorems regarding the fractional derivative and fractional power series. In Sect. 3, the solution to a general time-invariant fractional IVP is constructed with some related convergence results and error analysis. In Sect., our approach is applied to various nonlinear time-fractional dispersive equations with graphical and numerical analysis to illustrate the adequacy of the proposed approach. Finally, some concluding remarks are given in Sect. 5. Preliminaries Many definitions and studies of fractional derivatives have been proposed in the literature. Probably this is due to the fact that no harmonious definition preserves all properties of the classical integer-order derivative. These definitions include the Grunwald Letnikov, Riemann Liouville, Weyl, Riesz and Caputo versions. However, in the Caputo case, the derivative of a constant function is zero and one can properly define the initial conditions for the fractional differential equations which can be handled by using an analogy with the classical integer-order case. For these reasons, we adopt the Caputo fractional derivative definition in this work. Definition.1 Arealfunctionu(t, t >0isinthespaceCλ R 1 if there exists a real number a > λ such that u(t=t a v(twherev(t C[0,, and it is in the space Cλ n if u(n Cλ 1, n N.

3 Jaradat et al. Advances in Difference Equations ( :13 Page 3 of 1 Definition. The Riemann Liouville fractional integral operator of order >0associates with a real function u(t Cλ 1 1 is defined as Jt [u(t] = 1 t Ɣ( 0 (t τ 1 u(τ dτ, and Jt 0 is an identity operator. Definition.3 The Caputo time-fractional derivative of order >0ofu(t C 1 n, n N is defined as D n t [u(t] = Jt [u (n (t] if n 1< < n and D t [u(t] = u(n (tif = n. Similarly, for n being the smallest integer that exceeds, the Caputo time-fractional derivative operator of order is given as D n t [u(x, t] = Jt [ n u(x,t t n ]ifn 1< < n and D t [u(x, t] = n u(x,t t n if = n. Remark 1 A direct implementation of the Caputo derivative yields D t [tp ]= Ɣ(p+1 Ɣ(p +1 tp for p >0andD t [c]=0wherecisaconstant. Also, it is easy to see that the Caputo derivative is a left inverse of the Riemann Liouville integral but not a right inverse. Specifically, for n 1< n, n N,andu(t Cλ 1 n we have D t J t [u(t] = u(t andjt D t [u(t] = u(t n 1 u(k (0 + tk k!,wheret >0. It should be noted here that it suffices to consider the Caputo fractional derivative of order 0 < 1sinceD t [u(t] = D (n 1 t [u (n 1 (t] for arbitrary order n 1< n,where (n 1 (0, 1]. Definition. ([5] A fractional power series (FPS expansion is an infinite series about t = t 0 of the form c k(t t 0 k where 0 n 1< n, t t 0. Theorem.5 [5] Suppose that u(t has a FPS expansion about t 0 as above for t 0 t t 0 + R. If the D k t [u(t] are continuous on (t 0, t 0 + R for k N, then c k = Dk t [u(t 0 ] where D k Ɣ(k+1 t is the k-fold Caputo derivative and R is the radius of convergence. Definition.6 ([] A power series of the form f k (xt k, (.1 where x I = I 1 I m R m,0< 1, and t 0 is called a multi-fractional power series about t =0. Theorem.7 ([] Suppose that u(x, t has a multi-fractional power series representation about t =0as above for x Iand0 t R. If D k t [u(x, t] are continuous on I (0, R for each k N, then f k (x= Dk t [u(x,0] where R is the radius of convergence. Ɣ(k+1 3 Analytic solution of homogeneous time-invariant fractional IVP As our approach depends mainly on constructing an analytical solution of the timefractional differential equation under consideration, we first present, in a similar fashion to the classical power series, some essential convergence theorems pertaining to our proposed solution. Theorem 3.1 Let {f k (x} be a sequence of functions f k : I R. If (.1 is convergent for some t = t 0 >0,then it is convergent for all t (0, t 0.

4 Jaradat et al. Advances in Difference Equations ( :13 Page of 1 Proof Assume that (.1 is convergent for t = t 0 >0.Then,forfixedɛ 0 >0,thereexistsN N such that f k (xt0 k < ɛ 0 for all k N. It follows that if k N,wehave f k (xt k < ɛ 0 ( t t 0 k for all x I and t (0, t 0, which shows that f k(xt k is absolutely convergent (and so convergent. We remark here that if f 0 (x is a bounded function on I, then the convergence of (.1at some t = t 0 > 0 implies the convergence on [0, t 0. Corollary 3. Let {f k (x} be a sequence of functions f k : I R. If (.1 is divergent for some t = t 0 >0,then it is divergent for all t > t 0. Proof Suppose not. That is, (.1 is convergent for some t > t 0. Then from Theorem 3.1,it converges on (0, tandthusconvergesatt 0, which is a contradiction. Corollary 3.3 Let {f k (x} be a sequence of functions f k : I R. Then one of the following casesistrueasregards(.1: p 1 : The series converges only at t =0; p : the series converges for all t 0; p 3 : there exists R >0(called the radius of convergence such that (.1 converges for all t (0, R and diverges for all t > R. Proof Suppose it is the case that both p 1 and p are not valid. Then there exist, β R + such that (.1 convergesatt = and diverges at t = β. Therefore, the set T = {t > 0: f k(xt k converges} is nonempty and T (0, bytheorem3.1. ThusR := supt exists.now,ift > R, then(.1 is divergent and if 0 < t < R, then, by the definition of the supremum there exists t < t 0 T such that (.1isconvergentatt 0 and so by Theorem 3.1 convergent on (0, t 0. The other cases can be handled easily. Now, consider the following general homogeneous time-invariant fractional initial value problem: D [ ] ( t u(x, t = F u(x, t, 0 t < R, u(x,0=f (x, (3.1 where D t is the Caputo fractional operator with (0, 1], u(x, t is an unknown function, F is an analytic differential operator in the variables x = x 1,...,x m that involves both linear and nonlinear terms, R R, andf (x C (R m. In our next theorem, we exhibit a parallel scheme of the Taylor series method to solve problem (3.1. The method gives an analytical solution in the form of convergent multi-fractional power series without the need for linearization, perturbation, or discretization of the variables. Instead of equating terms with the same degree of homogeneity, our approach depends recursively on timefractional differentiation to obtain the unknown series coefficients. Notation 3. We denote the coefficient extraction operator for a multi-fractional power series G(x, t, which extracts a constant multiple of the coefficient of t n in G, by[t n ] G.

5 Jaradat et al. Advances in Difference Equations ( :13 Page 5 of 1 More precisely, for n 1 [ t n ] G = [ t n] g k (xt k = Ɣ(n +1g n (x. (3. Note that, for a multi-fractional power series representation, G(x, t = g k(xt k,we have D n [ ] t G(x, t t=0 = Ɣ(n +1g n (x= [ t n] G, (3.3 where D n t = D t D t (n times. Theorem 3.5 Suppose that the solution u(x, t of (3.1 has a convergent multi-fractional power series representation of the form (.1 for x I R m, and 0 t R. If D n t [u(x, t] are continuous on I (0, R for n N, then the solution of (3.1 is given analytically by u(x, t=f (x+ n=1 [t (n 1 ] F Ɣ(n +1 tn. (3. Proof As u(x, t satisfies the initial condition, we should have f 0 (x =f (x. Applying the operator D t to Eq. (3.1(n 1 times, using the linearity property of the Caputo operator, and Remark 1 we have for n 1 D (n 1 t [ ( ] F u(x, t = D (n 1[ [ ]] t D t u(x, t [ ] f k (xt k = D n t = k=n Ɣ(k +1f k (x Ɣ((k n +1 t(k n, (3.5 for all x I and 0 t R.Inparticularatt =0,weobtain f n (xɣ(n +1=D (n 1 [ ( ] t F u(x, t t=0 (3.6 for x I, and hence f n (x= D(n 1 t [F(u(x, t] t=0 Ɣ(n +1 as required. = [t(n 1 ] F Ɣ(n +1 (3.7 Asan immediate consequenceof Theorems.7 and 3.5, we obtain the following generalized Taylor formula.

6 Jaradat et al. Advances in Difference Equations ( :13 Page 6 of 1 Corollary 3.6 Let u(x, t be the solution of (3.1 under the same hypotheses of Theorem 3.5. Then u(x, t is analytic at t =0in the sense of fractional power series. I.e. u(x, t= n=0 D n t [u(x,0] Ɣ(n +1 tn. (3.8 As an immediate specialcase of Theorem 3.5, we have the following explicit description of the solution coefficients in terms of the previous coefficient. Corollary 3.7 If F is a linear differential operator with constant coefficients, i.e. F(u = m k j=1 i=0 a ij i u where a x i ij s are constants, then for n 1 j f n (x= 1 Ɣ(n +1 m j=1 k i=0 i a ij x i f n 1 (x. (3.9 j Remark In practical terms to find f n (x forn 1, it is sufficient to substitute the (n 1-truncated series n 1 f k(xt k of u(x, t ineq.(3.7 since the remaining terms contain higher powers O(t n andd (n 1 t [O(t n ] = O(t=0whent =0. Due to the complexity related to this kind of equations, it is not always possible to find a general expression for the series coefficients. In such a case, the solution can be found approximately as a partial sum of the series, u n (x, t = n f k(xt k in some reasonable interval of t, and thus the overall errors can be made smaller by adding more new terms as shown in the following case. Corollary 3.8 Let {f k (x} be a uniformly bounded sequence of functions f k : I R. Then (.1 is uniformly convergent for 0 t R <1.Moreover, u(x, t u n (x, t 0 as n. Proof By the uniform boundedness of {f k (x},thereexistsm >0suchthat f k (x M for all k N and x I.Thus f k (xt k MR k.since MRk is convergent geometric series (with ratio R < 1, by the comparison test, we see that f k(xt k is uniformly convergent as desired. Moreover, u(x, t un (x, t k=n+1 MR k = MR(n+1 1 R n 0. (3.10 Applications and discussions The goal of this section is to verify the applicability of our proposed approach, derived from Theorem 3.5, in studying the memory effects due to the time-fractional derivative on various nonlinear dispersive equations. It should be noted here that all necessary calculations and graphical analysis have been done by using Mathematica 10 software. Example 1In[39], Rosenau and Hyman introduced a class of solitary waves with compact support (called compactons to understand the role of nonlinear dispersion on pattern formation in liquid drops. These compactons are solitary waves free of exponential tails and they re-emerge with exactly the same coherent shape after a collision. Motivated by the

7 Jaradat et al. Advances in Difference Equations ( :13 Page 7 of 1 importance of taking into account the memory effects due to time-fractional derivative, we first consider the time-fractional version K (, of the third-order Rosenau Hyman equation, which reads D t u + a( u x + ( u = 0, (.1 xxx subject to the initial condition with compact support c 3a u(x,0= cos ( ax, x π a, 0, otherwise, (. where t 0, a R +,and0< 1. In this case we have F(u = (a(u x +(u xxx. In accordance with Theorem 3.5, the proposed fractional power series solution of Eqs. (.1 and (.hastheform u(x, t= c 3a cos ( ax + f 1 (xt + f (xt + f 3 (xt 3 +. (.3 Following Eq. (3.7 and taking into account Remark, we successively obtain the coefficients f n (x as follows: c ( ax f 1 (x=+ 3 aɣ( +1 sin c 3 ( ax, f (x= 6Ɣ( +1 cos, c ( a ax f 3 (x= 1Ɣ(3 +1 sin c 5 ( a ax, f (x=+ Ɣ( +1 cos, c 6 ( a f 5 (x=+ 3 ax 8Ɣ(5 +1 sin c 7 a ( ax, f 6 (x= 96Ɣ(6 +1 cos, c 8 ( a f 7 (x= 5 ax 19Ɣ(7 +1 sin c 9 a 3 ( ax, f 8 (x=+ 38Ɣ(8 +1 cos,.. and so on. In general, for n 1wehave c 3a Ɣ(n +1f n (x= ( 1k ( c a c 3a ( 1k ( c a k cos( ax k+1 sin( ax ifn =k, ifn =k +1. (. Therefore, in view of (3. we have the exact memory solution of K (, in closed form, [ u(x, t= c ( ax 1+cos 3a ( 1 k ( c at Ɣ(k +1 k ( ax ( 1 k ( c at ] k+1 + sin Ɣ((k = c [ 1+cos 3a ( ax ( E,1 c a t

8 Jaradat et al. Advances in Difference Equations ( :13 Page 8 of 1 + c at ( ax sin ( ] E,+1 c a t, (.5 where E,1 ( ande,+1 ( are the two-parameter Mittag Leffler functions. Similar versions of (.1 were handled using homotopy analysis/perturbation methods [0, 1] and the reduced differential transform approach [0], where the obtained solutions were consistent with (.5. Particularly when = 1, we have the exact solution for the classical K(, Rosenau Hyman equation [39] [ u(x, t= c ( ax ( 1 k ( c at k 1+cos 3a (k! ( ax ( 1 k ( ] c at k+1 + sin (k +1! = c [ ( ( ( ( ] ax c at ax c at 1+cos cos + sin sin 3a = c ( a(x ct 3a cos (.6 with compact support x ct π a. Figure 1 shows the effect of different values of the fractional derivative order 0 < 1on the cross section approximate solutions u (x, t= n=0 f n(xt n of the equation for K (, with fixed t =0.5andc = a = 1 in the compact support x 0.5 π. Clearly, u (x,0.5 for = 1 is in high agreement with the cross section of the exact solution for K(, on its compact support. Furthermore, u (x, 0.5 for different values of 0 < 1arecontinuously homotopic, as increases, to the corresponding one of the exact solution of K(,. Table 1 shows the consecutive errors con n (x, t= u n+1 (x, t u n (x, t for n =andn =10. It is clear that we have a remarkable accuracy even with few terms and we can have more Figure 1 The cross section t =0.5ofK (, approximate solutions for different values of with c = a =1 Table 1 The consecutive errors of K (, for different values 0 < < 1 at some points (x, t con (x, t con 10 (x, t =0.5 =0.50 =0.75 =0.5 =0.50 =0.75 ( π, E E E E E E 15 ( π, E E E E E E 11 ( π, E E E E E E 10 (π, E E E E E E 15 (π, E E E E E E 11 (π, E E E E E E 10

9 Jaradat et al. Advances in Difference Equations ( :13 Page 9 of 1 accuracy when n gets larger and/or t gets closer to zero. Also, we can see that the error approaches zero as the memory index approaches one. Example In[, 3], Zakharov and Kuznetsov introduced an isotropic two-dimensional equation ZK(, to describe the weakly nonlinear ion acoustic waves in a strongly magnetized lossless plasma. As our next example, we consider the time-fractional nonlinear dispersive equation ZK (,, D t u +( u x + ( u xxx + ( u = 0, (.7 yyx subject to the initial condition with compact support c 3 u(x, y,0= cos ( x+y, x + y π, 0, otherwise, (.8 where t 0and0< 1. From Theorem 3.5, the fractional power series solution for Eqs. (.7and(.8hastheform u(x, y, t= c 3 cos ( x + y + f 1 (x, yt + f (x, yt + f 3 (x, yt 3 +. (.9 Proceeding as before with straightforward calculations, the coefficients are obtained recursively as c ( x + y f 1 (x, y=+ 6Ɣ( +1 sin c 3 ( x + y, f (x, y= 1Ɣ( +1 cos c ( x + y f 3 (x, y= Ɣ(3 +1 sin c 5 ( x + y, f (x, y=+ 8Ɣ( +1 cos c 6 ( x + y f 5 (x, y=+ 96Ɣ(5 +1 sin c 7 ( x + y, f 6 (x, y= 19Ɣ(6 +1 cos c 8 ( x + y f 7 (x, y= 38Ɣ(7 +1 sin c 9 ( x + y, f 8 (x, y=+ 768Ɣ(8 +1 cos..,,,, and so on. For n 1, we have the following general form for the unknown coefficients: c 3 Ɣ(n +1f n (x, y= ( 1k ( c k cos( x+y ifn =k, c 3 ( 1k ( c k+1 sin( x+y ifn =k +1. (.10 Therefore, the memory solution in closed form is given by [ u(x, y, t= c ( x + y ( 1 k ( ct k 1+cos 3 Ɣ(k +1 ( x + y ( 1 k ( ct ] k+1 + sin Ɣ((k +1 +1

10 Jaradat et al. Advances in Difference Equations ( :13 Page 10 of 1 = c [ ( x + y 1+cos E,1 ( c 3 t ( + ct x + y E sin,+1 ( c ] t. (.11 In particular when = 1, we have the exact solution to the classical ZK(, equation [] [ u(x, y, t= c ( x + y 1+cos 3 ( 1 k (k! ( x + y ( 1 k + sin (k +1! = c ( x + y ct 3 cos ( ct k ( ] ct k+1. (.1 Figure shows the effect of the memory index 0 < 1 on the surface approximate solutions u (x, y, t ofzk (, with fixed t =0.5andc =1.Insomesense,u (x, y,0.5 for different values of are continuously communicated to reach the corresponding one when =1. Table shows the consecutive errors for n =andn = 10. Apparently, we have a smaller error when t is relatively small, n is large, or approaches one. Figure The approximate solutions u (x, y,0.5ofzk (, for different values of with c =1 Table The consecutive errors of K (, for different values 0 < < 1 at some points (x, y, t con (x, t con 10 (x, t =0.5 =0.50 =0.75 =0.5 =0.50 =0.75 ( π, π, E E E 06.01E E E 15 ( π, π, E E E E E E 1 ( π, π, E E E E E E 10 (π, π, E E E E E E 15 (π, π, E E E E E E 1 (π, π, E E E E E 07.6E 10

11 Jaradat et al. Advances in Difference Equations ( :13 Page 11 of 1 Example 3In [5], Rosenau proposed a dispersive dissipative entity DD(k, m, n which combines the interaction between convection, dispersion and dissipation. In this example, we consider the time-fractional version DD (1,,(theclassicaloneisknownasthe K(, -Burger equation, which reads D t u ( u x + ( u xxx ( u 1 = 0, (.13 xx subject to the initial condition u(x,0= (1+tanh( x 1 = (1+e x, where t 0and0< 1. Following Eq. (3.7, we recursively obtain the coefficients f n (xforn 1asfollows: 1 Ɣ(n +1f n (x= n e x : n even, 1 n e x : n odd. (.1 Therefore, in view of (3. we have the memory solution of DD (1,,inclosedform, u(x, t= 1 e x = 1+e x 1 Ɣ(k +1 ( t k + e x ( 1 t Ɣ((k k+1 [ t ( t ( t ] E,+1 E,1. (.15 Particularly with = 1, we have the exact solution u(x, t= (1+e x t = (1+tanh( t x 1 for the classical DD(1,,equation[6]. Figure 3 shows the behavior of the approximate solutions u (x,0.5 of DD (1,, for various values of 0 < 1. Evidently, u (x,0.5 when = 1 is in harmony with the exact solution in x [ π+1, π+1 ]. Moreover, u (x, 0.5 for different values of is continuously communicated until = 1 is reached. Thus a convenient solution is expected for various values of. Example In[7], Dey studied the role of the fifth-order dispersion term in the existence of the compacton solutions and as regards the soliton stability for the usual Korteweg de Vries (KdV equation. We next consider the time-fractional version K (,,1oftheform D t u + ( u x ( u xxx + ( u 1 5x = 0 (.16 Figure 3 The cross section t =0.5ofDD (1,, approximate solutions for different values of

12 Jaradat et al. Advances in Difference Equations ( :13 Page 1 of 1 subject to the initial condition u(x,0= 16c 1 cosh ( x, where t 0and0< 1. Following the formula (3.7, we recursively obtain the coefficients f n (x, n 1 as 1 follows: 16c 1 Ɣ(n +1f n (x= ( c n cosh( x :n even, 16c 1 ( c n sinh( x :n odd. (.17 Therefore in view of (3., the solution of K (,,1isgiveninclosedformby u(x, t= = [ ( 16c 1 x 1+cosh ( x ( 1 ct sinh Ɣ((k [ ( 16c 1 x 1+cosh ( ct x sinh ( 1 ct Ɣ(k +1 ] k+1 ( c E t,1 k ( c E t ],+1, (.18 which is identical to the solution obtained by using the homotopy perturbation method [8]. Particularly with = 1, we have the exact solution for the classical K(,,1equation[9], u(x, t= [ ( 16c 1 x 1+cosh ( x ( 1 ct sinh (k +1! ( 16c 1 x ct = cosh 1 ( 1 ct k (k! k+1 ] (.19 with compact support x ct π. Figure exhibitstheapproximatesolutionsu (x,0.5ofk (,, 1 for various values of 0< 1. Apparently, u (x,0.5 when = 1 is in harmony with the exact solution on its compact support. Further, u (x, 0.5 for different values of is continuously communicated until =1isreached. Figure The cross section t =0.5ofK (,,1 approximate solutions for different values of

13 Jaradat et al. Advances in Difference Equations ( :13 Page 13 of 1 5 Conclusion In this paper, a hybrid version of the power series method is presented to handle a general time-invariant fractional initial value problem of the form (3.1. The solution of (3.1 is given analytically in the form of a convergent multi-fractional power series (.1 without the need of any linearization, perturbation, or discretization methods. Several nonlinear dispersive examples were tested and the exact memory solutions were successfully obtained in closed form. In fact, the physical interpretation of these solutions is beyond the scope of this work. However, in some sense, the graphs of the n-term approximate memory solutions, labeled by the memory index 0 < 1, were continuously deformed as they represent the memory and hereditary effects. Acknowledgements The authors would like to thank the editor and the referees for their valuable time and consideration. Competing interests The authors declare that they have no competing interests. Authors contributions All authors contributed equally to the manuscript and typed, read and approved the final manuscript. Publisher s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Received: 6 January 018 Accepted: 1 April 018 References 1. Nigmatullin, R.R.: To the theoretical explanation of the Universal response. Phys. Status Solidi B 13, (198. Mainardi, F., Paradisi, P.: Fractional diffusive waves. J. Comput. Acoust. 9, ( Podlubny, I.: Fractional Differential Equations. Academic Press, San Diego (1999. Miller, K.S., Ross, B.: An Introduction to the Fractional Calculus and Fractional Differential Equations. Wiley, New York ( Samko, S.G., Kilbas, A.A., Marichev, O.I.: Fractional Integrals and Derivatives: Theory and Applications. Gordon and Breach, New York ( Oldham, K.B., Spanier, J.: The Fractional Calculus: Theory and Application of Differentiation and Integration to Arbitrary Order. Academic Press, New York ( Kilbas, A.A., Srivastava, H.M., Trujillo, J.J.: Theory and Applications of Fractional Differential Equations. Elsevier, Amsterdam ( Diethelm, K.: The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type. Springer, Berlin ( Jesus, I.S., Machado, J.A.T.: Fractional control of heat diffusion systems. Nonlinear Dyn. 5, 63 8 ( Chen, Y., Vinagre, B.M., Podlubny, I.: Continued fraction expansion to discretize fractional order derivatives an expository review. Nonlinear Dyn. 38, ( Shakeri, F., Dehghan, M.: Numerical solution of the Klein-Gordon equation via He s variational iteration method. Nonlinear Dyn. 51,89 97 ( Odibat, Z., Momani, S.: Application of variational iteration method to nonlinear differential equations of fractional order. Int. J. Nonlinear Sci. Numer. Simul. 7, 7 3 ( Jafari, H., Nazari, M., Baleanu, D., Khalique, C.M.: A new approach for solving a system of fractional partial differential equations. Comput. Math. Appl. 66, ( Ray, S.S., Bera, R.K.: Analytical solution of a fractional diffusion equation by Adomian decomposition method. Appl. Math. Comput. 17, ( Liu, J., Hou, G.: Numerical solutions of the space- and time-fractional coupled Burgers equations by generalized differential transform method. Appl. Math. Comput. 17, ( Odibat, Z., Kumar, S., Shawagfeh, N., Alsaedi, A., Hayat, T.: A study on the convergence conditions of generalized differential transform method. Math. Methods Appl. Sci. 0, 0 8 ( Dehghan, M., Manafian, J., Saadatmandi, A.: Solving nonlinear fractional partial differential equations using the homotopy analysis method. Numer. Methods Partial Differ. Equ. 6, 8 79 ( He, J.-H.: Homotopy perturbation method: a new nonlinear analytical technique. Appl. Math. Comput. 135, ( Kumar, A., Kumar, S.: A modified analytical approach for fractional discrete KdV equations arising in particle vibrations. Proc.Natl.Acad.Sci.IndiaSect.APhys.Sci.88, ( Bhrawy, A.H., Taha, T.M., Machado, J.A.T.: A review of operational matrices and spectral techniques for fractional calculus. Nonlinear Dyn. 81, ( Jafari, H., Kadkhoda, N., Baleanu, D.: Fractional lie group method of the time-fractional Boussinesq equation. Nonlinear Dyn. 81, (015

14 Jaradat et al. Advances in Difference Equations ( :13 Page 1 of 1. Alquran, M., Jaradat, I.: A novel scheme for solving Caputo time-fractional nonlinear equations: theory and application. Nonlinear Dyn. 91, ( Jaradat, I., Al-Dolat, M., Al-Zoubi, K., Alquran, M.: Theory and applications of a more general form for fractional power series expansion. Chaos Solitons Fractals 108, (018. El-Ajou, A., Abu-Arqub, O., Momani, S.: Approximate analytical solution of the nonlinear fractional KdV-Burgers equation: a new iterative algorithm. J. Comput. Phys. 93, ( El-Ajou, A., Abu-Arqub, O., Al-Zhour, Z., Momani, S.: New results on fractional power series: theories and applications. Entropy 15, ( Kumar, A., Kumar, S., Yan, S.-P.: Residual power series method for fractional diffusion equations. Fundam. Inform. 151, ( Zhang, Y., Kumar, A., Kumar, S., Baleanu, D., Yang, X.-J.: Residual power series method for time-fractional Schrödinger equations. J. Nonlinear Sci. Appl. 9, ( Alquran, M., Jaradat, H.M., Syam, M.I.: Analytical solution of the time-fractional Phi- equation by using modified residual power series method. Nonlinear Dyn. 90,55 59 ( Chen, P., Li, Y.: Existence of mild solutions for fractional evolution equations with mixed monotone nonlocal conditions. Z. Angew. Math. Phys. 65, ( Chen, P., Zhang, X., Li, Y.: Nonlocal problem for fractional stochastic evolution equations with solution operators. Fract. Calc.Appl.Anal.19, ( Chen, P., Zhang, X., Li, Y.: Approximation technique for fractional evolution equations with nonlocal integral conditions. Mediterr. J. Math. 1(1, 6 ( Chen, P., Zhang, X., Li, Y.: Study on fractional non-autonomous evolution equations with delay. Comput. Math. Appl. 73, ( Chen, P., Li, Y., Chen, Q., Feng, B.: On the initial value problem of fractional evolution equations with noncompact semigroup. Comput. Math. Appl. 67, (01 3. Wang, J., Zhou, Y.: A class of fractional evolution equations and optimal controls. Nonlinear Anal., Real World Appl. 1, 6 7 ( Ouyang, Z.: Existence and uniqueness of the solutions for a class of nonlinear fractional order partial differential equations with delay. Comput. Math. Appl. 61, ( Tatar, N.: Existence results for an evolution problem with fractional nonlocal conditions. Comput. Math. Appl. 60, ( Balachandran, K., Park, J.Y.: Nonlocal Cauchy problem for abstract fractional semilinear evolution equations. Nonlinear Anal. 71, ( Li, F., Liang, J., Xu, H.K.: Existence of mild solutions for fractional integrodifferential equations of Sobolev type with nonlocal conditions. J. Math. Anal. Appl. 391, ( Rosenau, P., Hyman, J.M.: Compactons: solitons with finite wavelengths. Phys. Rev. Lett. 70, ( Singh, J., Kumar, D., Swroop, R., Kumar, S.: An efficient computational approach for time-fractional Rosenau Hyman equation. Neural Comput. Appl. 017, 1 8 ( Odibat, Z.: Compact and noncompact structures for nonlinear fractional evolution equations. Phys. Lett. A 37, (008. Zakharov, V.E., Kuznetsov, E.A.: On three-dimensional solitons. Sov. Phys. JETP 39, ( Kivshar, Y.S., Pelinovsky, D.E.: Self-focusing and transverse instabilities of solitary waves. Phys. Rep. 331, (000. Wazwaz, A.M.: Nonlinear dispersive special type of the Zakharov Kuznetsov equation ZK(n, n with compact and noncompact structures. Appl. Math. Comput. 161, ( Rosenau, P.: On a class of nonlinear dispersive-dissipative interactions. Physica D 13,55 56 ( Wazwaz, A.M.: The tanh method for compact and noncompact solutions for variants of the KdV-Burger and the K(n, n-burger equations. Physica D 13, ( Dey, B.: Compacton solutions for a class of two parameter generalized odd-order Korteweg de Vries equations. Phys. Rev. E 57, ( Koçak, H., Özis, T., Yildirim, A.: Homotopy perturbation method for the nonlinear dispersive K(m, n, 1 equations with fractional time derivatives. Int. J. Numer. Methods Heat Fluid Flow 0, ( Zhu, Y., Tong, K., Chaolu, T.: New exact solitary-wave solutions for the K(,, 1 and K(3, 3, 1 equations. Chaos Solitons Fractals 33, (007