Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 982810, 7 pages

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

## Wavelet-Based Homotopy Analysis Method for Nonlinear Matrix System and Its Application in Burgers Equation

Hebei University of Engineering, Handan 056038, China

Received 17 March 2013; Accepted 14 June 2013

Academic Editor: Evangelos J. Sapountzakis

Copyright © 2013 Xing Ruyi. 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

To generalize the homotopy analysis method (HAM) to multidegree-of-freedom nonlinear system, the adaptive precise integration method (APIM) is introduced into the HAM, with which the almost exact value of the exponential matrix can be obtained. Combining the interval interpolation wavelet collocation method, HAM-based APIM can be employed to solve the nonlinear PDEs. As an example, Burgers equation is spatially discretized by the interval quasi-Shannon wavelet collocation method and solved by the proposed method to illustrate the effectiveness and great potential of the homotopy analysis method in nonlinear problems.

#### 1. Introduction

The homotopy analysis method (HAM) is developed in 1992 by Liao in [1–5]. This method has been successfully applied in various nonlinear problems in science and engineering in [6–11] and references therein. Based on homotopy of topology, the validity of the HAM is independent of whether or not there exist small parameters in the considered equation. Therefore, the HAM can overcome the foregoing restrictions and limitations of perturbation techniques [12, 13]. In fact, He [14, 15] proposed another asymptotic numerical method for solving nonlinear problems named homotopy perturbation method (HPM). The better improvement is adding an auxiliary parameter into the homotopy equation, which is helpful to eliminate the secular term in the perturbation solution. This can improve the rate of convergence greatly. Unlike analytical perturbation methods, HPM does not depend on small parameter which is difficult to find. Most researchers think of HPM as a special case of HAM, which is much simpler than classical HAM and has been widely used in solving various nonlinear problems [16]. However, the HAM applies mainly in single-degree-of-freedom systems; as for multidegree-of-freedom nonlinear differential equations, due to the correlation of the equations, the solving process is very tedious, and strong calculating technique is needed.

In an attempt to overcome the previous difficulties, the precise integration method (PIM) proposed by Zhong [17] in 1994 is employed, because it is effective in computation precision and stability for the solutions of linear differential equations [18, 19]. However, like other numerical methods, PIM has its own limitation. When it is applied to the problems of time-variant or nonlinear time integration, PIM usually expands the nonlinear term in Taylor series firstly, which is then solved by using an iteration method. In order to obtain high precision, a very small time step should be taken for PIM, which results in huge computational complexity and greatly restricts the application of PIM in engineering.

Burgers’ model of turbulence is a crucial fluid dynamic model, and the study of this model and theory of shock waves has been considered by many authors both for conceptual understanding of a class of physical flows and for testing various numerical methods. In this work, the Burgers equation is firstly transformed into a set of nonlinear matrix differential equations using wavelet collocation method [20–22]. Since the definition domain of wavelet transformation is an infinite interval, which can deduce the boundary effect, this can be eliminated by introducing an interval wavelet [23] instead of the common wavelet function.

Taking full advantage of the HAM and PIM, a novel approximate method for multidegree-of-freedom nonlinear matrix differential equation based on interval wavelet is presented in this paper.

#### 2. The HAM Based on the PIM

##### 2.1. Basic Idea of the Homotopy Analysis Method

Taking as an example, we consider the following nonlinear algebraic equation: First of all, we construct such a homotopy where is an initial guess of and is called embedding parameter. Obviously, at and , one has respectively. Thus, as increases from 0 to 1, varies continuously from to . Such kind of continuous variation is called deformation in topology [20]. Now, enforcing , that is, we have now a family of algebraic equations. Obviously, the solution of the previous family of algebraic equations is dependent upon the embedding parameter . So, the family of equations can be rewritten as

As the embedding parameter increases from 0 to 1, varies from the initial guess to the solution of . We call the family of equations like (5) the zeroth-order deformation equation

Because is now a function of the embedding parameter , we can expand it into Maclaurin series where is employed, and Here, the series (6) is called homotopy series, and is called the th-order homotopy derivative of . If the convergence radius of the previous Maclaurin series is no less than 1, using the relationship , one has

According to the fundamental theorem of calculus about Taylor series, the coefficient of the homotopy series (6) is unique.

Substituting (6) into (5) and rearranging based on powers of -terms, one has Whose solutions are

In this way, one obtains one by one in the order . Then one has 1st-order homotopy-series approximation 2nd-order homotopy-series approximation

##### 2.2. The Precise Integration Method for Nonlinear Matrix Differential Equations

Taking as an example, we consider the following matrix differential equations: where is an -dimensional unknown vector, and are matrix, and is a time variable. The general solution of (14) is To solve the integral part of (15), the linear interpolation approximation is employed where , , and is the th element of vector . Based on the previous derivations, the recurrence formula of (14) can be obtained as

Supposing the time step is constant, the matrix exponential function can be calculated accurately by using PIM. The calculation of the exponential function is based on the identity where is an integer. It is suggested to select , such as , .

