About this Journal Submit a Manuscript Table of Contents
International Journal of Differential Equations
Volume 2011 (2011), Article ID 545607, 15 pages
http://dx.doi.org/10.1155/2011/545607
Research Article

Solving Famous Nonlinear Coupled Equations with Parameters Derivative by Homotopy Analysis Method

Department of Applied Mathematics, Ferdowsi University of Mashhad, P.O. Box 1159, Mashhad 91779-48974, Iran

Received 15 May 2011; Accepted 4 July 2011

Academic Editor: Shaher M. Momani

Copyright © 2011 Sohrab Effati et al. 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

The homotopy analysis method (HAM) is employed to obtain symbolic approximate solutions for nonlinear coupled equations with parameters derivative. These nonlinear coupled equations with parameters derivative contain many important mathematical physics equations and reaction diffusion equations. By choosing different values of the parameters in general formal numerical solutions, as a result, a very rapidly convergent series solution is obtained. The efficiency and accuracy of the method are verified by using two famous examples: coupled Burgers and mKdV equations. The obtained results show that the homotopy perturbation method is a special case of homotopy analysis method.

1. Introduction

Fractional differential equations have gained importance and popularity during the past three decades or so, mainly due to its demonstrated applications in numerous seemingly diverse fields of science and engineering. For example, the nonlinear oscillation of earthquake can be modeled with fractional derivatives, and the fluid-dynamic traffic model with fractional derivatives can eliminate the deficiency arising from the assumption of continuum traffic flow. The differential equations with fractional order have recently proved to be valuable tools to the modeling of many physical phenomena [1, 2]. This is because of the fact that the realistic modeling of a physical phenomenon does not depend only on the instant time, but also on the history of the previous time which can also be successfully achieved by using fractional calculus. Most nonlinear fractional equations do not have exact analytic solutions, so approximation and numerical techniques must be used. The Adomain decomposition method [3], the homotopy perturbation method [4], the variational iteration method [5], and other methods have been used to provide analytical approximation to linear and nonlinear problems. However, the convergence region of the corresponding results is rather small. In 1992, Liao employed the basic ideas of the homotopy in topology to propose a general analytic method for nonlinear problems, namely, homotopy analysis method [610]. This method has been successfully applied to solve many types of nonlinear problems in science and engineering, such as the viscous flows of non-Newtonian fluids [11], the KdV-type equations [12], finance problems [13], fractional Lorenz system [14], and delay differential equation [15]. The HAM contains a certain auxiliary parameter which provides us with a simple way to adjust and control the convergence region and rate of convergence of the series solution.

The HAM offers certain advantages over routine numerical methods. Numerical methods use discretization which gives rise to rounding off errors causing loss of accuracy and requires large computer memory and time. This computational method yields analytical solutions and has certain advantages over standard numerical methods. The HAM method is better since it does not involve discretization of the variables and hence is free from rounding off errors and does not require large computer memory or time.

In this paper, we extend the application of HAM to discuss the explicit numerical solutions of a type of nonlinear-coupled equations with parameters derivative in this form: 𝜕𝛼𝑢𝜕𝑡𝛼=𝐿1(𝑢,𝑣)+𝑁1𝜕(𝑢,𝑣),𝑡>0,𝛽𝑣𝜕𝑡𝛽=𝐿2(𝑢,𝑣)+𝑁2(𝑢,𝑣),𝑡>0,(1.1) where 𝐿𝑖 and 𝑁𝑖 (𝑖=1,2) are the linear and nonlinear functions of 𝑢 and 𝑣, respectively, 𝛼 and 𝛽 are the parameters that describe the order of the derivative. Different nonlinear coupled systems can be obtained when one of the parameters 𝛼 or 𝛽 varies. The study of (1.1) is very necessary and significant because when 𝛼 and 𝛽 are integers, it contains many important mathematical physics equations.

The paper has been organized as follows. Notations and basic definitions are given in Section 2. In Section 3 the homotopy analysis method is described. In Section 4 applying HAM for two famous coupled examples: Burgers and mKdV equations. Discussion and conclusions are presented in Section 5.

2. Description on the Fractional Calculus

Definition 2.1. A real function 𝑓(𝑡),𝑡>0 is said to be in the space 𝐶𝜇, 𝜇𝑅 if there exists a real number 𝑝>𝜇, such that 𝑓(𝑡)=𝑡𝑝𝑓1(𝑡) where 𝑓1(0,), and it is said to be in the space 𝐶𝜇𝑛 l if and only if (𝑛)𝐶𝜇, 𝑛𝑁. Clearly 𝐶𝜇𝐶𝜈 if 𝜈𝜇.

