An improved homotopy analysis method with accelerated convergence for nonlinear problems

In this paper we propose an accelerated convergence method, which is combined with the homotopy analysis method (HAM), to solve nonlinear problems. The HAM is applied to obtain approximate expressions. According to the numbers of terms in the approximations, some ratio-control parameters are introduced in the solution expressions. By solving simultaneous algebraic equations, all artificial parameters can be optimally identified, including the so-called convergence-control parameter ħ. Twoexamples are given to illustrate the validity of the new method. Comparison with L-P perturbation method and Runge-Kutta method reveals that the improved HAM is better than the standard HAM and applies especially to the problems with complicated nonlinear terms.


Introduction
There are various types of nonlinear problems in the engineering field, most of which are strongly nonlinear and do not have an analytical solution.To date, many powerful methods have been developed to solve nonlinear differential equations, such as the Krylov-Bogoliubov-Mitropolsky (KBM) [1][2][3], the method of harmonic balance [4][5][6], Adomian's decomposition method (ADM) [7,8] and the small parameter method [9][10][11].However, these methods cannot provide us with a simple way to adjust and control the convergence region and rate of approximate solution series.The homotopy analysis method (HAM) [12][13][14][15] was proposed in 1992 by Liao [16].Based on the homotopy of topology, the advantage of the HAM is that it is independent of a small or large physical parameter and provides a way to control and ensure the convergence of approximate series.It has been proven that the HAM logically contains non-perturbations such as the Lyapunov artificial small parameter method, the ADM, the -expansion method, the variational iteration method and so on [17,18].Therefore, the HAM is more general.
However, to obtain a rather accurate homotopy approximation, the high-order deformation equation must be solved, which is usually complicated to compute.And in the standard HAM, the so called ℏ-curves must be plotted to find the optimal values of ℏ with which the approximations converge fastest.In this paper, we introduce an accelerated convergence method, by which the values of ℏ can be optimally recognized and the improved approximations converge faster than the standard homotopy approximations.

Basics of the HAM
One can transform a nonlinear problem into an infinite number of linear sub-problems and obtain series solutions via the homotopy analysis method, whether the nonlinear problem contains small parameters or not.Consider the following general nonlinear system: where is a linear operator, is a nonlinear operator, ( ) is a known function, ( ) is an unknown function to be determined and denotes temporal or spatial independent variables.Let ( ) denote an initial guess of the exact solution for ( ) , ℏ ≠ 0 denote an auxiliary convergence-control parameter, and ( ) ≠ 0 denote an auxiliary function.Then, the homotopy function is constructed as: where ∈ [0,1] is an embedding parameter, ℒ is an auxiliary linear operator with the property ℒ(0) = 0, and ℵ is a nonlinear operator defined as: Then, enforcing: the zero-order deformation equation would be constructed as: When = 0 and = 1, we have respectively: when = 0, and: when = 1.
It is easy to show that from Eq. ( 6) and Eq.(7): Thus, as the embedding parameter increases from 0 to 1, Φ( ; ) varies continuously from the initial guess ( ) to the exact solution ( ) of the original Eq. ( 1).Differentiating Eq. ( 5) times with respect to the embedding parameter and then setting = 0 and finally dividing them by !, one can obtain the so-called th-order deformation equation: where: for simplicity, and: Because Φ( ; ) depends on the embedding parameter ∈ [0,1], it can be expanded into a power series of by using the Tylor theorem as: where: ( ) can be found by solving Eq. (10), which is a linear equation.If the initial guess ( ), auxiliary linear operator ℒ , convergence-control parameter ℏ and auxiliary function ( ) are properly chosen, the series Eq. ( 14) is convergent to Φ( ; ) at = 1 , and one can have, according to Eq. ( 9) and Eq. ( 14), the homotopy-series solution: and the th-order homotopy-approximation:

