Research Article  Open Access
Solving Famous Nonlinear Coupled Equations with Parameters Derivative by Homotopy Analysis Method
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 fluiddynamic 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 [6โ10]. This method has been successfully applied to solve many types of nonlinear problems in science and engineering, such as the viscous flows of nonNewtonian fluids [11], the KdVtype 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 nonlinearcoupled equations with parameters derivative in this form: where and () 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 is said to be in the space , if there exists a real number , such that where , and it is said to be in the space l if and only if , . Clearly if .
Definition 2.2. The RiemannLiouville fractional integral operator () of order , of a function , , is defined as
is the wellknown Gamma function. Some of the properties of the operator , which we will need here, are as follows.
For , , and
Definition 2.3. For the concept of fractional derivative, there exist many mathematical definitions [1, 16โ19]. In this paper, the two most commonly used definitions: the Caputo derivative and its reverse operator RiemannLiouville 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 Here, we also need two basic properties about them:
Definition 2.4. The MittagLeffler function with is defined by the following series representation, valid in the whole complex plane:
3. Basic Idea of HAM
To describe the basic ideas of the HAM, we consider the operator form of (1.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 socalled zeroorder deformation equations: where is the embedding parameter, is a nonzero auxiliary parameter, is an auxiliary function, is an auxiliary linear operator, , are initial guesses of , and , are two unknown functions, respectively. It is important that one has great freedom to choose auxiliary things in HAM. Obviously, when and , the following holds: respectively. Thus, as increases from 0 to 1, the solution , varies from the initial guess , to the solution , . Expanding , in Taylor series with respect to , we have where 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 , then we have which must be one of solutions of original nonlinear equation, as proved by Liao [8]. As and , (3.2) and (3.3) become 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 zeroorder deformation equation (3.2). Define the vector Differentiating equations (3.2) and (3.3) times with respect to the embedding parameter and then setting and finally dividing them by , we have the socalled thorder deformation equation: where Applying the RiemannLiouville integral operator on both side of (3.10), we have It should be emphasized that , for 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 socalled 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 , 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 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: where , and the fractional operators and are defined as in (2.3). Assuming the initial value as The exact solutions of (4.1) for the special case: are For application of homotopy analysis method, in view of (4.1) and the initial condition given in (4.2), it is convenient to choose as the initial approximate of (4.1). We choose the linear operators with the property where is constant of integration. Furthermore, we define a system of nonlinear operators as We construct the zerothorder and the thorder deformation equations where We start with an initial approximation , , thus we can obtain directly the other components as The absolute error of the 6thorder HAM and exact solution with as shown in Figure 1. Also the absolute errors have been calculated in Table 1. Figure 2 shows the numerical solutions of the nonlinear coupled Burgers equations with parameters derivative with , . Figure 3 shows the explicit numerical solutions with , , and at .

(a)
(b)
(a)
(b)
(a)
(b)
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.
Remark 4.1. This example has been solved using homotopy perturbation method [21]. The graphs drawn and tables by 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: with the initial conditions, As we know, when (4.9) has the kinktype soliton solutions constructed by Fan [22], where , , and .
For application of homotopy analysis method, in view of (4.9) and the initial condition given in (4.10), it in convenient to choose as the initial approximate of (4.10). We choose the linear operators with the property where is constant of integration. Furthermore, we define a system of nonlinear operators as We construct the zerothorder and the thorder deformation equations where We start with an initial approximation , , with , ,, thus we can obtain directly the other components as follows:
The absolute error of the 6thorder HAM and exact solution with as shown in Figure 5. Also the absolute errors have been calculated for in Table 2. Figure 6 shows the numerical solutions of the nonlinear coupled Burgers equations with parameters derivative with , . Figure 7 shows the explicit numerical solutions with , , and at .

(a)
(b)
(a)
(b)
(a)
(b)
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.
Remark 4.2. This example has been solved using homotopy perturbation method [21]. The graphs drawn and tables by 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
 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 Site  Google Scholar
 K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974.
 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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Google Scholar
 S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problems, Ph.D. thesis, Shanghai Jiao Tong University, 1992.
 S. Liao, Beyond Perturbation: Introduction to the Homotopy Analysis Method, CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2003.
 S. Liao, โOn the homotopy anaylsis method for nonlinear problems,โ Applied Mathematics and Computation, vol. 147, pp. 499โ513, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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 Site  Google Scholar  Zentralblatt MATH
 T. Hayat, M. Khan, and M. Ayub, โOn nonlinear flows with slip boundary condition,โ Zeitschrift für Angewandte Mathematik und Physik, vol. 56, no. 6, pp. 1012โ1029, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 S. Abbasbandy and F. S. Zakaria, โSoliton solutions for the 5thorder KdV equation with the homotopy analysis method,โ Nonlinear Dynamics, vol. 51, no. 12, pp. 83โ87, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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 Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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 Site  Google Scholar
 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 Site  Google Scholar  Zentralblatt MATH
 B. J. West, M. Bologna, and P. Grigolini, Physics of Fractal Operators, Institute for Nonlinear Science, Springer, New York, NY, USA, 2003.
 K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, A WileyInterscience Publication, John Wiley & Sons, New York, NY, USA, 1993.
 S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives, Gordon and Breach Science Publishers, Yverdon, Switzerland, 1993.
 I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
 J. H. He, โHomotopy perturbation technique,โ Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 34, pp. 257โ262, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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 Site  Google Scholar  Zentralblatt MATH
 E. G. Fan, โSoliton solutions for a generalized HirotaSatsuma coupled KdV equation and a coupled MKdV equation,โ Physics Letters. A, vol. 282, no. 12, pp. 18โ22, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
Copyright
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.