Absolute/convective instabilities and the convective Mach number in a compressible mixing layer

Two aspects of the stability of a compressible mixing layer: Absolute/Convective instability and the convective Mach number were considered. It was shown that, for Mach numbers less than one, the compressible mixing layer is convectively unstable unless there is an appreciable amount of backflow. Also presented was a rigorous derivation of a convective Mach number based on linear stability theory for the flow of a multi-species gas in a mixing layer. The result is compared with the heuristic definitions of others and to selected experimental results

Viscous effects in the absolute-convective instability of the Batchelor vortex

International audienceThe effects of viscosity on the instability properties of the Batchelor vortex are investigated. The characteristics of spatially amplified branches are first documented in the convectively unstable regime for different values of the swirl parameter q and the co-flow parameter a at several Reynolds numbers Re. The absolute-convective instability transition curves, determined by the Briggs-Bers zero-group velocity criterion, are delineated in the (a,q)-parameter plane as a function of Re. The azimuthal wavenumber m of the critical transitional mode is found to depend on the magnitude of the swirl q and on the jet (a > -0.5) or wake (a < -0.5) nature of the axial flow. At large Reynolds numbers, the inviscid results of Olendraru et al. (1999) are recovered. As the Reynolds number decreases, the pocket of absolute instability in the (a,q)-plane is found to shrink gradually. At Re = 667, the critical transitional modes for swirling jets are m = -2 or m = -3 and absolute instability prevails at moderate swirl values even in the absence of counterflow. For higher swirl levels, the bending mode m = -1 becomes critical. The results are in good overall agreement with those obtained by Delbende et al. (1998) at the same Reynolds number. However, a bending (m = +1) viscous mode is found to partake in the outer absolute-convective instability transition for jets at very low positive levels of swirl. This asymmetric branch is the spatial counterpart of the temporal viscous mode isolated by Khorrami (1991) and Mayer and Powell (1992). At Re = 100, the critical transitional mode for swirling jets is m = -2 at moderate and high swirl values and, in order to trigger an absolute instability, a slight counterflow is always required. A bending (m = +1) viscous mode again becomes critical at very low swirl values. For wakes (a < -0.5) the critical transitional mode is always found to be the bending mode m = -1, whatever the Reynolds number. However, above q = 1.5, near-neutral centre modes are found to define a tongue of weak absolute instability in the (a,q)-plane. Such modes had been analytically predicted by Stewartson and Brown (1985) in a strictly temporal inviscid framework

Mode selection in swirling jet experiments: A linear stability analysis

International audienceThe primary goal of the study is to identify the selection mechanism responsible for the appearance of a double-helix structure in the pre-breakdown stage of so-called screened swirling jets for which the circulation vanishes away from the jet. The family of basic flows under consideration combines the azimuthal velocity profiles of Carton & McWilliams (1989) and the axial velocity profiles of Monkewitz (1988). This model satisfactorily represents the nozzle exit velocity distributions measured in the swirling jet experiment of Billant et al. (1998). Temporal and absolute/convective instability properties are directly retrieved from numerical simulations of the linear impulse response for different swirl parameter settings. A large range of negative helical modes, winding with the basic flow, are destabilized as swirl is increased, and their characteristics for large azimuthal wavenumbers are shown to agree with the asymptotic analysis of Leibovich & Stewartson (1983). However, the temporal study fails to yield a clear selection principle. The absolute/convective instability regions are mapped out in the plane of the external axial flow and swirl parameters. The absolutely unstable domain is enhanced by rotation and it remains open for arbitrarily large swirl. The swirling jet with zero external axial flow is found to first become absolutely unstable to a mode of azimuthal wavenumber m = -2, winding with the jet. It is suggested that this selection mechanism accounts for the experimental observation of a double-helix structure

Analysis of instability patterns in non-Boussinesq mixed convection using a direct numerical evaluation of disturbance integrals

