#### Abstract

A new iterative technique is employed to solve a system of nonlinear fractional partial differential equations. This new approach requires neither Lagrange multiplier like variational iteration method (VIM) nor polynomials like Adomian's decomposition method (ADM) so that can be more easily and effectively established for solving nonlinear fractional differential equations, and will overcome the limitations of these methods. The obtained numerical results show good agreement with those of analytical solutions. The fractional derivatives are described in Caputo sense.

#### 1. Introduction

In recent years, it has been turned out that fractional differential equations can be used successfully to model many phenomena in various fields such as fluid mechanics, viscoelasticity, physics, chemistry, and engineering. For instance, the fluid−dynamics traffic model with fractional derivatives [1] is able to eliminate the deficiency arising from the assumption of continuum traffic flow, and the nonlinear oscillation of earthquakes can be modeled by fractional derivatives [2]. Fractional differentiation and integration operators can also be used for extending the diffusion and wave equations [3]. Most of fractional differential equations do not have exact analytical solutions, hence considerable heed has been focused on the approximate and numerical solutions of these equations. Although variational iteration method [4−8] and Adomian’s decomposition method [9−14] are approaches that have been utilized extensively to provide analytical approximations of linear and nonlinear problems, they have limitations due to complicated algorithms of calculating Adomian polynomials for nonlinear fractional problems, and an inherent inaccuracy in determining the Lagrange multiplier for fractional equations. In this study, a new alternative procedure that needs no Lagrange multiplier or Adomian polynomials is used to obtain an analytical approximate solution of a system of nonlinear fractional partial differential equations (1.1) to illustrate the effectiveness, accuracy, and convenience of this method.

In this work, we consider the solution of generalized Hirota−Satsuma coupled KdV of time−fractional order which is presented by a system of nonlinear partial differential equations, of the form:

subject to the following initial conditions:

where and are arbitrary constants. The Hirota−Satsuma system of equations [15] was introduced to describe the interaction of two long waves with different dispersion relations. The case of in system (1.1) was solved by Wu et al. [16]. If then in one case the , and are travelling−wave solutions of system (1.1) when .

#### 2. Basic Definitions

In this section, there are some basic definitions and properties of the fractional calculus theory which are used in this paper.

*Definition 2.1. *A real function , is mentioned to be in the space if there exists a real number such that where and it is said to be in the space if .

*Definition 2.2. *The left−sided Riemann−Liouville fractional integral operator of order , of a function is defined as

The properties of the operator can be found in [1, 17], and we only mention the following (in this case, and ):

The Riemann−Liouville derivative has certain disadvantages, in trying way to model real−world phenomena with fractional differential equations. Therefore, we will employ a modification of fractional differential operator proposed by Caputo, in his work [18] on the theory of viscoelasticity.

*Definition 2.3. *The (left−sided) Caputo fractional derivative of is defined as
for

Also, we need two of its basic properties.

Lemma 2.4. *If and , then one has:
*

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

In this work, we consider the one−dimensional linear nonhomogeneous fractional partial differential equations in fluid mechanics, where the unknown function is assumed to be a causal function of time, that is, vanishing for .

*Definition 2.5. *For as the smallest integer that exceeds , the Caputo time−fractional derivative operator of order is defined as

For more information on the mathematical properties of fractional derivatives and integrals, one can consult the mentioned references.

#### 3. Basic Ideas of Fractional Iteration Method (FIM)

As pointed in [19], to illustrate fractional iteration method, we consider the following nonlinear fractional differential equation (more general form can be considered without loss of generality):

where the fractional differential operator is de*ﬁ*ned as in (2.3), , is a nonlinear function of , and is an unknown function to be determined later. We want to find a solution of (3.1) having the form

Let denote the so−called auxiliary function. Multiplying (3.1) by and then applying , the Riemann−Liouville fractional integral operator, of order defined by (2.1), on both sides of the resulted term yields

Let denote the so−called auxiliary parameter. Multiplying (3.3) by and then adding , the solution of (3.1), on both sides of the resulted term yields

