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.
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 Figure 1 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.
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 . 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 , 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 is the normal to the plane.
By conservation of mass principle with constant physical properties, the equation of continuity is We consider the two-dimensional reversed stagnation-point flow in unsteady state and the flow is bounded by an infinite plane , the fluid remains at rest when time . At , it starts impulsively in motion which is determined by the stream function
At large distances far above the planar boundary, the existence of the potential flow implies an inviscid boundary condition. It is given by where and are the components of flow velocity, is a constant proportional to , is the external flow velocity removing from the plane, and is the characteristic length. We have at and at , but the no-slip boundary at wall 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: 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 where and are the velocity components along and axes, and is the density.
Substituting and into the governing equations results in a simplified partial differential equation. From the definition of the stream function, we have
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:
The pressure gradient can be again reduced by a further differentiation equation (2.7b) with respect to . That is and (2.7a) reduces to The initial and boundary conditions are The last condition reduces the above differential equation (2.9) to the form with the boundary conditions Equation (2.11) is the similarity equation of the full Navier-Stokes equations at a two-dimensional reversed stagnation point. The coordinates and are replaced by a dimensionless variable . Under the boundary conditions , when the flow is in steady state such that , 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 . 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 Proudman and Johnson obtained that a similarity solution of (3.1) is in the form where and is a constant of integration. Robins and Howarth [6] estimated the value of 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 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 is relatively small, the solution (3.3) yields Substituting in (2.11) yields
This is exactly the classic stagnation-point problem (Hiemenz [1]) by changing the sign in . 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.
Although an asymptotic solution was obtained, it can easily been observed that this is not valid when the viscous term 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 time and is set in motion at such that at large distances far above the planar boundary the potential flow is a constant for all value of . Both Proudman and Johnson [5], and Robins and Howarth [6] have set and the corresponding boundary condition .
When the flow is in steady state such that , it was proven that the similarity velocity cannot ultimately approach to 1. The differential equation has no solution. If the potential flow is restricted not to be a constant, the boundary condition 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 where is a constant proportional to , is the external flow velocity removing from the plane, and is the characteristic length. These result in the governing equation (2.9) Integrating with respect to , we have Under the boundary conditions , the value of should be a constant and equal to . If the boundary condition is restricted not to be a constant, a particular time-dependence function may be expressed in the form where 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] Here is the time combined nondimensional variable and are the nondimensional velocity functions. Substitution of the similarity transformation yields an ordinary differential equation 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 and the equation becomes an autonomous differential equation In our analysis, is the dependent variable and is the independent variable. Equation (3.13) is reversed ranged as and the chain rule reduces equation (3.14) to a second-order ordinary differential equation Equation (3.15) is analytically solvable that the solution might be expressed as a low order polynomial. It is suggested that and substituting into equation (3.14) and comparing the coefficients in the powers of results in a system of linear algebraic equations Solving the related algebraic equation, we have Substituting the constant into (3.16) yields the first-order differential equation 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 The solution of Riccati equation can be obtained by a change of dependent variable, where the dependent variable is changed to by [11] By identifying , , and , the change of variables in (3.19) becomes so (3.19) becomes an second-order linear differential equation The general solution to this equation is where and are arbitrary constants. Applying this solution in (3.12) leads to the general solution of (3.11) Application of the impermeability condition leads to the determination of the constant , so the exact solution becomes Collecting results, the velocity functions become where is the nondimensional distance from the plate. In view of (3.27), the flow far from the boundary becomes
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 term . Another no-slip boundary condition is not satisfied completely near the wall region.
3.3. Numerical Solution
Since the analytical solution does not satisfy the no-slip condition , it is convenient to solve the similarity equation numerically. The similarity equation and the relevant boundary conditions are where 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 an th-order problem numerically is a common practice to reduce the equation to a system of first-order equations. Then, by defining , , , the ODE reduces to the form
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, because , anodesolver in MATLAB, can only solve the initial value problem. From (3.29), we gauss the value of such that .
The numerical solution for two-dimensional stagnation-point flow is shown in Figure 4.
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 term 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.
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 of , 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.
Acknowledgments
The research was partially supported by Research Committee of University of Macau through Grants no. RG079/09-10S//11T/SVK/FST and by Science and Technology Development Fund (FDCT) of Macao SAR through Grant no. 034/2009/A.