We investigate fully developed turbulent flow in curved channels to explore the interaction between turbulence and curvature-driven coherent structures. By focusing on two cases of mild and strong curvature, we examine systematically the effects of the Reynolds number through a campaign of direct numerical simulations, spanning flow regimes from laminar up to the moderately high Reynolds number – based on bulk velocity and channel height – of $87\,000$. Our analysis highlights the influence of curvature on the friction coefficient, showing that flow transition is anticipated by concave curvature and delayed by convex curvature. In the case of mild curvature, a frictional drag reduction compared with plane channel flow is found in the transitional regime. Spectral analysis reveals that the near-wall turbulence regeneration cycle is maintained in mildly curved channels, while it is absent or severely inhibited on the convex wall of strongly curved channels. Streamwise large-scale structures resembling Dean vortices are found to be weakly dependent on the Reynolds number and strongly affected by curvature: increasing curvature shifts these vortices towards the outer wall and reduces their size and coherence, limiting their contribution to streamwise velocity fluctuations and momentum transport. In the case of strong curvature, spanwise large-scale structures are also detected. These structures are associated with large pressure fluctuations and the suppression of turbulent stresses near the convex wall, where a region with negative turbulence production is observed and characterised via quadrant analysis.

We introduce a wall model for large-eddy simulation (WMLES) applicable to rough surfaces with Gaussian and non-Gaussian distributions for both the transitionally and fully rough regimes. The model is applicable to arbitrary complex geometries where roughness elements are assumed to be underresolved, i.e. subgrid-scale roughness. The wall model is implemented using a multi-hidden-layer feedforward neural network, with the mean geometric properties of the roughness topology and near-wall flow quantities serving as input. The optimal set of non-dimensional input features is identified using information theory, selecting variables that maximize information about the output while minimizing redundancy among inputs. The model also incorporates a confidence score based on Gaussian process modelling, enabling the detection of potentially low model performance for untrained rough surfaces. The model is trained using a direct numerical simulation (DNS) roughness database comprising approximately 200 cases. The roughness geometries for the database are selected from a large repository through active learning. This approach ensures that the rough surfaces incorporated into the database are the most informative, achieving higher model performance with fewer DNS cases compared with passive learning techniques. The performance of the model is evaluated both apriori and aposteriori in WMLES of turbulent channel flows with rough walls. Over 550 channel flow cases are considered, including untrained roughness geometries, roughness Reynolds numbers and grid resolutions for both transitionally and fully rough regimes. Our rough-wall model offers higher accuracy than existing models, generally predicting wall shear stress within an accuracy range of 1%–15 %. The performance of the model is also assessed on a high-pressure turbine blade with two different rough surfaces. We show that the new wall model predicts the skin friction and the mean velocity deficit induced by the rough surface on the blade within 1%–10 % accuracy except the region with transition or shock waves. This work extends the building-block flow wall model (BFWM) introduced by Lozano-Durán & Bae (2023. J. Fluid Mech. 963, A35) for smooth walls, expanding the BFWM framework to account for rough-wall scenarios.

Large-eddy simulations are analysed to determine the influence of suspended canopies, such as those formed in macroalgal farms, on ocean mixed layer (OML) deepening and internal wave generation. In the absence of a canopy, we show that Langmuir turbulence, when compared with wind-driven shear turbulence, results in a deeper OML and more pronounced internal waves beneath the OML. Subsequently, we examine simulations with suspended canopies of varying densities located in the OML, in the presence of a background geostrophic current. Intensified turbulence occurs in the shear layer at the canopy’s bottom edge, arising from the interaction between the geostrophic current and canopy drag. Structures resembling Kelvin–Helmholtz (KH) instability emerge as the canopy shear layer interacts with the underlying stratification, radiating internal waves beneath the OML. Both intensified turbulence and lower-frequency motions associated with KH-type structures are critical factors in enhancing mixing. Consequently, the OML depth increases by up to a factor of two compared with cases without a canopy. Denser canopies and stronger geostrophic currents lead to more pronounced KH-type structures and internal waves, stronger turbulence and greater OML deepening. Additionally, vertical nutrient transport is enhanced as the OML deepens due to the presence of the canopy. Considering that the canopy density investigated in this study closely represents offshore macroalgal farms, these findings suggest a mechanism for passive nutrient entrainment conducive to sustainable farming. Overall, this study reveals the intricate interactions between the suspended canopy, turbulent mixing and stratification, underscoring their importance in reshaping OML characteristics.

A collection of secondary instability calculations in streaky boundary layers is presented. The data are retrieved from well-resolved numerical simulations of boundary layers forced by free-stream turbulence (FST), considering different geometries and FST conditions. The stability calculations are performed before streak breakdown, taking place at various $Rey_x$ the Reynolds number based on the streamwise coordinate. Despite the rich streak population of various sizes, it is found that breaking streaks have similar aspect ratios, independently of the streamwise position where they appear. This suggests that wider streaks will break down further downstream than thinner ones, making the appearance of secondary instabilities somewhat independent of the streak’s wavelength. Moreover, the large difference in the integral length scale among the simulations suggests that this aspect ratio is also independent of the FST scales. An explanation for this behaviour is provided by showing that these breaking streaks are in the range of perturbations that can experience maximum transient growth according to optimal disturbance theory. This could explain why, at a given streamwise position, there is a narrow spanwise wavelength range where streak breakdown is more likely to occur.

Wave propagation in channels with area changes is a topic of significant practical interest that involves a rich set of coupled physics. While the acoustic wave problem has been studied extensively, the shock propagation problem has received less attention. In addition to its practical significance, this problem also introduces deep fundamental issues associated with how energy in propagating large-amplitude disturbances is redistributed upon interaction with inhomogeneities. This paper presents a study of shock scattering and entropy and vorticity coupling for shock wave propagation through discrete area changes. It compares results from computational fluid dynamics to those of one-dimensional quasi-steady calculations. The solution space is naturally divided into five ‘regimes’ based upon the incident shock strength and area ratio. This paper also presents perturbation methods to quantify the dimensionless scaling of physical effects associated with wave reflection/transmission and energy transfer to other disturbances. Finally, it presents an analysis of the ‘energetics’ of the interaction, quantifying how energy that initially resides in dilatational disturbances and propagates at the shock speed is redistributed into finite-amplitude reflected and transmitted waves as well as convecting vortical and entropy disturbances.

The formation process of a vortex pair generated by a two-dimensional starting jet has been investigated numerically over a range of Reynolds numbers from 500 to 2000. The effects of stroke ratio and nozzle configuration are examined. Only a single vortex pair can be observed in the vorticity field generated by small stroke ratios less than 10 while the leading vortex pair formed by larger stroke ratios eventually disconnects from the trailing jet. The formation numbers (13.6 and 9.3) for a straight nozzle and an orifice nozzle have been identified by the circulation criterion and they are further analysed by four other criteria. Using the contraction coefficient, formation numbers can be transformed into a universal value at about 16.5 for both nozzles. The effect of Reynolds number on the formation number is found to be within 12 % for parallel flow cases but it will increase up to 27 % for non-parallel flow cases due to shear-layer instability. A modified contraction-based slug model is proposed, and it can accurately predict the total invariants (e.g. circulation, hydrodynamic impulse and kinetic energy) shedding from the nozzle edge. Analytical estimation of the formation number is further conducted by matching the predicted total invariants to the Pierrehumbert model of steady vortex pairs. By assuming that pinch-off starts when the vortex pair achieves the steady state, two analytical models are proposed in terms of vortex impulse and translational velocity. The latter appears to be more appropriate to predict the formation number for two-dimensional flows.

Turbulent transonic buffet is an aerodynamic instability causing periodic (albeit, often irregular) oscillations of lift/drag in aerospace applications. Involving complex coupling between inviscid and viscous effects, buffet is characterised by shock wave oscillations and flow separation/reattachment. Previous studies have identified both two-dimensional (2-D) chordwise shock-oscillation and three-dimensional (3-D) buffet-/stall-cell modes. While the 2-D instability has been studied extensively, investigations of 3-D buffet have been limited to only low-fidelity simulations or experiments. Due to computational cost, almost all high-fidelity studies to date have been limited to narrow span-widths around 5 % of aerofoil chord length (aspect ratio, ), which is insufficiently wide to observe large-scale three-dimensionality. In this work, high-fidelity simulations are performed up to , on an infinite unswept NASA Common Research Model (CRM) wing profile at $Re=5\times 10^{5}$. At , intermittent 3-D separation bubbles are observed at buffet conditions. While previous Reynolds-averaged Navier–Stokes (RANS)/stability-based studies predict quasi-simultaneous onset of 2-D- and 3-D-buffet, a case that remains essentially 2-D is identified here. Strongest three-dimensionality was observed near low-lift phases of the buffet cycle at maximum flow separation, reverting to essentially 2-D behaviour during high-lift phases. Buffet was found to become 3-D when extensive mean flow separation was present. At , multiple 3-D separation bubbles form in a spanwise wavelength range $\lambda =1c$ to $1.5c$. Spectral proper orthogonal decomposition (SPOD) was applied to analyse the spatio/temporal structure of 3-D buffet-cells. In addition to the 2-D chordwise shock-oscillation mode (Strouhal number $St \approx 0.07-0.1$), 3-D modal structures were observed at the shock wave/boundary layer interaction at $St \approx 0.002-0.004$.