Definition 2.2. The Riemann-Liouville fractional integral operator (𝐽𝛼) of order 𝛼0, of a function 𝑓𝐶𝜇, 𝜇1, is defined as 𝐽𝛼1𝑓(𝑥)=Γ(𝛼)𝑥0(𝑥𝑡)𝛼1𝐽𝑓(𝑡)𝑑𝑡,𝑥>0.0𝑓(𝑥)=𝑓(𝑥).(2.1)Γ(𝛼) is the well-known Gamma function. Some of the properties of the operator 𝐽𝛼, which we will need here, are as follows.
For 𝑓𝐶𝜇, 𝜇1, 𝛼,𝛽0 and 𝛾1𝐽𝛼𝐽𝛽𝑓(𝑥)=𝐽𝛼+𝛽𝐽𝑓(𝑥),𝛼𝐽𝛽𝑓(𝑥)=𝐽𝛽𝐽𝛼𝐽𝑓(𝑥),𝛼𝑡𝛾=Γ(𝛾+1)𝑡Γ(𝛼+𝛾+1)𝛼+𝛾.(2.2)

Definition 2.3. For the concept of fractional derivative, there exist many mathematical definitions [1, 1619]. In this paper, the two most commonly used definitions: the Caputo derivative and its reverse operator Riemann-Liouville integral are adopted. That is because Caputo fractional derivative [1] allows the traditional assumption of initial and boundary conditions. The Caputo fractional derivative is defined as 𝐷𝛼𝑡𝜕𝑢(𝑥,𝑡)=𝛼𝑢(𝑥,𝑡)𝜕𝑡𝛼=1Γ(𝑛𝛼)𝑡0(𝑡𝜏)𝑛𝛼1𝜕𝑛𝑢(𝑥,𝑡)𝜕𝑡𝑛𝜕𝑑𝜏,𝑛1<𝛼<𝑛,𝑛𝑢(𝑥,𝑡)𝜕𝑡𝑛,𝛼=𝑛𝑁.(2.3) Here, we also need two basic properties about them: 𝐷𝛼𝐽𝛼𝐽𝑓(𝑥)=𝑓(𝑥),𝛼𝐷𝛼𝑓(𝑥)=𝑓(𝑥)𝑘=0𝑓(𝑘)0+𝑥𝑘𝑘!,𝑥>0.(2.4)

Definition 2.4. The Mittag-Leffler function 𝐸𝛼(𝑧) with 𝑎>0 is defined by the following series representation, valid in the whole complex plane: 𝐸𝛼(𝑧)=𝑛=0𝑧𝑛Γ(𝛼𝑛+1),𝛼>0,𝑧C.(2.5)

3. Basic Idea of HAM

To describe the basic ideas of the HAM, we consider the operator form of (1.1):𝑁𝐷𝛼𝑡𝑢𝑁𝐷(𝑥,𝑡)=0,𝑡>0,𝛽𝑡𝑣(𝑥,𝑡)=0,𝑡>0,(3.1) where 𝑁 is nonlinear operator, 𝐷𝛼𝑡 and 𝐷𝛽𝑡 stand for the fractional derivative and are defined as in (2.3), 𝑡 denotes an independent operator, and 𝑢(𝑥,𝑡), 𝑣(𝑥,𝑡) are unknown functions.

