scholarly journals Transient linear stability of pulsating Poiseuille flow using optimally time-dependent modes

2021 ◽  
Vol 927 ◽  
Author(s):  
J.S. Kern ◽  
M. Beneitez ◽  
A. Hanifi ◽  
D.S. Henningson

Time-dependent flows are notoriously challenging for classical linear stability analysis. Most progress in understanding the linear stability of these flows has been made for time-periodic flows via Floquet theory focusing on time-asymptotic stability. However, little attention has been given to the transient intracyclic linear stability of periodic flows since no general tools exist for its analysis. In this work, we explore the potential of using the recent framework of the optimally time-dependent (OTD) modes (Babaee & Sapsis, Proc. R. Soc. Lond. A, vol. 472, 2016, 20150779) to extract information about both the transient and the time-asymptotic linear stability of pulsating Poiseuille flow. The analysis of the instantaneous OTD modes in the limit cycle leads to the identification of the dominant instability mechanism of pulsating Poiseuille flow by comparing them with the spectrum and the eigenmodes of the Orr–Sommerfeld operator. In accordance with evidence from recent direct numerical simulations, it is found that structures akin to Tollmien–Schlichting waves are the dominant feature over a large range of pulsation amplitudes and frequencies but that for low pulsation frequencies these modes disappear during the damping phase of the pulsation cycle as the pulsation amplitude is increased beyond a threshold value. The maximum achievable non-normal growth rate during the limit cycle was found to be nearly identical to that in plane Poiseuille flow. The existence of subharmonic perturbation cycles compared with the base flow pulsation is documented for the first time in pulsating Poiseuille flow.

2010 ◽  
Vol 656 ◽  
pp. 417-447 ◽  
Author(s):  
ANIRBAN GUHA ◽  
IAN A. FRIGAARD

We present a detailed study of the linear stability of the plane Couette–Poiseuille flow in the presence of a crossflow. The base flow is characterized by the crossflow Reynolds number Rinj and the dimensionless wall velocity k. Squire's transformation may be applied to the linear stability equations and we therefore consider two-dimensional (spanwise-independent) perturbations. Corresponding to each dimensionless wall velocity, k ∈ [0, 1], two ranges of Rinj exist where unconditional stability is observed. In the lower range of Rinj, for modest k we have a stabilization of long wavelengths leading to a cutoff Rinj. This lower cutoff results from skewing of the velocity profile away from a Poiseuille profile, shifting of the critical layers and the gradual decrease of energy production. Crossflow stabilization and Couette stabilization appear to act via very similar mechanisms in this range, leading to the potential for a robust compensatory design of flow stabilization using either mechanism. As Rinj is increased, we see first destabilization and then stabilization at very large Rinj. The instability is again a long-wavelength mechanism. An analysis of the eigenspectrum suggests the cause of instability is due to resonant interactions of Tollmien–Schlichting waves. A linear energy analysis reveals that in this range the Reynolds stress becomes amplified, the critical layer is irrelevant and viscous dissipation is completely dominated by the energy production/negation, which approximately balances at criticality. The stabilization at very large Rinj appears to be due to decay in energy production, which diminishes like Rinj−1. Our study is limited to two-dimensional, spanwise-independent perturbations.


Author(s):  
Anirban Guha ◽  
Ian A. Frigaard

We have investigated the linear stability of plane Couette-Poiseuille flow in the presence of a cross-flow. The base flow is characterised by the cross flow Reynolds number, Ri and the dimensionless wall velocity, k. Corresponding to each k ∈ [0,1], we have observed two ranges of Ri for which the flow is unconditionally linearly stable. In the lower range, we have a stabilisation of long wavelengths leading to a cut-off Ri. In this range, cross-flow stabilisation and Couette stabilisation appear to act via very similar mechanisms in this range, leading to the potential for robust compensatory design of flow stabilisation using either mechanism. As Ri is increased, we see first destabilisation and then stabilisation at very large Ri. The instability is again a long wavelength mechanism. A linear energy analysis reveals that in this range the Reynolds stress becomes amplified, the critical layer is irrelevant and viscous dissipation is completely dominated by the energy production/negation, which approximately balances at criticality.


Author(s):  
H. Lüdeke ◽  
R. von Soldenhoff