The Fourier integrals representing linearised disturbances arising from an initially localised source are evaluated numerically for natural and mixed convection flows between two differentially heated plates. The corresponding spatio-temporal instability patterns are obtained for strongly non-Boussinesq high-temperature convection of air and are contrasted to their Boussinesq counterparts. A drastic change in disturbance evolution scenarios is found when a large cross-channel temperature gradient leads to an essentially nonlinear variation of the fluid's transport properties and density. In particular, it is shown that non-Boussinesq natural convection flows are convectively unstable while forced convection flows can be absolutely unstable. These scenarios are opposite to the ones detected in classical Boussinesq convection. It is found that the competition between two physically distinct instability mechanisms which are due to the action of the shear and the buoyancy are responsible for such a drastic change in spatio-temporal characteristics of instabilities. The obtained numerical results confirm and complement semi-analytical conclusions of Suslov 2007 on the absolute/convective instability transition in non-Boussinesq mixed convection. Generic features of the chosen numerical approach are discussed and its advantages and shortcomings are reported

Condition for convective instability of dark solitons

Simple derivation of the condition for the transition point from absolute instability of plane dark solitons to their convective instability is suggested. It is shown that unstable wave packet expands with velocity equal to the minimal group velocity of the disturbance waves propagating along a dark soliton. The growth rate of the length of dark solitons generated by the flow of Bose-Einstein condensate past an obstacle is estimated. Analytical theory is confirmed by the results of numerical simulations

Inflow/outflow boundary conditions and global dynamics of spatial mixing layers

The numerical simulation of incompressible spatially-developing shear flows poses a special challenge to computational fluid dynamicists. The Navier-Stokes equations are elliptic and boundary equations need to be specified at the inflow and outflow boundaries in order to compute the fluid properties within the region of interest. It is, however, difficult to choose inflow and outflow conditions corresponding to a given experimental situation. Furthermore the effects that changes in the boundary conditions or in the size of the computational domain may induce on the global dynamics of the flow are presently unknown. These issues are examined in light of recent developments in hydrodynamic stability theory. The particular flow considered is the spatial mixing layer but it was expected that similar phenomena were bound to occur in other cases such as channel flow, the boundary layer, etc. A short summary of local/global and absolute/convective instability concepts is given. The results of numerical simulations are presented which strongly suggest that global resonances may be triggered in domains of finite streamwise extent although the evolution of the perturbation vorticity field is everywhere locally convective. A relationship between finite domains and pressure sources which might help in devising a scheme to eliminate these difficulties is discussed

On the convectively unstable nature of optimal streaks in boundary layers

International audienceThe objective of the study is to determine the absolute/convective nature of the secondary instability experienced by finite-amplitude streaks in the flat-plate boundary layer. A family of parallel streaky base flows is defined by extracting velocity profiles from direct numerical simulations of nonlinearly saturated optimal streaks. The computed impulse response of the streaky base flows is then determined as a function of streak amplitude and streamwise station. Both the temporal and spatio-temporal instability properties are directly retrieved from the impulse response wave packet, without solving the dispersion relation or applying the pinching point criterion in the complex wavenumber plane. The instability of optimal streaks is found to be unambiguously convective for all streak amplitudes and streamwise stations. It is more convective than the Blasius boundary layer in the absence of streaks; the trailing edge-velocity of a Tollmien-Schlichting wave packet in the Blasius boundary layer is around 35% of the free-stream velocity, while that of the wave packet riding on the streaky base flow is around 70%. This is because the streak instability is primarily induced by the spanwise shear and the associated Reynolds stress production term is located further away from the wall, in a larger velocity region, than for the Tollmien-Schlichting instability. The streak impulse response consists of the sinuous mode of instability triggered by the spanwise wake-like profile, as confirmed by comparing the numerical results with the absolute/convective instability properties of the family of two-dimensional wakes introduced by Monkewitz (1988). The convective nature of the secondary streak instability implies that the type of bypass transition studied here involves streaks that behave as amplifiers of external noise