Because is a small time interval, is an extremely small time interval. Hence for the interval, the truncated Taylor expansion is applied with high precision where

Because is very small, the first five-term series expansion should be enough.

For computing the matrix , (17) could be factored as Such factorization should be iterated times. Noted that it is equivalent to do the following instruction: after the execution of the above instruction, .

Because the nonlinear term is expanded in Taylor series and the PIM uses an iterative method, the time step must be very small.

#### 3. Coupling Technique of the HAM and the PIM for Burgers Equation Based on Interval Wavelet

##### 3.1. Spatial Discretization of the Burgers Equation

Consider that the one-dimensional Burgers equation has the form [21–23] subject to the initial and boundary conditions where represents the time, while denotes the Reynolds number.

According to the basic idea of wavelet collocation method, the solution could be approximately represented as where are uniform discrete points and is the scaling function. Here we take the quasi-Shannon wavelet [24, 25] as the scaling function where is the window size parameter.

In order to take the quasi-Shannon scaling function as the basis function, is discretized in its definition domain And thereby the basis function is obtained where .

It is easy to see that has the interpolation property

Using its interpolation property, the first derivative and second derivative of can be obtained as follows:

Substituting (4) into (1), a series of nonlinear ODEs can be obtained

For simplicity, we can write (10) as matrix differential equation where

##### 3.2. Construction of Quasi-Shannon Interval Wavelet

The Lagrange polynomials , are introduced The following two weight coefficients are defined Then the interpolation basic function can be expressed as where is the number of he external points and is the support domain of the wavelet function.

It is easy to know that the quasi-Shannon interval wavelet is a linear combination of the quasi-Shannon scaling function . Therefore, the quasi-Shannon interval wavelet function possesses all the properties of the quasi-Shannon scaling function.

##### 3.3. The PIM for Matrix Differential Equations

We choose the time step as , assuming that the nonlinear term in (31) is linear in the relatively short time segment , the nonlinear term of (31) can be expressed in Taylor series, and the first two terms are taken. Therefore, (31) could be represented as where

Based on the previous derivation, the iterative solution of (38) can be obtained where the matrix exponential function can be calculated accurately by using PIM.

##### 3.4. The Homotopy Analysis Method Based on the Precise Integration Method

There are various ways to construct a homotopy function. For (11), a linear homotopy can be constructed as where is the initial approximation of and is the embedding parameter. As increases from 0 to 1, the solution varies from the initial guess to the solution , and it can be expanded in Taylor series about the embedding parameter where

that the approximate solutions are generally not valid for large . In order to improve the computation accuracy, we divide the time interval into a number of sections uniformly by setting the time step as , and use the iterative method. In such a way, a series of time nodes can be expressed as

Since the only requirement for is that from to , should change continuously from the initial guess to . Therefore, it is feasible to improve the computation accuracy by assuming that Then can be expanded in the form

Substituting (46) into (41) and rearranging based on powers of -terms, the following are obtained: where , and .

Substituting (47) into (46) and assuming , one has where , , . The matrix exponential function can be calculated by using PIM.

For simplicity, only the first two terms are taken as the approximate solution. Accordingly, the approximate solution is obtained

#### 4. Numerical Results and Discussion

The Burgers equation has analytical solution as follows:where is the th order modified Bessel function of the first kind. Figure 1 shows the analytical solutions with Re = 100 at , 0.4, 0.6, respectively. It is easy to see that the analytical solution evolves into a shock wave near , and the gradient becomes larger and larger with the increasing of Re.

In the analysis of the numerical result, we choose Reynolds number , the scaling parameter , the artificial parameter , and the number of the external points .

The comparison of the classical PIM and the HAM is shown in Figures 2 and 3 and Tables 1 and 2.

We can see clearly that when the time step is relatively large ( or ), the classical PIM fails actually it fails when is 0.04. However, the HAM has enough precision even when the time step increases to 0.1 and is more than 0.6.

As the time step becomes smaller, both the classical PIM and the HAM have higher precision; however, we can see clearly in Figure 3 that the HAM is still better than the classical PIM. Furthermore, from the comparison of (49) and (40), it can be seen that the computation cost of the proposed method is less than the classical PIM.

#### 5. Conclusions

In this paper, it is shown that the proposed method can solve the Burgers equation successfully. The numerical comparisons presented demonstrate that the solutions obtained by the classical PIM are not valid for large time step; the method proposed in this paper can obtain higher accuracy than the classical PIM even if the time step is small. Although the method proposed aims at Burgers equation, it is a promising tool for other nonlinear problems.

#### Acknowledgment

This work is supported by the National Key Technologies R & D Program of China under Grant no. 2012BAD35B02 and the National Natural Science Foundation of China under Grant no. 41171337.

#### References

- S. J. Liao,
*On the proposed homotopy analysis technique for nonlinear problems and its applications [Ph.D. dissertation]*, Shanghai Jiao Tong University, 1992. - S. J. Liao, “An approximate solution technique not depending on small parameters: a special example,”
*International Journal of Non-Linear Mechanics*, vol. 30, no. 3, pp. 371–380, 1995. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S. J. Liao, “A kind of approximate solution technique which does not depend upon small parameters—II. An application in fluid mechanics,”
*International Journal of Non-Linear Mechanics*, vol. 32, no. 5, pp. 815–822, 1997. View at Google Scholar · View at Scopus - S. J. Liao, “Homotopy analysis method: a kind of nonlinear analytical technique not depending on small parameters,”
*Shanghai J. Mech*, vol. 18, no. 3, pp. 196–200, 1997. View at Google Scholar - S. J. Liao,
*Beyond Perturbation: Introduction to the HomoTopy Analysis Method*, Champan and Hall/CRC Press, Boca Raton, Fla, USA, 2003. - M. Sajid and T. Hayat, “The application of homotopy analysis method to thin film flows of a third order fluid,”
*Chaos, Solitons & Fractals*, vol. 38, no. 2, pp. 506–515, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - M. Sajid, T. Hayat, and S. Asghar, “Comparison between the HAM and HPM solutions of thin film flows of non-Newtonian fluids on a moving belt,”
*Nonlinear Dynamics*, vol. 50, no. 1-2, pp. 27–35, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - S. J. Liao, “On the homotopy analysis method for nonlinear problems,”
*Applied Mathematics and Computation*, vol. 147, no. 2, pp. 499–513, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - S. Abbasbandy, Y. Tan, and S. J. Liao, “Newton-homotopy analysis method for nonlinear equations,”
*Applied Mathematics and Computation*, vol. 188, no. 2, pp. 1794–1800, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - T. Hayat, M. Khan, and M. Ayub, “On the explicit analytic solutions of an Oldroyd 6-constant fluid,”
*International Journal of Engineering Science*, vol. 42, no. 2, pp. 123–135, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S. Abbasbandy, “The application of homotopy analysis method to nonlinear equations arising in heat transfer,”
*Physics Letters A*, vol. 360, no. 1, pp. 109–113, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - A. H. Nayfeh,
*Perturbation Methods*, John Wiley & sons, New York, NY, USA, 2000. - J. A. Murdock,
*Perturbation Methods*, Cambridge Texts in Applied Mathematics, Cambridge University Press, Cambridge, UK, 1991. - J.-H. He, “Asymptotic methods for solitary solutions and compactons,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 916793, 130 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - J.-H. He, “Addendum: new interpretation of homotopy perturbation method,”
*International Journal of Modern Physics B*, vol. 20, no. 18, pp. 2561–2568, 2006. View at Publisher · View at Google Scholar · View at Scopus - S.-L. Mei, “Construction of target controllable image segmentation model based on homotopy perturbation technology,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 131207, 8 pages, 2013. View at Publisher · View at Google Scholar - W. X. Zhong, “On precise integration method,”
*Journal of Computational and Applied Mathematics*, vol. 163, no. 1, pp. 59–78, 2004. View at Publisher · View at Google Scholar · View at Scopus - W. X. Zhong, “Precise time-integration method for structural dynamics,”
*Journal of Dalian University of Technology*, vol. 34, no. 2, pp. 131–136, 1994. View at Google Scholar · View at Zentralblatt MATH · View at Scopus - W. X. Zhong, “Combined method for the solution of asymmetric Riccati differential equations,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 191, no. 1-2, pp. 93–102, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S. Bertoluzza and G. Naldi, “A wavelet collocation method for the numerical solution of partial differential equations,”
*Applied and Computational Harmonic Analysis*, vol. 3, no. 1, pp. 1–9, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S.-L. Mei, Q.-S. Lu, and S.-W. Zhang, “Adaptive wavelet precise integration method for partial differential equations,”
*Chinese Journal of Computational Physics*, vol. 21, no. 6, pp. 523–530, 2004. View at Google Scholar · View at Scopus - S.-L. Mei, Q.-S. Lu, S.-W. Zhang, and L. Jin, “Adaptive interval wavelet precise integration method for partial differential equations,”
*Applied Mathematics and Mechanics*, vol. 26, no. 3, pp. 364–371, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S.-L. Mei, H.-L. Lv, and Q. Ma, “Construction of interval wavelet based on restricted variational principle and its application for solving differential equations,”
*Mathematical Problems in Engineering*, vol. 2008, Article ID 629253, 14 pages, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - D.-C. Wan and G.-W. Wei, “The study of quasi wavelets based numerical method applied to burgers' equations,”
*Applied Mathematics and Mechanics*, vol. 21, no. 10, pp. 1099–1110, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - G. W. Wei, “Quasi wavelets and quasi interpolating wavelets,”
*Chemical Physics Letters*, vol. 296, no. 3-4, pp. 253–258, 1998. View at Google Scholar · View at Scopus