# Unsteady Reversed Stagnation-Point Flow over a Flat Plate

Authors: Vai Kuong Sin; Chon Kit Chio
Journal: ISRN Applied Mathematics
 (2012)

---

## Abstract

This paper investigates the nature of the development of two-dimensional laminar flow of an incompressible fluid at the reversed stagnation-point. Proudman and Johnson (1962) first studied the flow and obtained an asymptotic solution by neglecting the viscous terms. Robins and Howarth (1972) stated that this is not true in neglecting the viscous terms within the total flow field. Viscous terms in this analysis are now included, and a similarity solution of two-dimensional reversed stagnation-point flow is investigated by solving the full Navier-Stokes equations.


---

## Body

## 1. Introduction

The full Navier-Stokes equations are difficult or impossible to obtain an exact solution in almost every real situation because of the analytic difficulties associated with the nonlinearity due to convective acceleration. The existence of exact solutions is fundamental in their own right not only as solutions of particular flows, but also as agreeable in accuracy checks for numerical solutions.In some simplified cases, such as a fluid traveling through a rigid body (e.g., missile, sports ball, automobile, spaceflight vehicle), or in oil recovery industry, crude oil that can be extracted from an oil field is achieved by gas injection, as shown in Figure1 or, equivalently, an external flow impinges on a stationary point called stagnation-point that is on the surface of a submerged body in a flow, of which the velocity at the surface of the submerged object is zero. A stagnation-point flow develops and the flow in the vicinity of this stagnation point is governed by Navier-Stokes equations. The classic problems of two-dimensional stagnation-point flows can be analyzed exactly by Hiemenz [1]. The result is an exact solution for flow directed perpendicular to an infinite flat plate. Howarth [2] and Davey [3] extended the two-dimensional and axisymmetric flows to three dimensions, which are exact solutions of the Navier-Stokes equations, and Wang [4] studied the case for obliquely-impacting jets.Figure 1
Oil recovery industry.On the contrary, a solution against an infinite flat plane does not exist if the potential flow at infinity is reversed. The flow in the vicinity of this reversed stagnation point is governed by boundary-layer separation and vorticity generation and the reversed stagnation-point flow develops. Reversed stagnation-point flow is a flow in which the component of velocity normal to a wall is outward the wall everywhere in the region concerned, so that the vorticity created at the wall will be convected toward the wall, in opposition to viscous diffusion away from it.Proudman and Johnson [5] suggested that the convection terms dominate in considering the inviscid equation in the body of the fluid. By introducing a very simple function of a particular similarity variable and neglecting the viscous forces in their analytic result for region sufficient far from the wall, they obtained an asymptotic solution in reversed stagnation-point flow, describing the development of the region of separated flow for large time t. Robins and Howarth [6] have recently extended the asymptotic solution, finding the higher order terms by singular perturbation methods. They indicated that the viscous forces cannot be ignored in the governing equation because of a consistent asymptotic expansion in both this outer inviscid region and also in the inner region near the plane. Smith [7] generalized the solution of Proudman and Johnson with both viscous and convection terms in balance by considering the monotonic potential flow when the time is relatively large. Shapiro [8] obtained a solution for unsteady backward stagnation-point flow with injection or suction.These unsteady flows fit within a class of similarity transformations originally identified by Birkhoff using a group-theoretic approach.Numerical simulation of reversed stagnation-point flow with full Navier-Stokes equations has been studied in [9]. In the present study, the unsteady reversed stagnation-point flow is investigated. The flow is started impulsively in motion with a constant velocity away from near the stagnation point. A similarity solution of full Navier-Stokes equations is solved by applying numerical method.


## 2. Flow Analysis Model

