JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 192519 10.1155/2014/192519 192519 Research Article Approximate Solution of Fractional Nonlinear Partial Differential Equations by the Legendre Multiwavelet Galerkin Method Mohamed M. A. 1 Torky M. Sh. 2 Ezzat Magdy A. 1 Faculty of Science, Suez Canal University Ismailia, Ismailia Egypt scuegypt.edu.eg 2 The High Institute of Administration and Computer, Port Said University, Port Said Egypt psu.edu.eg 2014 12 2 2014 2014 25 05 2013 03 01 2014 04 01 2014 3 3 2014 2014 Copyright © 2014 M. A. Mohamed and M. Sh. Torky. This 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 Legendre multiwavelet Galerkin method is adopted to give the approximate solution for the nonlinear fractional partial differential equations (NFPDEs). The Legendre multiwavelet properties are presented. The main characteristic of this approach is using these properties together with the Galerkin method to reduce the NFPDEs to the solution of nonlinear system of algebraic equations. We presented the numerical results and a comparison with the exact solution in the cases when we have an exact solution to demonstrate the applicability and efficiency of the method. The fractional derivative is described in the Caputo sense.

1. Introduction

Nowadays, fractional differential equations have garnered a great deal of attention and appreciation recently due to its ability to provide an accurate description of different nonlinear phenomena. The process of development of models based on fractional order differential systems has lately gained popularity in the investigation of dynamical systems. The advantage of fractional order systems is that they allow greater degrees of freedom in the model. The field of chaos has also snatched the attention of the researchers and this contributes to a large amount of the current research these days.

In recent decades, fractional calculus has found diverse applications in different scientific and technological fields , such as thermal engineering, acoustics, electromagnetism, control, robotics, viscoelasticity, diffusion, edge detection, turbulence, signal processing, information sciences, communications, and many other physical processes and also in medical sciences. Fractional differential equations (FDEs) have also been applied in modeling many physical and engineering problems and fractional differential equations in nonlinear dynamic [6, 7]. The importance of getting approximate and exact solutions of nonlinear fractional differential equations in mathematics and physics remains an important problem that requires to discover new methods of approximate and exact solutions. However, finding exact solutions to these nonlinear fractional differential equations are difficult to obtain it . Therefore, the numerical methods used to deal with these equations  and they have largely been using some semianalytical techniques to solve these equations such as, differential transform method , Adomian decomposition method , Laplace decomposition method , homotopy perturbation method , and variational iteration method . The majority of these methods have shortcomings inbuilt such as calculating Adomian’s polynomials, the Lagrange multiplier, mixed results, and the large computational work.

The aim of this paper is to expand the application of Legendre multiwavelet Galerkin method to provide approximate solutions for initial value problems of fractional nonlinear partial differential equations and to make comparison with that obtained by other numerical methods.

2. Preliminaries and Notations

In this section, we give the definition of the Riemann-Liouville fractional derivative and fractional integral with some basic properties.

Definition 1.

