A Novel Version of HPM Coupled with the PSEM Method for Solving the Blasius Problem

Filobello-Nino, Uriel; Vazquez-Leal, Hector; Huerta-Chua, Jesus; Jimenez-Fernandez, Victor Manuel; Sandoval-Hernandez, Mario Alberto; Delgado-Alvarado, Enrique; Tlapa-Carrera, Victor Manuel

doi:https://doi.org/10.1155/2021/5909174

Discrete Dynamics in Nature and Society

On this page

Abstract Introduction Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 5909174 | https://doi.org/10.1155/2021/5909174

A Novel Version of HPM Coupled with the PSEM Method for Solving the Blasius Problem

Uriel Filobello-Nino,¹Hector Vazquez-Leal,^1,2Jesus Huerta-Chua,³Victor Manuel Jimenez-Fernandez,¹Mario Alberto Sandoval-Hernandez,⁴Enrique Delgado-Alvarado,³and Victor Manuel Tlapa-Carrera³

Academic Editor: Marcio Eisencraft

Received20 May 2021

Revised02 Jul 2021

Accepted20 Jul 2021

Published06 Aug 2021

Abstract

This work studies the nonlinear differential equation that models the Blasius problem (BP) which is of great importance in fluid dynamics. The aim is to obtain an approximate analytical expression that adequately describes the phenomenon considered. To find such approximation, we propose a new method denominated powered homotopy perturbation (PHPM). Unlike HPM algorithm, the successive integration process generated by PHPM will consider zero the constants of integration in each approximation, except the last one. In the same way, PHPM will propose an adequate initial trial function provided of some unknown parameters in such a way that it will not evaluate the initial conditions in the iterations of the process; therefore, this set of parameters will be employed with the purpose of adjusting in the best accurate way the proposed approximation at the final part of the process. As a matter of fact, we will note from this analysis that the proposed solution is compact and easy to evaluate and involves a sum of five exponential functions plus a linear part of two terms, which is ideal for practical applications. With the purpose to get a better approximation, we find useful to combine PHPM with the power series extender method (PSEM) which implies to add to the PHPM solution one rational function with parameters to adjust. From this proposal, we find an approximate solution competitive with others from the literature.

1. Introduction

As it is well known, the fluid motion is divided into two clearly defined regions; the first one concerns the region near the object where the friction effect is appreciable and is known as boundary layer. On the other hand, the second region is that where these effects can be neglected [1,2]. The boundary layer is characterized for being the region where the fluid velocity parallel to the surface is less 99% of the free stream velocity. Nevertheless, even in the case of a laminar flow, the solution of the equations describing the laminar boundary layer is complicated [1,2]. This work will search for a solution for the BP, which describes the two-dimensional steady state viscous flow over a semi-infinite flat plate. Specifically, the aim is to get an analytical approximate solution that adequately describes this relevant fluid mechanics problem. The difficulties of these scientific problems give rise to the proposal of new methods with the end to obtain solutions to the differential equations that govern nonlinear problems such as the one proposed here; however, the search for such solutions is not an easy task and justifies all the research efforts carried out on this topic. In effect, an exact solution to the proposed nonlinear problem can rarely be obtained [3], and for the same reason, several methods have been added to the best known classical methods. Next, we provide a list with some of most employed analytical methods in accordance with the literature: variational approaches [4–6], the tanh method [7], exp-function [8,9], Adomian’s decomposition method [10–15], parameter expansion [16], the homotopy perturbation method (HPM) [1, 6, 17–33], the perturbation method [34–36], the modified Taylor series method [37], the Picard method [38], the PSEM [39–42], the homotopy analysis method [25, 43], the variational iteration method [44], and the differential transform method [45], among others. Numerical methods play an important role with the purpose to provide algorithms to get approximations of differential equations [46–49]. This work introduces two novelties; on one side, it presents an effective modification of the standard HPM, the power homotopy perturbation method (PHPM), with the purpose of finding an analytical approximate solution for the BP. As it is well known, the HPM has been one of the most efficient methods proposed to get analytical approximate solutions for nonlinear problems. The HPM was first proposed by a Chinese mathematician Ji-Huan He to investigate a wide variety of nonlinear problems which arise in science and engineering. The HPM successfully couples the homotopy theory with the perturbation theory (but it is not restricted to small parameters as with traditional perturbation methods). The main idea is that a complicated problem is continuously deformed into a simple problem which is easy to solve in order to get an analytic exact or approximate solution [20]. On the other hand, we will solve the same problem, but now, the PHPM will provide a solution just for a part of the whole interval of the BP. Instead of looking for a solution for the other part of the interval and introducing a piece-wise kind solution [41] for the whole interval and even proposing, for instance, to express the final result in a compact expression by using the unit step function [50], we will employ the PSEM to obtain an approximate solution valid for all the definition domains of the proposed problem. As a matter of fact, after comparing the approximate solution obtained for this work with other approximations from the literature, we will see that the combination of the PHPM and PSEM offers, indeed, potential for the approximate solution of nonlinear problems. The rest of this work is as follows: Section 2 introduces the basic idea of the standard HPM. Section 3 explains thePHPM. Additionally, Section 4 explains the basic concept of the PSEM. On the other hand, Section 5 presents, with some details, the necessary elements for the nonlinear problem to solve. Besides, Section 6 presents the application of the proposed methods, in the search for an approximate solution for the third-order nonlinear ordinary differential equation which describes the BP. Section 7 offers a whole discussion about the proposed solutions and their comparison with other approximations found in the literature. Finally, a brief conclusion of the relevant aspects of this work is given in Section 8.