The viscous fluid flows in a rectangular Cartesian coordinates(x,y,z), Figure 2, illustrating the motion of external flow directly moving perpendicular out of an infinite flat plane wall. The origin is the so-called stagnation point and z is the normal to the plane.Figure 2
Coordinate system of reversed stagnation-point flow.By conservation of mass principle with constant physical properties, the equation of continuity is(2.1)∂u∂x+∂v∂y=0.
We consider the two-dimensional reversed stagnation-point flow in unsteady state and the flow is bounded by an infinite plane y=0, the fluid remains at rest when time t<0. At t=0, it starts impulsively in motion which is determined by the stream function
(2.2)ψ=-αxy.At large distances far above the planar boundary, the existence of the potential flow implies an inviscid boundary condition. It is given by(2.3)u=-αx.v=V0,
where u and v are the components of flow velocity, A is a constant proportional to V0/L, V0 is the external flow velocity removing from the plane, and L is the characteristic length. We have u=0 at x=0 and v=0 at y=0, but the no-slip boundary at wall (y=0) cannot be satisfied.Since for a viscous fluid the flow motion is determined by only two factors, the kinematic viscositiesν and α, we consider the following modified stream function:
(2.4)ψ=-Aνxf(η,τ),η=Aνy,τ=At,
where η is the nondimensional distance from wall and τ is the nondimensional time. Note that the stream function automatically satisfies equation of continuity (2.1). The Navier-Stokes equations [10] governing the unsteady flow with constant physical properties are
(2.5)∂u∂t+u∂u∂x+v∂u∂y=-1ρ∂p∂x+ν(∂2u∂x2+∂2u∂y2),∂v∂t+u∂v∂x+v∂v∂y=-1ρ∂p∂y+ν(∂2v∂x2+∂2v∂y2),
where u and v are the velocity components along x and y axes, and ρ is the density.Substitutingu and v into the governing equations results in a simplified partial differential equation. From the definition of the stream function, we have
(2.6)u=∂ψ∂y=-Axfη,v=-∂ψ∂x=Aνf.The governing equations can be simplified by a similarity transformation when several independent variables appear in specific combinations, in flow geometries involving infinite or semi-infinite surfaces. This leads to rescaling, or the introduction of dimensionless variables, converting the original system of partial differential equations into the following pair of partial differential equations:(2.7a)-A2xfητ+A2x(fη)2-A2xffηη=-1ρ∂p∂x-A2xfηηη,(2.7b)AAνfτ+AAνffη=-1ρ∂p∂y+AAνfηη.The pressure gradient can be again reduced by a further differentiation equation (2.7b) with respect to x. That is
(2.8)∂2p∂x∂y=0
and (2.7a) reduces to
(2.9)[fητ-(fη)2+ffηη-fηηη]η=0.
The initial and boundary conditions are
(2.10)f(η,0)≡η(η≠0),f(0,τ)=fη(0,τ)=0(t≠0),f(∞,τ)~η.
The last condition reduces the above differential equation (2.9) to the form
(2.11)fητ-(fη)2+ffηη-fηηη=-1,
with the boundary conditions
(2.12)f(0,τ)=fη(0,τ)=0,fη(∞,τ)=1.
Equation (2.11) is the similarity equation of the full Navier-Stokes equations at a two-dimensional reversed stagnation point. The coordinates x and y are replaced by a dimensionless variable η. Under the boundary conditions fη(∞,τ)=1, when the flow is in steady state such that fητ≡0, the differential equation has no solution.


## 3. Similarity Analysis

### 3.1. Asymptotic Solution

Whenτ is relatively small, Proudman and Johnson [5] first considered the early stages of the diffusion of the initial vortex sheet at y=0. They suggested that when the flow is near the wall region, the viscous forces are dominant, and the viscous term in the governing Navier-Stokes equations is important only near the boundary.On the contrary, the viscous forces were neglected away from the wall. The convection terms dominate the motion of external flow in considering the inviscid equation in the fluid. They considered the similarity of the inviscid equation(3.1)fητ-(fη)2+ffηη+1=0.
Proudman and Johnson obtained that a similarity solution of (3.1) is in the form
(3.2)f(η,τ)=eτF(γ),(3.3)F(γ)=γ-2c(1-e-cγ),
where γ=ηe-τ and c is a constant of integration. Robins and Howarth [6] estimated the value of c to be 3.51. This solution describes the flow in the outer region, moving away from the plane with a constant velocity. It can be checked that the viscous term fηηη is still small compared to the convective terms, so that their assumption of neglecting the viscous term is still valid.In the inner region, the viscous term cannot be neglected and the no-slip condition must be satisfied on the wall. Whenτ→∞ and η/eτ is relatively small, the solution (3.3) yields
(3.4)F~-γ=-ηe-τ,f=-η,f′=-1.
Substituting in (2.11) yields
(3.5)f′′′-ff′′+(f′)2-1=0,f(0)=f′(0)=0,f′(∞)=-1.This is exactly the classic stagnation-point problem (Hiemenz [1]) by changing the sign in f. It is a third-order nonlinear ordinary differential equation and does not have an analytic solution, and thus it is necessary to solve it numerically. The general features of the predicted streamline are sketched in Figure 3.Figure 3
Streamlines of reversed stagnation-point flow.Although an asymptotic solution was obtained, it can easily been observed that this is not valid when the viscous termfηηη is neglected within the total flow field. Robins and Howarth [6] indicated that there is a consistent asymptotic expansion in both outer inviscid region and also in the inner region that must exist close to the wall where the viscous forces need to be included. No exact solutions in both outer and inner regions were discovered.