By means of generalizing the traditional homotopy method, Liao [6] constructs the so-called zero-order deformation equations:𝜙(1𝑞)𝐿1(𝑥,𝑡,𝑞)𝑢0𝐷(𝑥,𝑡)=𝑞𝐻(𝑡)𝑁𝛼𝑡𝜙1𝜙(𝑥,𝑡,𝑞),(3.2)(1𝑞)𝐿2(𝑥,𝑡,𝑞)𝑣0𝐷(𝑥,𝑡)=𝑞𝐻(𝑡)𝑁𝛽𝑡𝜙2(𝑥,𝑡,𝑞),(3.3) where 𝑞[0,1] is the embedding parameter, 0 is a non-zero auxiliary parameter, 𝐻(𝑡)0 is an auxiliary function, 𝐿 is an auxiliary linear operator, 𝑢0(𝑥,𝑡), 𝑣0(𝑥,𝑡) are initial guesses of 𝑢(𝑥,𝑡), 𝑣(𝑥,𝑡) and 𝜙1(𝑥,𝑡,𝑞), 𝜙2(𝑥,𝑡,𝑞) are two unknown functions, respectively. It is important that one has great freedom to choose auxiliary things in HAM. Obviously, when 𝑞=0 and 𝑞=1, the following holds: 𝜙1(𝑥,𝑡,0)=𝑢0(𝑥,𝑡),𝜙1𝜙(𝑥,𝑡,1)=𝑢(𝑥,𝑡),2(𝑥,𝑡,0)=𝑣0(𝑥,𝑡),𝜙2(𝑥,𝑡,1)=𝑣(𝑥,𝑡),(3.4) respectively. Thus, as 𝑞 increases from 0 to 1, the solution 𝜙1(𝑥,𝑡,𝑞), 𝜙2(𝑥,𝑡,𝑞) varies from the initial guess 𝑢0(𝑥,𝑡), 𝑣0(𝑥,𝑡) to the solution 𝑢(𝑥,𝑡), 𝑣(𝑥,𝑡). Expanding 𝜙1(𝑥,𝑡,𝑞), 𝜙2(𝑥,𝑡,𝑞) in Taylor series with respect to 𝑞, we have𝜙1(𝑥,𝑡,𝑞)=𝑢0(𝑥,𝑡)++𝑚=1𝑢𝑚(𝑥,𝑡)𝑞𝑚,𝜙2(𝑥,𝑡,𝑞)=𝑣0(𝑥,𝑡)++𝑚=1𝑣𝑚(𝑥,𝑡)𝑞𝑚,(3.5) where𝑢𝑚1(𝑥,𝑡)=𝜕𝑚!𝑚𝜙1(𝑥,𝑡,𝑞)𝜕𝑞𝑚||𝑞=0,𝑣𝑚1(𝑥,𝑡)=𝜕𝑚!𝑚𝜙2(𝑥,𝑡,𝑞)𝜕𝑞𝑚||𝑞=0.(3.6) If the auxiliary linear operator, the initial guess, the auxiliary parameter , and the auxiliary function are so properly chosen, the series (3.5) converges at 𝑞=1, then we have𝑢(𝑥,𝑡)=𝑢0(𝑥,𝑡)++𝑚=1𝑢𝑚(𝑥,𝑡),𝑣(𝑥,𝑡)=𝑣0(𝑥,𝑡)++𝑚=1𝑣𝑚(𝑥,𝑡),(3.7) which must be one of solutions of original nonlinear equation, as proved by Liao [8]. As =1 and 𝐻(t)=1, (3.2) and (3.3) become𝜙(1𝑞)𝐿1(𝑥,𝑡,𝑞)𝑢0𝜙(𝑥,𝑡)+𝑞𝑁1(𝜙(𝑥,𝑡,𝑞)=0,1𝑞)𝐿2(𝑥,𝑡,𝑞)𝑢0(𝜙𝑥,𝑡)+𝑞𝑁2(𝑥,𝑡,;𝑞)=0,(3.8) which is used mostly in the homotopy perturbation method [20], where as the solution obtained directly, without using Taylor series. According to the definition (3.6), the governing equation can be deduced from the zero-order deformation equation (3.2). Define the vector 𝑢𝑛=𝑢0(𝑥,𝑡),𝑢1(𝑥,𝑡),,𝑢𝑛,𝑣(𝑥,𝑡)𝑛=𝑣0(𝑥,𝑡),𝑣1(𝑥,𝑡),,𝑣𝑛(𝑥,𝑡).(3.9) Differentiating equations (3.2) and (3.3) 𝑚 times with respect to the embedding parameter 𝑞 and then setting 𝑞=0 and finally dividing them by 𝑚!, we have the so-called 𝑚th-order deformation equation: 𝐿𝑢𝑚(𝑥,𝑡)𝜒𝑚𝑢𝑚1(𝑥,𝑡)=𝐻(𝑡)𝑅1,𝑚𝑢𝑚1,𝐿𝑣𝑚(𝑥,𝑡)𝜒𝑚𝑣𝑚1(𝑥,𝑡)=𝐻(𝑡)𝑅2,𝑚𝑣𝑚1,(3.10) where 𝑅1,𝑚𝑢𝑚1=1𝜕(𝑚1)!𝑚1𝐷𝛼𝑡𝜙1(𝑥,𝑡,𝑞)𝜕𝑞𝑚1||𝑞=0,𝑅2,𝑚𝑣𝑚1=1𝜕(𝑚1)!𝑚1𝐷𝛽𝑡𝜙2(𝑥,𝑡,𝑞)𝜕𝑞𝑚1||𝑞=0,𝜒𝑚=0,𝑚1,1,𝑚>1.(3.11) Applying the Riemann-Liouville integral operator 𝐽𝛼,𝐽𝛽 on both side of (3.10), we have𝑢𝑚(𝑥,𝑡)=𝜒𝑚𝑢𝑚1(𝑥,𝑡)𝜒𝑚𝑛1𝑖=0𝑢𝑖𝑚10+𝑡𝑖𝑖!+𝐻(𝑡)𝐽𝛼𝑅1,𝑚𝑢𝑚1,𝑣𝑚(𝑥,𝑡)=𝜒𝑚𝑣𝑚1(𝑥,𝑡)𝜒𝑚𝑛1𝑖=0𝑣𝑖𝑚10+𝑡𝑖𝑖!+𝐻(𝑡)𝐽𝛽𝑅2,𝑚𝑣𝑚1.(3.12) It should be emphasized that 𝑢𝑚(𝑥,𝑡), 𝑣𝑚(𝑥,𝑡) for 𝑚1 is governed by the linear equation (3.10), under the linear boundary conditions that come from original problem, which can be easily solved by symbolic computation software such as MATLAB. For the convergence of the above method we refer the reader to Liao's work. Liao [7] proved that, as long as a series solution given by the homotopy analysis method converges, it must be one of exact solutions. So, it is important to ensure that the solution series is convergent. Note that the solution series contain the auxiliary parameter , which we can choose properly by plotting the so-called -curves to ensure solution series converge.

