Multifidelity models (MFMs) can be used to construct predictive models for flow quantities of interest (QoIs) over the space of uncertain/design parameters, with the purpose of uncertainty quantification, data fusion and optimization. For numerical simulation of turbulence, there is a hierarchy of methodologies ranked by accuracy and cost, where each methodology may have several numerical/modelling parameters that control the predictive accuracy and robustness of its resulting outputs. Compatible with these specifications, the present hierarchical MFM strategy allows for simultaneous calibration of the fidelity-specific parameters in a Bayesian framework as developed by Goh et al. (Technometrics, vol. 55, no. 4, 2013, pp. 501–512). The purpose of the MFM is to provide an improved prediction, mainly interpolation over the range covered by training data, by combining lower- and higher-fidelity data in an optimal way for any number of fidelity levels; even providing confidence intervals for the resulting QoI. The capabilities of the MFM are first demonstrated on an illustrative toy problem, and it is then applied to three realistic cases relevant to engineering turbulent flows. The latter include the prediction of friction at different Reynolds numbers in turbulent channel flow, the prediction of aerodynamic coefficients for a range of angles of attack of a standard airfoil and the uncertainty propagation and sensitivity analysis of the separation bubble in the turbulent flow over periodic hills subject to geometrical uncertainties. In all cases, based on only a few high-fidelity data samples, the MFM leads to accurate predictions of the QoIs. The result of the uncertainty quantification and sensitivity analyses are also found to be accurate compared with the ground truth in each case.

The flow physics of turbulent boundary layers is investigated using spectral analysis based on the spanwise scale decomposition of the Reynolds stress transport equation, with data obtained from a direct numerical simulation of the turbulent boundary layer at $Re_\tau \simeq 2020$. Here, we extend the framework of Kawata & Alfredsson (Phys. Rev. Lett., vol. 120, 2018, p. 244501) for plane Couette flows to zero-pressure-gradient boundary layers. The equation contains three fundamental fluxes, which govern the Reynolds stress transport: (i) a scale flux of the interaction between small-scale and large-scale structures, and two spatial fluxes dominated by (ii) pressure and (iii) turbulent transport along the wall-normal direction. The scale flux reveals evidence of the inverse turbulent kinetic energy transfer, from small to large scales, occurring at the near-wall region, whereby for the scale flux of the Reynolds shear stress transport, the inverse transfer extends across the entire boundary layer. The wall-normal fluxes reveal the interactions occurring between scales at the buffer and logarithmic regions. In addition, there is interaction between the large-scale structures and the free stream flow occurring at the edge of the boundary layer, which was not observed in the Couette flow. Flow structures associated with inverse interscale transport of Reynolds shear stress are identified by applying conditional analysis to the spectrally decomposed velocity fields. While the inverse transport is interpreted as the net energy transfer from small-scale ejections ($Q2$) and sweeps ($Q4$) to the large-scale counterparts, conditional time estimates of the direct and inverse interscale transport reveal that both processes play a substantial role across a broad range of scales.

