Instability of optimal non-axisymmetric base-flow deviations in pipe Poiseuille flow

2007 ◽  
Vol 588 ◽  
pp. 189-215 ◽  
Author(s):  
GUY BEN-DOV ◽  
JACOB COHEN

The stability of pipe flow when mildly deviating from the developed Poiseuille profile by a non-axisymmetric azimuthally periodic distortion is examined. The motivation for this is to consider deviations, the origin of which may be attributed to small-amplitude disturbances having sinusoidal periodicity along the azimuthal coordinate, which are known to be the ones most amplified by the transient growth linear mechanism. A mathematical technique for finding the minimum energy density of azimuthally periodic deviations triggering exponential instability is presented. The results show that owing to bifurcations multiple solutions of optimal deviations exist. As the Reynolds number is increased additional bifurcations appear and create more distinct solutions. The different solutions correspond to different radial distributions of the deviations, and at Reynolds numbers of about 2000 they are distributed over less than a half of the pipe radius. It is found that the dependence of the optimal deviation velocity leading to instability on the Reynolds number Re is approximately 20/Re. A comparison to axisymmetric base-flow deviations shows that the minimum energy required for an azimuthally periodic deviation to trigger instability is almost twice that for the axisymmetric flow. However, azimuthally periodic deviations, which are shown to have a streaky pattern, may have a role in the self-sustaining process. They may be formed as a result of a transient growth amplification of initial streamwise rolls and can produce, via self-interactions between the resulting growing waves, patterns of streamwise rolls as well.

Author(s):  
Peter W Duck

The stability of developing entry flow in both two-dimensional channels and circular pipes is investigated for large Reynolds numbers. The basic flow is generated by uniform flow entering a channel/pipe, which then provokes the growth of boundary layers on the walls, until (far downstream) fully developed flow is attained; the length for this development is well known to be (Reynolds number)×the channel/pipe width/diameter. This enables the use of high-Reynolds-number theory, leading to boundary-layer-type equations which govern the flow; as such, there is no need to impose heuristic parallel-flow approximations. The resulting base flow is shown to be susceptible to significant, three-dimensional, transient (initially algebraic) growth in the streamwise direction, and, consequently, large amplifications to flow disturbances are possible (followed by ultimate decay far downstream). It is suggested that this initial amplification of disturbances is a possible and alternative mechanism for flow transition.


1998 ◽  
Vol 358 ◽  
pp. 357-378 ◽  
Author(s):  
M. NAGATA

The stability of nonlinear tertiary solutions in rotating plane Couette flow is examined numerically. It is found that the tertiary flows, which bifurcate from two-dimensional streamwise vortex flows, are stable within a certain range of the rotation rate when the Reynolds number is relatively small. The stability boundary is determined by perturbations which are subharmonic in the streamwise direction. As the Reynolds number is increased, the rotation range for the stable tertiary motions is destroyed gradually by oscillatory instabilities. We expect that the tertiary flow is overtaken by time-dependent motions for large Reynolds numbers. The results are compared with the recent experimental observation by Tillmark & Alfredsson (1996).


Author(s):  
Nariman Ashrafi

The nonlinear stability and bifurcation of the one-dimensional channel (Poiseuille) flow is examined for a Johnson-Segalman fluid. The velocity and stress are represented by orthonormal functions in the transverse direction to the flow. The flow field is obtained from the conservation and constitutive equations using the Galerkin projection method. Both inertia and normal stress effects are included. The stability picture is dramatically influenced by the viscosity ratio. The range of shear rate or Weissenberg number for which the base flow is unstable increases from zero as the fluid deviates from the Newtonian limit as decreases. Typically, two turning points are observed near the critical Weissenberg numbers. The transient response is heavily influenced by the level of inertia. It is found that the flow responds oscillatorily. When the Reynolds number is small, and monotonically at large Reynolds number when elastic effects are dominated by inertia.


1999 ◽  
Vol 395 ◽  
pp. 211-236 ◽  
Author(s):  
V. SHANKAR ◽  
V. KUMARAN