Remark 3.1. The parameters 𝛼 and 𝛽 can be arbitrarily chosen as, integer or fraction, bigger or smaller than 1. When the parameters are bigger than 1, we will need more initial and boundary conditions such as 𝑢0(𝑥,0), 𝑢0(𝑥,0), and the calculations will become more complicated correspondingly. In order to illustrate the problem and make it convenient for the readers, we only confine the parameters to [0,1] to discuss.

4. Application

4.1. The Nonlinear Coupled Burgers Equations with Parameters Derivative

In order to illustrate the method discussed above, we consider the nonlinear coupled Burgers equations with parameters derivative in an operator form:𝐷𝛼𝑡𝑢𝐿𝑥𝑥𝑢2𝑢𝐿𝑥𝑢+𝐿𝑥𝐷𝑢𝑣=0,(0<𝛼1),𝛽𝑡𝑣𝐿𝑥𝑥𝑣2𝑣𝐿𝑥𝑣+𝐿𝑥𝑢𝑣=0,(0<𝛽1),(4.1) where 𝑡>0, 𝐿𝑥=𝜕/𝜕𝑥 and the fractional operators 𝐷𝛼𝑡 and 𝐷𝛽𝑡 are defined as in (2.3). Assuming the initial value as𝑢(𝑥,0)=sin𝑥,𝑣(𝑥,0)=sin𝑥.(4.2) The exact solutions of (4.1) for the special case: 𝛼=𝛽=1 are𝑢(𝑥,𝑡)=𝑒𝑡sin𝑥,𝑣(𝑥,𝑡)=𝑒𝑡sin𝑥.(4.3) For application of homotopy analysis method, in view of (4.1) and the initial condition given in (4.2), it is convenient to choose𝑢0(𝑥,𝑡)=sin𝑥,𝑣0(𝑥,𝑡)=sin𝑥,(4.4) as the initial approximate of (4.1). We choose the linear operators𝐿1𝜙1(𝑥,𝑡,𝑞)=𝐷𝛼𝑡𝜙1,𝐿(𝑥,𝑡,𝑞)2𝜙2(𝑥,𝑡,𝑞)=𝐷𝛽𝑡𝜙2,(𝑥,𝑡,𝑞)(4.5) with the property 𝐿(𝑐)=0 where 𝑐 is constant of integration. Furthermore, we define a system of nonlinear operators as 𝑁1𝜙𝑖(𝑥,𝑡,𝑞)=𝐷𝛼𝑡𝜙1𝜕(𝑥,𝑡,𝑞)2𝜙1(𝑥,𝑡,𝑞)𝜕𝑥22𝜙1(𝑥,𝑡,𝑞)𝜕𝜙1(𝑥,𝑡,𝑞)+𝜕𝜙𝜕𝑥1(𝑥,𝑡,𝑞)𝜙2(𝑥,𝑡,𝑞)𝑁𝜕𝑥2𝜙𝑖(𝑥,𝑡,𝑞)=𝐷𝛽𝑡𝜙2𝜕(𝑥,𝑡,𝑞)2𝜙2(𝑥,𝑡,𝑞)𝜕𝑥22𝜙2(𝑥,𝑡,𝑞)𝜕𝜙2(𝑥,𝑡,𝑞)+𝜕𝜙𝜕𝑥1(𝑥,𝑡,𝑞)𝜙2(𝑥,𝑡,𝑞).𝜕𝑥(4.6) We construct the zeroth-order and the 𝑚th-order deformation equations where 𝑅1,𝑚𝑢𝑚1=𝐷𝛼𝑡𝑢𝑚1𝑢𝑚1𝑥𝑥2𝑚1𝑘=0𝑢𝑘𝑢𝑚1𝑘𝑥+𝑚1𝑘=0𝑢𝑘𝑣𝑚𝑘1𝑥,𝑅2,𝑚𝑣𝑚1=𝐷𝛽𝑡𝑣𝑚1𝑣𝑚1𝑥𝑥2𝑚1𝑘=0𝑣𝑘𝑣𝑚1𝑘𝑥+𝑚1𝑘=0𝑢𝑘𝑣𝑚𝑘1𝑥.(4.7) We start with an initial approximation 𝑢(𝑥,0)=sin(𝑥), 𝑣(𝑥,0)=sin(𝑥), thus we can obtain directly the other components as 𝑢1=𝑡𝑎sin(𝑥),𝑢Γ(𝑎+1)2=sin𝑥𝑎2Γ(𝑏+1)Γ(𝑎)2×𝑎Γ(𝑎+(1/2))2Γ(𝑏+1)Γ(𝑎)2Γ1𝑎+2𝑎2Γ(𝑏+1)Γ(𝑎)2Γ1𝑎+2+𝑎Γ(𝑎)Γ(2𝑎+1)𝑡𝑎Γ(𝑏+1)+21𝑎Γ(𝑎)Γ𝑎+2𝑡𝑎Γ(𝑏+1)+22Γ1𝑎Γ(𝑎)𝑎+2𝑡(𝑏+𝑎)cos(𝑥)22𝑡(2𝑎)1cos(𝑥)Γ(𝑏+1)Γ𝑎+2+2𝑡(2𝑎)𝑎2Γ(𝑏+1)Γ(𝑎)2𝑣1=𝑡𝑏sin(𝑥),𝑣Γ(𝑏+1)2=sin𝑥𝑏2Γ(𝑎+1)Γ(𝑏)2×𝑏Γ(𝑏+(1/2))2Γ(𝑎+1)Γ(𝑏)2Γ(2𝑏+1)+𝑏Γ(𝑏),Γ(2𝑏+1)𝑡𝑏Γ(𝑎+1)+2𝑏Γ(𝑏)Γ(2𝑏+1)𝑡𝑏Γ(𝑏+1)+22𝑡𝑏Γ(𝑏)Γ(2𝑏+1),(𝑎+𝑏)cos(𝑥)22𝑡(2𝑏)cos(𝑥)Γ(𝑎+1)Γ(2𝑏+1)+21𝑡(𝑏)𝑡𝑏𝑏2Γ(𝑎+1),Γ(𝑏)2(4.8) The absolute error of the 6th-order HAM and exact solution with =1 as shown in Figure 1. Also the absolute errors |𝑢(𝑡)𝜙6(𝑡)| have been calculated in Table 1. Figure 2 shows the numerical solutions of the nonlinear coupled Burgers equations with parameters derivative with =1, 𝛼=𝛽=1. Figure 3 shows the explicit numerical solutions with =1, 𝛼=1/4, and 𝛽=1/3 at 𝑡=0.02.