The left sided Riemann-Liouville fractional integral of order μ 0 ,  of a function f C α , α - 1 , is defined as (1) I μ f ( t ) = { 1 Γ ( μ ) 0 t f ( τ ) ( t - τ ) 1 - μ d τ , μ > 0 , t > 0 , f ( t )    μ = 0 .

Definition 2.

The (left sided) Caputo fractional derivative of f , f C - 1 m , m I N { 0 } , is defined as  (2) D * μ f ( t ) = μ f ( t ) t μ = { I m - μ [ m f ( t ) t m ] ,    μ > 0 , t > 0 , m f ( t ) t m μ = m . Note that 

(3) I t μ f ( x , t ) = 1 Γ ( μ ) 0 t f ( x , s ) ( t - s ) 1 - π d τ , μ > 0 , t > 0 ,

(4) D * μ f ( x , t ) = I t m - μ [ m f ( x , t ) t m ] , m - 1 < μ m ,

(5) I μ t γ = Γ ( γ + 1 ) Γ ( γ + μ + 1 ) t γ + μ ,    μ > 0 ,    γ > - 1 , t > 0 .

3. Properties of Legendre Multiwavelets 3.1. Wavelets [<xref ref-type="bibr" rid="B35">35</xref>, <xref ref-type="bibr" rid="B36">36</xref>]

Wavelets constitute a family of functions constructed from dilation and translation of a single function called the mother wavelet. When the dilation parameter a and the translation parameter b vary continuously we have the following family of continuous wavelets: (6) ψ a , b ( t ) = | a | - 1 / 2 ψ ( t - b a ) , a , b , a 0 . If we restrict the parameters a and b to discrete values as a = a 0 - k , b = n b 0 a 0 - k , a 0 > 1 , b 0 > 0 , and n , k we have the following family of discrete wavelets: (7) ψ k , n ( t ) = | a | k / 2 ψ ( a 0 k t - n b 0 ) , where ψ k , n ( t ) form a wavelet basis for L 2 ( R ) . In particular, when a 0 = 2 and b 0 = 1 then ψ k , n ( t ) forms an orthonormal basis .

3.2. Legendre Multiwavelets [<xref ref-type="bibr" rid="B36">36</xref>]

Legendre multiwavelets ψ n m ( t ) = ψ ( k , n , m , t ) have four arguments; n , n = 0 , 1 , 2 , , 2 k - 1 , k can assume any positive integer; m is the order for Legendre polynomials and t is the normalized time. They are defined on the interval [ 0 , 1 ] (8) ψ n m ( t ) = { 2 m + 1    2 k / 2 P m ( 2 k t - n ) , for    n 2 k t    n + 1 2 k 0 , otherwise , where m = 0 , 1 , , M - 1 , n = 0 , 1 , 2 , , 2 k - 1 . The coefficient 2 m + 1 is for orthonormality; P m ( t ) are the well-known shifted Legendre polynomials of order m which are defined on the interval [ 0 , 1 ] and can be determined with the aid of the following recurrence formula: (9) P 0 ( t ) = 1 , P 1 ( t ) = 2 t - 1 , P m + 1 ( t ) = ( 2 m + 1 m + 1 ) ( 2 t - 1 ) P m ( t ) - ( m m + 1 ) P m - 1 ( t ) , m = 1,2 , 3 , . Also the two-dimensional Legendre multiwavelet is defined as (10) ψ n 1 m 1 n 2 m 2 ( x , t ) = { A P m 1 ( 2 k 1 x - n 1 ) P m 2 ( 2 k 2 t - n 2 ) , for    n 1 2 k 1 x    n 1 + 1 2 k 1 n 2 2 k 2 t    n 2 + 1 2 k 2 0 , otherwise , where       A = ( 2 m 1 + 1 ) ( 2 m 2 + 1 )    2 ( k 1 + k 2 ) / 2 , n 1 and n 2 are defined similarly to n , k 1 and k 2 can assume any positive integer, m 1 and m 2 are the order for the Legendre polynomials, and ψ n 1 m 1 n 2 m 2 ( x , t ) forms a basis for L 2 ( [ 0 , 1 ] × [ 0 , 1 ] ) .

3.3. Function Approximation

A function f ( x , t ) defined over [ 0 , 1 ] × [ 0 , 1 ] can expand as (11) f ( x , t ) = n = 1 i = 0 l = 1 j = 0 c n , i , l , j ψ n , i ( x ) ψ l , j ( t ) . If the infinite series in (11) is truncated, it can be written as (12) f ( x , t ) = n = 1 2 k 1 - 1 i = 0 N l = 1 2 k 2 - 1 j = 0 M c n , i , l , j ψ n , i ( x ) ψ l , j ( t ) = Ψ T ( x ) F    Ψ ( t ) , where Ψ ( x ) and Ψ ( t ) are 2 k 1 ( M 1 + 1 ) × 1 and 2 k 2 ( M 2 + 1 ) × 1 matrices, respectively, given by (13) Ψ ( x ) = [ ψ 10 ( x ) , , ψ 1 M 1 ( x ) , , ψ 20 ( x ) , , ψ 2 M 1 ( x ) , , ψ ( 2 k 1 - 1 ) 0 ( x ) , , ψ ( 2 k 1 - 1 ) M 1 ( x ) ] , Ψ ( t ) = [ ψ 10 ( t ) , , ψ 1 M 1 ( t ) , , ψ 20 ( t ) , , ψ 2 M 1 ( t ) , , ψ ( 2 k 1 - 1 ) 0 ( t ) , , ψ ( 2 k 1 - 1 ) M 1 ( t ) ] . In addition, F is a 2 k 1 ( M 1 + 1 ) × 2 k 2 ( M 2 + 1 ) matrix whose elements can be calculated from (14) 0 1 ψ n i ( x ) ψ l j ( t ) f ( x , t ) d t d x , with n = 0 , 1 , , 2 k 1 - 1 , i = 0 , , M 1 , l = 0 , 1 , , 2 k 2 - 1 , j = 0 , , M 2 .

4. Solution of Nonlinear Fractional Partial Differential Equations

Consider the nonlinear fractional partial differential equation (15) D t α u = N ( u ) + g ( x , t ) ,          m < α < m + 1 , m    0 , with initial condition u ( x , 0 ) = f ( x ) .

Let (16) F ( u ) = D t α u - N ( u ) - g ( x , t ) .

A Galerkin approximation to (16) is constructed as follows. The approximation u N M is sought in the form of the truncated series (17) u N M ( x , t ) = { n = 1 2 k 1 i = 0 N l = 1 2 k 2 j = 0 M t c n , i , l , j ψ n , i ( x ) ψ l , j ( t ) + u ( x , 0 ) , for    m = 0 , n = 1 2 k 1 i = 0 N l = 1 2 k 2 j = 0 M t 2 c n , i , l , j ψ n , i ( x ) ψ l , j ( t ) + u ( x , 0 ) + t u t ( x , 0 ) , for    m = 1 , where ψ i j are the Legendre multiwavelet basis.

The expansion coefficients c n , i , l , j are determined by Galerkin equations: (18) F ( u N M ) , ψ n , i ψ l , j = 0 , where · denotes inner product defined as (19) F ( u N M ) , ψ n , i ψ l , j = 0 1 F ( u N M ) ( x , t ) ψ n i ( x ) ψ l j ( t ) d t d x . Galerkin equations (18) give a system of 2 k 1 - 1 ( N + 1 ) × 2 k 2 - 1 ( M + 1 ) equations that can be solved for the elements of (20) a n , i , l , j , i = 0,1 , , N , j = 0,1 , , M , n = 1,2 , , 2 k 1 , l = 1,2 , , 2 k 2 .

5. Illustrative Example

To demonstrate the effectiveness of the method, here we consider some linear fractional partial differential equations. The Legendre wavelets are defined only for t [ 0 , 1 ] ; we take a = 0 , b = 1 . The computations associated with the examples were performed using Mathematica and Maple.

Example 1.

Consider the nonlinear time-fractional diffusion equation in absence of both external force and reaction term  (21) α u t α - u x ( u u x ) = 0 , 0 < α 1 , t > 0 , with initial condition u ( x , 0 ) = x .

We applied the method presented in this paper for k 1 = k 2 = 0 and M = N = 1 ; from (17) we have (22) u N M ( x , t ) = i = 1 1 j = 0 1 t c 0 , i , 0 , j ψ 0 , i ( x ) ψ 0 , j ( t ) + x . Substituting (22) into (21) and using (18) and (19) we obtained the solution of (21) for different values of α = { 1 / 3 , 1 / 2 , 2 / 3,1 } . Table 1 shows the approximate solutions for (21) obtained for different values of α using the Legendre multiwavelet method and the Homotopy perturbation method . The values of α = 1 are the only case for which we know the exact solution u ( x , t ) = x + t and our approximate solution using Legendre multiwavelet method coincides with the approximate solution obtained using the Homotopy perturbation method . It is noted that only two bases of Legendre multiwavelet and fourth-order term of Homotopy perturbation method were used in evaluating the approximate solution of Table 1. Figure 1 shows a plot of u ( x , t ) with respect to x and t for different values of α = { 1 / 3 , 1 / 2 , 2 / 3,1 } . Figure 2 shows a plot of u ( x , t ) for various values of t and different values of α = { 1 / 3 , 1 / 2 , 2 / 3,1 } , x = 1 .

Approximate solutions for Example 1.

t x α = 1 / 3 α = 1 / 2 α = 2 / 3 α = 1
u HPM u LEG u HPM u LEG u HPM u LEG u HPM u LEG u EX
0.0 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25 0.25
0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50 0.50
0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75 0.75
1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00

0.25 0.25 0.95546 0.93263 0.81419 0.78225 0.68960 0.66479 0.50 0.50 0.50
0.50 1.20546 1.18263 1.06419 1.03225 0.93960 0.91479 0.75 0.75 0.75
0.75 1.45546 1.43263 1.31419 1.28225 1.18960 1.16479 1.00 1.00 1.00
1.00 1.70546 1.68263 1.56419 1.53225 1.43960 1.41479 1.25 1.25 1.25

0.50 0.25 1.13882 1.30877 1.04788 1.13276 0.94783 0.98312 0.75 0.75 0.75
0.50 1.38882 1.55877 1.29788 1.38276 1.19783 1.23312 1.00 1.00 1.00
0.75 1.63882 1.80877 1.54788 1.63276 1.44783 1.48312 1.25 1.25 1.25
1.00 1.88882 2.05877 1.79788 1.88276 1.69783 1.73312 1.50 1.50 1.50

0.75 0.25 1.26745 1.37843 1.22721 1.30152 1.16441 1.20499 1.00 1.00 1.00
0.50 1.51745 1.62843 1.47721 1.55152 1.41441 1.45499 1.25 1.25 1.25
0.75 1.76745 1.87843 1.72721 1.80152 1.66441 1.70499 1.50 1.50 1.50
1.00 2.01745 2.12843 1.97721 2.05152 1.91441 1.95499 1.75 1.75 1.75

1.00 0.25 1.36985 1.14160 1.37838 1.28854 1.35773 1.33039 1.25 1.25 1.25
0.50 1.61985 1.39160 1.62838 1.53854 1.60773 1.58039 1.50 1.50 1.50
0.75 1.86985 1.64160 1.87838 1.78854 1.85773 1.83039 1.75 1.75 1.75
1.00 2.11985 1.89160 2.12838 2.03854 2.10773 2.08039 2.00 2.00 2.00

(a) Plot of u ( x , t ) with respect to x and t at α = 1 / 3 . (b) Plot of u ( x , t ) with respect to x and t at α = 1 / 2 . (c) Plot of u ( x , t ) with respect to x and t at α = 2 / 3 . (d) Plot of u ( x , t ) with respect to x and t at α = 1 .

Plot of u ( x , t ) for various values of t and different values of α , x = 1 .

Example 2.

Consider the fractional nonlinear Klein-Gordon equation  (23) α u t α - x ( u x ) + u 2 = 0 , 1 α 2 , t > 0 , with initial condition u ( x , 0 ) = 1 + sin x ,    u t ( x , 0 ) = 0 .

We applied the method presented in this paper for k 1 = k 2 = 0 and M = N = 1 ; from (17) we have (24) u N M ( x , t ) = i = 1 1 j = 0 1 t 2 c 0 , i , 0 , j ψ 0 , i ( x ) ψ 0 , j ( t ) + 1 + sin x . Substituting (24) into (23) and using (17) we obtained the solution of (23) for different values of α = { 1 , 1.5 , 1.75 , 2 } . Table 2 shows the approximate solutions for (23) obtained for different values of α using the Legendre multiwavelet method. We have given the solution simulations in Figure 3 according to different values of α . Figure 4 shows a plot of u ( x , t ) for various values of t and different values of α , x = 1 .

Approximate solutions for Example 2.

t x α = 1 α = 1.5 1.75 α = 2
u LEG u LEG u LEG u LEG
0.25 0.25 0.918665978 1.132388298 1.164570723 1.188858256
0.50 1.023903100 1.316468982 1.361000156 1.394880289
0.75 1.099331862 1.470741307 1.527621230 1.571093962
1.00 1.132379628 1.582632636 1.651861308 1.704926639

0.50 0.25 0.670529477 0.862842686 0.954935119 1.028617931
0.50 0.683586761 0.936088923 1.063084555 1.165329668
0.75 0.666835687 0.979526801 1.141425631 1.272233047
1.00 0.607703616 0.980583683 1.177385711 1.336755430

0.75 0.25 0.502994454 0.552019181 0.676793302 0.7897781568
0.50 0.458476525 0.501019780 0.671719555 0.8269013752
0.75 0.384150235 0.420212023 0.636837450 0.8342162344
1.00 0.267442949 0.297023267 0.559574347 0.799150097

1.00 0.25 0.416060911 0.313169842 0.388441429 0.4954341090
0.50 0.348572389 0.173995974 0.272945982 0.415723105
0.75 0.251275507 0.005013748 0.127642174 0.306203742
1.00 0.111597629 −0.206349475 −0.060042629 0.154303383

(a) Plot of u ( x , t ) with respect to x and t at α = 1 . (b) Plot of u ( x , t ) with respect to x and t at α = 1.5 . (c) Plot of u ( x , t ) with respect to x and t at α = 1.75 . (d) Plot of u ( x , t ) with respect to x and t at α = 2 .

Plot of u ( x , t ) for various values of t and different values of α , x = 1 .

Example 3.

Consider the following nonlinear time-fractional equation [28, 39, 40]: (25) α u t α + u u x = x + x t 2 , 0 < α 1 , t > 0 , with initial condition u ( x , 0 ) = 0 .

We applied the method presented in this paper for k 1 = k 2 = 0 and M = N = 1 ; from (17) we have (26) u N M ( x , t ) = i = 1 1 j = 0 1 t 2 c 0 , i , 0 , j ψ 0 , i ( x ) ψ 0 , j ( t ) . Substituting (26) into (25) and using (17) we obtained the solution of (25) for different values of α = { 0.5 , 0.75 , 1 } . We have given the solution simulations in Figure 5 according to different values of α . Figure 6 shows a plot of u ( x , t ) for various values of t and different values of    α , x = 1 . Tables 3, 4, and 5 show the approximate solutions for (25) obtained using the Legendre multiwavelet method, the decomposition method , the variational iteration method , Homotopy perturbation method , and generalized differential transform method (GDTM)  for different values of α = { 0.5 , 0.75 , 1 } . The values of α = 1 are the only case for which we know the exact solution u ( x , t ) = x t and Table 4 provides that our approximate solution using Legendre multiwavelet is more accurate than the approximate solution obtained using the decomposition method, Homotopy perturbation method, and the variational iteration method. In addition to that, our approximate solution using Legendre multiwavelet as the approximate solution was obtained using GDTM and exact solution.

Numerical values when α = 0.5 .

t x u GDTM u HPM u VIM u ADM u Leg
0.2 0.25 0.123635 0.104573 0.103750 0.112844 0.092812
0.50 0.247270 0.209146 0.207499 0.225688 0.185624
0.75 0.370905 0.313720 0.311249 0.311249 0.278436
1.00 0.494540 0.418293 0.414999 0.451375 0.371248

0.4 0.25 0.177148 0.177229 0.172012 0.164004 0.161241
0.50 0.354295 0.354458 0.344025 0.328008 0.322483
0.75 0.531443 0.531686 0.516037 0.492011 0.483725
1.00 0.708590 0.708915 0.688050 0.656015 0.644967

0.6 0.25 0.226965 0.230500 0.215641 0.243862 0.205288
0.50 0.453931 0.461000 0.431283 0.487721 0.410577
0.75 0.680896 0.691499 0.646924 0.731581 0.615866
1.00 0.907861 0.921999 0.862566 0.975441 0.821155

Numerical values when α = 0.75 .

t x u GDTM u HPM u VIM u ADM u Leg
0.2 0.25 0.080835 0.078306 0.077933 0.078787 0.070083
0.50 0.161671 0.156612 0.155865 0.157574 0.140166
0.75 0.242596 0.234919 0.233798 0.236361 0.210249
1.00 0.323342 0.313225 0.311730 0.315148 0.280332

0.4 0.25 0.135446 0.136806 0.134855 0.128941 0.129648
0.50 0.270892 0.273612 0.269710 0.257881 0.259297
0.75 0.406338 0.410418 0.404565 0.386821 0.388945
1.00 0.541784 0.547225 0.539420 0.515762 0.518594

0.6 0.25 0.185529 0.185146 0.179990 0.177238 0.178696
0.50 0.371057 0.370292 0.359979 0.354477 0.357393
0.75 0.556586 0.555437 0.539969 0.531715 0.536089
1.00 0.742114 0.740583 0.719958 0.7089541 0.714786

Numerical values when α = 1 .

t x u GDTM u HPM u VIM u ADM u Leg u EX
0.2 0.25 0.05 0.049989 0.050309 0.050000 0.05 0.05
0.50 0.10 0.099978 0.100619 0.100000 0.10 0.10
0.75 0.15 0.149968 0.150928 0.150001 0.15 0.15
1.00 0.20 0.199957 0.201237 0.200001 0.20 0.20

0.4 0.25 0.10 0.099645 0.101894 0.100023 0.10 0.10
0.50 0.20 0.199290 0.203787 0.200046 0.20 0.20
0.75 0.30 0.298935 0.305681 0.300069 0.30 0.30
1.00 0.40 0.398580 0.407575 0.400092 0.40 0.40

0.6 0.25 0.15 0.147158 0.153094 0.150411 0.15 0.15
0.50 0.30 0.294317 0.306188 0.300823 0.30 0.30
0.75 0.45 0.441475 0.459282 0.451234 0.45 0.45
1.00 0.60 0.588634 0.612376 0.601646 0.60 0.60

Plot of u ( x , t ) with respect to x and t at (a) α = 0.5 , (b) α = 0.75 , and (c) α = 1 .

Plot of u ( x , t ) for various values of t and different values of α , x = 1 .

Example 4.

Consider the following time fractional advection nonhomogeneous equation [41, 42]: (27) α u t α + u u x = 2 t + x + t 3 + x t 2 , 0 < α 1 , t > 0 , with initial condition u ( x , 0 ) = 0 .

We applied the method presented in this paper for k 1 = k 2 = 0 and M = N = 1 ; from (17) we have (28) u N M ( x , t ) = i = 1 1 j = 0 1 t 2 c 0 , i , 0 , j ψ 0 , i ( x ) ψ 0 , j ( t ) .

Substituting (28) into (27) and using (17) we obtained the solution of (27) for different values of α = { 0.5 , 0.75 , 1 } . Table 6 shows the approximate solutions for (27) obtained for different values of α using the Legendre multiwavelet method and the Homotopy perturbation method . The values of α = 1 are the only case for which we know the exact solution u ( x , t ) = t 2 + x t and our approximate solution using Legendre multiwavelet method coincides with the approximate solution obtained using the Homotopy perturbation method [41, 42] and the approximate solution obtained using Adomian decomposition method and variational iteration method . We have given the solution simulations in Figure 7 according to different values of α = { 0.5 , 0.75 , 1 } . Figure 8 shows a plot of u ( x , t ) for various values of t and α = { 0.5,0.75,1 } , x = 1 .

Numerical values when α = { 0.5,0 . 75,1 } .

t x α = 0.5 α = 0.75 α = 1
u HPM u Leg u HPM u Leg u HPM u Leg u EX
0.2 0.25 0.243843 0.216416 0.153820 0.150051 0.09 0.09 0.09
0.50 0.351274 0.309228 0.232395 0.220134 0.14 0.14 0.14
0.75 0.458705 0.402040 0.310971 0.290218 0.19 0.19 0.19
1.00 0.566136 0.494853 0.389547 0.360301 0.24 0.24 0.24

0.4 0.25 0.534552 0.474657 0.388100 0.363967 0.26 0.26 0.26
0.50 0.667613 0.635899 0.514376 0.493616 0.36 0.36 0.36
0.75 0.800673 0.797140 0.640652 0.623264 0.46 0.46 0.46
1.00 0.933734 0.958382 0.766928 0.752913 0.56 0.56 0.56

0.6 0.25 0.934713 0.774722 0.726195 0.641747 0.51 0.51 0.51
0.50 1.083891 0.980011 0.890481 0.820443 0.66 0.66 0.66
0.75 1.233068 1.185300 1.054767 0.999140 0.81 0.81 0.81
1.00 1.382245 1.390589 1.219052 1.177837 0.96 0.96 0.96

Plot of u ( x , t ) with respect to x and t at (a) α = 0.5 , (b) α = 0.75 , and (c) α = 1 .

Plot of u ( x , t ) for various values of t and different values of α , x = 1 .

6. Conclusion

In this study, it is shown how Legendre multiwavelet can be applied to provide approximate solutions for initial value problems of fractional nonlinear partial differential equations. The Legendre multiwavelet properties are presented. The main characteristic of this approach is using these properties together with the Galerkin method to reduce the NFPDEs to the solution of nonlinear system of algebraic equations. In addition, we compered our results with that obtained by other numerical methods. The results show that the Legendre multiwavelet is a powerful mathematical tool for fractional nonlinear partial differential equations. We used Mathematica and Maple programs for computations in this paper.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Oldham K. B. Spanier J. The Fractional Calculus 1974 New York, NY, USA Academic Press MR0361633 Mainardi F. Fractional calculus: some basic problems in continuum and statistical mechanics Fractals and Fractional Calculus in Continuum Mechanics 1997 378 New York, NY, USA Springer 291 348 MR1611587 Podlubny I. Fractional Differential Equations 1999 198 San Diego, Calif, USA Academic Press MR1658022 Hilfer R. Applications of Fractional Calculus in Physics 1999 Orlando, Fla, USA Academic Press Dalir M. Bashour M. Applications of fractional calculus Applied Mathematical Sciences 2010 4 21–24 1021 1032 MR2607355 ZBL1195.26011 Jafari H. Daftardar-Gejji V. Solving a system of nonlinear fractional differential equations using Adomian decomposition Journal of Computational and Applied Mathematics 2006 196 2 644 651 10.1016/j.cam.2005.10.017 MR2249452 ZBL1099.65137 Lu J. Chen G. A note on the fractional-order Chen system Chaos, Solitons and Fractals 2006 27 3 685 688 Liu Y. Q. Ma J. H. Exact solutions of a generalized multi-fractional nonlinear diffusion equation in radical symmetry Communications in Theoretical Physics 2009 52 5 857 861 Anastassi Z. A. Simos T. E. Numerical multistep methods for the efficient solution of quantum mechanics and related problems Physics Reports 2009 482-483 1 240 10.1016/j.physrep.2009.07.005 MR2552375 Ertürk V. S. Momani S. Solving systems of fractional differential equations using differential transform method Journal of Computational and Applied Mathematics 2008 215 1 142 151 10.1016/j.cam.2007.03.029 MR2400624 ZBL1141.65088 Erturk V. S. Momani S. Odibat Z. Application of generalized differential transform method to multi-order fractional differential equations Communications in Nonlinear Science and Numerical Simulation 2008 13 8 1642 1654 10.1016/j.cnsns.2007.02.006 MR2398696 ZBL1221.34022 Odibat Z. Momani S. A generalized differential transform method for linear partial differential equations of fractional order Applied Mathematics Letters 2008 21 2 194 199 10.1016/j.aml.2007.02.022 MR2426978 ZBL1132.35302 Odibat Z. Momani S. Erturk V. S. Generalized differential transform method: application to differential equations of fractional order Applied Mathematics and Computation 2008 197 2 467 477 10.1016/j.amc.2007.07.068 MR2400669 ZBL1141.65092 Kurulay M. Bayram M. Approximate analytical solution for the fractional modified KdV by differential transform method Communications in Nonlinear Science and Numerical Simulation 2010 15 7 1777 1782 10.1016/j.cnsns.2009.07.014 MR2606417 ZBL1222.35172 Alomari A. K. A new analytic solution for fractional chaotic dynamical systems using the differential transform method Computers & Mathematics with Applications 2011 61 9 2528 2534 10.1016/j.camwa.2011.02.043 MR2795001 ZBL1221.65191 Liu J. Hou G. Numerical solutions of the space- and time-fractional coupled Burgers equations by generalized differential transform method Applied Mathematics and Computation 2011 217 16 7001 7008 10.1016/j.amc.2011.01.111 MR2775690 ZBL1213.65131 Çetinkaya A. Kiymaz O. The solution of the time-fractional diffusion equation by the generalized differential transform method Mathematical and Computer Modelling 2013 57 2349 2354 El-Sayed A. M. A. Gaber M. The Adomian decomposition method for solving partial differential equations of fractal order in finite domains Physics Letters A 2006 359 3 175 182 10.1016/j.physleta.2006.06.024 MR2253516 ZBL1236.35003 Jafari H. Daftardar-Gejji V. Solving linear and nonlinear fractional diffusion and wave equations by Adomian decomposition Applied Mathematics and Computation 2006 180 2 488 497 10.1016/j.amc.2005.12.031 MR2270026 ZBL1102.65135 Wang Q. Numerical solutions for fractional KdV-Burgers equation by Adomian decomposition method Applied Mathematics and Computation 2006 182 2 1048 1055 10.1016/j.amc.2006.05.004 MR2282548 ZBL1107.65124 Safari M. Ganji D. D. Moslemi M. Application of He's variational iteration method and Adomian's decomposition method to the fractional KdV-Burgers-Kuramoto equation Computers & Mathematics with Applications 2009 58 11-12 2091 2097 10.1016/j.camwa.2009.03.043 MR2557332 ZBL1189.65255 Yusufoğlu E. Numerical solution of Duffing equation by the Laplace decomposition algorithm Applied Mathematics and Computation 2006 177 2 572 580 10.1016/j.amc.2005.07.072 MR2291982 ZBL1096.65067 Khan Y. An effective modification of the Laplace decomposition method for nonlinear equations International Journal of Nonlinear Sciences and Numerical Simulation 2009 10 1373 1376 Jafari H. Khalique C. M. Nazari M. Application of the Laplace decomposition method for solving linear and nonlinear fractional diffusion-wave equations Applied Mathematics Letters 2011 24 11 1799 1805 10.1016/j.aml.2011.04.037 MR2812215 ZBL1231.65179 Gepreel K. A. Nofal T. A. Alotaibi F. M. Numerical solutions for the time and space fractional nonlinear partial differential equations Journal of Applied Mathematics 2013 2013 12 482419 10.1155/2013/482419 Gupta P. K. Singh M. Homotopy perturbation method for fractional Fornberg-Whitham equation Computers & Mathematics with Applications 2011 61 2 250 254 10.1016/j.camwa.2010.10.045 MR2754133 ZBL1211.65138 Wang Q. Homotopy perturbation method for fractional KdV equation Applied Mathematics and Computation 2007 190 2 1795 1802 10.1016/j.amc.2007.02.065 MR2339772 ZBL1122.65397 Momani S. Odibat Z. Homotopy perturbation method for nonlinear partial differential equations of fractional order Physics Letters A 2007 365 5-6 345 350 10.1016/j.physleta.2007.01.046 MR2308776 ZBL1203.65212 Yıldırım A. Koçak H. Homotopy perturbation method for solving the space-time fractional advection-dispersion equation Advances in Water Resources 2009 32 1711 1716 R Y. M. Noorani M. S. M. Hashim I. Variational iteration method for fractional heat- and wave-like equations Nonlinear Analysis: Real World Applications 2009 10 3 1854 1869 10.1016/j.nonrwa.2008.02.026 MR2502991 ZBL1172.35302 Odibat Z. Momani S. The variational iteration method: an efficient scheme for handling fractional partial differential equations in fluid mechanics Computers & Mathematics with Applications 2009 58 11-12 2199 2208 10.1016/j.camwa.2009.03.009 MR2557350 ZBL1189.65254 Guo S. Mei L. Li Y. Fractional variational homotopy perturbation iteration method and its application to a fractional diffusion equation Applied Mathematics and Computation 2013 219 11 5909 5917 10.1016/j.amc.2012.12.003 MR3018433 ZBL1273.65157 Luchko Y. Gorenflo R. An operational method for solving fractional differential equations with the Caputo derivatives Acta Mathematica Vietnamica 1999 24 2 207 233 MR1710779 ZBL0931.44003 Lukashchuk S. Y. An approximate solution method for ordinary fractional differential equations with the Riemann-Liouville fractional derivatives Communications in Nonlinear Science and Numerical Simulation 2014 19 2 390 400 MR3142482 10.1016/j.cnsns.2013.06.021 Gu J. S. Jiang W. S. The Haar wavelets operational matrix of integration International Journal of Systems Science 1996 27 623 628 Yousefi S. A. Legendre multiwavelet Galerkin method for solving the hyperbolic telegraph equation Numerical Methods for Partial Differential Equations 2010 26 3 535 543 10.1002/num.20445 MR2656299 ZBL1189.65231 Das S. Vishal K. Gupta P. K. Solution of the nonlinear fractional diffusion equation with absorbent term and external force Applied Mathematical Modelling 2011 35 8 3970 3979 10.1016/j.apm.2011.02.003 MR2793679 ZBL1221.35437 Kurulay M. Solving the fractional nonlinear Klein-Gordon equation by means of the homotopy analysis method Advances in Difference Equations 2012 2012, article 187 8 10.1186/1687-1847-2012-187 MR3016697 Odibat Z. Momani S. Numerical methods for nonlinear partial differential equations of fractional order Applied Mathematical Modelling 2008 32 28 39 Momani S. Odibat Z. A novel method for nonlinear fractional partial differential equations: combination of DTM and generalized Taylor's formula Journal of Computational and Applied Mathematics 2008 220 1-2 85 95 10.1016/j.cam.2007.07.033 MR2444156 ZBL1148.65099 Liu Y. Approximate solutions of fractional nonlinear equations using homotopy perturbation transformation method Abstract and Applied Analysis 2012 2012 14 752869 10.1155/2012/752869 MR2910714 ZBL1242.65136 Khan Y. Wu Q. Homotopy perturbation transform method for nonlinear equations using He's polynomials Computers & Mathematics with Applications 2011 61 8 1963 1967 10.1016/j.camwa.2010.08.022 MR2785523 ZBL1219.65119 Wazwaz A.-M. A comparison between the variational iteration method and Adomian decomposition method Journal of Computational and Applied Mathematics 2007 207 1 129 136 10.1016/j.cam.2006.07.018 MR2332954 ZBL1119.65103