Stability analysis of Boundary Layer in Poiseuille Flow Through A Modified Orr-Sommerfeld Equation
For applications regarding transition prediction, wing design and control of boundary layers, the fundamental understanding of disturbance growth in the flat-plate boundary layer is an important issue. In the present work we investigate the stability of boundary layer in Poiseuille flow. We normalize pressure and time by inertial and viscous effects. The disturbances are taken to be periodic in the spanwise direction and time. We present a set of linear governing equations for the parabolic evolution of wavelike disturbances. Then, we derive modified Orr-Sommerfeld equations that can be applied in the layer. Contrary to what one might think, we find that Squire’s theorem is not applicable for the boundary layer. We find also that normalization by inertial or viscous effects leads to the same order of stability or instability. For the disturbances flow (), we found the same critical Reynolds number for our two normalizations. This value coincides with the one we know for neutral stability of the known Orr-Sommerfeld equation. We noticed also that for all overs values of in the case correspond the same values of at whatever the normalization. We therefore conclude that in the boundary layer with a 2D-disturbance, we have the same neutral stability curve whatever the normalization. We find also that for a flow with hight hydrodynamic Reynolds number, the neu- tral disturbances in the boundary layer are two-dimensional. At last, we find that transition from stability to instability or the opposite can occur according to the Reynolds number and the wave number.
Institut de Mathématiques et de Sciences Physiques, BP: 613 Porto Novo, Bénin
The Abdus Salam International Centre for Theoretical Physics, Trieste, Italy
Boundary-layer theory is crucial in understanding why certain phenomena occur. It is well known that the instability of boundary layers is sensitive to the mean velocity profile, so that a small distortion to the basic flow may have a detrimental effect on its stability. Prandtl (1904) proposed that viscous effects would be confined to thin layers adjacent to boundaries in the case of the motion of fluids with very little viscosity i.e. in the case of flows for which the characteristic Reynolds number, , is large. In a more general sense we will use boundary-layer theory (BLT) to refer to any large-Reynolds-number. Ho and Denn studied low Reynolds number stability for plane Poiseuille flow by using a numerical scheme based on the shooting method. They found that at low Reynolds numbers no instabilities occur, but the numerical method led to artificial instabilities.Lee and Finlayson used a similar numerical method to study both Poiseuille and Couette flow, and confirmed the absence of instabilities at low Reynolds number. R. N. RAY et al  investigated the linear stability of plane Poiseuille flow at small Reynolds number of a conducting Oldroyd fluid in the presence of magnetic field. They found that viscoelastic parameters have destabilizing effect and magnetic field has a stabilizing effect in the field of flow but no instabilities are found.
In this paper, we study the linear stability of boundary layer in a plane
Poiseuille flow. For this, we derive two fourth-order equations name modified fourth-order Orr-Sommerfeld equations governing the stability analysis in boundary layer for the flow. The first is obtained by making dimensionless quantities by the inertial effects. The second takes into account the form adopted by the rheologists i.e. make the quantities dimensionless by normalizing by the viscous effects. This allowed us to see the effect of each type of normalization on the stability in the boundary layer. So, we solve numerically the corresponding eigenvalues problems. We employ Matlab in all our numerical computations to find eigenvalues.
The paper is organized as follows. In the second section the boundary layer theory is presented. In the third section we present the general formulation, highlighting the fundamental equations that model the flat-plate boundary layer flow according to the normalization by inertial and viscous effects. In the fourth section the modified Orr-Sommerfeld equations governing the stability analysis in boundary layer are checked and in the fifth section, analysis of the stability is investigated. The conclusions and perspectives are presented in the final section.
2 Boundary Layer Theory
When applying the theory of complex potential around an airfoil in considering the model of inviscid incompressible irrotational plan flow, we know that the model allows to deduce the lift but the drag is zero. This is contrary to experimental observations which show that the drag affects all flows of real fluids. These are viscous. They adhere to the walls and the tangential component of the velocity is zero if the wall is fixed. The latter condition can be satisfied by the perfect fluid. Moreover, the irrotational condition is far from reality as we know that the production of vorticity occurs at the walls. To remedy the deficiencies of the theory of perfect fluid, it must appeal to the theory of the boundary layer which is a necessary correction for flows with high Reynolds numbers. Theoris the boundary layer is due to L. Prandtl. The boundary layer is the area of the flow which is close to the wall or of an obstacle present in a uniform flow at the upstream infinity or on the confining walls of internal flow. Within the boundary layer is a thin zone, it is estimated that viscous effects are of the same magnitude as the inertial effects. The boundary layer is the place of intense generation of vorticity which will not diffuse into the area outside thereof. This leads to a very modern concept of comprehensive approach to the problem by breaking it down into two areas: firstly the boundary layer where we will consider the viscous effects in a model of simplified Navier-Stokes and other from the outer area where we will use the complex potential theory in the inviscid incompressible flow. This outer zone has speeds which are of the same order of magnitude as that of the incident flow.
The boundary layer along an obstacle is therefore thin since the fluid travel great distances downstream of the leading edge during the time interval during which the vortex diffuse only a small distance from the wall. The creation of vorticity in the boundary layer allows the physical realization of the fluid flow around the profile. This movement gives rise to a wake in the area near the trailing edge. The importance of the wake depend on the shape of the obstacle and the angle of incidence of the upstream flow at the leading edge.
We consider incompressible flow of a fluid with constant density and dynamic viscosity , past a body with typical length . We assume that a typical velocity scale is , and the Reynolds number is given by
For simplicity we will, for the most part, consider two-dimensional incompressible flows, although many of our statements can be generalised to three-dimensional flowsand/or compressible flows.
Boundary Layer Theory applies to flows where there are extensive inviscid regions separated by thin shear layers, say, of typical width . For one such shear layer take local dimensional Cartesian coordinates and along and across the shear layer respectively. Denote the corresponding velocity components by and respectively, pressure by and time by . On the basis of scaling arguments it then follows that
Further, it can also be deduced that the key approximations in classical Boundary Layer Theory are that the pressure is constant across the shear layer i.e.
and that streamwise diffusion is negligible, i.e. if represents any variable
The former approximation is more significant dynamically.
Now, using the transformations
and taking the limit , the Boundary Layer Theory equations can be deduced from the Navier-Stokes equations:
For flow past a rigid body the appropriate boundary conditions are
where is the inviscid slip velocity past the body. Further, from (6) evaluated at the edge of the boundary layer
We define the viscous blowing velocity out of the boundary layer to be
indicates the strength of blowing, or suction, at the edge of the boundary layer induced by viscous effects. It is good diagnostic for dynamically significant effects within the boundary layer much better than, say, the wall shear which can remain regular while becomes unbounded.
3 General formulation
Consider an incompressible boundary layer over a flat plate as illustrated in the following Figure.
The streamwise coordinate is scaled with the length scale , which is a fixed distance from the leading edge. The wall-normal and spanwise coordinates and , respectively, are scaled with the boundary-layer parameter , where is the kinematic viscosity of the fluid and is the streamwise freestream velocity at the distance from the leading edge. The streamwise velocity is scaled with , while the wall- normal and spanwise velocities and , respectively, are scaled with . The pressure is scaled with and the time is scaled with . The Reynolds numbers used here are defined as and . It is useful to note the relations . We want to study the linear stability of a high Reynolds number flow. The non-dimensional Navier–Stokes equations for an incompressible flow where pressure and time are normalized by inertial effects
are linearized around a two-dimensional, steady base flow to obtain the stability equations for the spatial evolution of three-dimensional, time-dependent disturbances . The base flow and the disturbances are scaled in the same way. The disturbances are taken to be periodic in the spanwise direction and time, which allows us to assume solutions of the form
where represents either one of the disturbances , , or . The complex streamwise wave number captures the fast wavelike variation of the modes and is therefore scaled with . itself is assumed to vary slowly with . Since is scaled with , the factor appears in front of the integral. The -dependence in the amplitude function includes the weak variation of the disturbances. The real spanwise wave number and angular frequency are scaled in a consistent way with and , respectively. Introducing the assump- tion (11) in the linearized Navier–Stokes equations and neglecting all third order terms in or higher, we arrive at the parabolized stability equations in boundary-layer scalings
If we normalize pressure and time by viscous effects i.e. scaled with and scaled with the Navier-Stokes equations take the following forms
The linearized Navier-Stokes equations with the previous disturbances under the same considerations become
4 Modified Orr-Sommerfeld Equation
Considering temporal the stability problem with and real, we will simplify the problem. Our strategy will be first to eliminate , , to leave a single equation in . This can be used finally to determine the linear stability (or instability) in the boundary layer oy the base flow. Remember that and if , the flow is stable; , the flow is unstable and we have neutral stability if .
Taking + we get
Using the continuity equation (16), (28) becomes
Operating across (29) with , using the assumptions in boundary layer and the equation (18) we get the modified Orr-Sommerfeld equation for the boundary layer
Considering normalization with viscous effects we would have
5 Stability Analysis
We put equations (30)-(31) respectively in the form
In order to investigate the application of the Squire theorem, consider first the normalization with inertial effects (32). Often, however, we are only interested in the instability that appears first as the control parameter is increased. In this case, Squire’s theorem tells us that we need only consider disturbances.
Consider a base state . Imagine a growing disturbance to this base state at Reynolds number , with wavenumbers , , and . This corresponds to a solution , with ( of the modified Orr-Sommerfeld equation
Now consider a disturbance at a Reynolds number . This has , and must satisfy the modified Orr-Sommerfeld equation
For values , and , this modified Orr-Sommerfeld equation has the form
which is exactly the same as 34. It must therefore have the same growing solution with .
Therefore, corresponding to the growing disturbance at with , , and , there exists a growing disturbance at with and . These conditions lead to and we get that the disturbances are two-dimensional. Finally, we get , not . And so, we can not applied the Squire’s theorem in the boundary layer of the flow.
Using the same assumptions in (33), we get first with i.e. and secondly with i.e. . So we have first and secondly because . We see therefore that we must take and so . We find the same result that the disturbances are two-dimensional and then Squire’s theorem can’t be applied.
Finally, we still consider our three-dimensional disturbances so without using the theorem of Squire. Thus we write with . This will allow us to deduce numerically whether the application of Squire’s theorem in the boundary layer. Indeed corresponds to a two dimensional disturbance i.e .
We employ Matlab (Windows Version) in all our numerical computations to find eigenvalues. A Poiseuille flow with the basic profile
The eigenvalue problems (32)-(33) are solved numerically with the suitable boundary conditions. The solutions are found in a layer bounded at with . The results of calculations are presented in the following figures. In each group of figures, the first one is figure , the second is figure and the third is figure . First we present the figures relatives to the eigenvalue problem (32).
For a fixed , we get figure 1 of vs for sequential values of in which shows the entire graph. and are the magnified versions of .
For a fixed , we get figure 2 of vs for sequential values of in which shows the entire graph. and are the magnified versions of .
For a fixed , we get figure 3 of vs for sequential values of in which shows the entire graph. and are the magnified versions of .
Secondly we present the figures relatives to the eigenvalue problem (33).
For a fixed , we get figure 4 of vs for sequential values of in which shows the entire graph. and are the magnified versions of .
For a fixed , we get figure 5 of vs for sequential values of in which shows the entire graph. and are the magnified versions of .
Through the figures and , it is easy to see that if we take a curve with in the instability area , we don’t have necessary in the two normalization . We also see through the figures that the normalization of time and pressure by inertial effects or viscous effects lead to the same order of stability/instability in the boundary layer. This confirms that in the boundary layer, viscous forces are on the same magnitude as the inertial forces i.e. the local Reynolds number is on unity order. We see also through the figures that at low Reynolds number the flow is stable but if the Reynolds number increase, unstability appears. The increase of the wave number induce also the stability of the flow. By figures and we find the first value of and for which the first transition from stability to instability occurs i.e. in the case . We find the values which corresponds exactly to the one we know as critical value of the neutral stability of Poiseuille flow. Note that this value is the same in figure which corresponds to normalization by inertial effects and in figure which corresponds to normalization by viscous effects. We noticed also that for all overs values of in the case correspond the same values of at whatever the normalization. We therefore conclude that in the boundary layer with a 2D-disturbance, we have the same neutral stability curve whatever the normalization.
In this paper, we have investigated the stability of boundary layer in Poiseuille flow. We have shown that the instability of the perturbed flow is governed by a remarkably equation named modified Orr-Sommerfeld equation. Contrary to what one might think, we find that Squire’s theorem is not applicable for the boundary layer. We find also that normalization by inertial or viscous effects leads to the same order of stability or instability. For the disturbances flow (), we found the same critical Reynolds number for our two normalizations. This value coincides with the one we know for neutral stability of the known Orr-Sommerfeld equation. We noticed also that for all overs values of in the case correspond the same values of at whatever the normalization. We therefore conclude that in the boundary layer with a 2D-disturbance, we have the same neutral stability curve whatever the normalization.
The authors thank IMSP-UAC for financial support.
-  J. G. Oldroyd. On the formulation of rheological equations of state. Proc. Roy. Soc. London Ser. A. 200 (1950), 523–541. MR 11,703a.
-  T. C. Ho and M. M. Denn. Stability of plane Poiseuille flow of a highly elastic liquid. J. Non- Newtonian Fluid Mech. 3 (1977), 179–195. Zbl 414.76008.
-  Bertolotti, F. P., Herbert, T. and Spalart.P. R Linear and nonlinear stability of the Blasius boundary layer.J. Fluid Mech, 1992, 242, 441–474.
-  K. C. Lee and B. A. Finlayson. Stability of plane Poiseuille and Couette flow of a Maxwell fluid. J. Non-Newtonain Fluid Mech. 21 (1986), 65–78. Zbl 587.76059.
-  R. C. Lock. The stability of the flow of an electrically conducting fluid between parallel planes under a transverse magnetic field. Proc. Roy. Soc. London Ser. A. 233 (1955), 105–125. MR 18,700a.
-  Abu-Ghannam, B. and Shaw, R. Natural transition of boundary layers - the effects of turbulence, pressure gradient, and flow history. J. Mech. Eng. Sci.1980- 22 (5), 213–228.
-  Andersson, P., Berggren, M. and Henningson, D. Optimal disturbances and bypass transition in boundary layers. Phys.Fluids 11 (1),1999, 134–150.
-  Andersson, P., Henningson, D. and Hanifi, A. On a stabilization procedure for the parabolic stability equations. J. Eng. Math.1998, 33, 311–332.
-  Landau, L. & Lifchitz, E. M. Fluid Mechanics. Butterworth Heinemann, Oxford, 1997.
-  Ray, Samad, and Chaudhury: Department of Mathematics, University of Burdwan, Burdwan-713104, West Bengal, India.