Stationary velocity-perturbation streaks have recently been identified in laminar swept-wing boundary-layer flow interacting with a surface forward-facing step. Streaky structures at the step promote early laminar-turbulent transition under certain conditions. This work utilizes direct numerical simulations to explore the mechanisms of growth of stationary streaks at the step and provides insight into their origin, nature, and spatial organization. The analysis is mainly focused on, but not restricted to, incoming perturbations in the form of stationary crossflow instability. Stationary streaky structures are found to be universal to swept forward-facing-step flow subjected to three-dimensional perturbations in the incoming boundary layer. The streaks at the step are primarily ascribed to the lift-up effect. They appear as a linear perturbation response of the highly sheared step flow to the cross-stream pattern of incoming perturbations. A mechanism of base-flow deceleration additionally contributes to feeding growth to the streaks. Linear stability analysis carried out through the harmonic Navier-Stokes method confirms that the streaks are a linear perturbation phenomenon. Effects of spanwise perturbation wavelength and effective sweep angle on the mechanisms of the streaks are also assessed.

Linear stability analyses are performed to study the dynamics of linear convective instability mechanisms in a laminar shock-wave/boundary-layer interaction at Mach 1.7. In order to account for all two-dimensional gradients elliptically, we introduce perturbations into an initial-value problem that are found as solutions to an eigenvalue problem formulated in a moving frame of reference. We demonstrate that this methodology provides results that are independent of the numerical setup, frame speed, and type of eigensolutions used as initial conditions. The obtained time-integrated wave packets are then Fourier-transformed to recover individual-frequency amplification curves. This allows us to determine the dominant spanwise wavenumber and frequency yielding the largest amplification of perturbations in the shock-induced recirculation bubble. By decomposing the temporal wave-packet growth rate into the physical energy-production processes, we provide an in-depth characterization of the convective instability mechanisms in the shock-wave/boundary-layer interaction. For the particular case studied, the largest growth rate is achieved in the near-vicinity of the bubble apex due to the wall-normal (productive) and streamwise (destructive) Reynolds-stress energy-production terms. We also observe that the Reynolds heat-flux effects are similar but contribute to a smaller extent.

Micro-ramps are popular passive flow control devices which can delay flow separation by re-energising the lower portion of the boundary layer. We compute the laminar base flow, the instantaneous transitional flow, and the mean flow around a micro-ramp immersed in a quasi-incompressible boundary layer at supercritical roughness Reynolds number. Results of our Direct Numerical Simulations (DNS) are compared with results of BiLocal stability analysis on the DNS base flow and independent tomographic Particle Image Velocimetry (tomo-PIV) experiments. We analyse relevant flow structures developing in the micro-ramp wake and assess their role in the micro-ramp functionality, i.e., in increasing the near-wall momentum. The main flow feature of the base flow is a pair of streamwise counter-rotating vortices induced by the micro-ramp, the so-called primary vortex pair. In the instantaneous transitional flow, the primary vortex pair breaks up into large-scale hairpin vortices, which arise due to linear varicose instability of the base flow, and unsteady secondary vortices develop. Instantaneous vortical structures obtained by DNS and experiments are in good agreement. Matching linear disturbance growth rates from DNS and linear stability analysis are obtained until eight micro-ramp heights downstream of the micro-ramp. For the setup considered in this article, we show that the working principle of the micro-ramp is different from that of classical vortex generators; we find that transitional perturbations are more efficient in increasing the near-wall momentum in the mean flow than the laminar primary vortices in the base flow.

A combined experimental and numerical approach to the analysis of the secondary stability of realistic swept-wing boundary layers is presented. Global linear stability theory is applied to experimentally measured base flows. These base flows are three-dimensional laminar boundary layers subject to spanwise distortion due to the presence of primary stationary crossflow vortices. A full three-dimensional description of these flows is accessed through the use of tomographic particle image velocimetry (PIV). The stability analysis solves for the secondary high-frequency modes of type I and type II, ultimately responsible for turbulent breakdown. Several pertinent parameters arising from the application of the proposed methodology are investigated, including the mean flow ensemble size and the measurement domain extent. Extensive use is made of the decomposition of the eigensolutions into the terms of the Reynolds-Orr equation, allowing insight into the production and/or destruction of perturbations from various base flow features. Stability results demonstrate satisfactory convergence with respect to the mean flow ensemble size and are independent of the handling of the exterior of the measurement domain. The Reynolds-Orr analysis reveals a close relationship between the type I and type II instability modes with spanwise and wall-normal gradients of the base flow, respectively. The structural role of the in-plane velocity components in the perturbation growth, topology and sensitivity is identified. Using the developed framework, further insight is gained into the linear growth mechanisms and later stages of transition via the primary and secondary crossflow instabilities. Furthermore, the proposed methodology enables the extension and enhancement of the experimental measurement data to the pertinent instability eigenmodes. The present work is the first demonstration of the use of a measured base flow for stability analysis applied to the swept-wing boundary layer, directly avoiding the modelling of the primary vortices receptivity processes.

Selective Frequency Damping (SFD) is a popular method for the computation of globally unstable steady-state solutions in fluid dynamics. The approach has two model parameters whose selection is generally unclear. In this article, a detailed analysis of the influence of these parameters is presented, answering several open questions with regard to the effectiveness, optimum efficiency and limitations of the method. In particular, we show that SFD is always capable of stabilising a globally unstable systems ruled by one unsteady unstable eigenmode and derive analytical formulas for optimum parameter values. We show that the numerical feasibility of the approach depends on the complex phase angle of the most unstable eigenvalue. A numerical technique for characterising the pertinent eigenmodes is presented. In combination with analytical expressions, this technique allows finding optimal parameters that minimise the spectral radius of a simulation, without having to perform an independent stability analysis. An extension to multiple unstable eigenmodes is derived. As computational example, a two-dimensional cylinder flow case is optimally stabilised using this method. We provide a physical interpretation of the stabilisation mechanism based on, but not limited to, this Navier–Stokes example.

As the community investigates more complex flows with stronger streamwise variations and uses more physically inclusive stability techniques, such as BiGlobal theory, there is a perceived need for more accuracy in the base flows. To this end, the implication is that using these more advanced techniques, we are now including previously neglected terms of O(Re²). Two corresponding questions follow: (1) how much accuracy can one reasonably achieve from a given set of basic-state equations and (2) how much accuracy does one need to converge more advanced stability techniques? The purpose of this paper is to generate base flow solutions to successively higher levels of accuracy and assess how inaccuracies ultimately affect the stability results. Basic states are obtained from solving the self-similar boundary-layer equations, and stability analyses with LST, which both share O(1/Re) accuracy. This is the first step toward tackling the same problem for more complex basic states and more advanced stability theories. Detailed convergence analyses are performed, allowing to conclude on how numerical inaccuracies from the basic state ultimately propagate into the stability results for different numerical schemes and instability mechanisms at different Mach numbers.