Abstract
We apply He's polynomials coupled with the diagonal Padé approximants for solving various singular and nonsingular boundary value problems which arise in engineering and applied sciences. The diagonal Padé approximants prove to be very useful for the understanding of physical behavior of the solution. Numerical results reveal the complete reliability of the proposed combination.
1. Introduction
With the rapid development of nonlinear sciences, many analytical and numerical techniques have been developed by various scientists for solving singular and nonsingular initial and boundary value problems which arise in the mathematical modeling of diversified physical problems related to engineering and applied sciences. The application of these problems involves physics, astrophysics, experimental and mathematical physics, nuclear charge in heavy atoms, thermal behavior of a spherical cloud of gas, thermodynamics, population models, chemical kinetics, and fluid mechanics see [1–68] and the references therein. Several techniques [1–68] including decomposition, variational iteration, finite difference, polynomial spline, differential transform, exp-function and homotopy perturbation have been developed for solving such problems. Most of these methods have their inbuilt deficiencies coupled with the major drawback of huge computational work. He [19–24] developed the homotopy perturbation method (HPM) for solving linear, nonlinear, initial and boundary value problems. The homotopy perturbation method was formulated by merging the standard homotopy with perturbation. Recently, Ghorbani and Saberi-Nadjafi [15, 16] introduced He’s polynomials 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. The basic motivation of this paper is to apply He’s polynomials coupled with the diagonal Padé approximants for solving singular and nonsingular boundary value problems. The Padé approximants are applied in order to make the work more concise and for the better understanding of the solution behavior. The use of Padé approximants shows real promise in solving boundary value problems in an infinite domain; see [42, 50, 56–59]. It is well known in the literature that polynomials are used to approximate the truncated power series. It was observed [42, 50, 56–59] that polynomials tend to exhibit oscillations that may give an approximation error bounds. Moreover, polynomials can never blow up in a finite plane and this makes the singularities not apparent. To overcome these difficulties, the obtained series is best manipulated by Padé approximants for numerical approximations. Using the power series, isolated from other concepts, is not always useful because the radius of convergence of the series may not contain the two boundaries. It is now well known that Padé approximants [42, 50, 56–59] have the advantage of manipulating the polynomial approximation into rational functions of polynomials. By this manipulation, we gain more information about the mathematical behavior of the solution. In addition, the power series are not useful for large values of. It is an established fact that power series in isolation are not useful to handle boundary value problems. This can be attributed to the possibility that the radius of convergence may not be sufficiently large to contain the boundaries of the domain. It is therefore essential to combine the series solution with the Padé approximants to provide an effective tool to handle boundary value problems on an infinite or semi-infinite domain. We apply this powerful combination of series solution and Padé approximants for solving a variety of boundary value problems. Precisely the proposed combination is applied on boundary layer problem, unsteady flow of gas through a porous medium, Thomas-Fermi equation, Flierl-Petviashivili (FP) equation, and Blasius problem. It is worth mentioning that Flierl-Petviashivili equation has singularity behavior at which is a difficult element in this type of equations. We transform the FP equation to a first-order initial value problem and He’s polynomials are applied to the reformulated first-order initial value problem which leads the solution in terms of transformed variable. The desired series solution is obtained by implementing the inverse transformation. The fact that the proposed algorithm solves nonlinear problems without using Adomian’s polynomials is a clear advantage of this technique over the decomposition method.
2. Homotopy Perturbation Method 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 [15, 16, 19–24, 41–50, 60, 63–68]. The homotopy perturbation method uses the homotopy parameter p as an expanding parameter [19–24] to obtain if 1, then (2.5) corresponds to (2.2) and becomes the approximate solution of the form It is well known that series (2.6) is convergent for most of the cases and also the rate of convergence is dependent on ; see [19–24]. We assume that (3.2) has a unique solution. The comparisons of like powers of give solutions of various orders. In sum, according to [15, 16], He’s HPM considers the nonlinear term as
where ’s are the so-called He’s polynomials [15, 16], which can be calculated by using the formula
of various orders.
3. Padé Approximants
A Padé approximant is the ratio of two polynomials constructed from the coefficients of the Taylor series expansion of a function . The Padé approximants to a function are given by [42, 50, 56–59]
where is polynomial of degree at most and is a polynomial of degree at most . The formal power series determine the coefficients of by the equation. Since we can clearly multiply the numerator and denominator by a constant and leave unchanged, we imposed the normalization condition
Finally, we require that have noncommon factors. If we write the coefficient of as
then by (3.6) and (3.7), we may multiply (3.3) by, which linearizes the coefficient equations. We can write out (3.5) in more details as
To solve these equations, we start with (3.6), which is a set of linear equations for all the unknown q’s. Once the q’s are known, then (3.7) gives and explicit formula for the unknown p’s, which complete the solution. If (3.6) and (3.7) are nonsingular, then we can solve them directly and obtain (3.8) [42, 50, 56–59], where (3.8) holds, and if the lower index on a sum exceeds the upper, the sum is replaced by zero:
To obtain diagonal Padé approximants of different order such as , or , we can use the symbolic calculus software Maple.
4. Numerical Applications
In this section, we apply He’s polynomials for solving boundary layer problem, unsteady flow of gas through a porous medium, Thomas-Fermi equation, Flierl-Petviashivili equation, and Blasius problem. The powerful Padé approximants are applied for making the work more concise and to get the better understanding of solution behavior.
Example 4.1 (see [51, 59]). Consider the following nonlinear third-order boundary layer problem which appears mostly in the mathematical modeling of physical phenomena in fluid mechanics [51, 59] with boundary conditions By applying the convex homotopy, we have comparing the co-efficient of like powers of , following approximants are made where are He’s polynomials. The series solution is given as
Example 4.2 (see [51, 57]). Consider the following nonlinear differential equation which governsthe unsteady flow of gas through a porous medium with the following boundary conditions: By applying the convex homotopy method we have By comparing the coefficient of like powers of , the following approximants are obtained: where are He’s polynomials. The series solution is given as The diagonal Padé approximants [51, 57] can be applied to analyze the physical behavior. Based on this, the Padé approximants produced the slope to be and by using Padé approximants we find where Using (4.11)–(4.13) gives the values of the initial slope listed in Table 3. The formulas (4.11) and (4.12) suggest that the initial slope depends mainly on the parameter, where. Table 3 exhibits the initial slopes for various values of. Table 4 exhibits the values of for for.
Example 4.3 (see [56]). Consider the following Thomas-Fermi (T-F) equation [6–13, 17, 31, 33, 34, 54] which arises in the mathematical modeling of various models in physics, astrophysics, solid state physics, nuclear charge in heavy atoms, and applied sciences: with boundary conditions By applying the convex homotopy, Now, we apply a slight modification in the conventional initial value and take instead of where By comparing the coefficient of like powers of , the following approximants are obtained The series solution is given as Setting the series solution is obtained as The diagonal Padé approximants can be applied [56] in order to study the mathematical behavior of the potential and to determine the initial slope of the potential.
Example 4.4 (see [42]). Consider the generalized variant of the Flierl-Petviashivili equation [37] with boundary conditions Using the transformation the generalized FP equation can be converted to the following first-order initial value problem: with initial conditions By applying the convex homotopy, we have The series solution after four iterations is given by where denote and the inverse transformation will yield where denote Diagonal Padé approximants can be applied [42] to find the roots of the FP monopole for Table 9 shows that the roots of the monopole converge to −1 as increases.
Example 4.5 (see [58, 59]). Consider the two-dimensional nonlinear inhomogeneous initial boundary value problem for the integro-differential equation related to the Blasius problem with boundary conditions where the constant is positive and defined by By applying the convex homotopy, we have Proceeding as before, the series solution is given as and consequently Now, we apply the diagonal Padé approximants to determine a numerical value for the constant by using the given condition. Padé approximant of usually converges on the entire real axis [58, 59]. Moreover, is free of singularities on the real axis. Substituting the boundary conditions in each Padé approximant which vanishes if the coefficient of with the highest power in the numerator vanishes. By solving the resulting polynomials of these coefficients, we obtain the values of α listed in Table 10 [58, 59].
5. Conclusion
In this paper, we applied a reliable combination of He’s polynomials and the diagonal Padé approximants for obtaining approximate solutions of various singular and nonsingular boundary value problems of diversified physical nature. The proposed algorithm is employed without using linearization, discretization, transformation, or restrictive assumptions. The fact that the suggested technique solves nonlinear problems without using Adomian’s polynomials is a clear advantage of this technique over the decomposition method.
Acknowledgment
The author is highly grateful to Brig (R) Qamar Zaman, the Vice Chancellor of HITEC University, Taxila Cantt, Pakistan for providing excellent research environment and facilities.