tab1
Table 1: The comparison of the results of the HAM (=1) and exact solution for the 𝑢(𝑥,𝑡), 𝛼=𝛽=1.
fig1
Figure 1: The comparison of the 6th-order HAM and exact solution with =1,𝛼=𝛽=1.
fig2
Figure 2: Explicit numerical solutions with=1, 𝛼=𝛽=1.
fig3
Figure 3: Explicit numerical solutions with =1, 𝛼=1/4, and𝛽=1/3.

As suggested by Liao [7], the appropriate region for is a horizontal line segment. We can investigate the influence of on the convergence of the solution series gevin by the HAM, by plotting its curve versus h, as shown in Figure 4.

545607.fig.004
Figure 4: The -curves obtained from the 5-order HAM approximate solution.

Remark 4.1. This example has been solved using homotopy perturbation method [21]. The graphs drawn and tables by =1 are in excellent agreement with that graphs drawn with HPM.

4.2. The Nonlinear Coupled mKdV Equations with Parameters Derivative

In order to illustrate the method discussed above, we consider the nonlinear coupled mKdV equations with parameters derivative in an operator form:𝐷𝛼𝑡1𝑢2𝑢𝑥𝑥𝑥+3𝑢2𝑢𝑥32𝑣𝑥𝑥3(𝑢𝑣)𝑥+3𝜆𝑢𝑥𝐷=0,𝛽𝑡𝑣+𝑣𝑥𝑥𝑥+3𝑣𝑣𝑥+3𝑢𝑥𝑣𝑥3𝑢2𝑣𝑥3𝜆𝑣𝑥=0,(4.9) with the initial conditions,𝑏𝑢(𝑥,0)=𝜆2𝑘+𝑘tanh(𝑘𝑥),𝑣(𝑥,0)=2𝑘1+𝑏+𝑏tanh(𝑘𝑥).(4.10) As we know, when 𝛼=𝛽=1 (4.9) has the kink-type soliton solutions𝑏𝑢(𝑥,𝑡)=𝜆2𝑘+𝑘tanh(𝑘𝜉),𝑣(𝑥,𝑡)=2𝑘1+𝑏+𝑏tanh(𝑘𝜉),(4.11) constructed by Fan [22], where 𝜉=𝑥+(1/4)(4𝑘26𝜆+6𝑘𝜆/𝑏+3𝑏2/𝑘2)𝑡, 𝑘0, and 𝑏0.