### 3.2. Particular Solution

In our two-dimensional model, the fluid remains at rest when timet<0 and is set in motion at t>0 such that at large distances far above the planar boundary the potential flow is a constant V0 for all value of t. Both Proudman and Johnson [5], and Robins and Howarth [6] have set V0=1 and the corresponding boundary condition fη(∞,τ)=1.When the flow is in steady state such thatfητ≡0, it was proven that the similarity velocity fη(η) cannot ultimately approach to 1. The differential equation has no solution. If the potential flow V0 is restricted not to be a constant, the boundary condition fη(∞,τ) results in a time dependent function and then we obtain another approach of similarity solution in reversed stagnation-point flow. As with the governing equation of reversed stagnation-point flow, we can write the stream function as
(3.6)ψ=-Aνxf(η,τ),η=Aνy,τ=At,
where A is a constant proportional to V0(τ)/L, V0(τ) is the external flow velocity removing from the plane, and L is the characteristic length. These result in the governing equation (2.9)
(3.7)[fητ-(fη)2+ffηη-fηηη]η=0.
Integrating with respect to η, we have
(3.8)fητ-(fη)2+ffηη-fηηη=-C(τ).
Under the boundary conditions fη(∞,τ)=1, the value of C(τ) should be a constant and equal to 1. If the boundary condition fη(∞,τ) is restricted not to be a constant, a particular time-dependence function C(τ) may be expressed in the form
(3.9)C(τ)=cτ2,
where c is an arbitrary constant. The partial differential equation can be simplified by a similarity transformation when a new similarity variable is introduced. This converts the original partial differential equation into an ordinary differential equation. For a time dependent function, we introduce the diffusion variable transformation [7]
(3.10)ς=η1τ,f(η,τ)=1τF(ς).
Here ς is the time combined nondimensional variable and F(ς) are the nondimensional velocity functions. Substitution of the similarity transformation yields an ordinary differential equation
(3.11)-12ςF′′-F′-F′2+FF′′-F′′′=-c,
where the prime denotes the derivative with respect to the variable ς.Equation (3.11) is a third-order nonlinear ordinary differential equation and a key step in obtaining an analytical solution is to reversedrange the equation in an autonomous differential equation. In mathematics, an autonomous differential equation is a system of ordinary differential equations which does not explicitly depend on the independent variable.In order to omitting the variableς in the differential equation, it is recognized a change of variable
(3.12)Q=F-12ς
and the equation becomes an autonomous differential equation
(3.13)QQ′′-2Q′-Q′2-Q′′′=-c+34.
In our analysis, P=Q′ is the dependent variable and Q is the independent variable. Equation (3.13) is reversed ranged as
(3.14)QP′-2P-P2-P′′=-c+34
and the chain rule reduces equation (3.14) to a second-order ordinary differential equation
(3.15)QPdPdQ-2P-P2-PddQ(PdPdQ)=-c+34.
Equation (3.15) is analytically solvable that the solution might be expressed as a low order polynomial. It is suggested that
(3.16)P=a+bQ+dQ2,
and substituting into equation (3.14) and comparing the coefficients in the powers of Q results in a system of linear algebraic equations
(3.17)2a2d+ab2+2a-a2=-c+34,8abd+b3+2b+ab=0,8ad2+(7b2+2)d=0,12bd2-bd=0,6d3-d2=0.
Solving the related algebraic equation, we have
(3.18)a=-32,b=0,c=34,d=16.
Substituting the constant into (3.16) yields the first-order differential equation
(3.19)Q′=-32+16Q2.
Equation (3.19) is Riccati equation, which is any ordinary differential equation that is quadratic in the unknown function. The standard form of Riccati equation is
(3.20)Q′=RQ2+SQ+T.
The solution of Riccati equation can be obtained by a change of dependent variable, where the dependent variable y is changed to q by [11]
(3.21)Q=-q′q1R.
By identifying R=1/6, S=0, and T=-3/2, the change of variables in (3.19) becomes
(3.22)Q=-q′(1/6)q=-6q′q
so (3.19) becomes an second-order linear differential equation
(3.23)q′′-14=0.
The general solution to this equation is
(3.24)q=Acosh⁡ς2+Bsinh⁡ς2,
where A and B are arbitrary constants. Applying this solution in (3.12) leads to the general solution of (3.11)
(3.25)F(ς)=ς2-3Asinh⁡ς/2+3Bcosh⁡ς/2Acosh⁡ς/2+Bsinh⁡ς/2.
Application of the impermeability condition F(0)=0 leads to the determination of the constant B=0, so the exact solution becomes
(3.26)F(ς)=ς2-3tanhς2.
Collecting results, the velocity functions become
(3.27)f(η,τ)=1τ(ς2-3tanhς2),
where ς=(A/ντ)y is the nondimensional distance from the plate. In view of (3.27), the flow far from the boundary (x,ς→∞) becomes
(3.28)f(η,τ)⟶1τ(ς2-3).We obtain a particular solution of the unsteady reversed stagnation-point flow. The above solution is obtained in the similarity framework for unsteady viscous flows. The appearance of this positive factor in the first terms of (3.28) shows that this remote flow is directed toward the axis of symmetry and away from the plate. The second term in (3.28) describes a uniform velocity directed toward the plate. An adverse pressure gradient near the wall region leads to boundary-layer separation and associated flow reversal.The particular solution is noteworthy in that it is completely analytical, but it is limited to the region far away from the plate in the presence of nonzero termF′(0)=-1. Another no-slip boundary condition F′(0)=0 is not satisfied completely near the wall region.


