Abstract and Applied Analysis

Volume 2013 (2013), Article ID 561980, 11 pages

http://dx.doi.org/10.1155/2013/561980

Research Article

Approximate Analytic Solutions of Time-Fractional Hirota-Satsuma Coupled KdV Equation and Coupled MKdV Equation

Jincun Liu and Hong Li

School of Mathematical Sciences, Inner Mongolia University, 235 West Daxue Road, Hohhot 010021, China

Received 29 October 2012; Accepted 28 January 2013

Academic Editor: Abdel-Maksoud A. Soliman

Copyright © 2013 Jincun Liu and Hong Li. 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.

Abstract.

By introducing the fractional derivative in the sense of Caputo and combining the pretreatment technique to deal with long nonlinear items, the generalized two-dimensional differential transform method is proposed for solving the time-fractional Hirota-Satsuma coupled KdV equation and coupled MKdV equation. The presented method is a numerical method based on the generalized Taylor series expansion which constructs an analytical solution in the form of a polynomial. The numerical results show that the generalized two-dimensional differential transform method is very effective for the fractional coupled equations.

1. Introduction

In the last past decade, the fractional differential equations have been widely used in various fields of physics and engineering. The analytical approximation of such problems has attracted great attention and became a considbased onerable interest in mathematical physics. Some powerful methods including the homotopy perturbation method [1], Adomian decomposition method [2, 3], variational iteration method [4, 5], homotopy analysis method [6, 7], fractional complex transform method [8], and generalized differential transform method [9] have been developed to obtain exact and approximate analytic solutions. These solution techniques are more clear and realistic methods for fractional differential equations, because they give the approximate solutions of the considered problems without any linearization or discretization.

The variational iteration method and the homotopy perturbation method were first proposed by Professor He in [10, 11], respectively. The idea of the variational iteration method is to construct correction functionals using general Lagrange multipliers identified optimally via the variational theory, and the initial approximations can be freely chosen with unknown constants. Recently, Wu and Lee proposed a fractional variational iteration method for fractional differential equation based on the modified Riemann CLiouville derivative, which is more effective to solve fractional differential equation [12]. This method has been developed by many authors, see [13–16] and the references cited therein. The homotopy perturbation method, which does not require a small parameter in an equation, has a significant advantage that it provides an analytical approximate solution to a wide range of nonlinear problems in applied sciences. Recently, the fractional complex transform is developed to convert the fractional differential equation to its differential partner and gave a geometrical explanation [8]. These methods are more effective for solving the linear and nonlinear fractional differential equations.

The differential transform method was used firstly by Zhou in 1986 to study electric circuits [17]. The differential transform is an iterative procedure based on the Taylor series expansion which constructs an analytic solution in the form of a polynomial. The method is well addressed in [18–23]. Recently, the generalized differential transform method is developed for obtaining approximate analytic solutions for some linear and nonlinear differential equations of fractional order [24, 25].

The aim of this paper is to directly extend the generalized two-dimensional differential transform method to obtain the approximate analytic solutions of a time-fractional Hirota-Satsuma coupled KdV equation,

and a time-fractional coupled MKdV equation,

where is a constant and and are parameters describing the order of the time-fractional derivatives of , , and , respectively. The fractional derivatives are considered in the Caputo sense. In the case of and , (1) and (2) reduce to the classical Hirota-Satsuma coupled KdV equation and coupled MKdV equation [26], respectively.

The Caputo fractional derivative is considered here because it allows traditional initial and boundary conditions to be included in the formulation of the problem.

Definition 1 (see [27]). The fractional derivative of in the Caputo sense is defined as

for and . Here, is the Riemann-Liouville integral operator of order , defined by

Definition 2 (see [28]). For to be the smallest integer that exceeds , the Caputo time-fractional derivative of order is defined as

and the space-fractional derivatives of Caputo type can be defined analogously.

Lemma 3 (see [29]). The Caputo fractional derivative of the power function satisfies

2. Generalized Two-Dimensional Differential Transform Method (GDTM)

Consider a function of two variables and suppose that it can be represented as a product of two single variable functions, that is, . Based on the properties of generalized two-dimensional differential transform, the function can be represented as

where , is called the spectrum of . The generalized two-dimensional differential transform of the function is given by

where , is times. In the case of and , the generalized two-dimensional differential transform (8) reduces to two-dimensional differential transform.

The fundamental theorems of the generalized two-dimensional differential transform are as follows.

Theorem 4. Suppose that , and are the differential transformations of the functions , and , respectively:

(a)if , then .(b)if , then .(c)if , then .(d)if , then .(e)if , then .

Theorem 5. If , function , where and has the generalized Taylor series expansion ,

(i) and arbitrary, or(ii), arbitrary, and for , where ,

then the generalized differential transform (8) becomes

If , and , then

Some details of the aformentioned theorems can be found in [30].

3. Applications of GDTM

3.1. Fractional Hirota-Satsuma Coupled KdV Equation

Consider the following time-fractional Hirota-Satsuma coupled KdV equation:

subject to the initial conditions

where , , , and are arbitrary constants.

The exact solutions of (11) and (12), for the special case of , given in [26], are

Suppose that the solutions , , and can be represented as the products of single-valued functions, respectively. Applying the generalized two-dimensional differential transform to both sides of (11) and using the related theorems, we have

The generalized two-dimensional differential transforms of the initial conditions can be obtained as follows:

Utilizing the recurrence relations (14) and the transformed initial conditions, we can obtain all the , , and with the help of Mathematica. Moreover, substituting all into (7), we obtain the series form solution

where

The closed forms of , , and are

Similarly, substituting all and all into (7), respectively, we obtain the series form solutions

The closed forms of , and are

The closed forms of , , and are

The approximate solutions of (11) in finite series forms are given by

In order to verify whether the approximate solutions of (22)–(24) lead to higher accuracy, we draw the figures of the approximate solutions of (22)–(24) with , as well as the exact solutions (13) when , , and . It can be seen from Figures 1(a) to 3(b) that the solutions obtained by the presented method is nearly identical with the exact solutions. So, we conclude that a good approximation is achieved by using the GDTM.

(a)

(b)

Figure 1: The surface shows the solution of (11): (a) exact solution; (b) approximate solution of (22) when , , and .

In the following, we will construct an approximate solution of (11) with the new initial conditions,

where , and are arbitrary constants.

The generalized two-dimensional differential transforms of the initial conditions of (25) are given by

Utilizing the recurrence relations in (14) and the transformed initial conditions, we can obtain the following approximate solutions:

where the closed forms of , and are

3.2. Fractional Coupled MKdV Equation

Consider the following time-fractional coupled MKdV equation:

subject to the initial conditions

where is an arbitrary constant.

The exact solutions of (29) and (30), for the special case of , given in [26], are

Equation system of (29) is more complex. In order to obtain an explicit iteration scheme, we follow the pretreatment technique introduced recently by Chang [20, 21] to deal with long nonlinear items. Firstly, we suppose that

According to (32), (29) can be equivalently written as the following form:

Suppose that the solutions and can be represented as the products of single-valued functions, respectively. Applying the generalized two-dimensional differential transform to both sides of (33) and using the related theorems, we have

Herein , and denote the differential transformations of the functions , and , respectively.

The generalized two-dimensional differential transforms of the initial conditions of (30) can be obtained as follows:

Utilizing the recurrence relations of (34) and the transformed initial conditions, we can obtain all the and with the help of Mathematica. Through the complex calculation which is similar to the solving process in Section 3.1, we have the approximate solutions of (29) in finite series,

where the closed forms of and are given by

The effectiveness and accuracy of the approximate solutions can be seen from the comparison figures.

Figures 4 and 5 show the approximate solutions of (36) and (37) and the exact ones of (31) with and , respectively. Comparing Figures 1(a) and 1(b) with Figures 2(a) and 2(b), we can see that the solutions obtained by different methods are nearly identical. From these figures, we can know that the series solutions converge rapidly, so a good approximation has been achieved.

(a)

(b)

Figure 2: The surface shows the solution of (11): (a) exact solution; (b) approximate solution of (23) when , , and .

(a)

(b)

Figure 3: The surface shows the solution of (11): (a) exact solution; (b) approximate solution (24) when , and .

(a)

(b)

Figure 4: The surface shows the exact solutions (31): (a) ; (b) when and .

(a)

(b)

Figure 5: The surface shows the approximate solutions of (29): (a) of (36); (b) of (37) when and .

4. Summary and Discussion

In this paper, combining the Caputo fractional derivative, the GDTM was applied to derive approximate analytical solutions of the time-fractional Hirota-Satsuma coupled KdV equation and coupled MKdV equation with initial conditions. The numerical solutions obtained from the GDTM are shown graphically. The obtained results demonstrate the reliability of the algorithm and its wider applicability to nonlinear fractional coupled partial differential equations.

In [8] and the references cited therein, the so-called fractional complex transform (FCT) is suggested to convert a fractional differential equation with Jumarie's modification of Riemann-Liouville derivative into its classical differential partner. According to the idea of FCT, for some fractional differential equations with Caputo time-fractional derivative , we assume that

consequently, from Lemma 3,

and it can be seen easily that

We can point out that the approximate solutions of (11) and (29) by GDTM can be derived by the similar method compared with FCT. Without loss of generality, we only consider (29) with the initial conditions of (30). In fact, suppose that

From (33), (40), and (41), we have the following iteration formulae:

for .

Setting and , and using (43), we can obtain the 3-order approximate solutions as (36) and (37) by using the symbol computational software Mathematica.

Acknowledgments

The project is supported by the NNSF of China (11061021), the NSF-IMU of China (2012MS0105, 2012MS0106), and the YSF-IMU of China (ND0811).

References

S. Momani and Z. Odibat, “Homotopy perturbation method for nonlinear partial differential equations of fractional order,” Physics Letters A, vol. 365, no. 5-6, pp. 345–350, 2007.

C. Li and Y. Wang, “Numerical algorithm based on Adomian decomposition for fractional differential equations,” Computers & Mathematics with Applications, vol. 57, no. 10, pp. 1672–1681, 2009.

S. R. Santanu, “Analytical solution for the space fractional diffusion equation by two-step Adomian decomposition method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1295–1306, 2009.

M. R. Yulita, M. S. M. Noorani, and I. Hashim, “Variational iteration method for fractional heat- and wave-like equations,” Nonlinear Analysis. Real World Applications, vol. 10, no. 3, pp. 1854–1869, 2009.

J. H. He, “Approximate analytical solution for seepage flow with fractional derivatives in porous media,” Computer Methods in Applied Mechanics and Engineering, vol. 167, no. 1-2, pp. 57–68, 1998.

H. Jafari and S. Seifi, “Solving a system of nonlinear fractional partial differential equations using homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 1962–1969, 2009.

H. Xu, S. J. Liao, and X.-C. You, “Analysis of nonlinear fractional partial differential equations with the homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1152–1156, 2009.

J. H. He, S. K. Elagan, and Z. B. Li, “Geometrical explanation of the fractional complex transform and derivative chain rule for fractional calculus,” Physics Letters A, vol. 376, no. 4, pp. 257–259, 2012.

Z. Odibat and S. Momani, “A generalized differential transform method for linear partial differential equations of fractional order,” Applied Mathematics Letters, vol. 21, no. 2, pp. 194–199, 2008.

J. H. He, “Approximate solution of nonlinear differnetial equations with convolution product nonlinearities,” Computer Methods in Applied Mechanics and Engineering, vol. 167, pp. 69–73, 1998.

J. H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.

G. Wu and E. W. M. Lee, “Fractional variational iteration method and its application,” Physics Letters A, vol. 374, no. 25, pp. 2506–2509, 2010.

G. Wu, “A fractional variational iteration method for solving fractional nonlinear differential equations,” Computers & Mathematics with Applications, vol. 61, no. 8, pp. 2186–2190, 2011.

J. H. He, “A short remark on fractional variational iteration method,” Physics Letters A, vol. 375, no. 38, pp. 3362–3364, 2011.

J. Ji, J. Zhang, and Y. Dong, “The fractional variational iteration method improved with the Adomian series,” Applied Mathematics Letters, vol. 25, no. 12, pp. 2223–2226, 2012.

J. H. He and X.-H. Wu, “Variational iteration method: new development and applications,” Computers & Mathematics with Applications, vol. 54, no. 7-8, pp. 881–894, 2007.

J. K. Zhou, Differential TransFormation and Its Applications for Electric CircuIts, Wuhan University Press, Wuhan, China, 1986, (in Chinese).

A. S. V. Ravi Kanth and K. Aruna, “Solution of singular two-point boundary value problems using differential transformation method,” Physics Letters A, vol. 372, no. 26, pp. 4671–4673, 2008.

A. Fatma, “Solutions of the system of differential equations by differential transform method,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 547–567, 2004.

S. H. Chang and I. L. Chang, “A new algorithm for calculating one-dimensional differential transform of nonlinear functions,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 799–805, 2008.

S. H. Chang and I.-L. Chang, “A new algorithm for calculating two-dimensional differential transform of nonlinear functions,” Applied Mathematics and Computation, vol. 215, no. 7, pp. 2486–2494, 2009.

K. C. Chen and S. H. Ho, “Solving partial differential equations by two-dimensional differential transform method,” Applied Mathematics and Computation, vol. 106, no. 2-3, pp. 171–179, 1999.

A. S. V. Ravi Kanth and K. Aruna, “Differential transform method for solving the linear and nonlinear Klein-Gordon equation,” Computer Physics Communications, vol. 180, no. 5, pp. 708–711, 2009.

S. Momani, Z. Odibat, and V. S. Erturk, “Generalized differential transform method for solving a space- and time-fractional diffusion-wave equation,” Physics Letters A, vol. 370, no. 5-6, pp. 379–387, 2007.

M. Kurulay and M. Bayram, “Approximate analytical solution for the fractional modified KdV by differential transform method,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 7, pp. 1777–1782, 2010.

E. G. Fan, “Soliton solutions for a generalized Hirota-Satsuma coupled KdV equation and a coupled MKdV equation,” Physics Letters A, vol. 282, no. 1-2, pp. 18–22, 2001.

I. Podlubny, Fractional Differential Equations, Academic Press, New York, NY, USA, 1999.

M. Caputo, “Linear models of dissipation whose Q is almost frequency independent. Part II,” A Journal of the Royal Statistical Society, vol. 13, pp. 529–539, 1967.

M. Ishteva, R. Scherer, and L. Boyadjiev, “On the caputo operator of fractional calculus and c-Laguerrefunctions,” http://www.cc.gatech.edu/~mishteva/paperMSRJ.pdf

S. Momani and Z. Odibat, “A novel method for nonlinear fractional partial differential equations: combination of DTM and generalized Taylor's formula,” Journal of Computational and Applied Mathematics, vol. 220, no. 1-2, pp. 85–95, 2008.