For application of homotopy analysis method, in view of (4.9) and the initial condition given in (4.10), it in convenient to choose𝑏𝑢(𝑥,0)=𝜆2𝑘+𝑘tanh(𝑘𝑥),𝑣(𝑥,0)=2𝑘1+𝑏+𝑏tanh(𝑘𝑥),(4.12) as the initial approximate of (4.10). We choose the linear operators𝐿1𝜙1(𝑥,𝑡,𝑞)=𝐷𝛼𝑡𝜙1,𝐿(𝑥,𝑡,𝑞)2𝜙2(𝑥,𝑡,𝑞)=𝐷𝛽𝑡𝜙2,(𝑥,𝑡,𝑞)(4.13) with the property 𝐿(𝑐)=0 where 𝑐 is constant of integration. Furthermore, we define a system of nonlinear operators as 𝑁1𝜙𝑖(𝑥,𝑡,𝑞)=𝐷𝛼𝑡𝜙11(𝑥,𝑡,𝑞)2𝜕3𝜙1(𝑥,𝑡,𝑞)𝜕𝑥3+3𝜙1(𝑥,𝑡,𝑞)2𝜕𝜙1(𝑥,𝑡,𝑞),3𝜕𝑥2𝜕2𝜙2(𝑥,𝑡,𝑞)𝜕𝑥2𝜕𝜙31(𝑥,𝑡,𝑞)𝜙2(𝑥,𝑡,𝑞)𝜕𝑥+3𝜆𝜕𝜙1(𝑥,𝑡,𝑞),𝑁𝜕𝑥2𝜙𝑖(𝑥,𝑡,𝑞)=𝐷𝛽𝑡𝜙2+𝜕(𝑥,𝑡,𝑞)3𝜙2(𝑥,𝑡,𝑞)𝜕𝑥3+3𝜙2(𝑥,𝑡,𝑞)𝜕𝜙2(𝑥,𝑡,𝑞),𝜕𝑥+3𝜕𝜙1(𝑥,𝑡,𝑞)𝜕𝑥𝜕𝜙2(𝑥,𝑡,𝑞)𝜕𝑥3𝜙1(𝑥,𝑡,𝑞)2𝜕𝜙2(𝑥,𝑡,𝑞)𝜕𝑥3𝜆𝜕𝜙2(𝑥,𝑡,𝑞).𝜕𝑥(4.14) We construct the zeroth-order and the 𝑚th-order deformation equations where 𝑅1,𝑚𝑢𝑚1=𝐷𝛼𝑡𝑢𝑚112𝑢𝑚1𝑥𝑥𝑥+3𝑚1𝑖=0𝑢𝑖𝑚1𝑖𝑘=0𝑢𝑘𝑣𝑚1𝑖𝑘𝑥,32𝑣𝑚1𝑥𝑥3𝑚1𝑘=0𝑢𝑘𝑣𝑚𝑘1𝑥𝑢+3𝜆𝑚1𝑥,𝑅2,𝑚𝑣𝑚1=𝐷𝛽𝑡𝑣𝑚1+𝑣𝑚1𝑥𝑥𝑥+3𝑚1𝑘=0𝑣𝑘𝑣𝑚1𝑘𝑥+3𝑚1𝑘=0𝑢𝑘𝑥𝑣𝑚𝑘1𝑥,3𝑚1𝑖=0𝑢𝑖𝑚1𝑖𝑘=0𝑢𝑘𝑣𝑚1𝑖𝑘𝑥𝑣3𝜆𝑚1𝑥.(4.15) We start with an initial approximation 𝑢(𝑥,0)=(𝑏/2𝑘)+𝑘tanh(𝑘𝑥), 𝑣(𝑥,0)=(𝜆/2)(1+(𝑘/𝑏))+𝑏tanh(𝑘𝑥), with 𝑘=0.1, 𝑏=1,𝑘=1/3, thus we can obtain directly the other components as follows: 𝑢1=11771620𝑡𝑎1+tanh((1/3)𝑥)2,𝑢Γ(𝑎+1)2=1437400656100Γ(𝑎+𝑏+1)𝑎2Γ(𝑏+1)Γ(𝑎)2Γ(2𝑎+1)145800,Γ(𝑎+𝑏+1)𝑎2Γ(𝑏+1)Γ(𝑎)21Γ(2a+1)tanh3𝑥656100Γ(𝑎+𝑏+1)𝑎2,Γ(𝑏+1)Γ(𝑎)2Γ(2𝑎+1)145800Γ(𝑎+𝑏+1)𝑎2Γ(𝑏+1)Γ(𝑎)21Γ(2𝑎+1),tanh3𝑥3275102𝑡(𝑏+𝑎)𝑎2Γ(𝑏+1)Γ(𝑎)2Γ(2𝑎+1)+13100402𝑡(𝑏+𝑎),𝑎2Γ(𝑏+1)Γ(𝑎)21Γ(2𝑎+1)tanh3𝑥29825302𝑡(𝑏+𝑎)𝑎2Γ(𝑏+1)Γ(𝑎)2,1Γ(2𝑎+1)tanh3𝑥4317790Γ(𝑎+𝑏+1)𝑎Γ(𝑎)Γ(2𝑎+1)𝑡𝑎𝑣Γ(𝑏+1),+1=1213540𝑡𝑏1+tanh((1/3)𝑥)2,𝑣Γ(𝑏+1)2=1145800145800𝑏2Γ(𝑎+1)Γ(𝑏)21Γ(2𝑏+1)tanh3𝑥145800𝑏2,Γ(𝑎+1)Γ(𝑏)21Γ(2𝑏+1)tanh3𝑥9720𝑏2Γ(𝑎+1)Γ(𝑏)2Γ(2𝑏+1),327510𝑏Γ(𝑏)Γ(2𝑏+1)𝑡𝑏Γ(𝑎+1)3275102𝑏Γ(𝑏)Γ(2𝑏+1)𝑡𝑏Γ(𝑎+1),+3177902𝑏Γ(𝑏)Γ(2𝑏+1)𝑡(𝑏+𝑎)6355802𝑏Γ(𝑏)Γ(2𝑏+1)𝑡(𝑏+𝑎),1tanh3𝑥22824802𝑏Γ(𝑏)Γ(2𝑏+1)𝑡(𝑏+𝑎)1tanh3𝑥3+14124021𝑏,tanh3𝑥4+(4.16)

