Numerical Solution of Blasius Equation through Neural Networks Algorithm

In this paper mathematical techniques have been used for the solution of Blasius differential equation. The method uses optimized artificial neural networks approximation with Sequential Quadratic Programming algorithm and hybrid AST-INP techniques. Numerical treatment of this problem reported in the literature is based on Shooting and Finite Differences Method, while our mathematical approach is very simple. Numerical testing showed that solutions obtained by using the proposed methods are better in accuracy than those reported in literature. Statistical analysis provided the convergence of the proposed model.


Introduction
Blasius differential equation is the mother of all boundary layer equations in fluid mechanics.It is almost a hundred-year-old differential equation and still an active topic among the all time researchers.Blasius derived the famous Blasius equation by using transform technique.Equation was discussed in many articles by analytical and numerical ways.Many scientists investigated interesting results of Blasius equation either from mathematical point of view or from engineering prospective.Howart numerically solved this differential equation and com-pared the results with [1].Further, analytical solutions which are uniformly valid over the whole domain do not exist till 1999.Liao in his best paper gave the result by Homotopy Analysis Method (HPM) [2] [3].G.I. Shishkin [4] showed asymptotic behavior of differential and difference solutions to get different scheme with a finite number of nodes for enough long interval.The standard non-homogenous Blasius equation is The Blasius equation describes the velocity profile of fluid in a boundary layer.It is a basic equation in the fluid mechanics which appears in the study of flow of an incompressible viscous fluid over a semi-infinite plane.Blasius equation is basically derived from classical Navier Stock equation [5]- [10].This well-known equation is investigated by many researchers to find its solution.There are many analytical and numerical methods, like Decomposition Method (DM), Iteration Method (IM), Homotopy Analysis Method (HAM) and Parameter Iteration Method (PIM), applied for solutions of differential equations.
We use Sequential Quadratic Programming (SQP) and Active set-Interior point technique (AST-INP) which is hybrid technique as optimization tools in MATLAB to solve the Blasius differential equation.In order to increase the accuracy, we repeated this proposed algorithm for several time for random selection of variable values with a moderate number of weights or variables.Moreover, a detailed statistical analysis was provided in both cases for the validity of this proposed method.In this paper we presented the numerical analysis of said equation on the bases of experiments and showed that the present solutions are highly accurate as compared to other methods.
The rest of the paper is organized as follows.In Section 2, importance of artificial neural networks and application is presented.In Section 3, we formulate the Blasius problem and propose a mathematical model for the numerical treatment of it with the help of activation function called Log-Sigmoid basis on logarithmic function.Numerical results and their graphical details are presented in next section.A brief statistical analysis and summery is presented in last section.

Artificial Neural Networks
Artificial intelligence technique is suitable in order to solve different types of differential equation.Lee and kang work with parallel processor computers to solve ordering differential equation by using Hopfield neural network model.Meads and Fernandez used B1 splines and feed forward neural network architecture to solve nonlinear and linear ordinary differential equation.Hybrid artificial neural network-nelder-mead method used by malek and shekari to solve higher order linear differential equation.A new bilaterally approach was introduced to find upper and lower bound of blasius equation by Lee.But bilaterally approach could not satisfy the boundary condition.But present study introduced the method to solve blasius differential equation which satisfies the boundary condition.
Computational models of biological brain are the example of neural networks.As the brain works, neural networks comprises large number of interconnected neurons.Each neuron has the ability to perform simple computation.As compared to biological neuron, an artificial neuron is much simpler.The construction of Artificial neural networks (ANN) is hidden in one or more layers where the factual processing is performed through weighted connections.Each neurons in the covered layer have connections to all neurons in uncovered (output) layer as shown in Figure 1.Application of such system is wide range.An artificial neural networks can learn to perform complex tasks like system identification, function approximation, trend prediction, pattern recognition and process control.Neurons in the input region only behave as buffer for dividing the input signals to neurons in the covered region.

Mathematical Formulation
The Blasius equation is given by ( ) with boundary conditions ( ) ( ) Further Equation (1) satisfies the asymptotic condition ( ) We construct a mathematical model based on active set with fitness function.The system is based on numerical computation through which we obtain optimum variables of the proposed model [11]- [13].
The solution ( ) u η of the differential Equation ( 5) along with it's derivatives, can be approximated by the following continuous relations as in neural network methodology and the activation function Log-Sigmoid is defined as ( ) Blasius equation mathematical model with the help of above activation function was developed to approximate the solution of Equations (5-8) with 2 nd and 3 rd order derivatives.Using this methodology the solution ( ) u η can be approximated by ( ) where m is the number of neurons, δ , ω , and β are real-valued bounded adaptive parameters or weights can be expressed an array W as ( ) , , , , , , , , , , , The second and third order derivatives of solution ( ) u η can be approximation by the continuous relations where ( ) The mathematical model for Equation ( 5) can be formulated by a linear combination of networks Equations (10-12) is called a differential equation neural networks.
The fitness function for proposed model ∈ has been formulated for the Equations (5-8) using Mathematical model by defining the error as the sum of mean squared errors:
where ( ) And 2 ∈ for initial values can be defined as For 1