2. Standard HPM

The standard HPM was conceived with the end to approach various kinds of nonlinear problems [19,20]. The HPM is a combination of the classical perturbation technique and the homotopy (whose origin lies in the branch of mathematics denominated topology), but it is not restricted to small parameters as with traditional perturbation methods (to say PM). For example, the HPM requires neither small parameter nor linearization and usually requires few iterations to obtain solutions with good precision [19,20].

To conceive how the HPM works, we consider a general nonlinear differential equation, which can be expressed as follows:with the following boundary conditions:where is a general differential operator, is a boundary operator, is a given analytical function, and is the domain boundary for . Next, assuming that can be divided into two operators and , where is linear and nonlinear; (1) can be expressed as

An homotopy can be constructed in accordance with [19,20].orwhere is a homotopy parameter, whose values belong to the range of ; is the first approximation for the solution of (3) that satisfies the boundary conditions.

Next, we will assume that the solution for (4) or (5) can be written as a power series of in the formand after substituting (6) into (5) and equating identical powers of terms, there can be found values for the sequence , , , after solving a coupled system of differential equations.

Taking the limit , an approximate solution for (1) is obtained as follows:where denotes the starting point for the initial conditions and , , , are the initial conditions of the problem.

3. The Proposed Method: Powered Homotopy Perturbation Method (PHPM)

From the above, we note that the standard HPM establishes that the values of higher-order approximations , , have to be zero if they are evaluated in the initial condition value . The proposed method PHPM essentially follows the same steps of the HPM, but unlike it, the successive integration process of the coupled system (7)–(10) is performed without evaluating the initial conditions mentioned by the HPM. In a sequence, the process of successive approximations (7)–(10) will be affected considering zero the constants of integration in each approximation, excepting the last one, which will keep all their integration constants that are originated from the process of solution. In this stage, one approximate solution for the problem to solve is obtained again considering the limit (11). Nevertheless, we note that this solution possesses some integration constants to evaluate, even other parameters that could rise of the freedom of the homotopy technique [22,23]. These parameters are evaluated from the initial conditions applied to (11) and the knowledge of other mathematical and physical properties of the nonlinear problem. Therefore, the solution for a problem is expressed in terms of solving an algebraic system of equations.

The scheme for the PHPM, analogous to (7)–(10), isassuming that is the last iteration considered in order to provide an approximate solution for (1), and then, it is possible to rewrite (15) symbolically aswhere and are inverse operators and , , , are constants to determine.

Taking the limit , it is possible to obtain an approximate solution for (1) as follows:

The constants , , , are found from the conditions

We will see that the Blasius problem is a good candidate to apply this scheme.

4. Basic Concept of the PSEM

In accordance with [39], we begin assuming a nonlinear differential equation which is expressed as follows:

The boundary conditions of (19) can be expressed aswhere is a boundary operator; is a given analytical function; is the boundary for domain ; and and represent linear and nonlinear operators, respectively. Finally, expresses differentiation along the normal drawn outwards from . The following step consists in assuming that the solution of (19) can be expressed as a power series:where represents the coefficients of the power series.

In accordance with PSEM flexibility, series (21) can be obtained by some approximate methods from the literature: the Taylor series method [37], power series method [39], homotopy perturbation method [1,17–23], perturbation method [34–36], homotopy analysis method [43], variational iteration method [44], differential transform method [45], and Adomian decomposition method [10–15], among others.

Next, we assume that the solution for differential equation (19) can be expressed as a finite sum of certain functions in accordance with [39].orwhere are, in principle, unknown constants to be determined by the PSEM, are arbitrary trial functions, and and represent the orders of approximation of (22) and (23), respectively. As a matter of fact, from now on, we will designate (22) and (23) as the PSEM trial function (TF).

The next step consists in obtaining the Taylor series of (22) or (23) so that we get the power series of the formwhere Taylor coefficients are expressed in terms of parameters . In particular, are the coefficients of the order zero Taylor expansion applied to the finite sum of functions (22) and rational function (23), respectively. Finally, we equate the coefficients of the corresponding power series (24) or (25) with (21) in order to get the values . After substituting them into (22) or (23), we finally obtain the PSEM approximation.

5. Problem Formulation

This section will provide the sufficient elements required for the rest of this work. As an application of the PHPM, we will get an analytical approximate solution for the nonlinear third-order Blasius ordinary differential equation (BODE) which can be written in the following terms:where denotes the dimensionless stream function and is a dimensionless coordinate proportional to the perpendicular distance from the front of the plate, inversely proportional to the thickness of the boundary layer. As it is well known, the boundary layer is defined as the adjacent region to a solid surface where viscous forces are relevant and its thickness is defined where the derivative is . From [51], it results that the values of and defining the boundary layer are related by (see Section 6).

On the other hand, Blasius found a series solution for system (26) and (27) [52] which is given byand was numerically determined with the value 0.332057336 [53,54].

Such as it occurs for every mathematical series, (28) converges only for a finite interval (see [51]), where is a number such that or ; therefore, this series does not totally describe the Blasius problem, given that it requires a solution for the whole interval .