The characteristics of turbulent boundary layers (TBLs) subjected to adverse pressure gradients are analysed through well-resolved large-eddy simulations. The geometries under study are the NACA0012 and NACA4412 wing sections, at $0^{\circ }$ and $5^{\circ }$ angle of attack, respectively, both of them at a Reynolds number based on inflow velocity and chord length of $Re_{c}=400\,000$. The turbulence statistics show that adverse pressure gradients (APGs) have a significant effect on the mean velocity, velocity fluctuations and turbulent kinetic energy budget, and this effect is more prominent on the outer region of the boundary layer. Furthermore, the effect of flow history is assessed by means of an integrated Clauser pressure-gradient parameter $\overline{\unicode[STIX]{x1D6FD}}$ (Vinuesa et al., Flow Turbul. Combust., vol. 99, 2017, pp. 565–587), through the study of cases with matching local values of $\unicode[STIX]{x1D6FD}$ and the friction Reynolds number $Re_{\unicode[STIX]{x1D70F}}$ to isolate this effect. Our results show a noticeable effect of the flow history on the outer region, however the differences in the near-wall peak of the tangential velocity fluctuations appear to be mostly produced by the local APG magnitude. The one-dimensional power-spectral density shows energetic small scales in the outer region of APG TBLs, whereas these energetic scales do not appear in zero-pressure-gradient (ZPG) TBLs, suggesting that small scales near the wall are advected towards the outer layer by the APG. Moreover, the linear coherence spectra show that the spectral outer peak of high-Reynolds-number ZPG TBLs is highly correlated with the near-wall region (Baars et al., J. Fluid Mech., vol. 823, 2017, R2), unlike APG TBLs which do not show such a correlation. This result, together with the different two-dimensional spectra of APG and high-Reynolds-number ZPG TBLs, suggests different energisation mechanisms due to APG and increase in Reynolds number. To the authors’ knowledge, this is the first in-depth analysis of the TBL characteristics over wings, including detailed single-point statistics, spectra and coherence.

Direct numerical simulations (DNSs) are performed to analyse the secondary flow of Prandtl’s second kind in fully developed spanwise-periodic channels with in-plane sinusoidal walls. The secondary flow is characterized for different combinations of wave parameters defining the wall geometry at $Re_{h}=2500$ and 5000, where $h$ is the half-height of the channel. The total cross-flow rate in the channel $Q_{yz}$ is defined along with a theoretical model to predict its behaviour. Interaction between the secondary flows from opposite walls is observed if $\unicode[STIX]{x1D706}\simeq h\simeq A$, where $A$ and $\unicode[STIX]{x1D706}$ are the amplitude and wavelength of the sinusoidal function defining the wall geometry. As the outer-scaled wavelength ($\unicode[STIX]{x1D706}/h$) is reduced, the secondary vortices become smaller and faster, increasing the total cross-flow rate per wall. However, if the inner-scaled wavelength ($\unicode[STIX]{x1D706}^{+}$) is below 130 viscous units, the cross-flow decays for smaller wavelengths. By analysing cases in which the wavelength of the wall is much smaller than the half-height of the channel $\unicode[STIX]{x1D706}\ll h$, we show that the cross-flow distribution depends almost entirely on the separation between the scales of the instantaneous vortices, where the upper and lower bounds are determined by $\unicode[STIX]{x1D706}/h$ and $\unicode[STIX]{x1D706}^{+}$, respectively. Therefore, the distribution of the secondary flow relative to the size of the wave at a given $Re_{h}$ can be replicated at higher $Re_{h}$ by decreasing $\unicode[STIX]{x1D706}/h$ and keeping $\unicode[STIX]{x1D706}^{+}$ constant. The mechanisms that contribute to the mean cross-flow are analysed in terms of the Reynolds stresses and using quadrant analysis to evaluate the probability density function of the bursting events. These events are further classified with respect to the sign of their instantaneous spanwise velocities. Sweeping events and ejections are preferentially located in the valleys and peaks of the wall, respectively. The sweeps direct the instantaneous cross-flow from the core of the channel towards the wall, turning in the wall-tangent direction towards the peaks. The ejections drive the instantaneous cross-flow from the near-wall region towards the core. This preferential behaviour is identified as one of the main contributors to the secondary flow.