### 3.3. Numerical Solution

Since the analytical solution does not satisfy the no-slip conditionF′(0)=0, it is convenient to solve the similarity equation numerically. The similarity equation and the relevant boundary conditions are
(3.29)-12ςF′′-F′-F′2+FF′′-F′′′=-c,F(0)=F′(0)=0,F′(∞)=12,
where c=3/4 in order to satisfy the unsteady viscous flows in the outer region.This equation is a third-order nonlinear ordinary differential equation. It is convenient when solving an ODE system numerically to describe the problem in terms of a system of first-order equations in MATLAB[12].For example when solving annth-order problem numerically is a common practice to reduce the equation to a system of n first-order equations. Then, by defining y1=F, y2=F′, y3=F′′, the ODE reduces to the form
(3.30)dydx=[y2y3c-12ςy3-y2-y22+y1*y3].The first task is to reduce the equation above to a system of first-order equations and define in MATLAB a function to return these. Later, we need to change the boundary value into initial value, becauseode45, anodesolver in MATLAB, can only solve the initial value problem. From (3.29), we gauss the value of F′′(0)=-1.54306 such that F′(∞)=1/2.The numerical solution for two-dimensional stagnation-point flow is shown in Figure4.Figure 4
Numerical solutions of reversed stagnation-point flow.This solution is a similarity solution of the reversed stagnation-point flow over a flat plate. It describes an unsteady viscous flow in both outer and inner regions. A single dividing streamline plane separates streamlines approaching the plate from external flow streamlines. Though the viscous termF′′′ is still small compared to the convective terms in the outer region, this is not true in neglecting the viscous terms within the total flow field. The similarity velocity field is shown in Figure 5.Figure 5
Similarity velocity field as a function ofς.


## 3.1. Asymptotic Solution

Whenτ is relatively small, Proudman and Johnson [5] first considered the early stages of the diffusion of the initial vortex sheet at y=0. They suggested that when the flow is near the wall region, the viscous forces are dominant, and the viscous term in the governing Navier-Stokes equations is important only near the boundary.On the contrary, the viscous forces were neglected away from the wall. The convection terms dominate the motion of external flow in considering the inviscid equation in the fluid. They considered the similarity of the inviscid equation(3.1)fητ-(fη)2+ffηη+1=0.
Proudman and Johnson obtained that a similarity solution of (3.1) is in the form
(3.2)f(η,τ)=eτF(γ),(3.3)F(γ)=γ-2c(1-e-cγ),
where γ=ηe-τ and c is a constant of integration. Robins and Howarth [6] estimated the value of c to be 3.51. This solution describes the flow in the outer region, moving away from the plane with a constant velocity. It can be checked that the viscous term fηηη is still small compared to the convective terms, so that their assumption of neglecting the viscous term is still valid.In the inner region, the viscous term cannot be neglected and the no-slip condition must be satisfied on the wall. Whenτ→∞ and η/eτ is relatively small, the solution (3.3) yields
(3.4)F~-γ=-ηe-τ,f=-η,f′=-1.
Substituting in (2.11) yields
(3.5)f′′′-ff′′+(f′)2-1=0,f(0)=f′(0)=0,f′(∞)=-1.This is exactly the classic stagnation-point problem (Hiemenz [1]) by changing the sign in f. It is a third-order nonlinear ordinary differential equation and does not have an analytic solution, and thus it is necessary to solve it numerically. The general features of the predicted streamline are sketched in Figure 3.Figure 3
Streamlines of reversed stagnation-point flow.Although an asymptotic solution was obtained, it can easily been observed that this is not valid when the viscous termfηηη is neglected within the total flow field. Robins and Howarth [6] indicated that there is a consistent asymptotic expansion in both outer inviscid region and also in the inner region that must exist close to the wall where the viscous forces need to be included. No exact solutions in both outer and inner regions were discovered.


