About this Journal Submit a Manuscript Table of Contents
Abstract and Applied Analysis
Volume 2014 (2014), Article ID 495734, 6 pages
http://dx.doi.org/10.1155/2014/495734
Research Article

Constructing Uniform Approximate Analytical Solutions for the Blasius Problem

Department of Statistics and Computer Science, Kunsan National University, Gunsan 573-701, Republic of Korea

Received 17 November 2013; Revised 5 February 2014; Accepted 6 February 2014; Published 12 March 2014

Academic Editor: Lucas Jodar

Copyright © 2014 Beong In Yun. 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

We propose a simple constructive method which assures uniform accuracy of the approximate analytical solutions for the Blasius problem on the semi-infinite interval . The method is based on a weight function having an S-shape to reflect a series solution near the origin and a reference solution far from the origin. Numerical results show the efficiency of the proposed method.

1. Introduction

For the Blasius problem subject to the boundary conditions we recall the well-known properties [13] of the so-called Blasius function as follows:(i) with (ii) with .

Though the Blasius problem looks simple, search for an approximate analytical solution is known to be quite difficult. Until now, in the literature [422], lots of analytical methods have been proposed. Recently, in approximation of the solutions of nonlinear differential equations in unbounded domain, several efficient spectral methods [2327] have been proposed. These methods reduce solving the nonlinear equation to solving a system of nonlinear algebraic equations.

In this paper, we introduce a weight function in (8) whose values cluster to for and to for when is large enough. Then, employing a series approximate solutions for the Blasius function near the origin and a reference solution away from the origin, we propose a weighted averaging method (11) based on the function . The presented analytical solution , a smooth function on interval , highly reflects the near origin solution for and the faraway solution for . Furthermore, the solution can be continuously extended to the semi-infinite interval . For practical performance, a procedure to choose appropriate parameters in is included. In addition, to improve the accuracy of , we propose a corrected approximation formula including an auxiliary term which properly reflects the behavior of the deviation . Results of numerical experiments, compared with the aforementioned existing method [27], illustrate availability of the proposed method.

2. Series Solutions and a Reference Solution

For simplicity we consider the case of and . The power series of the Blasius stream function for this case is known as where and the coefficients are computed from the recurrence [1] In fact, the series becomes This series, however, converges for . In this paper, we will use a partial sum with an integer , for an approximate solution to the Blasius function near the origin.

On the other hand, for a reference solution approximating far from the origin, we consider the following linear function: based on the property (ii) in the previous section.

Figure 1(a) illustrates graphs of the series approximate solutions with on the interval . Therein, the dotted line indicates the numerical solution for the Blasius function . It is observed that overshoots when is even and undershoots when is odd. In addition, Figure 1(b) shows the graph of the reference solution which undershoots . To illustrate motivation of the main idea proposed in the next section, graphs of the differences with and are included in Figure 2, where is replaced by the numerical solution.

fig1
Figure 1: Approximations of the series solutions for each in (a) and the reference solution in (b).
495734.fig.002
Figure 2: Differences with and indicated by the thin lines and the thick line, respectively.

3. Uniform Approximate Analytical Solutions

For some and we introduce a weight function defined as It should be noted that and it is strictly increasing on the interval with for any . In addition, for a large it follows that This implies that the value of goes close to 0 for and to 1 for as increases. Figure 3 shows the graphs of with and , for example.

495734.fig.003
Figure 3: Behavior of the weight function with for each .

Moreover, we can find that the inverse function of takes a form of

In order to improve the accuracy of the approximate solutions for the Blasius function, we propose a weighted average of the series solution and the reference solution as Therein, for given and , we may take the optimal value of , denoted by , which minimizes the -norm of the residual function defined as

From the property (9) of the weight function , it follows that for large enough with This implies that the point is a threshold between the near origin series solution and the faraway reference solution .

We now summarize the procedure to choose the parameters , , and in the proposed solution in (11) as follows.(S1)Considering the undershoot of the reference solution , take an even integer in the series solution which overshoots the Blasius function (see Figure 1).(S2)Choose a length of the interval for some satisfying or .(S3)Find the optimal exponent of which minimizes defined in (12), that is, satisfies As a result, we may expect that the presented approximate solution with the parameters determined by the procedure (S1)–(S3) will become a corrected approximate solution which improves accuracy of both the series solution and the reference solution over the interval .