The absolute error of the 6th-order HAM and exact solution with =1 as shown in Figure 5. Also the absolute errors |𝑢(𝑡)𝜙6(𝑡)| have been calculated for in Table 2. Figure 6 shows the numerical solutions of the nonlinear coupled Burgers equations with parameters derivative with =1, 𝛼=𝛽=1. Figure 7 shows the explicit numerical solutions with =1, 𝛼=1/2, and 𝛽=2/3 at 𝑡=0.002.

tab2
Table 2: The comparison of the results of the HAM (=1) and exact solution for the 𝑢(𝑥,𝑡), 𝛼=𝛽=1.
fig5
Figure 5: The comparison of the 6th-order HAM and exact solution with=1, 𝛼=𝛽=1, 𝜆=0.1, 𝑏=1, and 𝑘=1/3.
fig6
Figure 6: Explicit numerical solutions with=1, 𝛼=𝛽=1, 𝜆=0.1, 𝑏=1, and 𝑘=1/3.
fig7
Figure 7: Explicit numerical solutions with =1, 𝛼=1/2, 𝛽=2/3, 𝜆=0.1, 𝑏=1, and 𝑘=1/3.

As suggested by Liao [7], the appropriate region for is a horizontal line segment. We can investigate the influence of on the convergence of the solution series gevin by the HAM, by plotting its curve versus h, as shown in Figure 8.

545607.fig.008
Figure 8: The -curves obtained from the 5-order HAM approximate solution.

Remark 4.2. This example has been solved using homotopy perturbation method [21]. The graphs drawn and tables by =1 are in excellent agreement with those graphs drawn with HPM.

5. Conclusion

In this paper, based on the symbolic computation MATLAB, the HAM is directly extended to derive explicit and numerical solutions of the nonlinear coupled equations with parameters derivative. HAM provides us with a convenient way to control the convergence of approximation series by adapting , which is a fundamental qualitative difference in analysis between HAM and other methods. The obtained results demonstrate the reliability of the HAM and its wider applicability to fractional differential equation. It, therefore, provides more realistic series solutions that generally converge very rapidly in real physical problems. MATLAB has been used for computations in this paper.

References

  1. M. Caputo, “Linear models of dissipation whose Q is almost frequency independent, part II,” Geophysical Journal of the Royal Astronomical Society, vol. 13, no. 5, pp. 529–539, 1967. View at Publisher · View at Google Scholar
  2. K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974.
  3. S. Momani and N. Shawagfeh, “Decomposition method for solving fractional Riccati differential equations,” Applied Mathematics and Computation, vol. 182, no. 2, pp. 1083–1092, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  4. Z. Odibat and S. Momani, “Modified homotopy perturbation method: application to quadratic Riccati differential equation of fractional order,” Chaos, Solitons and Fractals, vol. 36, no. 1, pp. 167–174, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  5. Z. Odibat and S. Momani, “Application of variation iteration method to nonlinear differential equations of fractional order,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 1, no. 7, pp. 15–27, 2006.
  6. S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems, Ph.D. thesis, Shanghai Jiao Tong University, 1992.
  7. S. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2003.
  8. S. Liao, “On the homotopy anaylsis method for nonlinear problems,” Applied Mathematics and Computation, vol. 147, pp. 499–513, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. S. Liao, “Comparison between the homotopy analysis method and homotopy perturbation method,” Applied Mathematics and Computation, vol. 169, no. 2, pp. 1186–1194, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  10. S. Liao, “Homotopy analysis method: a new analytical technique for nonlinear problems,” Journal of Communications in Nonlinear Science and Numerical Simulation, vol. 2, no. 2, pp. 95–100, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  11. T. Hayat, M. Khan, and M. Ayub, “On non-linear flows with slip boundary condition,” Zeitschrift für Angewandte Mathematik und Physik, vol. 56, no. 6, pp. 1012–1029, 2005. View at Publisher · View at Google Scholar · View at MathSciNet
  12. S. Abbasbandy and F. S. Zakaria, “Soliton solutions for the 5th-order KdV equation with the homotopy analysis method,” Nonlinear Dynamics, vol. 51, no. 1-2, pp. 83–87, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  13. S. P. Zhu, “An exact and explicit solution for the valuation of American put options,” Quantitative Finance, vol. 6, no. 3, pp. 229–242, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  14. A. K. Alomari, M. S. M. Noorani, R. Nazar, and C. P. Li, “Homotopy analysis method for solving fractional Lorenz system,” Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 7, pp. 1864–1872, 2010. View at Publisher · View at Google Scholar
  15. A. K. Alomari, M. S. M. Noorani, and R. Nazar, “Solution of delay differential equation by means of homotopy analysis method,” Acta Applicandae Mathematicae, vol. 108, no. 2, pp. 395–412, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  16. B. J. West, M. Bologna, and P. Grigolini, Physics of Fractal Operators, Institute for Nonlinear Science, Springer, New York, NY, USA, 2003.
  17. K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, A Wiley-Interscience Publication, John Wiley & Sons, New York, NY, USA, 1993.
  18. S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives, Gordon and Breach Science Publishers, Yverdon, Switzerland, 1993.
  19. I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
  20. J. H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  21. Y. Chen and H. An, “Homotopy perturbation method for a type of nonlinear coupled equations with parameters derivative,” Applied Mathematics and Computation, vol. 204, no. 2, pp. 764–772, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  22. 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet