The recirculating flow at the rear of a flat-base three-dimensional body with ground proximity is investigated for different body attitudes defined by the pitch angle varying in the range $-1.5^\circ \lt \alpha \lt +2.6^\circ$ and the yaw angle in the range $0^\circ \lt \beta \lt +12^\circ$. Experiments measuring the three components of the mean velocity field in two perpendicular planes intersecting the recirculation area as well as the base pressure distribution are conducted for 50 different attitudes. They provide a clear correlation between the orientation of the spatially averaged reversed flow and the gradient at the centre of the base pressure distribution. Both vectors are found to be in the same so-called w-plane, that is perpendicular to the base of which the azimuthal position changes with the body attitude due to either the flow orientation at the base separation or sometimes to a ground separation for large nose-up pitch. Numerical simulations of the same geometry realised for 10 attitudes show satisfactory agreement with the force coefficients measured in the experiment. Base flow variations induced by attitude changes are also well captured, particularly that of the w-plane. The full three-dimensional simulation data are used to show that the inner structure of the separation bubble is always a tilted recirculation torus, where the tilt orientation is given by the base pressure gradient. At the bubble closure, a pair of longitudinal vortices symmetrically located on both sides of the w-plane are permanently observed with circulations consistent with the circulation of the dividing streamline separation in the w-plane.

The inviscid mechanism, driving flow instabilities in a $1:3$, planar and symmetric sudden expansion, is discerned through a sensitivity-based protocol, also referred to as inviscid structural sensitivity analysis, with a specific focus on the onset and nature of the secondary instability. The fundamental idea of this methodology is to change the contribution of viscosity solely in the global stability equations, while freezing the base-flow field at the critical conditions. This is practically implemented by decoupling the Reynolds number that serves as the control parameter for determining the steady base flow from that governing the disturbance evolution, in order to repeat the structural sensitivity analysis while progressively increasing the Reynolds number in the linearised equations only. Accordingly, the sequence of structural sensitivity maps enables us to highlight the flow regions where the inviscid instability mechanism acts. The numerical results reveal that the classical structural sensitivity analysis accurately locates the wavemaker region within the primary recirculation zone, but only its inviscid limit can unveil that the core of the instability coincides with the centre of the primary vortex: a hallmark of an elliptic instability. To validate the global findings, the results of the inviscid structural sensitivity analysis are compared with those obtained from geometric optics. The agreement of the two approaches confirms the inviscid character of the instability, thereby providing a complete picture of the nature of the secondary bifurcation.

We discuss flow-induced vibrations of an equilateral triangular prism confined to travel on a circular path when placed in the concave or convex orientations with respect to the flow. In each orientation, we consider three different initial angles for the prism. In Case 1, one side of the prism sees the flow first; in Case 2, one sharp edge sees the flow first; and in Case 3, one side of the prism is parallel to the incoming flow. We show that the response of the structure as well as the observed wake depend heavily on both the orientation and the initial angle of the prism. Case 1 exhibits vortex-induced vibration (VIV) in the concave orientation and galloping in the convex orientation. Case 2 does not oscillate in the concave orientation; however, oscillates about a mean deflection after a critical reduced velocity in the convex orientation. Case 3 exhibits small-amplitude oscillations in the concave orientation about a mean deflection, while in the convex orientation, exhibits VIV at low reduced velocities, followed by an asymmetric response with VIV features in a half-cycle and galloping features in the other half, and divergence at higher reduced velocities. These different types of responses are accompanied by a myriad of vortex patterns in the wake, from two single vortices shed in the wake in each cycle of oscillations to two vortex pairs, two sets of co-rotating vortices, and a combination of single vortices and vortex pairs depending on the prism’s orientation and its initial angle.