The idea of the accelerated convergence method
Only the series solutions can be obtained for most nonlinear problems.It makes sense to ensure and accelerate the convergence of the series.Once the base functions that express the solution are chosen, the proportion of the terms of the solution is pivotal to the convergence rate.We try to introduce some ratio-control parameters in the approximate solution expressions, combined with the original governing equations and the given conditions, to control the proportion between the terms.If there are + 1 non-constant terms in the approximate expressions, ratio-control parameters should be introduced as multipliers of the first terms.Assume that the th-order homotopy approximation of Eq. ( 1) can be expressed as: Then we introduce ratio-control parameters in Eq. ( 18) as follows: For initial value or boundary value problems, the values of ( ), ′( ), …, and ( ) ( ) can be obtained by continuous derivation of Eq. ( 1) if the given conditions is at = , then an equation can be written as follows: , ( ) + , ( ) + ⋯ + , ( ) + , ( ) + = ( ), , ( ) + , ( ) + ⋯ + , ( ) + , ( ) = ( ), ⋯ ( ) , ( ) + ( ) , ( ) + ⋯ + ( ) , ( ) + ( ) , ( ) = ( ) ( ). ( 2233.AN IMPROVED HOMOTOPY ANALYSIS METHOD WITH ACCELERATED CONVERGENCE FOR NONLINEAR PROBLEMS.HONG-WEI LI, JUN WANG, LI-XIN LU, ZHI-WEI WANG Note that the numbers of equations should be equal to the numbers of unknowns in Eq. ( 20) including the convergence-control parameter ℏ .Sometimes Eq. ( 20) consists of nonlinear algebraic equations, thus there may be multiple solutions.Therefore, to choose the proper values of parameters, we can use the squared residual: , , , , ⋯ , , , ℏ = ; , , , , ⋯ , , , ℏ , where, denotes the number of the solutions of Eq. ( 20) and: Then, one can obtain the optimal values of parameters for which ( , , , , ⋯ , , , ℏ ) is the minimum.
Two examples are evaluated in the following to illustrate the specific method.

Duffing equation
First, consider the duffing equation: where the dot denotes the derivative with respect to .The base functions and the initial approximation of the solution are chosen respectively as: where is the frequency to be determined and: The auxiliary linear operator can be chosen as: and the non-zero auxiliary function ( ) as 1 for simplicity.By using the homotopy analysis method, one can easily obtain the second-order homotopy approximations of Eq. ( 23) as: where: There are three non-constant terms in the second-order approximation ( ) .Thus, we introduce two so-called ratio-control parameters and into Eq.(28) as: One can easily obtain the expressions of ( ), ( ) ( ) and ( ) ( ) by differentiating Eq. ( 29) continuously with respect to , respectively and the exact values of (0), (0), ( ) (0) and ( ) (0) from Eq. ( 23).Then, simultaneous equations can be written as follows: In order to choose the proper values of unknowns, for the sake of computational efficiency, we calculate the discrete squared residual: where denotes the order of approximations, is an integer and: Obviously, Eq. ( 31) is a good approximation of Eq. ( 21) for large enough .In this example, we use = 100.The values of the unknows , , and ℏ can be obtained by solving Eq. (30) and using Eq.(31).Then the improved frequency can be obtained from Eq. ( 27).The second order standard homotopy approximations can be easily obtained by plotting the so called ℏ-curves for different cases.Table 1 gives the values of parameter ℏ which minimizes the discrete squared residual in the standard HAM (SHAM).To illustrate the validity of the improved homotopy analysis method (IHAM) for this example, comparison of frequency is given in Table 2.

Simple pendulum attached to a rotating rigid frame[12]
The governing equation is: subject to the initial conditions: To apply the HAM to solve the above problem, the following approximation is used: as a result, we obtain the following approximate differential equation: Obviously, the more terms of the Maclaurin series considered, the more efficient the solution series obtained but the more complicated the problem will be.To illustrate the validity and simplicity of the accelerated convergence method, we only use the first two terms of the Maclaurin series to approximate the original Eq. ( 34) as a duffing-like Eq. ( 37).The base functions and the initial approximation of the solution are chosen as Eq. ( 24) and Eq. ( 25), respectively.The second order approximate expressions of Eq. ( 37) can be easily obtained by the HAM.There are three non-constant terms in the second-order approximate expression ( ) .Thus, the improved expression ( ) should contain two ratio-control parameters like Eq. ( 29).
The second order homotopy approximate frequency of Eq. ( 37) reads: Note that, to establish the simultaneous equations like Eq. (30) in this example, one should obtain the exact values of (0) , (0) , ( ) (0) and ( ) (0) by differentiating the original governing Eq. (34) continuously with respect to and using the initial conditions.
We also obtain the values of convergence-control parameter ℏ which minimizes the discrete squared residual in the SHAM by plotting the ℏ-curves and compare the minimum of discrete squared residual in the SHAM and the IHAM, as shown in Table 3.To illustrate the validity of the accelerated convergence method for this example, comparison of frequency is given in Table 4.

Comparison with L-P perturbation method
Lindstedt-Poincaré perturbation method [19,20] (L-P perturbation method) is a powerful analytical technique to weak nonlinear problems.Eq. (34) have a periodic solution.In the idea of L-P perturbation, one should express the solution ( ) and the frequency squared in power series with respect to : where, | | ≪ 1.By substituting Eq. (39) and Eq.(40) into Eq.( 23), one obtain: By equating coefficients of like powers of , one can obtain the set of recursive linear differential equations: and so on.By solving above equations and eliminating the secular terms, we can obtain three terms order L-P approximation as: Fig. 3 and Fig. 4 give the solution comparisons of L-P perturbation method with SHAM and IHAM in different weak nonlinear cases.HONG-WEI LI, JUN WANG, LI-XIN LU, ZHI-WEI WANG

Results and discussion
The standard homotopy analysis method (SHAM) gives us a way to control the ratio and region of approximations for nonlinear problems by choosing the values of ℏ, thus in general it can give better solutions than other analytical methods.Without plotting ℏ-curves, the same and even more accurate results can be obtained by improved homotopy analysis method (IHAM) in various cases, as shown in Table 1.Table 2 shows that the accelerated convergence method can give more accurate values of frequencies than ℏ-curves.We also applied the IHAM to a problem with complicated nonlinear terms and obtained better results than ones by plotting the ℏ-curves in different cases, as shown in Table 4, Table 5, Fig. 1 and Fig. 2.This successful application means that the problems with complicated nonlinear terms can be reduced to simple ones.One can only analyze simplified duffing-like equations by the HAM and obtain greatly accurate results of the original problems by the accelerated convergence method.For weak nonlinear problems, L-P perturbation can give accurate enough approximate solution.Comparison with L-P perturbation reveals that IHAM is of high accuracy for various cases, as shown in Fig. 3 and Fig. 4.

Conclusions
In this paper, we proposed an accelerated convergence method and improved the HAM.Faster convergent results can be obtained by our method than plotting the ℏ-curves.The idea of the method is to control the proportion between the homotopy approximation terms by introducing some ratio-control parameters in the approximate solution expressions.All artificial parameters including the so called convergence-control parameter ℏ can be optimally identified by solving algebraic equations.The current work has elucidated the validity and justifications of this combination.

Fig. 1 .Fig. 2 .
Fig. 1.Comparison of the exact solution with the analytic approximations in the case of = 1, = 0.5, = 0.5π for example 2. Symbols: numerical solution; solid line: improved homotopy approximation by accelerated convergence; dashed line: standard homotopy approximation by plotting ℏ-curves IMPROVED HOMOTOPY ANALYSIS METHOD WITH ACCELERATED CONVERGENCE FOR NONLINEAR PROBLEMS.HONG-WEI LI, JUN WANG, LI-XIN LU, ZHI-WEI WANG

Table 2 .
Comparison of corresponding to various parameters of system for example 1