Energy budget analysis and neutral curve characteristics for the linear instability of Couette–Poiseuille flow

2021 ◽  
Vol 33 (3) ◽  
pp. 034102
Author(s):  
Srinivas Kirthy K. ◽  
Sourabh S. Diwan
1976 ◽  
Vol 73 (1) ◽  
pp. 153-164 ◽  
Author(s):  
P.-A. Mackrodt

The linear stability of Hagen-Poiseuille flow (Poiseuille pipe flow) with superimposed rigid rotation against small three-dimensional disturbances is examined at finite and infinite axial Reynolds numbers. The neutral curve, which is obtained by numerical solution of the system of perturbation equations (derived from the Navier-Stokes equations), has been confirmed for finite axial Reynolds numbers by a few simple experiments. The results suggest that, at high axial Reynolds numbers, the amount of rotation required for destabilization could be small enough to have escaped notice in experiments on the transition to turbulence in (nominally) non-rotating pipe flow.


When two parallel plates move normal to each other with a slow time-dependent speed, the velocity field developed in the intervening film of fluid is approximately that of plane Poiseuille flow, except that the magnitude of the velocity is dependent on time and on the coordinate parallel to the planes. This fact is intrinsic to Reynolds’ lubrication theory, and can be shown to follow from the Navier-Stokes equations when both the modified Reynolds number ( Re M ) and an aspect ratio ( δ ) are small. The modified Reynolds number is the product of δ and an actual Reynolds number ( Re ), which is based on the gap between the planes and on a characteristic velocity. The occurrence of flow instability and of turbulence in the film depend on Re . Typical values of Re , which are known to be required for the linear instability of plane Poiseuille flow, are of order 6000. This condition can be achieved, even if Re M is of order 1, provided that δ is of order 10 -4 . Such parameter values are typical of lubrication problems. The Orr-Sommerfeld equation governing flow instability is derived in this paper by use of the WKBJ technique, δ being the approximate small parameter to represent the small length-scale of the disturbance oscillations compared with the larger scale of the basic laminar flow. However, the coefficients in the Orr-Sommerfeld equation depend on slow space and time variables. Consequently the eigenrelation, derivable from the Orr-Sommerfeld equation and the associated boundary conditions, constitutes a nonlinear first-order partial differential equation for a phase function. This equation is solved by use of Charpit’s method for certain special forms of the time-dependent gap between the planes, followed by detailed numerical calculations. The relation between time-dependence and flow instability is delineated by the calculated results. In detail the nature of the instability can be described as follows. We consider a disturbance wave at or near a particular station, the initial distribution of amplitude being gaussian in the slow coordinate parallel to the planes. In the context of the Orr-Sommerfeld equation and its eigenrelation, the particular station implies an equivalent Reynolds number, while the initial distribution of the disturbance wave implies an equivalent wavenumber. As time increases, the disturbance wave can be considered to move in the instability diagram of equivalent wavenumber against Reynolds number, in the sense that these parameters are time- and space-dependent for the evolution of the disturbance-wave system. For our detailed calculations we use a quadratic approximation to the eigenrelation, an approximation which is quite accurate. If the initial distribution implies a point within the neutral curve, when the plates are squeezed together the equivalent wavenumber falls while the equivalent Reynolds number rises, and amplification takes place until the lower branch of the neutral curve is nearly crossed. If the plates are pulled apart (dilatation) the equivalent wavenumber rises, while the Reynolds number drops, and amplification takes place until the upper branch of the neutral curve has been just crossed. In the case of dilatation the transition from amplification to damping takes place more quickly than for the case of squeezing, in part due to the geometry of the neutral curve.


1979 ◽  
Vol 15 (6) ◽  
pp. 1589-1592 ◽  
Author(s):  
L. W. Gay ◽  
L. J. Fritschen

1981 ◽  
Vol 23 ◽  
pp. 61-76 ◽  
Author(s):  
C.H.M. Van Bavel ◽  
J. Damagnez ◽  
E.J. Sadler

2002 ◽  
Vol 451 ◽  
pp. 35-97 ◽  
Author(s):  
S. JONATHAN CHAPMAN

Certain laminar flows are known to be linearly stable at all Reynolds numbers, R, although in practice they always become turbulent for sufficiently large R. Other flows typically become turbulent well before the critical Reynolds number of linear instability. One resolution of these paradoxes is that the domain of attraction for the laminar state shrinks for large R (as Rγ say, with γ < 0), so that small but finite perturbations lead to transition. Trefethen et al. (1993) conjectured that in fact γ <−1. Subsequent numerical experiments by Lundbladh, Henningson & Reddy (1994) indicated that for streamwise initial perturbations γ =−1 and −7/4 for plane Couette and plane Poiseuille flow respectively (using subcritical Reynolds numbers for plane Poiseuille flow), while for oblique initial perturbations γ =−5/4 and −7/4 Here, through a formal asymptotic analysis of the Navier–Stokes equations, it is found that for streamwise initial perturbations γ =−1 and −3/2 for plane Couette and plane Poiseuille flow respectively (factoring out the unstable modes for plane Poiseuille flow), while for oblique initial perturbations γ =−1 and −5/4. Furthermore it is shown why the numerically determined threshold exponents are not the true asymptotic values.


Cardiometry ◽  
2014 ◽  
pp. 17-34
Author(s):  
Sergey Chefranov ◽  
◽  
Alexander Chefranov

1994 ◽  
Vol 263 ◽  
pp. 133-150 ◽  
Author(s):  
I. A. Frigaard ◽  
S. D. Howison ◽  
I. J. Sobey

The stability to linearized two-dimensional disturbances of plane Poiseuille flow of a Bingham fluid is considered. Bingham fluids exhibit a yield stress in addition to a plastic viscosity and this description is typically applied to drilling muds. A non-zero yield stress results in an additional parameter, a Bingham number, and it is found that the minimum Reynolds number for linear instability increases almost linearly with increasing Bingham number.


Sign in / Sign up

Export Citation Format

Share Document