## 3.2. Particular Solution

In our two-dimensional model, the fluid remains at rest when timet<0 and is set in motion at t>0 such that at large distances far above the planar boundary the potential flow is a constant V0 for all value of t. Both Proudman and Johnson [5], and Robins and Howarth [6] have set V0=1 and the corresponding boundary condition fη(∞,τ)=1.When the flow is in steady state such thatfητ≡0, it was proven that the similarity velocity fη(η) cannot ultimately approach to 1. The differential equation has no solution. If the potential flow V0 is restricted not to be a constant, the boundary condition fη(∞,τ) results in a time dependent function and then we obtain another approach of similarity solution in reversed stagnation-point flow. As with the governing equation of reversed stagnation-point flow, we can write the stream function as
(3.6)ψ=-Aνxf(η,τ),η=Aνy,τ=At,
where A is a constant proportional to V0(τ)/L, V0(τ) is the external flow velocity removing from the plane, and L is the characteristic length. These result in the governing equation (2.9)
(3.7)[fητ-(fη)2+ffηη-fηηη]η=0.
Integrating with respect to η, we have
(3.8)fητ-(fη)2+ffηη-fηηη=-C(τ).
Under the boundary conditions fη(∞,τ)=1, the value of C(τ) should be a constant and equal to 1. If the boundary condition fη(∞,τ) is restricted not to be a constant, a particular time-dependence function C(τ) may be expressed in the form
(3.9)C(τ)=cτ2,
where c is an arbitrary constant. The partial differential equation can be simplified by a similarity transformation when a new similarity variable is introduced. This converts the original partial differential equation into an ordinary differential equation. For a time dependent function, we introduce the diffusion variable transformation [7]
(3.10)ς=η1τ,f(η,τ)=1τF(ς).
Here ς is the time combined nondimensional variable and F(ς) are the nondimensional velocity functions. Substitution of the similarity transformation yields an ordinary differential equation
(3.11)-12ςF′′-F′-F′2+FF′′-F′′′=-c,
where the prime denotes the derivative with respect to the variable ς.Equation (3.11) is a third-order nonlinear ordinary differential equation and a key step in obtaining an analytical solution is to reversedrange the equation in an autonomous differential equation. In mathematics, an autonomous differential equation is a system of ordinary differential equations which does not explicitly depend on the independent variable.In order to omitting the variableς in the differential equation, it is recognized a change of variable
(3.12)Q=F-12ς
and the equation becomes an autonomous differential equation
(3.13)QQ′′-2Q′-Q′2-Q′′′=-c+34.
In our analysis, P=Q′ is the dependent variable and Q is the independent variable. Equation (3.13) is reversed ranged as
(3.14)QP′-2P-P2-P′′=-c+34
and the chain rule reduces equation (3.14) to a second-order ordinary differential equation
(3.15)QPdPdQ-2P-P2-PddQ(PdPdQ)=-c+34.
Equation (3.15) is analytically solvable that the solution might be expressed as a low order polynomial. It is suggested that
(3.16)P=a+bQ+dQ2,
and substituting into equation (3.14) and comparing the coefficients in the powers of Q results in a system of linear algebraic equations
(3.17)2a2d+ab2+2a-a2=-c+34,8abd+b3+2b+ab=0,8ad2+(7b2+2)d=0,12bd2-bd=0,6d3-d2=0.
Solving the related algebraic equation, we have
(3.18)a=-32,b=0,c=34,d=16.
Substituting the constant into (3.16) yields the first-order differential equation
(3.19)Q′=-32+16Q2.
Equation (3.19) is Riccati equation, which is any ordinary differential equation that is quadratic in the unknown function. The standard form of Riccati equation is
(3.20)Q′=RQ2+SQ+T.
The solution of Riccati equation can be obtained by a change of dependent variable, where the dependent variable y is changed to q by [11]
(3.21)Q=-q′q1R.
By identifying R=1/6, S=0, and T=-3/2, the change of variables in (3.19) becomes
(3.22)Q=-q′(1/6)q=-6q′q
so (3.19) becomes an second-order linear differential equation
(3.23)q′′-14=0.
The general solution to this equation is
(3.24)q=Acosh⁡ς2+Bsinh⁡ς2,
where A and B are arbitrary constants. Applying this solution in (3.12) leads to the general solution of (3.11)
(3.25)F(ς)=ς2-3Asinh⁡ς/2+3Bcosh⁡ς/2Acosh⁡ς/2+Bsinh⁡ς/2.
Application of the impermeability condition F(0)=0 leads to the determination of the constant B=0, so the exact solution becomes
(3.26)F(ς)=ς2-3tanhς2.
Collecting results, the velocity functions become
(3.27)f(η,τ)=1τ(ς2-3tanhς2),
where ς=(A/ντ)y is the nondimensional distance from the plate. In view of (3.27), the flow far from the boundary (x,ς→∞) becomes
(3.28)f(η,τ)⟶1τ(ς2-3).We obtain a particular solution of the unsteady reversed stagnation-point flow. The above solution is obtained in the similarity framework for unsteady viscous flows. The appearance of this positive factor in the first terms of (3.28) shows that this remote flow is directed toward the axis of symmetry and away from the plate. The second term in (3.28) describes a uniform velocity directed toward the plate. An adverse pressure gradient near the wall region leads to boundary-layer separation and associated flow reversal.The particular solution is noteworthy in that it is completely analytical, but it is limited to the region far away from the plate in the presence of nonzero termF′(0)=-1. Another no-slip boundary condition F′(0)=0 is not satisfied completely near the wall region.