Transonic buffet is a complex and strongly nonlinear unstable flow sensitive to variations in the incoming flow state. This poses great challenges for establishing accurate-enough reduced-order models, limiting the application of model-based control strategies in transonic buffet control problems. To address these challenges, this paper presents a time-variant modelling approach that incorporates rolling sampling, recursive parameter updating and inner iteration strategies under dynamic incoming flow conditions. The results demonstrate that this method successfully overcomes the difficulty in designing appropriate training signals and obtaining unstable steady base flow. Additionally, it improves the global predictive capability and identification efficiency of linear models for nonlinear flow-system responses by more than one order of magnitude. Furthermore, two adaptive control strategies – minimum variance control and generalised predictive control – are validated as effective based on the time-variant reduced-order model through numerical simulations of the transonic buffet flow over the NACA 0012 aerofoil. The adaptive controllers effectively regulate the unstable eigenvalues of the flow system, achieving the desired control outcomes. They ensure that the shock wave buffet phenomenon does not recur after control is applied, and that the actuator deflection, specifically the trailing-edge flap, returns to zero. Moreover, the control results further confirm the global instability essence of transonic buffet flow from a control perspective, thereby deepening the cognition of this nonlinear unstable flow.

The effects of Reynolds number across ${\textit{Re}}=1000$, $2500$, $5000$ and $10\,000$ on separated flow over a two-dimensional NACA0012 airfoil at an angle of attack of $\alpha =14^\circ$ are investigated through biglobal resolvent analysis. We identify modal structures and energy amplifications over a range of frequencies, spanwise wavenumbers, and values of the discount parameter, providing insights across various time scales. Using temporal discounting, we find that the shear-layer dynamics dominates over short time horizons, while the wake dynamics becomes the primary amplification mechanism over long time horizons. Spanwise effects also appear over long time horizons, sustained by low frequencies. The low-frequency and high-wavenumber structures are found to be dominated by elliptic mechanisms within the recirculation region. At a fixed angle of attack and across the Reynolds numbers, the response modes shift from wake-dominated structures at low frequencies to shear-layer-dominated structures at higher frequencies. The frequency at which the dominant mechanism changes is independent of the Reynolds number. Comparisons at a different angle of attack ($\alpha =9^\circ$) show that the transition from wake to shear-layer dynamics with increasing frequency only occurs if the unsteady flow is three-dimensional. We also study the dominant frequencies associated with wake and shear-layer dynamics across the angles of attack and Reynolds numbers, and confirm characteristic scaling laws from the literature.

This numerical investigation focuses on the mechanisms, flow topology and onset of Kelvin–Helmholtz instabilities (KHIs), that drive the leading-edge shear-layer destabilisation in the wake of wall-mounted long prisms. Large-eddy simulations are performed at ${\textit{Re}} = 2.5\times 10^3, 5\times 10^3$ and $1\times 10^4$ for prisms with a range of aspect ratio (AR, height-to-width) between $0.25$ and $1.5$, and depth ratios (DR, length-to-width) of $1{-}4$. Results show that shear-layer instabilities enhance flow irregularity and modulate spanwise vortex structures. The onset of KHI is strongly influenced by depth ratio, such that long prisms (${\textit{DR}}= 4$) experience earlier initiation compared with shorter ones (${\textit{DR}}= 1$). At higher Reynolds numbers, the onset of KHI shifts upstream towards the leading-edge, intensifying turbulence kinetic energy and increasing flow irregularity, especially for long prisms. The results further show that in this configuration, energy transfer from the secondary recirculation region contributes to the destabilisation of the leading-edge shear layer by reinforcing low-frequency modes. A feedback mechanism is identified wherein energetic flow structures propagate upstream through reverse boundary-layer flow, re-energising the leading-edge shear layer. Quantification using probability density functions reveals rare, intense upstream energy convection events, driven by this feedback mechanism. These facilitate the destabilisation process regardless of Reynolds number. This study provides a comprehensive understanding of the destabilisation mechanisms for leading-edge shear layers in the wake of wall-mounted long prisms.

The turbulent transport of momentum, energy and passive scalar is investigated in the flow around a rectangular cylinder of aspect ratio 5 : 1 – a geometry representative of separating and reattaching flows from sharp-edged bodies. The study is based on direct numerical simulation (DNS) conducted at Reynolds numbers up to ${\textit{Re}} = 14\,000$, based on the cylinder thickness, with Schmidt number fixed at ${\textit{Sc}} = 0.71$. At this Reynolds number, the flow exhibits features of asymptotic high-${\textit{Re}}$ behaviour. Budgets of mean momentum, Reynolds stresses, mean scalar and scalar fluxes provide a detailed view of the underlying transport mechanisms. The mean momentum balance elucidates the role of turbulence in entraining free stream fluid, promoting shear-layer reattachment, sustaining backflow in the recirculation region and regulating wake dynamics through large-scale vortex shedding. The leading-edge shear layer is the main site of turbulence production, with energy injected into streamwise fluctuations and redistributed to cross-flow components by pressure–strain interactions. As ${\textit{Re}}$ increases, vertical fluctuations increasingly return energy to the mean upward flow, stabilising the separation bubble height. Turbulent transport dominates scalar redistribution. Scalar fluxes are primarily generated by interactions between Reynolds stresses and scalar gradient, and modulated by pressure-scalar gradient effects. An a priori evaluation of eddy-viscosity and diffusivity models quantifies the misalignment between modelled and DNS-resolved stress and flux tensors, as well as the inhomogeneity of eddy transport coefficients. This analysis deepens the understanding of transport phenomena in bluff-body flows approaching the asymptotic regime, and underpins the validation and improvement of turbulence models for separating and reattaching flows.

Natural fliers and marine swimmers twist and turn their lifting or control surfaces to manipulate the unsteady forces experienced in air and water. The passive deformation of such surfaces has been investigated by several researchers, but the aspect of controlled deformation has received comparatively less attention. In this paper, we experimentally measure the forces and the flow fields of a flat-plate wing (aspect ratio (AR) = 3), translating at a constant Reynolds number (Re) of 10 000, with a dynamically twisting span. We show that the unsteady forces can be dependably estimated by a three-dimensional discrete vortex model. In this model, we account for the leading-edge separation with the help of the leading-edge-suction parameter. Experiments are conducted for two angles of attack (AoAs), $5^\circ$ and $15^\circ$. In addition, two rates of twisting are implemented where part of the leading edge, closer to the tip region, is twisted away from the incoming flow, increasing the effective AoA. The results show that twisting away from the flow augments the lift forces in all cases, although the rate of increase of lift is higher for the highest twist rate. The act of twisting causes an increase in effective AoA beyond the static stall angle in the AoA $=15^\circ$ case. This is highlighted by a distinct dip in the force data following the initial rise after twisting is activated. The increase in effective AoA from the reference case (without twisting) causes separation of the flow below the mid-span. This, in turn, creates higher levels of vorticity in those regions and results in a leading-edge vortex with increased cross-section and strength when compared with the reference case without twisting. Finally, we apply force partitioning and reveal that dynamic twisting leads to a localised increase in vorticity-induced forces along the twisted part of the span, which is approximately twice that of the untwisted case.

Over-expansion flow can generate asymmetric shock wave interactions, which lead to significant lateral forces on a nozzle. However, there is still a lack of a suitable theory to explain the phenomenon of asymmetry. The current work carefully investigates the configurations of shock wave interactions in a planar nozzle, and proposes a theoretical method to analyse the asymmetry of over-expansion flows. First, various possible flow patterns of over-expansion flows are discussed, including regular and Mach reflections. Second, the free interaction theory and the minimum entropy production principle are used to analyse the boundary layer flow and main shock wave interactions, establish the relationship between the separation shock strength and separation position, and predict asymmetric configurations. Finally, experiments are conducted to validate the theoretical method, and similar experiments from other studies are discussed to demonstrate the effectiveness of the proposed method. Results demonstrate that the direction of asymmetric over-expansion flow is random, and the separated flow strives to adopt a pattern with minimal total pressure loss. Asymmetric interaction is a mechanism through which the flow can achieve a more efficient thermodynamic balance by minimising entropy production.

We analyse the small-scale characteristics, such as enstrophy, total strain and normality/non-normality, in the three-dimensional, separated flow around a NACA 0018 wing using direct numerical simulations. The angle of attack is $10^\circ$ and the Reynolds number (based on the chord length) is $Re_c=5000$. The role of non-normality is investigated by performing Schur decomposition of the velocity gradient tensor. We also apply the Schur decomposition to derive new expressions for the production of enstrophy and total strain arising from the mean flow inhomogeneity. We focus on two sections of the flow, across the recirculating zone and along the transitioning shear layer, and compare our results with homogeneous isotropic turbulence (HIT). Within the recirculating region, the non-normality index is approximately 0 (and close to the HIT value), indicating almost equal normal and non-normal contributions. However, in the separating layer non-normal effects strongly dominate, especially in the region of kinetic energy growth. Only in the decay region do the values of the non-normality index gradually approximate HIT values. The production of enstrophy due to vortex stretching is dominated by the mixed (interaction) term, where normal strain stretches non-normal vorticity. The same component also dominates the strain self-amplification term. The contributions of different QR regions to the production terms are also examined. Production due to mean strain rate is triggered upstream compared with production due to fluctuating strain fields.