In addition, we may extend to the semi-infinite interval continuously by setting for all , which assures sufficient accuracy over the interval for as can be observed in Figures 1(b) and 2.

For example, when we take , from Figure 2, we can find and thus we may set . The optimal exponent is which is obtained by the software, Mathematica V.9. By the similar way, we can choose the values of and for other cases of . Table 1 includes the results for the some small values, , where indicates the nearest integer to the optimal exponent .

tab1
Table 1: Values of , , and obtained by ()–().

Figure 4 illustrates the availability of the presented approximate solution with given in Table 1. Additionally, numerical results for the -norm errors of the approximate solution and its derivatives are given in Table 2.

tab2
Table 2: -norm errors of and its derivatives and .
fig4
Figure 4: Approximation of the weighted average with the parameters and its error in (a) and those of the related velocity profile in (b).

4. Further Improvement of the Approximate Solution

In a particular case of , observing the behavior of the difference error , we propose a correction formula by adding an auxiliary term to the formula as follows: where is the maximum of the absolute error at the point . Values of and are numerically evaluated as Numerical implementation for results in the errors Comparing the results with those in Table 2, one can find that the corrected approximation and its first derivative reasonably improve the accuracy of and .

For comparison with the existing approximation method, we consider the modified generalized Laguerre function Tau method introduced in the literature [27] such as based on the generalized Laguerre polynomials for and a scaling parameter . For the unknown coefficients ’s the Tau method [28, 29] is used, which generates a nonlinear system of algebraic equations. Thus a Newton-like iterative method is required to determine the coefficients ’s as a result.

Table 3 includes numerical results of the relative errors , and , with the parameters , for the Blasius function . Additionally, numerical results of and , and for the first derivative are given in Table 4. In the tables, the relative errors are defined as for an approximation to the Blasius function . Therein, and are replaced by the numerical solutions for a set of nodes . From Tables 3 and 4 we can see that the presented approximations and are less accurate than and on the region , and vice versa outside the region. However, it is also noticed that the inferiority of the presented approximations is quite overcome by the corrected approximation and .

tab3
Table 3: Relative errors for the Blasius function .
tab4
Table 4: Relative errors for the first derivative .

5. Conclusions

For the Blasius problem on the semi-infinite interval we proposed a uniformly accurate approximation formula in (11). The proposed method employs the weight function in (8) to combine a near origin series solution and a faraway reference solution.

To improve the accuracy, we introduced a correction method in (17) which includes an additional term reflecting the deviation of from the Blasius function . As a result we can observe that the presented method is available for approximation to and while the approximation to the second derivative is not so effective.

Comparing the presented solutions and with the existing solution , a solution from the generalized Laguerre spectral approach [27] based on Tau method, we summarize advantages of the presented method with discussions as follows.(i)The presented solution is composed of simple forms of known solutions, that is, a series solution and a reference solution + , while the spectral method requires solving a nonlinear system of algebraic equations. This implies that the presented method will save number of evaluations in numerical implementation.(ii)The corrected solution highly improves accuracy of with a small number of terms , and numerical results show that it is comparable to the spectral solution with .(iii)There is a room for further improvement of the present method, for example, by replacing the weight function by some more appropriate one or employing other partial solutions instead of or .

To conclude, though the presented method is limitedly applicable to the Blasius problem unlike the spectral methods, we may expect to develop an extensive method for other nonlinear differential equations motivated by the advantages above.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF-2013R1A1A4A03005079).