Since the aim of this work is to provide an analytical approximate solution valid for the entire domain and not only for small , it will result useful to express condition (27) in the asymptotic form:for some constant value of . As a matter of fact, we will see that the geometric determination of this parameter will be relevant for the solution of the asymptotic part of this problem.

6. Application of the PHPM with the End to Obtain Two Approximate Solutions for the Nonlinear Blasius Problem

In accordance with the previous section, we have to find an approximate solution for BODE.

From Section 3, we establish the following homotopy, by using the flexibility of PHPM:or

Following PHPM algorithm, we will assume that (33) admits a series solution of the formin such a way that, after substituting (34) into (33), we get the following system of linear differential equations:and so on.

Next, we propose a solution depending of two parameters, from the asymptotic characteristics of the Blasius problem, and let it bewhere positive quantities and will be determined from the properties of the problem.

Assuming that we keep until the third iteration of (34), we will solve (35)–(38).

Thus, by substituting (38) into (35), we getand after integrating three times, without keeping the integration constants, we obtain

In the same way, substituting (38) and (40) into (36), we get

Integrating three times, without keeping the integration constants, we obtain

The third iteration results from replacing (38), (40), and (42) into (37) so that we get

It is possible to integrate (43) in the same way three times; assuming that it is the last iteration, we have to add tree integration constants produced by tree integrations in such a way that we write

We note that, from asymptotic condition (30), it is required that the last term of (44) has to be zero; besides, from the same condition, it is clear that . Therefore, (44) has to adopt the form

After taking the limit , we obtain an approximate solution for (26) with boundary conditions (27) as follows:see (38), (40), (42), and (45).

From the above, it is clear that (46) is an analytical approximate solution for (31) which satisfies the correct asymptotic condition (30).

In accordance with (46), we have obtained an approximate solution with the schematic formwhere represents all the terms of (46) which contain only negative exponentials.

As it was already mentioned, the boundary layer is defined as the adjacent region to a solid surface where viscous forces are relevant and its thickness is defined where the derivative is .

Also, it was mentioned that this point corresponds to the coordinates , so that, for , it is expected that (47) adopts the asymptotic form (30) with good precision. Thus, and

Geometrically, (48) is the equation of the tangent line to the point and can be determined by substituting the coordinates into (48) to obtain the equation

Thus, (47) can be written aswhich ensures the asymptotic form of the solution, even from the point for which .

With the purpose of ensuring that (50) satisfies the initial conditions of the problem, it is necessary to adequately evaluate the values of the constants and . For it, we will substitute the first condition and the second condition of (27) into (50) in order to get an algebraic system for the unknown quantities.

The solution for nonlinear system (51) is given by the values

Thus, substituting these values into (50), we get an approximate solution for (26) with boundary conditions (27).

From Table 1, we observe that the proposed approximate solution has, in general terms, an acceptable precision. The comparison between our approximate solution and the exact solution shows a bad performance at the beginning of the interval and a gradual improvement; as a matter of fact, from onwards, the absolute error committed is less than 0.2 (e.g., for , the committed relative error is about ). Better results are obtained for values relatively big of the dimensionless coordinate. Although our approximation provides a good idea of the problem under study, next we propose other approximate solutions starting from the following idea.

From iterations (38), (40), (42), and (45), we note that as higher-order iterations are considered, a bigger denominator appears in the successive iterations. Given that we have considered until the third-order approximation, to considerably improve, more iterations have to be included, and for the same reason, we would have to consider a long cumbersome final approximation. In order to avoid that mentioned above, we will introduce two approximations, the first one for the interval and the second for the infinite interval .

The solution for will consist in recalculating the parameters (52) to obtain an analytical approximate solution for (31) in the mentioned interval. From series solution (21), it is possible to get the following four terms:

As we commented in Section 4, series (28) converges for values in the following interval , where is a number such that or . In a sequence, the four terms from (28) provide the main contribution in this interval; thus, with the purpose of recalculating the values of the parameters from (42), we propose substituting the following points of the solution curve into (50) in order to obtain an algebraic system for the unknown quantities and , and , which are obtained from (45).

After solving the resulting 2 × 2 nonlinear system, we get

Therefore, by substituting (54) into (50), we get an analytical approximate solution for the Blasius problem valid for .

On the other hand, in order to get a solution for , we will employ the PSEM. As explained in Section 4, next we assume the following rational function valid for the abovementioned interval:

The coefficients of (55) and are determined by using PSEM algorithm, and the presence of is with the end to contribute that (55) vanishes adequately.

In accordance with the PSEM, we propose the sum of (55) and (50) (with coefficients (54)) as a solution valid for the whole interval.

Next, we obtain a Taylor expansion around for (56), and this result is matched with series (53) (that is, we equate the coefficients of both power series) with the purpose to get an algebraic equation system with the end to determine the coefficients of (55).

As a consequence, we determine the following values:

Thus, substituting (57) into (56), we obtain an analytical approximate solution for the BP. From Table 1, we note that the accuracy of (56) with coefficients (54) and (57) is clearly better than the one obtained for the PHPM (comparing the first two columns with the exact solution).

7. Discussion

Next, we will present the main results obtained for this work. In particular, we will emphasize the advantages of the proposed PHPM for the case where there exists additional information about the problem to solve.