This study investigates the energy exchange between coherent structures in flows over four low-aspect-ratio (low-) plates using the tomographic particle image velocimetry dataset originally obtained by Zhu et al. (2024. J. Fluid Mech. 983, A35). The chord-based Reynolds number is $5400$, with fixed angle of attack $6 ^\circ$. In this study, multiscale proper orthogonal decomposition is applied to extract the coherent structures, including those associated with the vortex-shedding frequency $St_1$ and its subharmonic counterpart. Subsequently, the coherent kinetic energy budget is analysed with a focus on inter-scale energy transfer. This study demonstrates that the energy transfer between the scales centred at $St_1$ and $0.5\,St_1$ can exhibit a reverse or forward direction, depending on the transformation pattern of the leading-edge vortices (LEVs). Specifically, different triadic interactions are excited during the LEV transformation, and manifest themselves during the formation of hairpin vortices downstream. Understanding this nonlinear energy transfer is essential for elucidating mechanisms underlying the development of turbulence in three-dimensional flows over low- plates.

The flow over a cambered NACA 65(1)–412 airfoil at $Re\,=\,2\times 10^4$ is described based on a high-order direct numerical simulation. Simulations are run over a range of angles of attack, $\alpha$, where a number of instabilities in the unsteady, three-dimensional flow field are identified. The balance and competing effects of these instabilities are responsible for significant and abrupt (with respect to $\alpha$) changes in flow regime, with measurable consequences in time-averaged, integrated force coefficients, and in the far-wake footprint. At low $\alpha$, the flow is strongly influenced by vortex roll-up from the pressure side at the trailing edge. The interaction of this large-scale structure with shear and three-dimensional modal instabilities in the separated shear layer and associated wake region on the suction side, explains the transitions and bifurcations of the the flow states as $\alpha$ increases. The transition from a separation at low $\alpha$ to reattachment and establishment of a laminar separation bubble at the trailing edge at critical $\alpha$ is driven by instabilities within the separated shear layer that are absent at lower angles. Instabilities of different wavelengths are then shown to pave the path to turbulence in the near wake.

This paper presents a numerical investigation of the turbulence transition phenomenon in the wake of wall-mounted prisms. Large-eddy simulations are performed at $Re = 1\times 10^3 {-}5\times 10^3$ for prisms with a range of aspect ratio (height to width) from $0.25$ to $1.5$, and depth ratios (length to width) between $1$ and $4$. The results show that the wake irregularity is enhanced with increasing depth ratio, evidenced by higher turbulent kinetic energy (${\approx}90\,\%$) near the leading edge, and the onset of irregular, unsteady vortex shedding. This is attributed to interactions between Kelvin–Helmholtz instability (KHI) of the shear layer and large-scale vortex shedding, and it is induced by an unsteady shear layer, resembling flapping-like motion. These interactions elevate the flow momentum due to increased turbulence intensity and mixing, contributing to the wake transition phenomenon. To this end, this study defines the role of depth ratio in the transition phenomenon by showing that increasing depth ratio (e.g. from $1$ to $4$) leads to earlier onset of KHIs in the shear layer. These instabilities intensify with depth ratio, resulting in stronger interactions between shear layer and large-scale vortex shedding. Specifically, KHI-induced vortices interact more frequently with large-scale wake structures for higher depth ratio prisms, exciting larger flow fluctuations and irregular wake patterns. This interaction alters the frequency and coherence of vortex shedding, revealing a complex coupling mechanism that drives the transition to turbulence.

Flow field in the near wake of a short-finite circular cylinder at $L/D=1.0$ with an angle of attack between 0$^\circ$–15$^\circ$, where the transition from the non-reattaching flow to the reattaching flow appears, is investigated in wind tunnel tests with a supportless condition. Stereo particle image velocimetry measurements were applied to the experiments at the Reynolds number of $3.46\times 10^4$, and velocity fields in the near wake were obtained. The data was mainly analysed using spectral proper orthogonal decomposition. Characteristic large-scale wake structures of recirculation bubble pumping and large-scale vortex shedding were identified in the near wake of the cylinder regardless of the angle of attack. The phase difference of expansion and contraction of the recirculation flow appears in the recirculation bubble pumping at $\alpha \neq 0^\circ$. On the other hand, the eigenfunctions of velocity fluctuations at the vortex shedding frequency show a similar spatial pattern regardless of $\alpha$. Frequency analyses of wake position calculated from the reconstructed velocity field clarified that peak frequency is different between two in-plane directions when $\alpha \neq 0^\circ$. In addition, three vortex shedding patterns (anticlockwise/clockwise circular and flapping) are identified not only at $\alpha =0^\circ$ but also $\alpha \neq 0^\circ$. The feature of wake position in the radial direction for each pattern is observed regardless of the angle of attack. The relationship between the recirculation bubble pumping and the wake position in the radial direction is apparent in the non-reattaching flow but is weaker with $\alpha$ in the reattaching flow.

Direct numerical simulation is performed for flow separation over a bump in a turbulent channel. Comparisons are made between a smooth bump and one where the lee side is covered with replicas of shark denticles – dermal scales that consist of a slender base (the neck) and a wide top (the crown). As flow over the bump is under an adverse pressure gradient (APG), a reverse pore flow is formed in the porous cavity region underneath the crowns of the denticle array. Remarkable thrust is generated by the reverse pore flow as denticle necks accelerate the fluid passing between them in the upstream direction. Several geometrical features of shark denticles, including some that had not previously been considered hydrodynamically functional, are identified to form the two-layer denticle structure that enables and sustains the reverse pore flow and thrust generation. The reverse pore flow is activated by the APG before massive flow detachment. The results indicate a proactive, on-demand drag reduction mechanism that leverages and transforms the APG into a favourable outcome.

Low-frequency phenomena in an incompressible pressure-induced laminar separation bubble (LSB) on a flat plate is investigated using direct numerical simulation. The LSB configuration of Spalart and Strelets (J. Fluid Mech., vol. 403, 2000, pp. 329–349) is used. Wall pressure spectra indicate low-frequency-flapping $(St \sim 0.08)$ and high-frequency-shedding $(St \sim 1.52)$ regimes. Conditional velocity averages based on the fraction of reversed flow reveal the low frequency as an expansion/contraction of the LSB. While the high frequency only exhibits exponential growth within the LSB up to breakdown of the spanwise rollers, the low frequency and velocity fluctuations exhibit exponential growth upstream of separation. Instantaneous flow fields reveal large streamwise streaky structures forming within the LSB and extending past reattachment, much like high and low speed streaks in turbulent boundary layers. A predominance of sweep-like events ($Q4$) is observed during contraction and of ejection-like events ($Q2$) during expansion. These motions appear as dominant low-frequency modes in three-dimensional proper orthogonal and dynamic mode decompositions, exhibiting spatial amplification from separation to reattachment. The advection of a group of spanwise alternating streaky structures past the LSB results in an overall contraction after which the bubble expands to its ‘unforced’ state in the absence of the streaks. The low frequency then corresponds to the time it takes for streaks to form, amplify and advect past the LSB from separation to reattachment. This behaviour is linked to the mean flow deformation reported by Marxen and Rist (J. Fluid Mech., vol. 660, 2010, pp. 37–54), where the presence of streaks results in reduced mean bubble size. The formation of these streaky structures, in the absence of free stream turbulence, may be attributed to an absolute instability of the LSB due to the development of a secondary bubble within the primary.

The spatio-temporal scales, as well as a comprehensive self-sustained mechanism of the reattachment unsteadiness in shock wave/boundary layer interaction, are investigated in this study. Direct numerical simulations reveal that the reattachment unsteadiness of a Mach 7.7 laminar inflow causes over $26\,\%$ variation in wall friction and up to $20\,\%$ fluctuation in heat flux at the reattachment of the separation bubble. A statistical approach, based on the local reattachment upstream movement, is proposed to identify the spanwise and temporal scales of reattachment unsteadiness. It is found that two different types, i.e. self-induced and random processes, dominate different regions of reattachment. A self-sustained mechanism is proposed to comprehend the reattachment unsteadiness in the self-induced region. The intrinsic instability of the separation bubble transports vorticity downstream, resulting in an inhomogeneous reattachment line, which gives rise to baroclinic production of quasi-streamwise vortices. The pairing of these vortices initiates high-speed streaks and shifts the reattachment line upstream. Ultimately, viscosity dissipates the vortices, triggering instability and a new cycle of reattachment unsteadiness. The temporal scale and maximum vorticity are estimated with the self-sustained mechanism via order-of-magnitude analysis of the enstrophy. The advection speed of friction, derived from the assumption of coherent structures advecting with a Blasius-type boundary layer, aligns with the numerical findings.

We investigate effect of porous insert located upstream of the separation edge of a backward-facing step (BFS) in early transitional regime as a function of Reynolds number. This is an example of hydrodynamic system that is a combination of separated shear flow with large amplification potential and porous materials known for efficient flow destabilisation. Spectral analysis reveals that dynamics of BFS is dominated by spectral modes that remain globally coherent along the streamwise direction. We detect two branches of characteristic frequencies in the flow and with Hilbert transform we characterise their spatial support. For low Reynolds numbers, the dynamics of the flow is dominated by lower frequency, whereas for sufficiently large Reynolds numbers cross-over to higher frequencies is observed. Increasing permeability of the porous insert results in decrease in Reynolds number value, at which frequency cross-over occurs. By comparing normalised frequencies on each branch with local stability analysis, we attribute Kelvin–Helmholtz and Tollmien–Schlichting instabilities to upper and lower frequency branches, respectively. Finally, our results show that porous inserts enhance Kelvin–Helmholtz instability and promote transition to oscillator-type dynamics. Specifically, the amplitude of vortical (BFS) structures associated with higher-frequency branch follows Landau model prediction for all investigated porous inserts.

The turbulent wake behind a flat-back Ahmed body is investigated using stacked stereoscopic particle image velocimetry. The wake is disturbed by a steady jet from the centre of the base and the effects are quantified for key blowing rates. The unactuated wake exhibits bistable dynamics in the horizontal plane that are completely subdued for the optimal blowing case, yielding a base drag reduction of 9 %. The three-dimensional mean wake is reconstructed and used to evaluate the wake mass fluxes whose equilibrium determines the recirculation length. The results for the unactuated wake show that up to 80 % of replenishment fluid flux entering the recirculation bubble from the free-stream flow is provided through the low-pressure side of the base, where the symmetry-breaking shear layer roll-up occurs near the base. For the optimal blowing configuration, where the wake becomes symmetric, the flux of wake replenishment is severely reduced. This flow configuration results in elongated shear layers on all sides, which terminate the bubble with a roll-up of reduced intensity at a further downstream location. The dominant cause of bubble growth and the accompanying drag reduction is attributed to the momentum of the base blowing, and the new regime is referred to as the ‘favourable momentum regime’. Similar trends are observed when the model is at $5^{\circ }$ yaw where a reduction of drag and yaw-induced asymmetry are obtained. Proper orthogonal decomposition of the wake reveals the coherent structures related to the bistable flow and the symmetric wake under optimal blowing coefficient.

High angle of attack flows over swept three-dimensional wings based on the NACA 0015 profile are studied numerically at low Reynolds numbers. Linear stability analysis is used to compute instability and receptivity of the flow via the respective three-dimensional (triglobal) direct and adjoint eigenmodes. The magnitude of the adjoint eigenvectors is used to identify regions of maximum flow receptivity to momentum forcing. It is found that such regions are located above the primary three-dimensional separation line, their spanwise position varying with wing sweep. The wavemaker region corresponding to the leading global eigenmode is computed and found to lie inside the laminar separation bubble (LSB) at the spanwise location of peak recirculation. Increasing the Reynolds number leads to the wavemaker becoming more compact in the spanwise direction, and concentrated in the top and bottom shear layers of the LSB. As sweep is introduced, the wavemaker moves towards the wing tip, following the spanwise displacement of maximum recirculation. Flow modifications resulting from application of different types of forcing are studied by direct numerical simulation initialised with insights gained from stability analysis. Periodic forcing at the regions of maximum receptivity to momentum forcing results in greater departure from the baseline case compared to same (low, linear) amplitude forcing applied elsewhere, underlining the potential of linear stability analysis to identify optimal regions for actuator positioning.