Flows with velocity profiles very different from the parabolic velocity profile can occur in the entrance region of a tube as well as in tubes with converging/diverging cross-sections. In this paper, asymptotic and numerical studies are undertaken to analyse the temporal stability of such ‘non-parabolic’ flows in a flexible tube in the limit of high Reynolds numbers. Two specific cases are considered: (i) developing flow in a flexible tube; (ii) flow in a slightly converging flexible tube. Though the mean velocity profile contains both axial and radial components, the flow is assumed to be locally parallel in the stability analysis. The fluid is Newtonian and incompressible, while the flexible wall is modelled as a viscoelastic solid. A high Reynolds number asymptotic analysis shows that the non-parabolic velocity profiles can become unstable in the inviscid limit. This inviscid instability is qualitatively different from that observed in previous studies on the stability of parabolic flow in a flexible tube, and from the instability of developing flow in a rigid tube. The results of the asymptotic analysis are extended numerically to the moderate Reynolds number regime. The numerical results reveal that the developing flow could be unstable at much lower Reynolds numbers than the parabolic flow, and hence this instability can be important in destabilizing the fluid flow through flexible tubes at moderate and high Reynolds number. For flow in a slightly converging tube, even small deviations from the parabolic profile are found to be sufficient for the present instability mechanism to be operative. The dominant non-parallel effects are incorporated using an asymptotic analysis, and this indicates that non-parallel effects do not significantly affect the neutral stability curves. The viscosity of the wall medium is found to have a stabilizing effect on this instability.


2019 ◽  
Vol 864 ◽  
pp. 5-44 ◽  
Author(s):  
David Fabre ◽  
Raffaele Longobardi ◽  
Paul Bonnefis ◽  
Paolo Luchini

The unsteady axisymmetric flow through a circular aperture in a thin plate subjected to harmonic forcing (for instance under the effect of an incident acoustic wave) is a classical problem first considered by Howe (Proc. R. Soc. Lond. A, vol. 366, 1979, pp. 205–223), using an inviscid model. The purpose of this work is to reconsider this problem through a numerical resolution of the incompressible linearized Navier–Stokes equations (LNSE) in the laminar regime, corresponding to $Re=[500,5000]$. We first compute a steady base flow which allows us to describe the vena contracta phenomenon in agreement with experiments. We then solve a linear problem allowing us to characterize both the spatial amplification of the perturbations and the impedance (or equivalently the Rayleigh conductivity), which is a key quantity to investigate the response of the jet to acoustic forcing. Since the linear perturbation is characterized by a strong spatial amplification, the numerical resolution requires the use of a complex mapping of the axial coordinate in order to enlarge the range of Reynolds number investigated. The results show that the impedances computed with $Re\gtrsim 1500$ collapse onto a single curve, indicating that a large Reynolds number asymptotic regime is effectively reached. However, expressing the results in terms of conductivity leads to substantial deviation with respect to Howe’s model. Finally, we investigate the case of finite-amplitude perturbations through direct numerical simulations (DNS). We show that the impedance predicted by the linear approach remains valid for amplitudes up to order $10^{-1}$, despite the fact that the spatial evolution of the perturbations in the jet is strongly nonlinear.


2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.


1968 ◽  
Vol 90 (1) ◽  
pp. 109-114 ◽  
Author(s):  
Ahmed R. Wazzan ◽  
T. Okamura ◽  
A. M. O. Smith