This article presented a modification of the HPM, the modified HPM (PHPM), with the purpose of obtaining an analytical approximate solution for the relevant BP. Essentially, the procedure follows the standard HPM, but avoiding the integration constants in each iteration. Once one decides what is the last iteration to conclude the iterative procedure, the constants are added in this step, and the final approximate solution contains these parameters. We also mentioned that this scheme could result useful in the cases where, for example, the knowledge of certain mathematics and physical information could aid in knowing the value of some of these parameters (of course, some constants are determined from the initial and boundary conditions proposed by the problem). To begin, we were able to express the solution of this problem as the sum of an asymptotic part and other valid for relatively small values of . Along this work, we showed the manner to introduce some conditions with the end to determine the values of some constants. For example, the boundary condition implies that asymptotically, the derivative function tends to the constant value one and for the same reason, the dimensionless stream function has to assume, in this same regime, the rather simple form , .

Comparing this asymptotic form with (44), we noted that it is required that the last term of (44) has to be zero; besides, from the same condition, it is clear that . Therefore, the particular form of expressing the PHPM solution easily allows calculating these values. In the same way, we saw that the Blasius asymptotic solution admits a geometry interpretation which allowed knowing the value of . Therefore, we note the convenience of using the PHPM in order to determinate the asymptotic part of the proposed solution, taking advantage of the mathematics and physical information that concerns with the problem.

On the other hand, the PHPM easily allowed calculating the other unknown parameters and . This procedure was implemented in order to ensure that (50) satisfied the initial conditions of the problem. In accordance with PHPM algorithm, we substituted (50) into the first condition and the second condition of (27) with the end to obtain an algebraic system for the unknown quantities.

The analysis from Table 1 shows a moderate precision; as a matter of fact, we got the best performance for values next to and higher than this value. This result finds an explanation in the fact that the PHPM solution ensured, from the beginning, the correct behaviour of the asymptotic part (50). Besides, from the same table, we note that the absolute error for values between and is in the order of two to three hundredths, which represent an acceptable precision. It is expected that, for improving the precision of (50), higher-order iterations be required and/or an initial trial function that looks like (38) be provided but with more adjusting parameters (other trial functions could be tried with the correct asymptotic behaviour of the problem, for example, a sum of two terms like (38), and we will provide fourth adjusting parameters) in order to get a better result.

It is important to note the power and flexibility of the homotopy scheme proposed in this work. However, in future research, we propose exploring alternative ways to construct other homotopy schemes. For example, the homotopy perturbation method with two expanding parameters [28], homotopy perturbation method with three expansions [29], Li–He homotopy perturbation method [30], homotopy perturbation with an auxiliary parameter [31], and He Laplace method [32, 33].

7.1. Combination of the PHPM and PSEM

Next, in order to obtain a better approximation, it was proposed to recalculate the parameters and mentioned above, but changing the strategy. We conceived to divide the solution of the problem into two parts, the first one valid for the interval and the second for the infinite interval .

With respect to the solution in the infinite interval, the procedure consisted in recalculating the parameters (54) with the purpose of obtaining a new analytical approximate solution for (26)–(27) in the mentioned interval.

That is, we repeated the same steps that yield in (50) with parameters (52) but shortening the interval as it was indicated. We expected to improve the above approximation because the solution of the problem at the beginning of the interval is strongly affected, not only from the asymptotic term , but for , while as we consider bigger values of , the asymptotic part becomes the predominant, for which the solution for both regions is indeed different. Therefore, if we ask the PHPM to describe only partially the interval, it is expected to get a better approximation, which, in fact, was what occurred.

The manner to recalculate the values of the and of (50) was substituting some points of the solution curve into (50) in order to obtain an algebraic system for the mentioned unknown quantities. The points used for this purpose are , and , which were obtained from (53) (see also (28)). After solving the resulting 2 × 2 nonlinear system, we got the values (54). It is important to clarify why we did not employ the points correspondent, for instance, to and , instead of the points mentioned above which correspond to and . By using and , the PHPM would have obtained a solution only for the asymptotic part of the problem with a very good expectation to describe it with precision. We emphasize that two points are required for this purpose, which are obtained from the Blasius series. As explained in Section 5, the convergence of this series lies in the interval , where is a number such that .

Thus, to avoid calculating points from the Blasius series too close to its convergence upper end point, we preferred using values of , totally within the convergence interval.

In order to get a good approximation for the interval , a simple choice would have been to use (28). After all, we have seen that, for within its convergence interval, even four terms from (28) provide good results (53). In this way, we could provide a piece-wise kind of solution, and employing unit step functions , we may get a final result by using a compact expression. In a sequence, in theory, it would have been even possible to get a differentiable approximate expression by using differentiable and continuous analytical approximate functions for the step function. However, in general, many times, that mentioned earlier is not possible but only approximated because to match two approximate solutions in the join points is not an easy task.

With the purpose of getting a real continuous and differentiable solution for the whole interval of this problem, we proposed using the PSEM in order to get an accurate solution that includes the interval .

