Review Article  Open Access
Some Relatively New Techniques for Nonlinear Problems
Abstract
This paper outlines a detailed study of some relatively new techniques which are originated by He for solving diversified nonlinear problems of physical nature. In particular, we will focus on the variational iteration method (VIM) and its modifications, the homotopy perturbation method (HPM), the parameter expansion method, and expfunction method. These relatively new but very reliable techniques proved useful for solving a wide class of nonlinear problems and are capable to cope with the versatility of the physical problems. Several examples are given to reconfirm the efficiency of these algorithms. Some open problems are also suggested for future research work.
1. Introduction
With the rapid development of nonlinear sciences, many analytical and numerical techniques have been developed by various scientists. Most of the developed techniques have their limitations like limited convergence, divergent results, linearization, discretization, unrealistic assumptions and noncompatibility with the versatility of physical problems [1â€“100]. He [16â€“40] developed a number of efficient and reliable techniques for solving a wide class of nonlinear problems. These relatively new but very powerful methods proved to be fully synchronized with the complexities of the physical problems, see [1â€“7, 11â€“40, 49â€“80, 84â€“100] and the references therein. In the present study, we will focus our attention on Heâ€™s variational iteration (VIM), homotopy perturbation (HPM), modified variational iteration (MVIMS), parameter expansion, and expfunction methods. The variational iteration method (VIM) was suggested and proposed by He [17â€“24] in its preliminary form in 1999. The method has been used to solve nonlinear differential equations. In a subsequent work [23, 24] the VIM was developed into a full theory for solving diversified physical problems of versatile nature. It is to be highlighted that the variational iteration method (VIM) is also very effective for solving differentialdifference equations, see [49, 95] and the references therein. Moreover, He [17, 18, 27â€“38] introduced another wonderful technique, namely, homotopy perturbation (HPM) by merging the standard homotopy and perturbation. The HPM is independent of the drawbacks of the coupled techniques and absorbs all their positive features. It is to be noted that homotopy perturbation is a kind of perturbation method which can take full advantage of various perturbation methods, while using the homotopy technique to guarantee simple solution procedure. It is to worth mentioning that if the initial solution is suitably chosen, then only one or two iterations are enough to get the appropriate result, see [17, 18, 27â€“38]. The expfunction method was first proposed by He and Wu [39] in 2006. The method was originally suggested to search for solitary solutions and periodic solutions of nonlinear wave equations. It always leads to a generalized solution with free parameters which can be determined by using the initial/boundary conditions. The most interesting part is transformation between periodic and solitary solutions by using the socalled HeWu transformations. The method [17, 39, 40, 84â€“87, 97, 100] is always used as a tool to find exact solutions, but it can be utilized also for finding solutions approximately including the solutions for boundary value problems. It is to be highlighted that the present study would also outline the He parameterexpansion technique [18, 25â€“27]. The parameter expansion technique includes the modified LindstedtPioncare and book keeping parameter methods and previously called the parameterexpanding method. He [18] in his review article in 2006 also explained that the method does not require to construct a homotopy. These efficient techniques have been applied to a wide class of nonlinear problems, see [1â€“7, 11â€“40, 49â€“80, 84â€“100] and the references therein. With the passage of time some modifications in Heâ€™s variational iteration method (VIM) has been introduced by various authors. Abbasbandy [1, 2] made the coupling of Adomianâ€™s polynomials with the correction functional (VIMAP) of the VIM and applied this reliable version for solving Riccati differential and Klein Gordon equations. In a later work, Noor and MohyudDin [62, 64, 74] exploited this concept for solving various singular and nonsingular boundary and initial value problems. Recently, Ghorbani et al. [13, 14] introduced Heâ€™s polynomials (which are calculated from Heâ€™s homotopy perturbation method) by splitting the nonlinear term and also proved that Heâ€™s polynomials are fully compatible with Adomianâ€™s polynomials but are easier to calculate and are more user friendly. More recently, Noor and MohyudDin [60, 66â€“69, 72, 73] combined Heâ€™s polynomials and correction functional of the VIM and applied this reliable version (VIMHP) to a number of physical problems. It has been observed [60, 66â€“69, 72, 73] that the modification based on Heâ€™s polynomials (VIMHP) which was developed by Noor and MohyudDin is much easier to implement as compared to the one (VIMAP) where the socalled Adomianâ€™s polynomials along with their complexities are used. The basic motivation of the present study is the review of these very powerful and reliable techniques which have been originated by He for solving various nonlinear initial and boundary value problems of diversified physical nature. Several examples are given to reveal the efficiency and potential of these relatively new techniques. We have also pointed out that the techniques discussed in this paper can be extended for solving obstacle, free, moving, and contact problems, which arise in various fields of pure and engineering sciences. This is another aspect of future research work. The interested readers are advised to explore this avenue for innovative and novel applications of these techniques.
2. ExpFunction Method
Consider the general nonlinear partial differential equation of the type
Using a transformation
where and are constants, we can rewrite (2.1) in the following nonlinear ODE;
according to the expfunction method, which was developed by He and Wu [39], we assume that the wave solutions can be expressed in the following form
where , and are positive integers which are known to be further determined, and are unknown constants. We can rewrite (2.4) in the following equivalent form:
This equivalent formulation plays an important and fundamental part for finding the analytic solution of problems [5, 17, 39, 40, 53, 54, 58, 59, 75, 76, 84â€“87, 97â€“100]. To determine the value of and , we balance the linear term of highest order of (2.4) with the highest order nonlinear term. Similarly, to determine the value of â€‰and â€‰, we balance the linear term of lowest order of (2.3) with lowest order nonlinear term.
Example 2.1 (see [58]). Consider the ZKMEW (2.6) Introducing a transformation as we can covert (2.6) into ordinary differential equations where the prime denotes the derivative with respect to . The trial solution of the (2.7) can be expressed as follows, as shown in (2.6): To determine the value of and , we balance the linear term of highest order of (2.7) with the highest order nonlinear term, we obtain
Case 1. We can freely choose the values of but we will illustrate that the final solution does not strongly depend upon the choice of values of c and d . For simplicity, we set and , then the trial solution yields Substituting (2.10) into (2.7), we have where are constants obtained by Maple 11. Equating the coefficients of to be zero, we obtain Solution of (2.12) will yield We, therefore, obtained the following generalized solitary solution of (2.6) as follows: where , and â€‰are real numbers.Figure 1 depicts the soliton solutions of (2.6), when In case and are imaginary numbers, the obtained soliton solutions can be converted into periodic solutions or compactlike solutions. Therefore, we write and consequently, (2.14) becomes The above expression can be rewritten in expanded form: where . If we search for periodic solutions or compactlike solutions, the imaginary part in (2.16) must be zero, hence Figure 2 depicts the periodic solutions of (2.6) when .
Case 2. If and then the trial solution, (2.6) reduces to Proceeding as before, we obtain Hence, we get the generalized solitary solutions of (2.6) as follows: where , and k are real numbers.
Remark 2.2. It is worth mentioning that the transformation which is used to transform the solitary solutions to periodic or compactonlike solutions was first proposed by He and Wu [39] and is called the HeWu transformation. Moreover, the interpretation of this transformation is given by He [17].
3. Variational Iteration Method (VIM) and its Modifications
To illustrate the basic concept of the Heâ€™s VIM, we consider the following general differential equation:
where is a linear operator, a nonlinear operator and g (x ) is the inhomogeneous term. According to variational iteration method [1â€“7, 11, 17â€“24, 44, 49, 50, 60â€“64, 66â€“74, 77â€“80, 88, 90], we can construct a correction functional as follows: where â€‰is a Lagrange multiplier [17â€“24], which can be identified optimally via variational iteration method. The subscripts n denote the n th approximation, is considered as a restricted variation. That is, (3.2) is called a correction functional. The solution of the linear problems can be solved in a single iteration step due to the exact identification of the Lagrange multiplier. The principles of variational iteration method and its applicability for various kinds of differential equations are given in [17â€“24]. In this method, it is required first to determine the Lagrange multiplier optimally. The successive approximation of the solution u will be readily obtained upon using the determined Lagrange multiplier and any selective function consequently, the solution is given by We summarize some useful iteration formulae [23, 24] which would be used in the subsequent section:
3.1. Variational Iteration Method Using Heâ€™s Polynomials (VIMHP)
This modified version of variational iteration method [60, 66â€“69, 72, 73] is obtained by the elegant coupling of correction functional (2.7) of variational iteration method (VIM) with Heâ€™s polynomials and is given by
comparisons of like powers of give solutions of various orders.
3.2. Variational Iteration Method Using Adomianâ€™s Polynomials (VIMAP)
This modified version of VIM is obtained by the coupling of correction functional (2.3) of VIM with Adomianâ€™s polynomials [1, 2, 62, 64, 70, 74] and is given by
where are the socalled Adomianâ€™s polynomials and are calculated for various classes of nonlinearities by using the specific algorithm developed in [81â€“83].
Example 3.1 (see [61]). Consider the following singularly perturbed sixthorder Boussinesq equation with initial conditions The exact solution of the problem is given as The correction functional is given by Making the correction functional stationary, the Lagrange multiplier can easily be identified as we get the following iterative formula Consequently, following approximants are obtained The series solution is given by Table 1 exhibits the absolute error between the exact and the series solutions. Higher accuracy can be obtained by introducing some more components of the series solution.