References

  1. J. P. Boyd, “The Blasius function: computations before computers, the value of tricks, undergraduate projects, and open research problems,” SIAM Review, vol. 50, no. 4, pp. 791–804, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  2. S. Finch, “Prandtl-Blasius flow,” 2008, http://www.people.fas.harvard.edu/~sfinch/csolve/bla.pdf.
  3. H. Schlichting, Boundary-Layer Theory, McGraw-Hill, New York, NY, USA, 7th edition, 1979. View at MathSciNet
  4. 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  5. G. Adomian, “Solution of the Thomas-Fermi equation,” Applied Mathematics Letters, vol. 11, no. 3, pp. 131–133, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  6. F. M. Allan and M. I. Syam, “On the analytic solutions of the nonhomogeneous Blasius problem,” Journal of Computational and Applied Mathematics, vol. 182, no. 2, pp. 362–371, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  7. J. Biazar, M. Gholami Porshokuhi, and B. Ghanbari, “Extracting a general iterative method from an Adomian decomposition method and comparing it to the variational iteration method,” Computers & Mathematics with Applications, vol. 59, no. 2, pp. 622–628, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  8. J. P. Boyd, “The Blasius function in the complex plane,” Experimental Mathematics, vol. 8, no. 4, pp. 381–394, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. B. K. Datta, “Analytic solution for the Blasius equation,” Indian Journal of Pure and Applied Mathematics, vol. 34, no. 2, pp. 237–240, 2003. View at Zentralblatt MATH · View at MathSciNet
  10. J.-H. He, “Approximate analytical solution of Blasius' equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 3, no. 4, pp. 260–263, 1998. View at Publisher · View at Google Scholar · View at Scopus
  11. J.-H. He, “A review on some new recently developed nonlinear analytical techniques,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 1, no. 1, pp. 51–70, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  12. J.-H. He, “A simple perturbation approach to Blasius equation,” Applied Mathematics and Computation, vol. 140, no. 2-3, pp. 217–222, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  13. J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  14. J.-H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 1-2, pp. 87–88, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  15. S.-J. Liao, “A uniformly valid analytic solution of two-dimensional viscous flow over a semi-infinite flat plate,” Journal of Fluid Mechanics, vol. 385, pp. 101–128, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  16. S.-J. Liao, “An explicit, totally analytic approximate solution for Blasius' viscous flow problems,” International Journal of Non-Linear Mechanics, vol. 34, no. 4, pp. 759–778, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  17. J. Lin, “A new approximate iteration solution of Blasius equation,” Communications in Nonlinear Science & Numerical Simulation, vol. 4, no. 2, pp. 91–94, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  18. J. Parlange, R. D. Braddock, and G. Sander, “Analytical approximations to the solution of the Blasius equation,” Acta Mechanica, vol. 38, no. 1-2, pp. 119–125, 1981. View at Publisher · View at Google Scholar
  19. Ö. Savaş, “An approximate compact analytical expression for the Blasius velocity profile,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 10, pp. 3772–3775, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  20. A.-M. Wazwaz, “The variational iteration method for solving two forms of Blasius equation on a half-infinite domain,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 485–491, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  21. B. I. Yun, “Intuitive approach to the approximate analytical solution for the Blasius problem,” Applied Mathematics and Computation, vol. 215, no. 10, pp. 3489–3494, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  22. B. I. Yun, “Approximate analytical solutions using hyperbolic functions for the generalized Blasius problem,” Abstract and Applied Analysis, vol. 2012, Article ID 581453, 10 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  23. W. Bao and J. Shen, “A generalized-Laguerre-Hermite pseudospectral method for computing symmetric and central vortex states in Bose-Einstein condensates,” Journal of Computational Physics, vol. 227, no. 23, pp. 9778–9793, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  24. J. Shen and L.-L. Wang, “Some recent advances on spectral methods for unbounded domains,” Communications in Computational Physics, vol. 5, no. 2–4, pp. 195–241, 2009. View at MathSciNet
  25. K. Parand, M. Shahini, and M. Dehghan, “Rational Legendre pseudospectral approach for solving nonlinear differential equations of Lane-Emden type,” Journal of Computational Physics, vol. 228, no. 23, pp. 8830–8840, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  26. K. Parand, M. Dehghan, A. R. Rezaei, and S. M. Ghaderi, “An approximation algorithm for the solution of the nonlinear Lane-Emden type equations arising in astrophysics using Hermite functions collocation method,” Computer Physics Communications, vol. 181, no. 6, pp. 1096–1108, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  27. K. Parand, M. Dehghan, and A. Taghavi, “Modified generalized Laguerre function tau method for solving laminar viscous flow: the Blasius equation,” International Journal of Numerical Methods for Heat & Fluid Flow, vol. 20, no. 6-7, pp. 728–743, 2010. View at Publisher · View at Google Scholar · View at MathSciNet
  28. M. Dehghan and A. Saadatmandi, “A tau method for the one-dimensional parabolic inverse problem subject to temperature overspecification,” Computers & Mathematics with Applications, vol. 52, no. 6-7, pp. 933–940, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  29. A. Saadatmandi and M. Dehghan, “Numerical solution of a mathematical model for capillary formation in tumor angiogenesis via the tau method,” Communications in Numerical Methods in Engineering with Biomedical Applications, vol. 24, no. 11, pp. 1467–1474, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet