In this work, the stability and transition to turbulence over blunt flat plates with different leading-edge radii are investigated computationally. The benchmark experimental work for comparative studies is conducted by Borovoy et al. (AIAA J., vol. 60, 2022, pp. 497–507). The freestream Mach number is 5, the unit Reynolds number is $6\times 10^7$ m$^{-1}$, and the maximum nose-tip radius 3 mm exceeds the experimental reversal value. High-resolution numerical simulation and stability analysis are performed. Three-dimensional broadband perturbation is added on the far field boundary to initiate the transition. The highlight of this work is that the complete physical process is considered, including the three-dimensional receptivity, linear and nonlinear instabilities, and transition. The experimental reversal phenomenon is reproduced favourably in the numerical simulation for the first time. Linear stability analysis shows that unstable first and second modes are absent in the blunt-plate flows owing to the presence of the entropy layer, although these modes are evident in the sharp-leading-edge case. Therefore, the transition on the blunt plate is due to non-modal instabilities. Numerical results for all the blunt-plate cases reveal the formation of streamwise streaky structures downstream of the nose (stage I) and then the presence of intermittent turbulent spots in the transitional region (stage II). In stage I, a preferential spanwise wavelength approximately 0.9 mm is selected for all the nose-tip radii, and low-frequency components are dominant. In stage II, high-frequency secondary instabilities appear to grow, which participate in the eventual breakdown. By contrast, leading-edge streaks are not remarkable in the sharp-leading-edge case, where transition is induced by oblique first and Mack second modes. The transition reversal beyond the critical nose-tip radius arises from an increasing magnitude of the streaky response in the early stage, while the transition mechanism stays similar qualitatively.

Ekman pumping is a phenomenon induced by no-slip boundary conditions in rotating fluids. In the context of Rayleigh–Bénard convection, Ekman pumping causes a significant change in the linear stability of the system compared with when it is not present (that is, stress-free). Motivated by numerical solutions to the marginal stability problem of the incompressible Navier–Stokes equation (iNSE) system, we seek analytical asymptotic solutions which describe the departure of the no-slip solution from the stress-free one. The substitution of normal modes into a reduced asymptotic model yields a linear system for which we explore analytical solutions for various scalings of wavenumber. We find very good agreement between the analytical asymptotic solutions and the numerical solutions to the iNSE linear stability problem with no-slip boundary conditions.

The dominant mode instability in hypersonic boundary-layer transition is the so-called second-mode instability, which has a peculiar nature strongly coupled with thermoacoustic phenomena. In linear stability theory, the unstable wave is associated with one of the two eigenvalues that originate from the acoustic branches, referred to as slow and fast modes. Interestingly, the unstable mode (slow or fast) reaches its maximum amplification as the other mode (fast or slow) attains a minimum. The phase velocity of the two modes is then very close, and this phenomenon is called synchronization. The aim of the present study is to unravel the physical mechanism that explains the second-mode growth. To that aim, second-order nonlinear equations are written for the disturbances given by linear stability. In this framework, entropy, kinetic energy and temperature energy budgets are obtained up to second order. The budgets are scrutinized for various Mach numbers and for adiabatic and cold-wall thermal conditions. Perturbation entropy budgets clearly show the process is a reversible one. An energy exchange between kinetic energy and temperature energy of the weakly nonlinear modes is driven by pressure–dilatation terms. As underlined in previous studies, the unstable mode experiences an alternate heating and cooling near the wall, which is shown to be a rather nonlinear process. The change in fluctuating thermal energy in the form of a dilatational wave is sustained by pumping disturbance kinetic energy through the pressure–dilatation term, the direction of the conversion being driven by the relative phase between pressure and dilatation. This process is similar for the slow and fast modes, the unstable mode being amplified and the other being damped. No change in the process has been noted at the location of the synchronization, suggesting that the modes have the same nature but evolve independently.

In this paper, we study the effect of lateral wall vibrations on the excitation and evolution of non-modal perturbations in hypersonic boundary layers subject to low-frequency freestream vortical disturbances (FSVDs). A novel, high-efficiency numerical approach, combining the harmonic weakly nonlinear Navier–Stokes and nonlinear parabolised stability equation approaches, is developed, which is sufficient to accommodate both the rapid distortion of the perturbation in the leading-edge vicinity and the nonlinear development of finite-amplitude high-order harmonics in the downstream region. The boundary-layer response to low-frequency FSVDs shows a longitudinal streaky structure, for which the temperature perturbation shows much greater magnitude than the streamwise velocity perturbation. The lateral vibration induces a Stokes layer solution for the spanwise velocity perturbation, which interacts with the FSVD-induced perturbations and leads to a suppression of the non-modal perturbation and an enhancement of the downstream modal perturbation. The new perturbations excited by the FSVD–vibration interaction strengthen as the vibration intensifies, and they could become comparable with the FSVD-induced perturbations in downstream locations at a high vibration intensity, indicating a remarkable modification of the streaky structure and its instability property. Secondary instability (SI) analyses based on the streaky base flow indicate that the vibration could enhance or suppress the SI modes, depending on their initial phases over the vibration period. Overall, the average effect is that the low-frequency and high-frequency SI modes are stabilised and destabilised by the vibration, respectively. Since the high-frequency SI modes undergo higher amplifications, the subsequent bypass transition is likely to be promoted by relatively strong vibrations.

The boundary-layer stability on a section of a rotating wind turbine blade with an FFA-W3 series aerofoil at a chord Reynolds number of $3 \times 10^5$, with varying rotation and radii, is studied with direct numerical simulations and linear stability analyses. Low rotation does not significantly affect transition in the outboard blade region. The relative insensitivity to rotation is due to a laminar separation bubble near the leading edge, spanwise-deformed by a primary self-excited instability, promoting the secondary absolute instability of the Kelvin–Helmholtz (KH) vortices and rapid transition. Moderate increases in rotation, or moving inboard, stabilise the flow by accelerating the attached boundary layer and possibly inducing competition between cross-flow and KH modes. This delays separation and transition. Initially, for high rotation rates or radial locations close to the hub, transition is delayed. Nevertheless, strong stationary and travelling cross-flow modes are eventually triggered, spanwise modulating the KH rolls and shifting the transition line close to the leading edge. Cross-flow velocities as high as $56\,\%$ of the free stream velocity directed towards the blade tip are reached at the transition location. For radial locations farther from the hub, the effective angle of attack is decreased, and cross-flow transition occurs at lower rotation rates. The advance or delay of the transition line compared with a non-rotating configuration depends on the competing rotation effects of stabilising the attached boundary layer and triggering cross-flow modes in the separation flow region.

The development of thermal boundary layers and plume near a section-triangular roof under different isothermal heating conditions has been the focus of numerous numerical studies. However, flow transition in this type of flow has never been observed experimentally. Here, phase-shifting interferometry and thermistor measurements are employed to experimentally observe and quantify the flow transitions in a buoyancy-driven flow over an isothermal section-triangular roof. Visualisation of temperature contours is conducted across a wide range of Rayleigh numbers from laminar at 103 to chaotic state at 4 × 106. Power spectral density of the temperature measurements reveals the type of bifurcations developing as the Rayleigh number is increased. This flow transition is characterised as a complex bifurcation route with the presence of two fundamental frequencies, a low and a high frequency. We found that the thermal stratification in the environment plays a significant role in the flow transition. The spatial development of flow is also quantitatively and qualitatively described. In addition to clarifying flow transition in experiments, the work demonstrates the implementation of phase-shifting interferometry and punctual temperature measurements for characterisation of near-field flow over a heated surface.

We perform a detailed numerical study of modal and non-modal stability in oblique Couette–Poiseuille profiles, which are among the simplest examples of three-dimensional boundary layers. Through a comparison with the Orr–Sommerfeld operator for the aligned case, we show how an effective wall speed succinctly characterizes modal stability. Large-scale parameter sweeps reveal that the misalignment between the pressure gradient and wall motion is, in general, destabilizing. For flows that are sufficiently oblique, the instability is found to depend exclusively on the direction of wall motion and not on its speed, a conclusion supported, in part, by the perturbation energy budget and the evolution of the critical layers. Closed forms for the critical parameters in this regime are derived using a simple analysis. From a non-modal perspective, pseudoresonance is examined through the resolvent and the $\epsilon$-pseudospectra. An analysis of the unforced initial value problem shows that the maximum energy gain is highly dependent on both the magnitude and direction of the wall velocity. However, the strongest amplification is always achieved for configurations that are only weakly skewed. Finally, the optimal perturbations appear to develop via a lift-up effect enhanced by an Orr-like mechanism, the latter driven by cross-flow shear.

We investigated the effects of two-dimensional sharp-edged rectangular bumps on Tollmien–Schlichting (TS) wave evolution using direct numerical simulation. The bump height, $h$, ranged from 5 % to 40 % of the local displacement thickness, $\delta ^*$. Behind the bump, a recirculating flow region could be formed whose length increased nonlinearly with $h$. The bump height effect on the TS wave, which was the dominant, scaled super-exponentially with $h$. We also showed a substantial effect of the $\delta ^*$-based Reynolds number, ${\textit {Re}} _{\delta ^*}$. Firstly, the bump wake extended with ${\textit {Re}} _{\delta ^*}$, promoting larger TS wave growth rates. The second effect is related to proximity to the upper branch of the instability loop, accounting for the influence of the TS frequency, as well. It dictates the bump impact increases as it gets closer to transition, either by the bump moving downstream or the transition moving upstream. For a 40 % high bump, for example, changing the ${\textit {Re}} _{\delta ^*}$ at the bump location from 1500 to 2000 increased $\Delta N$ by a factor of 2 ($\Delta N$ represents a measure of a surface irregularity effect on the smooth plate N-factor). We also found that $(\Delta N)_{max}$ increases linearly with ${\textit {Re}} _{hh}$. Results in the subsonic regime showed that the bump impact attenuates with Mach number up to 0.7 but above it, stabilisation is surpassed by the destabilising effect caused by the recirculation lengthening. This is mostly associated with the bump wake that extends with the pressure gradient which increases substantially towards the sonic speed. This is enhanced if the surface is adiabatic rather than isothermal.

A previous paper of the authors (Duck & Stephen, J. Fluid Mech., vol. 917, 2021, A56) considered the effect of three-dimensional, temporally periodic, linear and incompressible disturbances on a Blasius boundary layer, in particular when the disturbance wavelength is both comparable to and longer than the boundary-layer thickness. This previous study revealed that, unlike the two-dimensional counterpart, a mode exists that exhibits regimes of downstream spatial growth. In this paper we extend the analysis to the compressible regime, based on the boundary-region equations methodology. The aforementioned unstable mode is seen to persist into the compressible regime, and is studied using a combination of numerical and asymptotic methods. The paper adopts several approaches. First is a numerical approach in which the spatial development of the disturbances is assessed. This then leads to a consideration of the far-downstream behaviour, using (several) asymptotic limits. Of some note, in addition to unstable modes found in the incompressible case, is the existence of a further class of instability, not found in the incompressible case (which is also analysed asymptotically), corresponding to what amounts to an inviscid instability. The far-downstream analysis enables a (sub-)classification into entropy and non-entropy modes. The former, according to this analysis, are spatially damped, with one caveat, as revealed by our marching procedure, which highlights how spatial development of disturbances can be important.

Numerous studies have indicated that turbulence typically initiates along the boundary layer of the stationary disk within a rotor–stator cavity. To describe the transition process to turbulence on the stationary side of a closed rotor–stator cavity, a comprehensive approach combining global linear stability analysis with direct numerical simulation was adopted in the present study. The proposed model aligns with that of Yim et al. (J. Fluid Mech., vol. 848, 2018, pp. 631–647), who investigated the stability characteristics of the rotating-disk boundary layer in a rotor–stator cavity. In order to achieve a stable inflow for the stationary-disk boundary layer, we rotate the shroud together with the rotating disk. Through careful global stability analysis, the predominant spiral mode exhibiting the highest instability in the boundary layer of the stationary disk was discerned, corroborating observations from simulations. Initially, the spiral mode undergoes linear amplification, reaches a state of linear saturation and enters the nonlinear regime. Following nonlinear saturation in the flow field, a circular wave mode arises due to the influence of mean flow distortion. As the Reynolds number attained a sufficiently high level, the interplay between the downstream-propagating circular mode and spiral mode amplified disturbances in the boundary layer of the stationary disk, ultimately leading to the development of localised turbulence at the mid-radius of the rotor–stator cavity. Notably, the present study is the first to elucidate the coexistence of laminar–transitional–turbulent flow states in the stationary-disk boundary layer through direct numerical simulations.

We investigated the stability of the bottom boundary layer (BBL) beneath periodic internal solitary waves (ISWs) of depression over a flat bottom through two-dimensional direct numerical simulations. We explored the convective versus absolute/global nature of the BBL instability in response to changes in Reynolds number, and the sensitivity of the instability to seeding noise in the front of the ISW – spanning laboratory to geophysical scales. The BBL was laminar at $Re_{ISW}=90$ and convectively unstable at $Re_{ISW}=300$. At laboratory-scale $Re_{ISW}=300$, the convective wave packet was periodically amplified by each successive ISW, until vortex shedding occurred. The associated noise-amplification behaviour potentially explains the discrepancies of the critical $Re_{ISW}$ between the lock–release laboratory experiments and our Dubreil–Jacotin–Long-initialized numerical simulations as the result of the difference in background noise. Instability energy decreased under the front shoulder of the ISW, analogous to flow relaminarization under a favourable pressure gradient. At geophysical-scale $Re_{ISW}=900$, the BBL was initially convectively unstable, and then the instability tracked with the ISW, appearing phenomenologically similar to a global instability. The simulated initial convective instability at both $Re_{ISW}=300$ and $Re_{ISW}=900$ is in agreement with local linear stability analysis which predicts that the instability group speed is always lower than the ISW celerity. Increased free stream perturbations in front of the ISW and larger $Re_{ISW}$ shift the location of vortex shedding (and enhanced bed shear stress) beneath the wave, closer to the ISW trough, thereby potentially changing the location of maximum sediment resuspension, in agreement with field observations at higher $Re_{ISW}$.

The non-Oberbeck–Boussinesq effects on the stability of a vertical natural convection boundary layer are investigated using the linearised disturbance equations for air flows up to a temperature difference of $\Delta T=100\,{\rm K}$. Based on the linear stability results, the neutral curve is shown to be sensitive to the choice of reference temperature. When evaluated using the film temperature $T_f$, a lower film Grashof number is required to trigger the linear instability for larger $\Delta T$. The relative contributions of shear and buoyant production to the perturbation kinetic energy budget reveals that the marginally unstable modes are amplified based on different mechanisms: for lower wavenumbers at relatively small Grashof number, the instability is driven by buoyancy; whereas for higher wavenumbers and larger Grashof number, the flow becomes unstable due to a shear instability. The use of reference temperature is found to scale the shear- and buoyant-driven instabilities differently so that no single reference temperature definition would collapse the neutral curves. The linear stability result further demonstrates that at a given Grashof number a higher temperature difference would give a larger amplification rate of the perturbation, which then leads to an earlier onset of the nonlinearities when evaluated at $T_f$. Finally, by comparing the amplification rates obtained from direct numerical simulation and the linear stability results, the extent of the linear regime is determined for $\Delta T = 100\,{\rm K}$.

A global transient linear stability analysis of the three-dimensional time-dependent flow around an aerofoil undergoing small-amplitude pitching motion is performed using the optimally time-dependent (OTD) framework. The most salient linear instabilities associated with the instantaneous basic state are computed and tracked over time. The resulting OTD modes reflect the variations in the basic state and can be used as predictors of its spatial and temporal evolution, including the formation of a laminar separation bubble and its gradual spanwise modulation via primary global instability, leading to secondary instability and finally rapid breakdown to turbulence. The study confirms and expands upon earlier stability analyses of the same case based on the local properties of spanwise averaged velocity profiles in the bubble that predicted the onset of absolute instability soon followed by rapid breakdown of the separation bubble. The three-dimensional structure of the most unstable OTD mode is extracted, which compares well with both the locally absolutely unstable mode and the evolution of the basic state itself.

We study the stability of a zero-pressure gradient boundary layer subjected to free-stream disturbances by means of local stability analysis. The dataset under study corresponds to a direct numerical simulation (DNS) of a flat plate with a sharp leading edge in realistic wind tunnel conditions, with a turbulence level of 3.45 % at the leading edge. We present a method to track the convective evolution of the secondary instabilities of streaks by performing sequential stability calculations following the wave packet, connecting successive unstable eigenfunctions. A scattered nature, in time and space, of secondary instabilities is seen in the stability calculations. These instabilities can be detected before they reach finite amplitude in the DNS, preceding the nucleation of turbulent spots, and whose appearance is well correlated to the transition onset. This represents further evidence regarding the relevance of secondary instabilities of streaks in the bypass transition in realistic flow conditions. Consistent with the spatio-temporal nature of this problem, our approach allows us to integrate directly the local growth rates to obtain the spatial amplification ratio of the individual instabilities, where it is shown that instabilities reaching an $N$-factor in the range [2.5,4] can be directly correlated to more than 65 % of the nucleation events. Interestingly, it is found that high amplification is not only attained by modes with high growth rates, but also by instabilities with sustained low growth rates for a long time.

This paper reports on the efficacy of the Görtler number in scaling the laminar-turbulent boundary-layer transition on rotating cones facing axial inflow. Depending on the half-cone angle $\psi$ and axial flow strength, the competing centrifugal and cross-flow instabilities dominate the transition. Traditionally, the flow is evaluated by using two parameters: the local meridional Reynolds number $Re_l$ comparing the inertial versus viscous effects and the local rotational speed ratio $S$ accounting for the boundary-layer skew. We focus on the centrifugal effects, and evaluate the flow fields and reported transition points using Görtler number based on the azimuthal momentum thickness of the similarity solution and local cone radius. The results show that Görtler number alone dominates the late stages of transition (maximum amplification and turbulence onset phases) for a wide range of investigated $S$ and half-cone angle ($15^{\circ } \leq \psi \leq 50^{\circ }$), although the early stage (critical phase) seems to be not determined by the Görtler number alone on the broader cones ($\psi =30^{\circ }$ and $50^{\circ }$) where the primary cross-flow instability dominates the flow. Overall, this indicates that the centrifugal effects play an important role in the boundary-layer transition on rotating cones in axial inflow.

The growth and characteristics of linear, oblique instabilities on a highly swept fin on a straight cone in Mach 6 flow are examined. Large streamwise pressure gradients cause doubly inflected cross-flow profiles and reversed flow near the wall, which necessitates using the harmonic linearized Navier–Stokes equations. The cross-flow instability is responsible for the most-amplified disturbances, however, not all disturbances show typical cross-flow characteristics. Distinct differences in perturbation structure are shown between small ($\sim$3–5 mm) and large ($\sim$10 mm) wavelength disturbances at the unit Reynolds number $Re' = 11 \times 10^6$ m$^{-1}$. As a result, amplification measurements based solely on wall quantities bias a most-amplified disturbance assessment towards larger wavelengths and lower frequencies than would otherwise be determined by an off-wall total-energy approach. A spatial-amplification energy-budget analysis demonstrates (i) that wall-normal Reynolds-flux terms dictate the local growth rate, despite other terms having a locally larger magnitude and (ii) that the Reynolds-stress terms are responsible for large-wavelength disturbances propagating closer to the wall compared with small-wavelength disturbances. Additionally, the effect of free-stream unit Reynolds number and small yaw angles on the perturbation amplification and energy budget is considered. At a higher Reynolds number ($Re' = 22 \times 10^6$ m$^{-1}$), the most-amplified wavelength shrinks. Perturbations do not behave self-similarly in the thinner boundary layer, and the shift in most-amplified wavelength is due to decreased dissipation relative to the lower-Reynolds-number case. Small yaw angles produce a streamwise shift in the boundary layer and disturbance amplification. The yaw results quantify a potential uncertainty source in experiments and flight.

This paper investigates the linear and nonlinear evolution of radiating modes in supersonic boundary layers under the influence of impinging sound waves. It is found that the ensuing boundary-layer response is extraordinarily large for a subset of the sound frequency and incident angle, and the resonant over-reflection, corresponding to the reflection coefficient becoming infinite, occurs at a particular pairing of frequency and incident angle. At this point, the reflected wave coincides with a locally neutral radiating mode, which emits spontaneously sound in the form of Mach waves. A fundamental resonance takes place between the incident wave and the radiating mode. Viewed in a developing boundary layer, the response is rendered finite by introducing non-parallelism and nonlinear effects near the neutral location of the radiating mode, where the sound wave directly excites the radiating mode and/or acts on the pre-existing radiating mode. Inhomogeneous amplitude equations are derived to describe the excitation as well as the nonlinear development of the radiating mode in the two regimes where non-equilibrium and non-parallelism play a leading-order role, respectively. A composite amplitude equation is then constructed to take into account both non-parallelism and non-equilibrium effects. This amplitude equation is, with an appropriate initial condition, solved to quantify the impact of the impinging sound wave on the linear and nonlinear instability characteristics of the radiating mode. The far-field analysis shows that the Mach wave field of the radiating mode is changed significantly due to the incident sound.

This paper investigates the linear and nonlinear evolution of radiating modes under the influence of the spontaneously emitted Mach waves in a simple set-up of the supersonic boundary layers that develop in the entry region of a channel formed by two parallel semi-infinite flat plates. Two scenarios are considered. The first occurs in the boundary layers having identical wall conditions, where the Mach wave emitted by a radiating mode in one boundary layer influences the instability in the other. The second scenario takes place when the wall temperatures are different, in which case the spontaneously radiated Mach wave is reflected by the other boundary layer back to act on the radiating mode. Appropriate amplitude equations with the acoustic feedback effect being accounted for are derived. In each case, the effect of the spontaneously emitted sound contributes a linear term of delay type to the respective amplitude equation. For the first scenario, analytical and numerical studies of the amplitude equations show that due to the back action of the spontaneously radiated Mach wave, the amplitude exhibits rapid oscillations, and in the case of enhanced feedback effects, its envelope experiences near extinction followed by resurrection. The study of the coupled equations shows that the two modes with different initial amplitudes either undergo oscillations before attenuating, or terminate a finite-distance singularity at different locations. For the second scenario, the acoustic feedback produces similar effects in a broad range of wall temperature. The effects become pronounced, and the dependence on the wall temperature becomes more sensitive when the latter approaches the value corresponding to the resonance. Estimates suggest that such acoustic feedback is likely to be present in typical wind tunnel experiments and models for scramjet combustors.

Boundary-layer disturbances are analysed on a $5^{\circ }$ half-angle blunted cone in Mach 5, high-enthalpy flow ($h_0 = 9\ {\rm MJ}\ {\rm kg}^{-1}$) with a low wall-to-edge temperature ratio, $T_w/T_e = 0.18$. Schlieren and focused laser differential interferometry (FLDI) are utilized to assess the structures and frequency content associated with disturbances. Wave packets are identified from bursts of modal content on time-resolved spectrograms. Bandpass filtering, proper orthogonal decomposition (POD) and space–time POD are then applied to the schlieren data. Bandpass filtering suggests the presence of wave packet dispersion and elongation indicative of slow-acoustic-wave synchronization. Modal reconstruction techniques indicate the radiation of content outside the boundary layer and distinct orientation changes within disturbances, potentially the first experimental evidence of the supersonic-mode instability in such a flow field. Cross-bicoherence computations are carried out for discrete time segments of data from both schlieren and FLDI data. They demonstrate that the most dominant nonlinear interactions are the fundamental–first-harmonic and the fundamental–low-frequency interactions.