Example 3.2 (see [72]). Consider the following Helmholtz equation: with initial conditions The correction functional is given as Making thecorrection functional stationary, the Lagrange multiplier can be identified as we obtained Applying the variational iteration method using Heâ€™s polynomials (VIMHP), we get Comparing the coefficient of like powers of p , following approximants are obtained: The solution is given as Table 2 exhibits the approximate solution obtained by using the HPM, ADM, and VIMHP. It is clear that the obtained results are in high agreement with the exact solutions. Higher accuracy can be obtained by using more terms.


Example 3.3 (see [62]). Consider the following nonlinear SchrÃ¶dinger equation with initial conditions The correction functional is given as Making the correction functional stationary, the Lagrange multipliers can be identified as consequently Applying the variational iteration method using Adomianâ€™s polynomials (VIMAP): where are the socalled Adomianâ€™s polynomials. First few Adomianâ€™s polynomials for nonlinear SchrÃ¶dinger equation are as follows: Employing these polynomials in the above iterative scheme, following approximants are obtained: The solution in a series form is given by and in a closed form by
Remark 3.4. It is worth mentioningthat although both the modified versions of variational iteration method (VIM) are compatible yet the modification based upon Heâ€™s polynomials (VIMHP) is much easier to implement and is more user friendly as compared to VIMAP where Adomianâ€™s polynomials along with their complexities are used.
4. Homotopy Perturbation Method (HPM) and Heâ€™s Polynomials
To explain the Heâ€™s homotopy perturbation method, we consider a general equation of the type,
where is any integral or differential operator. We define a convex homotopy by
where is a functional operator with known solutions , which can be obtained easily. It is clear that, for
we have
This shows that continuously traces an implicitly defined curve from a starting point to a solution function . The embedding parameter monotonically increases from zero to unit as the trivial problem continuously deforms the original problem The embedding parameter can be considered as an expanding parameter [13, 14, 17, 18, 26â€“38, 51, 52, 55â€“57, 60, 65â€“69, 72, 73, 89, 91â€“93, 96]. The homotopy perturbation method uses the homotopy parameter p as an expanding parameter [17, 18, 26â€“38] to obtain
if 1, then (4.5) corresponds to (4.2) and becomes the approximate solution of the form,
It is well known that series (4.5) is convergent for most of the cases and also the rate of convergence is dependent on (); see [17, 18, 26â€“38]. We assume that (4.6) has a unique solution. The comparisons of like powers of p give solutions of various orders. In sum, according to [13, 14], Heâ€™s HPM considers the nonlinear term as
where â€™s are the socalled Heâ€™s polynomials [13, 14], which can be calculated by using the formula
Example 4.1 (see [51]). Consider the following seventhorder generalized KdV (SOGKdV) equation where , with initial conditions where , and is an arbitrary parameter. Applying the convex homotopy, we get Comparing the coefficient of like powers of p where are the Heâ€™s polynomials. The series solution is given by