The theory of two-dimensional instability of laminar flow of water over solid surfaces is extended to include the effects of heat transfer. The equation that governs the stability of these flows to Tollmien-Schlichting disturbances is the Orr-Sommerfeld equation “modified” to include the effect of viscosity variation with temperature. Numerical solutions to this equation at high Reynolds numbers are obtained using a new method of integration. The method makes use of the Gram-Schmidt orthogonalization technique to obtain linearly independent solutions upon numerically integrating the “modified Orr-Sommerfeld” equation using single precision arithmetic. The method leads to satisfactory answers for Reynolds numbers as high as Rδ* = 100,000. The analysis is applied to the case of flow over both heated and cooled flat plates. The results indicate that heating and cooling of the wall have a large influence on the stability of boundary-layer flow in water. At a free-stream temperature of 60 deg F and wall temperatures of 60, 90, 120, 135, 150, 200, and 300deg F, the critical Reynolds numbers Rδ* are 520, 7200, 15200, 15600, 14800, 10250, and 4600, respectively. At a free-stream temperature of 200F and wall temperature of 60 deg F (cooled case), the critical Reynolds number is 151. Therefore, it is evident that a heated wall has a stabilizing effect, whereas a cooled wall has a destabilizing effect. These stability calculations show that heating increases the critical Reynolds number to a maximum value (Rδ* max = 15,700 at a temperature of TW = 130 deg F) but that further heating decreases the critical Reynolds number. In order to determine the influence of the viscosity derivatives upon the results, the critical Reynolds number for the heated case of T∞ = 40 and TW = 130 deg F was determined using (a) the Orr-Sommerfeld equation and (b) the present governing equation. The resulting critical Reynolds numbers are Rδ* = 140,000 and 16,200, respectively. Therefore, it is concluded that the terms pertaining to the first and second derivatives of the viscosity have a considerable destabilizing influence.


Author(s):  
Isabella M. Gavarini ◽  
Alessandro Bottaro ◽  
Frans T. M. Nieuwstadt

Transition in a cylindrical pipe flow still eludes thorough understanding. Most recent advances are based on the concept of transient growth of disturbances, but even this scenario is not fully confirmed by DNS and/or experiments. Based on the fact that even the most carefully conducted experiment is biased by uncertainties, we explore the spatial growth of disturbances developing on top of an almost ideal, axially invariant Poiseuille flow. The optimal deviation of the base flow from the ideal parabolic profile is computed by a variational tecnique, and unstable modes, driven by an inviscid mechanism, are found to exist for very small values of the norm of the deviation, at low Reynolds numbers.


Author(s):  
A. Inasawa ◽  
K. Toda ◽  
M. Asai

Disturbance growth in the wake of a circular cylinder moving at a constant acceleration is examined experimentally. The cylinder is installed on a carriage moving in the still air. The results show that the critical Reynolds number for the onset of the global instability leading to a self-sustained wake oscillation increases with the magnitude of acceleration, while the Strouhal number of the growing disturbance at the critical Reynolds number is not strongly dependent on the magnitude of acceleration. It is also found that with increasing the acceleration, the Ka´rma´n vortex street remains two-dimensional even at the Reynolds numbers around 200 where the three-dimensional instability occurs to lead to the vortex dislocation in the case of cylinder moving at constant velocity or in the case of cylinder wake in the steady oncoming flow.


2000 ◽  
Vol 123 (2) ◽  
pp. 347-358 ◽  
Author(s):  
P. Bagchi ◽  
M. Y. Ha ◽  
S. Balachandar

Direct numerical solution for flow and heat transfer past a sphere in a uniform flow is obtained using an accurate and efficient Fourier-Chebyshev spectral collocation method for Reynolds numbers up to 500. We investigate the flow and temperature fields over a range of Reynolds numbers, showing steady and axisymmetric flow when the Reynolds number is less than 210, steady and nonaxisymmetric flow without vortex shedding when the Reynolds number is between 210 and 270, and unsteady three-dimensional flow with vortex shedding when the Reynolds number is above 270. Results from three-dimensional simulation are compared with the corresponding axisymmetric simulations for Re>210 in order to see the effect of unsteadiness and three-dimensionality on heat transfer past a sphere. The local Nusselt number distribution obtained from the 3D simulation shows big differences in the wake region compared with axisymmetric one, when there exists strong vortex shedding in the wake. But the differences in surface-average Nusselt number between axisymmetric and three-dimensional simulations are small owing to the smaller surface area associated with the base region. The shedding process is observed to be dominantly one-sided and as a result axisymmetry of the surface heat transfer is broken even after a time-average. The one-sided shedding also results in a time-averaged mean lift force on the sphere.


Sign in / Sign up

Export Citation Format

Share Document