This paper proposes a resolution to the conundrum of the roles of convective and absolute instability in transition of the rotating-disk boundary layer. It also draws some comparison with swept-wing flows. Direct numerical simulations based on the incompressible Navier–Stokes equations of the flow over the surface of a rotating disk with modelled roughness elements are presented. The rotating-disk flow has been of particular interest for stability and transition research since the work by Lingwood (J. Fluid Mech., vol. 299, 1995, pp. 17–33) where an absolute instability was found. Here stationary disturbances develop from roughness elements on the disk and are followed from the linear stage, growing to saturation and finally transitioning to turbulence. Several simulations are presented with varying disturbance amplitudes. The lowest amplitude corresponds approximately to the experiment by Imayama et al. (J. Fluid Mech., vol. 745, 2014a, pp. 132–163). For all cases, the primary instability was found to be convectively unstable, and secondary modes were found to be triggered spontaneously while the flow was developing. The secondary modes further stayed within the domain, and an explanation for this is a proposed globally unstable secondary instability. For the low-amplitude roughness cases, the disturbances propagate beyond the threshold for secondary global instability before becoming turbulent, and for the high-amplitude roughness cases the transition scenario gives a turbulent flow directly at the critical Reynolds number for the secondary global instability. These results correspond to the theory of Pier (J. Engng Maths, vol. 57, 2007, pp. 237–251) predicting a secondary absolute instability. In our simulations, high temporal frequencies were found to grow with a large amplification rate where the secondary global instability occurred. For smaller radial positions, low-frequency secondary instabilities were observed, tripped by the global instability.

For problems governed by a non-normal operator, the leading eigenvalue of the operator is of limited interest and a more relevant measure of the stability is obtained by considering the harmonic forcing causing the largest system response. Various methods for determining this so-called optimal forcing exist, but they all suffer from great computational expense and are hence not practical for large-scale problems. In the present paper a new method is presented, which is applicable to problems of arbitrary size. The method does not rely on timestepping, but on the solution of linear systems, in which the inverse Laplacian acts as a preconditioner. By formulating the search for the optimal forcing as an eigenvalue problem based on the resolvent operator, repeated system solves amount to power iterations, in which the dominant eigenvalue is seen to correspond to the energy amplification in a system for a given frequency, and the eigenfunction to the corresponding forcing function. Implementation of the method requires only minor modifications of an existing timestepping code, and is applicable to any partial differential equation containing the Laplacian, such as the Navier-Stokes equations. We discuss the method, first, in the context of the linear Ginzburg-Landau equation and then, the two-dimensional lid-driven cavity flow governed by the Navier-Stokes equations. Most importantly, we demonstrate that for the lid-driven cavity, the optimal forcing can be computed using a factor of up to 500 times fewer operator evaluations than the standard method based on exponential timestepping.

This paper introduces a new method based on the diagnostic plot (Alfredsson et al., Phys. Fluids, vol. 23, 2011, 041702) to assess the convergence towards a well-behaved zero-pressure-gradient (ZPG) turbulent boundary layer (TBL). The most popular and well-understood methods to assess the convergence towards a well-behaved state rely on empirical skin-friction curves (requiring accurate skin-friction measurements), shape-factor curves (requiring full velocity profile measurements with an accurate wall position determination) or wake-parameter curves (requiring both of the previous quantities). On the other hand, the proposed diagnostic-plot method only needs measurements of mean and fluctuating velocities in the outer region of the boundary layer at arbitrary wall-normal positions. To test the method, six tripping configurations, including optimal set-ups as well as both under- and overtripped cases, are used to quantify the convergence of ZPG TBLs towards well-behaved conditions in the Reynolds-number range covered by recent high-fidelity direct numerical simulation data up to a Reynolds number based on the momentum thickness and free-stream velocity $Re_{\unicode[STIX]{x1D703}}$ of approximately 4000 (corresponding to 2.5 m from the leading edge) in a wind-tunnel experiment. Additionally, recent high-Reynolds-number data sets have been employed to validate the method. The results show that weak tripping configurations lead to deviations in the mean flow and the velocity fluctuations within the logarithmic region with respect to optimally tripped boundary layers. On the other hand, a strong trip leads to a more energized outer region, manifested in the emergence of an outer peak in the velocity-fluctuation profile and in a more prominent wake region. While established criteria based on skin-friction and shape-factor correlations yield generally equivalent results with the diagnostic-plot method in terms of convergence towards a well-behaved state, the proposed method has the advantage of being a practical surrogate that is a more efficient tool when designing the set-up for TBL experiments, since it diagnoses the state of the boundary layer without the need to perform extensive velocity profile measurements.

Turbulent boundary layers under adverse pressure gradients are studied using well-resolved large-eddy simulations (LES) with the goal of assessing the influence of the streamwise pressure-gradient development. Near-equilibrium boundary layers were characterized through the Clauser pressure-gradient parameter $\unicode[STIX]{x1D6FD}$. In order to fulfil the near-equilibrium conditions, the free stream velocity was prescribed such that it followed a power-law distribution. The turbulence statistics pertaining to cases with a constant value of $\unicode[STIX]{x1D6FD}$ (extending up to approximately 40 boundary-layer thicknesses) were compared with cases with non-constant $\unicode[STIX]{x1D6FD}$ distributions at matched values of $\unicode[STIX]{x1D6FD}$ and friction Reynolds number $Re_{\unicode[STIX]{x1D70F}}$. An additional case at matched Reynolds number based on displacement thickness $Re_{\unicode[STIX]{x1D6FF}^{\ast }}$ was also considered. It was noticed that non-constant $\unicode[STIX]{x1D6FD}$ cases appear to approach the conditions of equivalent constant $\unicode[STIX]{x1D6FD}$ cases after long streamwise distances (approximately 7 boundary-layer thicknesses). The relevance of the constant $\unicode[STIX]{x1D6FD}$ cases lies in the fact that they define a ‘canonical’ state of the boundary layer, uniquely characterized by $\unicode[STIX]{x1D6FD}$ and $Re$. The investigations on the flat plate were extended to the flow around a wing section overlapping in terms of $\unicode[STIX]{x1D6FD}$ and $Re$. Comparisons with the flat-plate cases at matched values of $\unicode[STIX]{x1D6FD}$ and $Re$ revealed that the different development history of the turbulent boundary layer on the wing section leads to a less pronounced wake in the mean velocity as well as a weaker second peak in the Reynolds stresses. This is due to the weaker accumulated effect of the $\unicode[STIX]{x1D6FD}$ history. Furthermore, a scaling law suggested by Kitsios et al. (Intl J. Heat Fluid Flow, vol. 61, 2016, pp. 129–136), proposing the edge velocity and the displacement thickness as scaling parameters, was tested on two constant-pressure-gradient parameter cases. The mean velocity and Reynolds-stress profiles were found to be dependent on the downstream development. The present work is the first step towards assessing history effects in adverse-pressure-gradient turbulent boundary layers and highlights the fact that the values of the Clauser pressure-gradient parameter and the Reynolds number are not sufficient to characterize the state of the boundary layer.

A numerical investigation of two locally applied drag-reducing control schemes is carried out in the configuration of a spatially developing turbulent boundary layer (TBL). One control is designed to damp near-wall turbulence and the other induces constant mass flux in the wall-normal direction. Both control schemes yield similar local drag reduction rates within the control region. However, the flow development downstream of the control significantly differs: persistent drag reduction is found for the uniform blowing case, whereas drag increase is found for the turbulence damping case. In order to account for this difference, the formulation of a global drag reduction rate is suggested. It represents the reduction of the streamwise force exerted by the fluid on a plate of finite length. Furthermore, it is shown that the far-downstream development of the TBL after the control region can be described by a single quantity, namely a streamwise shift of the uncontrolled boundary layer, i.e. a changed virtual origin. Based on this result, a simple model is developed that allows the local drag reduction rate to be related to the global one without the need to conduct expensive simulations or measurements far downstream of the control region.

Direct numerical simulations based on the incompressible nonlinear Navier–Stokes equations of the flow over the surface of a rotating disk have been conducted. An impulsive disturbance was introduced and its development as it travelled radially outwards and ultimately transitioned to turbulence has been analysed. Of particular interest was whether the nonlinear stability is related to the linear stability properties. Specifically three disk-edge conditions were considered; (i) a sponge region forcing the flow back to laminar flow, (ii) a disk edge, where the disk was assumed to be infinitely thin and (iii) a physically realistic disk edge of finite thickness. This work expands on the linear simulations presented by Appelquist et al. (J. Fluid. Mech., vol. 765, 2015, pp. 612–631), where, for case (i), this configuration was shown to be globally linearly unstable when the sponge region effectively models the influence of the turbulence on the flow field. In contrast, case (ii) was mentioned there to be linearly globally stable, and here, where nonlinearity is included, it is shown that both cases (ii) and (iii) are nonlinearly globally unstable. The simulations show that the flow can be globally linearly stable if the linear wavepacket has a positive front velocity. However, in the same flow field, a nonlinear global instability can emerge, which is shown to depend on the outer turbulent region generating a linear inward-travelling mode that sustains a transition front within the domain. The results show that the front position does not approach the critical Reynolds number for the local absolute instability, $R=507$. Instead, the front approaches $R=583$ and both the temporal frequency and spatial growth rate correspond to a global mode originating at this position.

The concept of edge states is investigated in the asymptotic suction boundary layer in relation to the receptivity process to noisy perturbations and the nucleation of turbulent spots. Edge tracking is first performed numerically, without imposing any discrete symmetry, in a large computational domain allowing for full spatial localisation of the perturbation velocity. The edge state is a three-dimensional localised structure recurrently characterised by a single low-speed streak that experiences erratic bursts and planar shifts. This recurrent streaky structure is then compared with predecessors of individual spot nucleation events, triggered by non-localised initial noise. The present results suggest a nonlinear picture, rooted in dynamical systems theory, of the nucleation process of turbulent spots in boundary-layer flows, in which the localised edge state plays the role of state-space mediator.

Turbulence in the asymptotic suction boundary layer is investigated numerically at the verge of laminarisation using direct numerical simulation. Following an adiabatic protocol, the Reynolds number $Re$ is decreased in small steps starting from a fully turbulent state until laminarisation is observed. Computations in a large numerical domain allow in principle for the possible coexistence of laminar and turbulent regions. However, contrary to other subcritical shear flows, no laminar–turbulent coexistence is observed, even near the onset of sustained turbulence. High-resolution computations suggest a critical Reynolds number $Re_{g}\approx 270$, below which turbulence collapses, based on observation times of $O(10^{5})$ inertial time units. During the laminarisation process, the turbulent flow fragments into a series of transient streamwise-elongated structures, whose interfaces do not display the characteristic obliqueness of classical laminar–turbulent patterns. The law of the wall, i.e. logarithmic scaling of the velocity profile, is retained down to $Re_{g}$, suggesting a large-scale wall-normal transport absent in internal shear flows close to the onset. In order to test the effect of these large-scale structures on the near-wall region, an artificial volume force is added to damp spanwise and wall-normal fluctuations above $y^{+}=100$, in viscous units. Once the largest eddies have been suppressed by the forcing, and thus turbulence is confined to the near-wall region, oblique laminar–turbulent interfaces do emerge as in other wall-bounded flows, however only transiently. These results suggest that oblique stripes at the onset are a prevalent feature of internal shear flows, but will not occur in canonical boundary layers, including the spatially growing ones.

Direct numerical simulations of the fully developed turbulent flow through a porous square duct are performed to study the effect of the permeable wall on the secondary cross-stream flow. The volume-averaged Navier–Stokes equations are used to describe the flow in the porous phase, a packed bed with porosity ${\it\varepsilon}_{c}=0.95$. The porous square duct is computed at $\mathit{Re}_{b}\simeq 5000$ and compared with the numerical simulations of a turbulent duct with four solid walls. The two boundary layers on the top wall and porous interface merge close to the centre of the duct, as opposed to the channel, because the sidewall boundary layers inhibit the growth of the shear layer over the porous interface. The most relevant feature in the porous duct is the enhanced magnitude of the secondary flow, which exceeds that of a regular duct by a factor of four. This is related to the increased vertical velocity, and the different interaction between the ejections from the sidewalls and the porous medium. We also report a significant decrease in the streamwise turbulence intensity over the porous wall of the duct (which is also observed in a porous channel), and the appearance of short spanwise rollers in the buffer layer, replacing the streaky structures of wall-bounded turbulence. These spanwise rollers most probably result from a Kelvin–Helmholtz type of instability, and their width is limited by the presence of the sidewalls.

Direct numerical simulations of turbulent channel flows at friction Reynolds numbers (Re) of 550, 1000 and 1500 are used to analyse the turbulent production, transfer and dissipation mechanisms in the compound space of scales and wall distances by means of the Kolmogorov equation generalized to inhomogeneous anisotropic flows. Two distinct peaks of scale-energy source are identified. The first, stronger one, belongs to the near-wall cycle. Its location in the space of scales and physical space is found to scale in viscous units, while its intensity grows slowly with $\mathit{Re}$, indicating a near-wall modulation. The second source peak is found further away from the wall in the putative overlap layer, and it is separated from the near-wall source by a layer of significant scale-energy sink. The dynamics of the second outer source appears to be strongly dependent on the Reynolds number. The detailed scale-by-scale analysis of this source highlights well-defined features that are used to make the properties of the outer turbulent source independent of Reynolds number and wall distance by rescaling the problem. Overall, the present results suggest a strong connection of the observed outer scale-energy source with the presence of an outer region of turbulence production whose mechanisms are well separated from the near-wall region and whose statistical features agree with the hypothesis of an overlap layer dominated by attached eddies. Inner–outer interactions between the near-wall and outer source region in terms of scale-energy fluxes are also analysed. It is conjectured that the near-wall modulation of the statistics at increasing Reynolds number can be related to a confinement of the near-wall turbulence production due to the presence of increasingly large production scales in the outer scale-energy source region.

Numerical simulations of the flow developing on the surface of a rotating disk are presented based on the linearized incompressible Navier–Stokes equations. The boundary-layer flow is perturbed by an impulsive disturbance within a linear global framework, and the effect of downstream turbulence is modelled by a damping region further downstream. In addition to the outward-travelling modes, inward-travelling disturbances excited at the radial end of the simulated linear region, $r_{end}$, by the modelled turbulence are included within the simulations, potentially allowing absolute instability to develop. During early times the flow shows traditional convective behaviour, with the total energy slowly decaying in time. However, after the disturbances have reached $r_{end}$, the energy evolution reaches a turning point and, if the location of $r_{end}$ is at a Reynolds number larger than approximately $R=594$ (radius non-dimensionalized by $\sqrt{{\it\nu}/{\rm\Omega}^{\ast }}$, where ${\it\nu}$ is the kinematic viscosity and ${\rm\Omega}^{\ast }$ is the rotation rate of the disk), there will be global temporal growth. The global frequency and mode shape are clearly imposed by the conditions at $r_{end}$. Our results suggest that the linearized Ginzburg–Landau model by Healey (J. Fluid Mech., vol. 663, 2010, pp. 148–159) captures the (linear) physics of the developing rotating-disk flow, showing that there is linear global instability provided the Reynolds number of $r_{end}$ is sufficiently larger than the critical Reynolds number for the onset of absolute instability.

We study the turbulent Ekman layer at moderately high Reynolds number, $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}1600 < \mathit{Re} = \delta _{E}G/\nu < 3000$, using direct numerical simulations (DNS). Here, $\delta _{E} = \sqrt{2\nu /f}$ is the laminar Ekman layer thickness, $G$ the geostrophic wind, $\nu $ the kinematic viscosity and $f$ is the Coriolis parameter. We present results for both neutrally, moderately and strongly stably stratified conditions. For unstratified cases, large-scale roll-like structures extending from the outer region down to the wall are observed. These structures have a clear dominant frequency and could be related to periodic oscillations or instabilities developing near the low-level jet. We discuss the effect of stratification and $\mathit{Re}$ on one-point and two-point statistics. In the strongly stratified Ekman layer we observe stable co-existing large-scale laminar and turbulent patches appearing in the form of inclined bands, similar to other wall-bounded flows. For weaker stratification, continuously sustained turbulence strongly affected by buoyancy is produced. We discuss the scaling of turbulent length scales, height of the Ekman layer, friction velocity, veering angle at the wall and heat flux. The boundary-layer thickness, the friction velocity and the veering angle depend on $Lf/u_\tau $, where $u_\tau $ is the friction velocity and $L$ the Obukhov length scale, whereas the heat fluxes appear to scale with $L^+=L u_\tau /\nu $.

Interactions between seabirds and commercial fishing activities have been well documented but little information is available regarding the impacts of more traditional fishing practices on seabird populations. We interviewed fishermen, administered questionnaires, and made field-based observations to determine the extent to which artisanal fisheries interact with and affect seabirds in the fjords and channels of the Chonos archipelago in southern Chile. Our surveys indicated a positive perception of seabirds as useful indicators of marine productivity and in their role scavenging fish waste and discards associated with fishing operations. However, the surveys also revealed that fishermen routinely establish seasonal camps for collecting seabird eggs and adults for food or bait and introduce feral predators to seabird breeding colonies on islands. Understanding the traditional practices of fishermen is critical for the future of community-based conservation of the region's marine resources and biodiversity.

The dynamics on the laminar–turbulent separatrix is investigated numerically for boundary-layer flows in the subcritical regime. Constant homogeneous suction is applied at the wall, resulting in a parallel asymptotic suction boundary layer (ASBL). When the numerical domain is sufficiently extended in the spanwise direction, the coherent structures found by edge tracking are invariably localized and their dynamics shows bursts that drive a remarkable regular or irregular spanwise dynamics. Depending on the parameters, the asymptotic dynamics on the edge can be either periodic in time or chaotic. A clear mechanism for the regeneration of streaks and streamwise vortices emerges in all cases and is investigated in detail.

Results are presented from a direct numerical simulation of a particle-laden spatially developing turbulent boundary layer up to . The peculiar feature of a boundary-layer flow seeded with heavy particles is the variation of the local dimensionless parameters defining the fluid–particle interactions along the streamwise direction. Two different Stokes numbers can be defined, one using inner flow units and the other with outer units. Since these two Stokes numbers exhibit different decay rates in the streamwise direction, we find a decoupled particle dynamics between the inner and the outer region of the boundary layer. Preferential near-wall particle accumulation is similar to that observed in turbulent channel flow, while different behaviour characterizes the outer region. Here the concentration and the streamwise velocity profiles are found to be self-similar and to depend only on the local value of the outer Stokes number and the rescaled wall-normal distance. These new results are powerful in view of engineering and environmental applications and corresponding flow modelling.

Direct numerical simulations of subcritical rotating, stratified and magneto-hydrodynamic wall-bounded flows are performed in large computational domains, focusing on parameters where laminar and turbulent flow can stably coexist. In most cases, a regime of large-scale oblique laminar-turbulent patterns is identified at the onset of transition, as in the case of pure shear flows. The current study indicates that this oblique regime can be shifted up to large values of the Reynolds number by increasing the damping by the Coriolis, buoyancy or Lorentz force. We show evidence for this phenomenon in three distinct flow cases: plane Couette flow with spanwise cyclonic rotation, plane magnetohydrodynamic channel flow with a spanwise or wall-normal magnetic field, and open channel flow under stable stratification. Near-wall turbulence structures inside the turbulent patterns are invariably found to scale in terms of viscous wall units as in the fully turbulent case, while the patterns themselves remain large-scale with a trend towards shorter wavelength for increasing . Two distinct regimes are identified: at low Reynolds numbers the patterns extend from one wall to the other, while at large Reynolds number they are confined to the near-wall regions and the patterns on both channel sides are uncorrelated, the core of the flow being highly turbulent without any dominant large-scale structure.