AbstractTo determine allowable tolerances between successive suction panels at hybrid laminar wings with suction surfaces, direct numerical simulations of Tollmien–Schlichting waves over different steps are carried out for realistic suction rates on a wind tunnel configuration. Simulations at given suction panel positions over forward and backward facing steps are carried out by the use of a high-order method for the direct simulation of Tollmien–Schlichting wave growth. Comparisons between high-fidelity direct numerical simulations and quick linear stability calculations have shown capabilities and limits of the well-validated linear stability theory design approach.


2019 ◽  
Vol 75 ◽  
pp. 29-43 ◽  
Author(s):  
Souvik Ghosh ◽  
Jean-Christophe Loiseau ◽  
Wim-Paul Breugem ◽  
Luca Brandt

Fluids ◽  
2021 ◽  
Vol 6 (12) ◽  
pp. 432
Author(s):  
Silvia C. Hirata ◽  
Mohamed Najib Ouarzazi

The onset of thermal instabilities in the plane Poiseuille flow of weakly elastic fluids is examined through a linear stability analysis by taking into account the effects of viscous dissipation. The destabilizing thermal gradients may come from the different temperatures imposed on the external boundaries and/or from the volumetric heating induced by viscous dissipation. The rheological properties of the viscoelastic fluid are modeled using the Oldroyd-B constitutive equation. As in the Newtonian fluid case, the most unstable structures are found to be stationary longitudinal rolls (modes with axes aligned along the streamwise direction). For such structures, it is shown that the viscoelastic contribution to viscous dissipation may be reduced to one unique parameter: γ=λ1(1−Γ), where λ1 and Γ represent the relaxation time and the viscosity ratio of the viscoelastic fluid, respectively. It is found that the influence of the elasticity parameter γ on the linear stability characteristics is non-monotonic. The fluid elasticity stabilizes (destabilizes) the basic Poiseuille flow if γ<γ* (γ>γ*) where γ* is a particular value of γ that we have determined. It is also shown that when the temperature gradient imposed on the external boundaries is zero, the critical Reynolds number for the onset of such viscous dissipation/viscoelastic-induced instability may be well below the one needed to trigger the pure hydrodynamic instability in weakly elastic solutions.


Fluids ◽  
2021 ◽  
Vol 6 (10) ◽  
pp. 342
Author(s):  
Maksims Zigunovs ◽  
Andrei Kolyshkin ◽  
Ilmars Iltins

Linear stability analysis of a steady convective flow in a tall vertical annulus caused by nonlinear heat sources is conducted in the paper. Heat sources are generated as a result of a chemical reaction. The effect of radial cross-flow through permeable porous walls of the annulus is analyzed. The problem is relevant to biomass thermal conversion. The base flow solution is obtained by solving nonlinear boundary value problem. Linear stability analysis is performed, using collocation method. The calculations show that radial inward or outward flow has a stabilizing effect on the flow, while the increase in the Frank–Kamenetskii parameter (proportional to the intensity of the chemical reaction) destabilizes the flow. The increase in the Reynolds number based on the radial velocity leads to the appearance of the second minimum on the marginal stability curves. The rate of increase in the critical Grashof number with respect to the Reynolds number is different for inward and outward radial flows.


Entropy ◽  
2020 ◽  
Vol 22 (12) ◽  
pp. 1353
Author(s):  
Hirotaka Morimatsu ◽  
Takahiro Tsukahara

Direct numerical simulations were carried out with an emphasis on the intermittency and localized turbulence structure occurring within the subcritical transitional regime of a concentric annular Couette–Poiseuille flow. In the annular system, the ratio of the inner to outer cylinder radius is an important geometrical parameter affecting the large-scale nature of the intermittency. We chose a low radius ratio of 0.1 and imposed a constant pressure gradient providing practically zero shear on the inner cylinder such that the base flow was approximated to that of a circular pipe flow. Localized turbulent puffs, that is, axial uni-directional intermittencies similar to those observed in the transitional circular pipe flow, were observed in the annular Couette–Poiseuille flow. Puff splitting events were clearly observed rather far from the global critical Reynolds number, near which given puffs survived without a splitting event throughout the observation period, which was as long as 104 outer time units. The characterization as a directed-percolation universal class was also discussed.


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.


Sign in / Sign up

Export Citation Format

Share Document