A Lie group integrator to solve the hydromagnetic stagnation point flow of a second grade fluid over a stretching sheet

1 Department of Mathematics, Basic Science Faculty, University of Bonab, P.O. Box 55513-95133, Bonab, Iran 2 Faculty of Engineering Technology, Amol University of Special Modern Technological, Amol, Iran 3 Department of Mathematics, College of Sciences, Jazan University, Jazan 45142, Saudi Arabia 4 Department of Basic Sciences, University of Engineering and Technology, Peshawar, Pakistan 5 Section of Mathematics, International Telematic University Uninettuno, Corso Vittorio Emanuele II, 39,00186 Roma, Italy


Introduction
Stagnation point flow reveals in practically all flow branches of engineering and science. Free stagnation point or a line exists interior of the fluid domain while in some situations flow is stagnated by a solid wall.
Hiemenz [1] demonstrated that Navier-Stokes equations are able to describe the stagnation point flow. Plane stagnation point flow which is called the Hiemenz flow, firstly discussed by Hiemenz as follows: f ′′′ + f f ′′ − ( f ′ ) 2 + 1 = 0, (1.1) with the boundaries: Gersten et al. [2] generalized these results to the stagnation point flow in three dimensions. Homotopy analysis method is used by Liao [3] to the 2D laminar viscous flow over a semi-infinite plate. Existence and uniqueness consequences for nonlinear third-order ODEs appeared in the stagnation point flow of a hydromagnetic fluid are discussed by Van Gorder et al. [4]. The flow characteristics in an electrically conducting second grade fluid pervaded by a uniform magnetic field over a stretching sheet is discussed by Vajravelu et al. [5]. Steady 2D stagnation-point flow of an incompressible viscous electrically conducting fluid over a flat deformable sheet is studied by Mahapatra et al. [6]. Van Gorder et al. [7] established the existence and uniqueness consequences for fourth-order ODEs over the semi-infinite interval, arising in the hydromagnetic stagnation point flow of a second grade fluid over a stretching sheet.
The concept of Lie groups plays an significant rule in establishing some intense numerical approaches to integrate ODEs, which conserve the geometrical properties. Preserving the geometry construction under discretization, is essential in the improvement of a qualitatively true performance in the minimization of approximation errors. So, along with the geometric structure and invariance of the ODEs, novel approaches may be invented, which are more stable and effective than the common approximation techniques.
A novel geometric method based on GL(6, R) is introduced for the magneto-hemodynamic flow in a semi-porous channel in [8]. Moreover, in [9], a nonlinear heat transfer equation is considered by S L(2, R). Order of differential equations and different initial or boundary conditions conclude different geometric algorithms for solving differential equations. That is, numerical algorithms based on e.g. GL(6, R) and GL(4, R) is completely different. Combination of radial basis functions with GPS are discussed in [10]. Some other recent papers of geometric numerical approach, Lie group analysis, and exact solution can be found in .
This paper is organized as follows. In the further section, we will touch only a few aspects of the Mathematical formulation for the hydromagnetic stagnation point flow. In Section 3, we derive an interesting algorithm to integrate the underlying nonlinear equation. Section 4 deals with the numerical results. Final remarks are provided in the conclusion section.

Mathematical formulation
Let us investigate the 2D steady stagnation-point flow of an electrically conducting fluid in the existence of a uniform transverse magnetic field towards a flat surface coinciding with the plane y = 0, the flow being enclosed to the region y > 0. Two equivalent and converse forces are considered in direction of the x−axis such that the wall is stretched keeping the origin fixed, and a uniform magnetic field B 0 is forced along the y−axis (Figure 1).