With this end, the PSEM proposed the rational function (55) and, therefore, the sum of (50) and (55) (with coefficients (54) for (50)) as a solution for the entire interval. Next, a Taylor expansion is obtained around for (56), and what results is matched, term by term, with series (53) up to the eleventh power in order to propose an algebraic system of equations to determine the coefficients of (55). This procedure is successful because it not only ensures the correctness of continuity and the differentiability for all the points of the whole interval but also takes advantage of both the accuracy of (53) in the initial part of the interval and the accuracy of the PHPM solution for infinite interval .

Finally, from Table 1, we compared the results obtained for this work with other approximations reported in the literature. In columns two and three appear the approximations discussed for this work.

As commented, in general, the PHPM with parameters (52) had an acceptable precision. Our approximation showed a bad approximation in the beginning of the interval, and after a gradual improvement, we note that, from onwards, the absolute error committed is less than 0.2; in particular, for relatively large values of , better results are obtained. Compared with the PHPM (54) and PSEM, the PHPM (52) obtained a lower performance, and the reasons were, in part, the subject of Section 6. Nevertheless, it is important to note that the performance of the PHPM is essentially the same for big values of .

On the other hand, the LTNHPM [57] showed a better performance than the PHPM (52), but it only provides results from to , whereby the proposed solution from the PHPM offered a better performance in general. In this work [57], a series with good results was found, but in general, the power series requires many terms in order to get an accurate approximate solution when it is applied to differential equations. As a matter of fact, the series that gave rise to values in Table 1 required 13 terms, and as we consider bigger values of , more and more terms of the series are needed, which is not convenient for practical applications. From the results of Table 1, it is clear that [1] presents a better precision than the PHPM (52) for values comprised between and , while from onwards, the PHPM (52) is better. It ratifies the correctness of the implemented strategy for the PHPM with the purpose of calculating the asymptotic part of the BP. Another point to consider is unlike the proposed method, the method in [1] calculated its unknown coefficients by using an optimization method, the least square method, which represents an advantage for the method in [1].

With respect to the comparison between the PHPM (52) and IPM [55], Table 1 shows that the PHPM provides a better accuracy, except for values comprised from to , which is very early in the interval. In general, the IPM offered the lower performance of the compared methods. On the contrary, Ahmad et al. [51] showed the best precision, which is deduced from Table 1. Nevertheless, Ahmad et al. [51] employed a Padé approximant of the derivative of (53) as the starting point and added corrective terms both in the numerator and denominator of the mentioned Padé to propose a new expression which yielded good precision. This proposal even successful for the BP is not a general procedure for finding solutions to nonlinear differential equations. Finally, the comparison of the proposed solution PHPM (52) with [56] shows that the proposal from the WRM is better for interval although the PHPM was competitive in that interval except at the early small values close to zero. Again, the PHPM (52) was better for the asymptotic part of the solution.

As explained, in order to improve the PHPM (52) approximate solution, we introduced the combination of the PHPM and PSEM. From Table 1, we noted that the obtained result was remarkably better. The combination of the PHPM (54) and PSEM is more precise than the PHPM (52), and it is better than the LTNHPM [57] although in interval , the performance of both was very similar. In the same way, the proposed solution was better than that in [1] although this employed an optimization procedure to calculate its coefficients. Likewise, we note that the PHPM (54) and PSEM were clearly better than the IPM and also overcome the WRM despite that this last is a method widely recognized. In this case, both methods achieved a good performance. Finally, as already said, the method in [51] obtained the best accuracy, but the difference of accuracy between the solutions obtained by the PHPM (54) and PSEM [51] was, most of the times, of only thousandths. An important point is that unlike the one proposed by this work, Ahmad et al. [51] just provided a particular procedure valid in principle for the Blasius problem. Finally, even we could present other HPM approximations for the Blasius equation [21], but this article only presented results for the short interval . Figure 1 shows the absolute error committed for the method solutions exposed earlier.

In brief, the results obtained for the PSEM combined with the PHPM showed that it is possible to couple two methods, which, in principle, are very different, but combining them can provide a total solution with a wider domain; as a matter of fact, the PHPM delivered the asymptotic part of the problem and the PSEM focused on solving the part closer to , in such a way that the coupling of both methods yielded a new function defined for the whole interval with good precision, as discussed, in relation with Table 1. The matter of improving the convergence of the PHPM and PHPM and PSEM combination could be considered in a future theme of investigation. As a matter of fact, methodologies such as those considered in [24, 30] could be useful with this purpose.

8. Conclusions

This work introduced a modification of the standard HPM, the PHPM, which is particularly useful when it is possible to dispose additional information, both mathematical and physical, for a proposed nonlinear problem. Unlike the HPM which establishes that the values of higher-order approximations have to be zero if they are evaluated for the initial condition, the PHPM proposed that the successive integrations process of the coupled system (7)–(10) is performed without evaluating the initial conditions mentioned by the HPM, considering zero the constants of integration in each iteration, except for the last one, which will keep all its integration constants that are originated by the process. The whole solution (11) will possess some integration constants to evaluate even other parameters that could rise from the freedom of the homotopy technique, and the PHPM evaluates these quantities from the initial conditions applied to (11) and the knowledge of other mathematical and physical properties of the nonlinear problem. The proposed PHPM expresses in a natural way the solution as a sum of an asymptotic part and other part that describes the solution for values closer to the origin. From the evidence presented for this article, it is expected that PHPM contributes to change the idea that an effective method has to be necessarily long and cumbersome as occurs with the HAM and Adomian decomposition method, among others the proposal of PHPM is that of a method that is both precise and easy to implement the proposed method allows incorporating the knowledge of the mathematical and physical properties that emanate from the problem to solve and not only the initial conditions

On the other hand, as occurred in this work, sometimes, it is expected that, for improving the precision of the proposed solutions for the PHPM it is required both to obtain higher order iterations and to introduce an initial trial function provided with several adjusting parameters, which could lead to a cumbersome procedure. To avoid that inconvenience, this work improved the results obtained by the PHPM by using a combination of the PHPM and PSEM which showed to have high potential for practical problems. The success of the combination is that the PSEM took advantage of the good performance of a known series solution for the initial part of the interval and the PHPM provided the asymptotic part of the solution with good accuracy. In summary, the PHPM and the union of the PSEM with PHPM demonstrated potential for future applications.

Nomenclature

BP:	Blasius problem
HPM:	Homotopy perturbation method
PHPM:	Powered homotopy perturbation method
PSEM:	Power series extender method
PM:	Perturbation method
BODE:	Blasius ordinary differential equation
LTNHPM:	Laplace transform coupled with the new homotopy perturbation method
IPM:	Iteration perturbation method
WRM:	Weighted residual method.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors would like to thank Rogelio Alejando Callejas-Molina and Roberto Ruiz Gomez for their contribution to this project.

References

U. Filobello-, H. Vazquez-Le, R. Castaneda- et al., “An approximate solution of blasius equation by using hpm method,” Asian Journal of Mathematics & Statistics, vol. 5, no. 2, pp. 50–59, 2012.
View at: Publisher Site | Google Scholar
L. D. Landau and E. M. Lifshitz, “Theoretical physics,” Fluid Mechanics, vol. 6, 1987.
View at: Google Scholar
U. Filobello-Nino, H. Vazquez-Leal, Y. Khan et al., “A handy exact solution for flow due to a stretching boundary with partial slip,” Revista mexicana de física E, vol. 59, no. 51–55, 2013.
View at: Google Scholar
L. M. B. Assas, “Approximate solutions for the generalized KdV-Burgers’ equation by He’s variational iteration method,” Physica Scripta, vol. 76, no. 2, pp. 161–164, 2007.
View at: Publisher Site | Google Scholar
M. Kazemnia, S. A. Zahedi, M. Vaezi, and N. Tolou, “Assessment of modified variational iteration method in bvps high-order differential equations,” Journal of Applied Sciences, vol. 8, no. 22, pp. 4192–4197, 2008.
View at: Publisher Site | Google Scholar
R. Noorzad, A. T. Poor, and M. Omidvar, “Variational iteration method and homotopy-perturbation method for solving burgers equation in fluid dynamics,” Journal of Applied Sciences, vol. 8, no. 2, pp. 369–373, 2008.
View at: Publisher Site | Google Scholar
D. J. Evans and K. R. Raslan, “The tanh function method for solving some important non-linear partial differential equations,” International Journal of Computer Mathematics, vol. 82, no. 7, pp. 897–905, 2005.
View at: Publisher Site | Google Scholar
F. Xu, “A generalized soliton solution of the konopelchenko-dubrovsky equation using He’s exp-function method,” Zeitschrift für Naturforschung A, vol. 62, no. 12, pp. 685–688, 2007.
View at: Publisher Site | Google Scholar
J. Mahmoudi, N. Tolou, I. Khatami, A. Barari, and D. D. Ganji, “Explicit solution of nonlinear zk-bbm wave equation using exp-function method,” Journal of Applied Sciences, vol. 8, no. 2, pp. 358–363, 2008.
View at: Publisher Site | Google Scholar
G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988.
View at: Publisher Site | Google Scholar
E. Babolian and J. Biazar, “On the order of convergence of adomian method,” Applied Mathematics and Computation, vol. 130, no. 2, pp. 383–387, 2002.
View at: Publisher Site | Google Scholar
K. Ali and M. Abadyan, “Efficiency of modified adomian decomposition for simulating the instability of nano-electromechanical switches: comparison with the conventional decomposition method,” Trends in Applied Sciences Research, vol. 7, no. 1, p. 57, 2012.
View at: Google Scholar
K. Ali and M. Abadyan, “Evaluating the ability of modified adomian decomposition method to simulate the instability of freestanding carbon nanotube: comparison with conventional decomposition method,” Journal of Applied Sciences, vol. 11, no. 19, pp. 3421–3428, 2011.
View at: Google Scholar
S. K. Vanani, S. Heidari, and M. Avaji, “A low-cost numerical algorithm for the solution of nonlinear delay boundary integral equations,” Journal of Applied Sciences, vol. 11, no. 20, pp. 3504–3509, 2011.
View at: Publisher Site | Google Scholar
S. Hossien Chowdhury, “A comparison between the modified homotopy perturbation method and adomian decomposition method for solving nonlinear heat transfer equations,” Journal of Applied Sciences, vol. 11, no. 7, pp. 1416–1420, 2011.
View at: Publisher Site | Google Scholar
L.-N. Zhang and L. Xu, “Determination of the limit cycle by he’s parameter-expansion for oscillators in a u3/(1 + u2) potential,” Zeitschrift für Naturforschung A, vol. 62, no. 7-8, pp. 396–398, 2007.
View at: Publisher Site | Google Scholar
H. Aminikhah and M. Hemmatnezhad, “A novel effective approach for solving nonlinear heat transfer equations,” Heat Transfer-Asian Research, vol. 41, no. 6, pp. 459–467, 2012.
View at: Publisher Site | Google Scholar
J.-H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 1, pp. 87-88, 2006.
View at: Publisher Site | Google Scholar
J. -H. He, “Recent development of the homotopy perturbation method,” Topological Methods in Nonlinear Analysis, vol. 31, no. 2, pp. 205–209, 2008.
View at: Google Scholar
J. -H. He, “A coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000.
View at: Publisher Site | Google Scholar
D. D. Ganji, H. Babazadeh, F. Noori, M. M. Pirouz, and M. Janipour, “An application of homotopy perturbation method for non-linear blasius equation to boundary layer flow over a flat plate,” International Journal of Nonlinear Science, vol. 7, no. 4, pp. 399–404, 2009.
View at: Google Scholar
U. Filobello-Nino, H. Vazquez-Leal, A. Sarmiento-Reyes et al., “Laplace transform-homotopy perturbation method with arbitrary initial approximation and residual error cancelation,” Applied Mathematical Modelling, vol. 41, pp. 180–194, 2017.
View at: Publisher Site | Google Scholar
U. Filobello-Nino, H. Vazquez-Leal, M. M. Rashidi et al., “Laplace transform homotopy perturbation method for the approximation of variational problems,” SpringerPlus, vol. 5, no. 1, pp. 1–33, 2016.
View at: Publisher Site | Google Scholar
X.-X. Li and C.-H. He, “Homotopy perturbation method coupled with the enhanced perturbation method,” Journal of Low Frequency Noise, Vibration and Active Control, vol. 38, no. 3-4, pp. 1399–1403, 2019.
View at: Publisher Site | Google Scholar
J.-H. He and Y. O. El-Dib, “Homotopy perturbation method for fangzhu oscillator,” Journal of Mathematical Chemistry, vol. 58, no. 10, pp. 2245–2253, 2020.
View at: Publisher Site | Google Scholar
Y. O. El-Dib, G. M. Moatimid, and N. S. Elgazery, “Stability analysis of a damped nonlinear wave equation,” Journal of Applied and Computational Mechanics, vol. 6, pp. 1394–1403, 2020.
View at: Google Scholar
N. Anjum and Q. Tul Ain, “Application of he’s fractional derivative and fractional complex transform for time fractional camassa-holm equation,” Thermal Science, vol. 24, no. 5, pp. 3023–3030, 2020.
View at: Publisher Site | Google Scholar
J. H. He, “Homotopy perturbation method with two expanding parameters,” Indian Journal of Physics, vol. 88, no. 2, pp. 193–196, 2014.
View at: Publisher Site | Google Scholar
J.-H. He and Y. O. El-Dib, “Homotopy perturbation method with three expansions,” Journal of Mathematical Chemistry, vol. 59, no. 4, pp. 1139–1150, 2021.
View at: Publisher Site | Google Scholar
N. Anjum, J.-H. He, Q. Tul Ain, and D. Tian, “Li-he’s modified homotopy perturbation method for doubly-clamped electrically actuated microbeams-based microelectromechanical system,” Facta Universitatis, Series: Mechanical Engineering, vol. 21, 2021.
View at: Google Scholar
D.-N. Yu, J.-H. He, and A. G. Garcıa, “Homotopy perturbation method with an auxiliary parameter for nonlinear oscillators,” Journal of Low Frequency Noise, Vibration and Active Control, vol. 38, no. 3-4, pp. 1540–1554, 2019.
View at: Publisher Site | Google Scholar
J.-H. He, G. M. Moatimid, and D. R. Mostapha, “Nonlinear instability of two streaming-superposed magnetic Reiner-Rivlin Fluids by He-Laplace method,” Journal of Electroanalytical Chemistry, vol. 895, Article ID 115388, 2021.
View at: Publisher Site | Google Scholar
N. Anjum, C.-H. He, and J. -H. He, Two-scale Fractal Theory for the Population Dynamics, Fractals, Singapore, 2021.
U. Filobello-Nino, H. Vazquez-Leal, A. Sarmiento-Reyes et al., “The study of heat transfer phenomena using pm for approximate solution with dirichlet and mixed boundary conditions,” Applied and Computational Mathematics, vol. 2, no. 6, pp. 143–148, 2013.
View at: Publisher Site | Google Scholar
U. Filobello-Nino, H. Vazquez-Leal, Y. Khan et al., “Using perturbation methods and Laplace-Padé approximation to solve nonlinear problems,” Miskolc Mathematical Notes, vol. 14, no. 1, pp. 89–101, 2013.
View at: Publisher Site | Google Scholar
U. Filobello-Nino, H. Vazquez-Leal, B. Benhammouda et al., “A handy approximation for a mediated bioelectrocatalysis process, related to michaelis-menten equation,” Springer Plus, vol. 3, no. 1, pp. 1–6, 2014.
View at: Publisher Site | Google Scholar
H. Vazquez-Leal, M. Sandoval-Hernandez, R. Castaneda-Sheissa, U. Filobello-Nino, and A. Sarmiento-Reyes, “Modified taylor solution of equation of oxygen diffusion in a spherical cell with michaelis-menten uptake kinetics,” International Journal of Applied Mathematical Research, vol. 4, no. 2, p. 253, 2015.
View at: Publisher Site | Google Scholar
U. Filobello-Nino, H. Vazquez-Leal, A. Perez-Sesma et al., “On a practical methodology for solving BVP problems by using a modified version of picard method,” Applied Mathematics & Information Sciences, vol. 10, no. 4, pp. 1355–1367, 2016.
View at: Publisher Site | Google Scholar
H. Vazquez-Leal and A. Sarmiento-Reyes, “Power series extender method for the solution of nonlinear differential equations,” Mathematical Problems in Engineering, vol. 2015, Article ID 717404, 7 pages, 2015.
View at: Publisher Site | Google Scholar
H. Vazquez-Leal, J. L. Garcia-Gervacio, U. A. Filobello-Nino, M. A. Sandoval-Hernandez, and A. L. Herrera-May, “Psem approximations for both branches of lambert w function with applications,” Discrete Dynamics in Nature and Society, vol. 2019, Article ID 8267951, 9 pages, 2019.
View at: Publisher Site | Google Scholar
U. Filobello-Nino, H. Vazquez-Leal, A. L. Herrera-May et al., “A handy, accurate, invertible and integrable expression for dawson’s function,” Acta Universitaria, vol. 29, Article ID e2124, 2020.
View at: Google Scholar
M. A. Sandoval-Hernandez, H. Vazquez-Leal, U. Filobello-Nino, and L. Hernandez-Martinez, “New handy and accurate approximation for the Gaussian integrals with applications to science and engineering,” Open Mathematics, vol. 17, no. 1, pp. 1774–1793, 2019.
View at: Publisher Site | Google Scholar
S. Liao, “Homotopy analysis method: a new analytic method for nonlinear problems,” Applied Mathematics and Mechanics, vol. 19, no. 10, pp. 957–962, 1998.
View at: Google Scholar
J. -H. He and X. -H. Wu, “Variational iteration method: new development and applications,” Computers & Mathematics with Applications, vol. 54, no. 7, pp. 881–894, 2007.
View at: Publisher Site | Google Scholar
J. M. W. Munganga, J. N. Mwambakana, R. Maritz, T. A. Batubenge, and G. M. Moremedi, “Introduction of the differential transform method to solve differential equations at undergraduate level,” International Journal of Mathematical Education in Science and Technology, vol. 45, no. 5, pp. 781–794, 2014.
View at: Publisher Site | Google Scholar
M. M. Peiravi, J. Alinejad, D. Ganji, and S. Maddah, “Numerical study of fins arrangement and nanofluids effects on three-dimensional natural convection in the cubical enclosure,” Transport Phenomena in Nano and Micro Scales, vol. 7, no. 2, pp. 97–112, 2019.
View at: Google Scholar
M. M. Peiravi and J. Alinejad, “Hybrid conduction, convection and radiation heat transfer simulation in a channel with rectangular cylinder,” Journal of Thermal Analysis and Calorimetry, vol. 140, no. 6, pp. 2733–2747, 2020.
View at: Publisher Site | Google Scholar
J. Alinejad and M. M. Peiravi, “Numerical analysis of secondary droplets characteristics due to drop impacting on 3d cylinders considering dynamic contact angle,” Meccanica, vol. 55, no. 10, pp. 1975–2002, 2020.
View at: Publisher Site | Google Scholar
M. M. Peiravi, J. Alinejad, D. D. Ganji, and S. Maddah., “3d optimization of baffle arrangement in a multi-phase nanofluid natural convection based on numerical simulation,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 30, no. 5, pp. 2583–2605, 2020.
View at: Google Scholar
D. G. Zill, A First Course in Differential Equations with Modeling Applications, Cengage Learning, Boston, MA, USA, 2012.
F. Ahmad, W. H. Al-Barakati, and A. - Barakati, “An approximate analytic solution of the blasius problem,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1021–1024, 2009.
View at: Publisher Site | Google Scholar
B. Heinrich, “Grenzschichten in flüssigkeiten mit kleiner reibung,” Zeitschrift für angewandte Mathematik und Physik, vol. 56, no. 1–37, 1908.
View at: Google Scholar
L. Howarth and L. Bairstow, “On the solution of the laminar boundary layer equations,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 164, no. 919, pp. 547–579, 1938.
View at: Publisher Site | Google Scholar
A. Asaithambi, “Solution of the Falkner-Skan equation by recursive evaluation of Taylor coefficients,” Journal of Computational and Applied Mathematics, vol. 176, no. 1, pp. 203–214, 2005.
View at: Publisher Site | Google Scholar
H. E. Ji-Huan, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006.
View at: Google Scholar
U. C. Karabulut and A. Kiliç, “Blasius denkleminin çözümü için çeşitli teknikler,” Balıkesir Üniversitesi Fen Bilimleri Enstitüsü Dergisi, vol. 20, pp. 129–142, 2018.
View at: Publisher Site | Google Scholar
H. Aminikhah, “Analytical approximation to the solution of nonlinear blasius’ viscous flow equation by ltnhpm,” ISRN Mathematical Analysis, vol. 2012, Article ID 957473, 9 pages, 2012.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Uriel Filobello-Nino 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.

PDF Download Citation

Download other formats

Order printed copies

Views

355

Downloads

450

Citations