Numerical Solution of the Blasius Viscous Flow Problem by Quartic B-Spline Method

where the prime denotes the derivatives with respect to x. In addition to the unknown function f, the solution of (1) and (2) is characterized by the value of α = f(0). Blasius [1] in 1908 found the exact solution of boundary layer equation over a flat plate. A highly accurate numerical solution of Blasius equation has been provided by Howarth [2], who obtained the initial slope α = f ex(0) = 0.332057. Liu and Chang [3] have developed a new numerical technique; they have transformed the governing equation into a nonlinear secondorder boundary value problem by a new transformation technique, and then they have solved it by the Lie group shooting method. He [4, 5] gave a solution in a family of power series with parameter p by means of the perturbation method for solving this equation. Bender et al. [6] proposed a simple approach using δ-expansion to obtain accurate totally analytical solution of viscous fluid flow over a flat plate. Aminikhah [7] used LTNHPM to obtain an analytical approximation to the solution of nonlinear Blasius viscous flow equation. Recently, the fixed point method (FPM) [8], which is based on the fixed point concept in functional analysis, is adopted to acquire the explicit approximate analytical solution to the nonlinear differential equation. Finally, efforts [9, 10] have been made to obtain the solution at the surface boundary and changed the problem from a boundary value differential equation into an initial value one. The solution in the entire domain, however, still requires computation. In this paper, the quartic B-spline approximations are employed to construct the numerical solution for solving the Blasius equation.Theunknowns are obtainedwith using optimization.Theproposedmethod is applied to the problem and the computed results are compared with those of Howarth’s method to demonstrate its efficiency.


Introduction
One of the well-known equations arising in fluid mechanics and boundary layer approach is Blasius differential equation.The classical Blasius [1] equation is a third-order nonlinear two-point boundary value problem, which describes twodimensional incompressible laminar flow over a semi-infinite flat plate at high Reynolds number, 2  () +  ()   () = 0,  ≥ 0, with boundary conditions where the prime denotes the derivatives with respect to .In addition to the unknown function , the solution of (1) and ( 2) is characterized by the value of  =   (0).Blasius [1] in 1908 found the exact solution of boundary layer equation over a flat plate.A highly accurate numerical solution of Blasius equation has been provided by Howarth [2], who obtained the initial slope  =   ex (0) = 0.332057.Liu and Chang [3] have developed a new numerical technique; they have transformed the governing equation into a nonlinear secondorder boundary value problem by a new transformation technique, and then they have solved it by the Lie group shooting method.He [4,5] gave a solution in a family of power series with parameter  by means of the perturbation method for solving this equation.Bender et al. [6] proposed a simple approach using -expansion to obtain accurate totally analytical solution of viscous fluid flow over a flat plate.Aminikhah [7] used LTNHPM to obtain an analytical approximation to the solution of nonlinear Blasius viscous flow equation.Recently, the fixed point method (FPM) [8], which is based on the fixed point concept in functional analysis, is adopted to acquire the explicit approximate analytical solution to the nonlinear differential equation.Finally, efforts [9,10] have been made to obtain the solution at the surface boundary and changed the problem from a boundary value differential equation into an initial value one.The solution in the entire domain, however, still requires computation.
In this paper, the quartic B-spline approximations are employed to construct the numerical solution for solving the Blasius equation.The unknowns are obtained with using optimization.The proposed method is applied to the problem and the computed results are compared with those of Howarth's method to demonstrate its efficiency.
Let () be the quartic B-spline function at the nodal points.Then approximate solution of (1) can be written as [13] where   () are the quartic B-spline functions and   are the unknown coefficients.Each B-spline covers the five elements so that an element is covered by five B-splines.The values of   () and its derivatives are tabulated in Table 1.
In order to solve system (8), we direct attention to that () and its derivatives satisfied in (1) and ( 2) and initial slope   (0) = 0.332057 approximately, and then we have Therefore, the error vector  of the approximation can be written as Now, we wish to minimize the error norm, the  2 -norm, such that By solving (11), we can get the values of   for  = −2, −1, . . .,  + 1.Now, with substituting the values  −2 ,  −1 ,  0 , . . .,  +1 into (7), the approximate value of () and its derivatives will be ensured.

Numerical Results
In this section, approximation values of (),   (), and   () based on proposed methods for some values of   are presented.These values for ℎ = 0.2, 0.1, 0.01 have been calculated and these results are compared with the Howarth results.Tables 2, 3, and 4 are made to compare between the present results and results given by Howarth for approximation values of (),   (), and   (), respectively.In Figures 1, 2, and 3, one can also see the comparison between our results and Howarth's results.

Conclusion
In this survey, the quartic B-spline approximations are used to solve the Blasius equation.This method led to a system of nonlinear equations.The unknowns are obtained by minimizing the error norm.The computed results are compared with those of DTM, LTNHPM, and Howarth's methods to demonstrate the validity and applicability of the technique.This method is simple in applicability and the results show that the solutions will become more accurate with reducing step size.The computations associated with the examples in this paper were performed using MATLAB R2015a.

2 International
, where ℎ =  +1 −   ,  = 0, 1, . . .,  − 1.The quartic B-splines are defined upon Journal of Engineering Mathematics an increasing set of  + 1 knots over the problem domain plus 8 additional knots outside the problem domain; the 8 additional knots are positioned as

Figure 1 :
Figure 1: The comparison of () between our results and Howarth's results with ℎ = 0.01.

Figure 2 :
Figure 2: The comparison of   () between our results and Howarth's results with ℎ = 0.01.

Figure 3 :
Figure 3: The comparison of   () between our results and Howarth's results with ℎ = 0.01.

Table 1 :
Coefficient of quartic B-splines and its derivatives at knots   .

Table 2 :
Comparison of () between some numerical results and our results using different values of ℎ.

Table 3 :
Comparison of   () between some numerical results and our results using different values of ℎ.

Table 4 :
Comparison of   () between some numerical results and our results using different values of ℎ.