However (3.4) can be solved iteratively as follows:

In (3.5), the subscript *n* denotes the *nth* iteration, and provided that the right hand of it, that is, , is a contractive mapping. The convergence of (3.5) is ensured by Banach’s fixed point theorem [20], as is shown in [19].

Now, we introduce a new convenient technique for controlling the convergence region and rate of solution series for this method. Assume that we gain a family of solution series in the auxiliary parameter by the means of fractional iteration method. Like HAM, by plotting the amount of the function or one of its derivatives at a particular point with respect to the auxiliary parameter (which is the so−called *−curve*), we can obtain a proper value of that ensures the convergence of the solution series. This proper value of corresponds to the curve segment nearly parallel to the horizontal axis in the *−curve* plot. Therefore, if we set any value in this region, which is so−called *the valid region of *, we are quite sure that the corresponding solution series converge.

Having freedom for choosing the auxiliary function , the auxiliary parameter , the initial approximation and the fractional integral order that is, fundamental to the validity and flexibility of the FIM, we can suppose that all of them are properly chosen, therefore, the iterative scheme (3.5) will converge to the exact solution. Accordingly, the successive approximations , , of the solution will be obtained by choosing that at least satis*ﬁ*es the initial and/or boundary conditions. Consequently, the exact solution may be obtained by using .

#### 4. Applications

In this section, we implement fractional iteration method to generalized Hirota−Satsuma coupled KdV of time−fractional order when . For convenience in applying FIM method, we choose the initial conditions given in (1.2) as the initial approximations:

Choosing and , we can construct the iterative scheme (3.5) for investigation of the traveling wave solution of (1.1) as follows:

Substituting the initial approximations, (4.1), into (4.2), for the case , yields

#### 5. Result and Discussion

In this section, four figures are presented corresponding to FIM results and exact solutions for the solitary wave solutions and with the initial conditions (1.2), when , , . Furthermore, numerical values for the case *, *and , are obtained for and .

Demonstrating the exactness of FIM, the numerical results are presented and only few iterations are required to achieve accurate solutions. The convergence of FIM for the generalized fractional−order Hirota−Satsuma−coupled KdV equation is controllable, using the so−called *−curves* presented in Figure 1 which are obtained based on the fourth−order FIM approximate solutions. In general, by the means of the so−called *−curve*, it is straight forward to choose a proper value of which ensures that the solution series is convergent. This proper value of corresponds to the curve segment nearly parallel to the horizontal axis. Both exact results and approximate solutions obtained for the first four approximations are plotted in Figures 2, 3, and 4. There are no visible differences in two solutions of each pair of diagrams.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

Tables 1, 2, and 3 show the numerical values by FIM when and , for and respectively.

#### 6. Conclusion

In this paper, the fractional iteration method (FIM) has been successfully applied to study Hirota−Satsuma−coupled KdV of time−fractional−order equation. FIM results are compared with the exact solutions and those obtained by Homotopy perturbation method [21].

The results show that fractional iteration method is a powerful and efficient technique in finding exact and approximate solutions for nonlinear partial differential equations of fractional order. The method provides the user with more realistic series solutions that converge very rapidly in real physical problems.

Compared with the ADM and VIM, the FIM has following advantages, [19].

(1)The auxiliary parameter provides us with a convenient way to modify and control the convergence region of the solution.(2)The solution of a given nonlinear problem can be expressed by an infinite number of solution series and thus can be more efficiently approximated by a better selection of the auxiliary parameter values.(3)Unlike the ADM, the FIM method is free from the need to use Adomian polynomials.(4)This method has no need for the Lagrange multiplier, correction functional, stationary conditions, the variational theory, and so forth, which eliminates the complications that exist in the VIM.(5)The fractional iteration method can be easily comprehended with only a basic knowledge of fractional calculus.(6)Compared to the ADM and VIM, the presented method proves simpler in its principles and more convenient for computer algorithms.In this work, we used Maple Package to calculate the series obtained by fractional iteration method.