## 3.3. Numerical Solution

Since the analytical solution does not satisfy the no-slip conditionF′(0)=0, it is convenient to solve the similarity equation numerically. The similarity equation and the relevant boundary conditions are
(3.29)-12ςF′′-F′-F′2+FF′′-F′′′=-c,F(0)=F′(0)=0,F′(∞)=12,
where c=3/4 in order to satisfy the unsteady viscous flows in the outer region.This equation is a third-order nonlinear ordinary differential equation. It is convenient when solving an ODE system numerically to describe the problem in terms of a system of first-order equations in MATLAB[12].For example when solving annth-order problem numerically is a common practice to reduce the equation to a system of n first-order equations. Then, by defining y1=F, y2=F′, y3=F′′, the ODE reduces to the form
(3.30)dydx=[y2y3c-12ςy3-y2-y22+y1*y3].The first task is to reduce the equation above to a system of first-order equations and define in MATLAB a function to return these. Later, we need to change the boundary value into initial value, becauseode45, anodesolver in MATLAB, can only solve the initial value problem. From (3.29), we gauss the value of F′′(0)=-1.54306 such that F′(∞)=1/2.The numerical solution for two-dimensional stagnation-point flow is shown in Figure4.Figure 4
Numerical solutions of reversed stagnation-point flow.This solution is a similarity solution of the reversed stagnation-point flow over a flat plate. It describes an unsteady viscous flow in both outer and inner regions. A single dividing streamline plane separates streamlines approaching the plate from external flow streamlines. Though the viscous termF′′′ is still small compared to the convective terms in the outer region, this is not true in neglecting the viscous terms within the total flow field. The similarity velocity field is shown in Figure 5.Figure 5
Similarity velocity field as a function ofς.


## 4. Conclusion

The foregoing study constitutes a similarity solution of the unsteady Navier-Stokes equations for reversed stagnation-point flow in the idealized case of an infinite plane boundary. In order to analyze the flow for nonzero values ofx, it is required to convert the full Navier-Stokes equations. This problem is now being studied by applying numerical method. The solution is obtained in the classical similarity framework for unsteady viscous flows.In the case of numerical methods, a brief analysis of the new solution is discussed. When the flow was near the plane wall region, the viscous forces were dominant, and the viscous term in the governing Navier-Stokes equations was important only near the boundary. On the contrary, the viscous forces were negligible when they were away from the wall.In this paper we have examined the case of two-dimensional reversed stagnation-point flow. However, with the establishment of this frame work, the more important practical properties in engineering and technology application, like the velocity of wall is function of time, the temperature of wall is function of time and distance from wall, can be investigated and they would be the next phase of this study.



---

*Source: 430432-2012-06-18.xml*