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.

The interaction between a forward-facing step (FFS) and single-frequency Tollmien–Schlichting (TS) waves is investigated with experiments and two-dimensional (2-D) direct numerical simulations (DNS). Dedicated hot-wire anemometry and particle image velocimetry measurements in the vicinity of the FFS provide characterisation of the perturbation field, as well as validation of the DNS results. Comparison between experiments, 2-D DNS, and linear parabolised stability equations confirm the 2-D nature of the flow and the linearity of the instability mechanisms around the FFS. Upstream of the step, TS waves are gradually amplified by the increasing adverse pressure gradient. In the step vicinity, both mean flow and perturbation field exhibit abrupt distortion, with decoupling of the base flow-oriented growth rate components indicating significant non-modal evolution. Downstream of the step, the mean flow recovers to baseline conditions, but the perturbation field remains highly distorted. Linear stability theory results suggest the presence of superimposed modes on the original TS mode in this region. Despite their decay in the streamwise direction, their presence imprints modifications in the TS wave growth and shape, manifested as the tilting of the perturbation structure in and against the mean flow shear direction. This initiates a reversed Orr mechanism, characterised by a region of stabilisation followed by destabilisation further downstream. Eventually, the TS waves realign to their asymptotic (modal) behaviour. Overall, the FFS destabilises the TS wave far downstream. However, the streamwise extent and magnitude of the stabilisation downstream of the FFS remain significant.

The Parabolized Stability Equations (PSE), Adaptive Harmonic Linearized Navier-Stokes (AHLNS) and Harmonic Navier-Stokes (HNS) solvers are used to analyze the linear and nonlinear stability of swept-wing boundary layers under the influence of smooth wall deformations of varying size and geometry. Special attention is given to the validity of the slowly varying flow assumption of PSE via a comparison with AHLNS and HNS results. The surface deformations analyzed in this work are found to affect the development of the primary stationary crossflow instability mode as well as higher harmonics. Analysis of the locally most amplified mode reveals successive modulation of the growth rate in the vicinity of the surface deformation. This process was found to be largely governed by linear terms and driven by the base flow modification due to the deformed wall. Similarly, the base flow modification causes higher harmonics to experience a significant destabilization. This is followed by stabilization as nonlinear interactions become dominant. The PSE methodology proved capable of predicting the stability response for small wall deformations with only minor amplitude discrepancies compared to HNS results. The main difference was found to occur in the wall-normal velocity profiles of the mean flow distortion mode. The deviations of the PSE results compared to harmonic stability methods increased as the protuberance was made steeper. Moreover, the PSE framework was not able to converge for all cases nonlinearly.