Figure 1.
A profile of the physical problem.
In the usual notation, the MHD equations for steady 2D stagnation-point flow in the boundary layer along the stretching surface are [6,34]: where ρ, σ, B 0 are density, electrical conductivity and magnetic field respectively. Also, (u, v) are the fluid velocity components toward x and y. U(x) stands for the stagnation-point velocity in the inviscid free stream. For this problem, the boundaries are Substituting (2.5) into (2.2), concludes The self-similar solutions will normally be made to meet Eq (2.6) with additional nonlinear terms in the governing differential equation, where K ≥ 0 is the viscoelastic parameter. By the additional condition f ′′ (∞) = 0, we investigate the following fourth order differential equation: where M = σ 0 B 2 0 /ρB ≥ 0 is the magnetic parameter, K = ωB/ν ≥ 0 is the viscoelastic parameter, s ≥ 0, C is the parameter due to the stagnation point flow, and χ is the stretching of the sheet.
† gl(4, R) is the Lie algebra associated with the Lie group GL(4, R). where Therefore, each of this Lie group elements may be erected by a single parameter G(r), with r ∈ [0, 1].
From (3.18) and (3.22) we have: Thus, for a given r, we may get the values y f 1 , y 0 3 , y 0 4 and y f 4 by the further explained process. Imposing the initial guesses of y 0 3 , y 0 4 , y f 1 and y f 4 , and the given boundary values of y 0 1 = s, y 0 2 = χ, y f 2 = C and y 0 3 = 0, with a specified r ∈ [0, 1], and from For a fixed value of r, we are able to compute y 0 3 and y 0 4 from Eq (3.25) and then approximately integrate Eq (3.7) by GPS from 0 to η ∞ , and compare the end value of y f 2 and y f 3 with the exact ones y 2 (η ∞ ) = C and y 3 (η ∞ ) = 0 respectively. That is, we have to find the minimum value of In the appendix, we give a brief algorithm of our proposed method to find the initial values of problem.
In the non-geometric numerical methods, investigation of the method stability is an important task which guarantees the validity of techniques. Without considering the convergence and error analysis, results may lead to solutions with no physical meaning or high-level of errors. This issue may arise from that, we do not know information about the geometric structure of original problem. However, in this paper we introduce a GL(4, R) method which conserve the geometrical properties of problem. Preserving the geometry construction under discretization, is essential in the improvement of a qualitatively true performance in the minimization of approximation errors.

Numerical results
We emphasize that the stepping approach developed for integrating the IVP needs the initial values (0) and f 4 (0) = f ′′′ (0) for the fourth-order ODEs. Both of the last initial values in Eq (2.8) are unknown. Therefore, approximating f ′′ (0) and f ′′′ (0) converts Eq (2.8) into an IVP. It is worth pointing out that every IVP can be integrated by GPS.
Indeed, our proposed method in this paper is split into two parts. In the first part, initial conditions of the problem are obtain by the GL(4, R) method. Our method is implemented by the convergence criterion of (3.26) with ϵ = 10 −8 . For this value of convergence criterion, only four or five steps of algorithm, mentioned in the Appendix, are needed. So, very fast convergence occur in this part of the proposed method Second part of this technique starts when we obtain the initial values of BVP from the first part. Second part which is nominated as the GPS, is an explicit iteration method which is wellknown as the weighted Euler method. Convergence of this part can be guaranteed from the geometric structure of method.
The results of this paper were announced with η ∞ = 7 and ∆ξ = 0.01 for the GPS procedure. There is no loss of generality in assuming y 0 3 = y 0 4 = y f 1 = y f 4 = 1 as initial guesses. Some mis-matching errors with various M, C, s, K and χ are plotted in Figure 2 It is observable that the profiles of f (ξ) and f ′ (ξ) decrease, when the magnetic parameter increases, while the reverse situation occur for f ′′ (ξ) and it decrease more rapidly toward zero as ξ tends to η ∞ .          The double mesh principle is applied to estimate the error e N (ξ) = | f N (ξ) − f 2N (ξ)|, (4.1) where N is the number of GPS iterations. Residual errors of the present method applied for Eq (2.8) with N = 10000 are plotted in Figure 8 for various parameters of equation. This figure shows that our method is reliable for current types of equations. Besides, obtained solutions and results are in good agreement with the results of [7].

Conclusions
In this paper, a novel iterative technique based on the Lie group of GL(4, R), as the first time, is proposed to solve the hydromagnetic stagnation point flow of a second grade fluid over a stretching sheet. Structures of the Lie group elements in GL(4, R) and corresponding Lie algebra elements in gl(4, R) are completely innovative for the underlying problem in this paper. Against the non-geometric numerical techniques, the proposed method uses the geometric structure of original problem and this guarantees that obtained results are reliable and the proposed method avoid the ghost solutions. Effects of equation parameters are graphically considered. Point-wise errors for various parameters show the power and effectiveness of proposed method.