Optimization Procedure for Numerical Solution
Furthermore, we provide some detail about the procedural steps for the optimization in MATLAB built-in function is given below.The generic flow diagram of the overall process is shown in Figure 2.
Step 1: Initialization: A vector with randomly generated bounded real values of length equal to the number of weights in each Mathematical model acts as the starting point for each solver: ( ) , , , , , , , , Here m represents the number of neurons.
Step 2: Fitness Evaluation: The MATLAB built-in function for constrained optimization problems is invoked for each model.
Step 3: Termination Criteria: Terminate the execution of the solver, if any of the following criteria is satisfied: • required level of predefined fitness achieved, i.e., 14 10 − ∈≤ .• total number of iterations executed, as listed in Table 1.
Step 4: Storage: Save the final optimal weights (variables) along with fitness values and computational time taken by the algorithm.
Step 5: Statistical Analysis: Repeat steps 1 to 4 for sufficiently large number of times to perform an effective and reliable statistical analysis.

Numerical Results
In this section, we show the output of the proposed method by the numerical results of Blasius equation.Further,  with hybrid AST-INP and SQP algorithms.Furthermore, we provided the statistical analysis with several time run of optimization tools for Mean, Median, STD and Variance tabulated in Table 2 and Table 3 for both cases.Which showed that the present solution is highly accurate as compared to others methods present in literature.

Statistical Analysis and Discussion
On the basis of the simulations and results obtained in the previous section, it can be concluded that Blasius differential equation can be solved by stochastic computational intelligence technique, like SQP optimization algorithm, AST-INP hybrid technique, supported with simulating annealing.The differential equation neural networks trained by SQP algorithm and AST-INP are better stochastic optimizers as compared to other algorithms.The statistical analysis for this case with AST-INP and SQP algorithm are tabulated in Table 2 and Table 3.These results showed the better accuracy of numerical data with reference solution.We presented minimum value, mean, median and STD for the accuracy of our solver for 300 time multi-runs with minimum time as shown in Table 2 and Table 3.Thus it can be stated that proposed computing approach is reliable, effective and easily applicable for complex differential equation of order three.In our future work, we intend to use other hybrid computational intelligence algorithms like GA-SQP, GA-INP and GA-AST to solve these problems using Bessel's polynomial as active function.

= 1 u
; with the initial and boundary condition ( ) ∞ = , where ( ) u x is dimension less stream function and is the similarity ordinate.

2 u and ( ) 3 u
represented 2nd and 3rd derivative with respect to η respectively.

Figure 2 .
Figure 2. Flow chart of Blasius equation model.

Furthermore, Figure 3 andFigure 3 .
Figure 3.A 3D view of neural network model strained with SQP for Blasius equation.

Figure 4 .
Figure 4.A 2D view of a set of variables of neural network model strained with SQP for Blasius equation.

Figure 5 .
Figure 5.A 2D view of neural network model strained with hybrid (AST-INP) for Blasius equation.

3 -
dimensional view of weights or variables obtained from AST-INP optimizer.In Figure 7, it is also shown that the confidence level of absolute error at 20 10 − is 90 percent and maximum values lies between 20 10 − and 18 10 − .We have shown the Chi-squares distribution fitness with normal form to the numerical data through several time simulation as shown in Figure 8. Furthermore, in Figure 9, we presented the comparison of reported results mean with exact solution (reference solution).

Figure 6 .
Figure 6.A 3D view of set of variables of neural network model strained with hybrid (AST-INP) for Blasius equation.

Figure 7 .
Figure 7. Chi Square curve of 300 results shows the confidence level of our data.

Figure 8 .
Figure 8. Fitness of Chi-square distribution with proposed data of Blasius equation.

Figure 9 .
Figure 9. Mean VS Exact value of Blasius equation.

Table 1 .
Parameter settings for the function "fmincon" in MATLAB simulations.

Table 2 .
Statistical analysis of solution of Blasius equation with hybrid (AST-INP).

Table 3 .
Statistical analysis of solution of Blasius equation with SQP.