Motivated by the need for compact descriptions of the evolution of non-classical wakes behind yawed wind turbines, we develop an analytical model to predict the shape of curled wakes. Interest in such modelling arises due to the potential of wake steering as a strategy for mitigating power reduction and unsteady loading of downstream turbines in wind farms. We first estimate the distribution of the shed vorticity at the wake edge due to both yaw offset and rotating blades. By considering the wake edge as an ideally thin vortex sheet, we describe its evolution in time moving with the flow. Vortex sheet equations are solved using a power series expansion method, and an approximate solution for the wake shape is obtained. The vortex sheet time evolution is then mapped into a spatial evolution by using a convection velocity. Apart from the wake shape, the lateral deflection of the wake including ground effects is modelled. Our results show that there exists a universal solution for the shape of curled wakes if suitable dimensionless variables are employed. For the case of turbulent boundary layer inflow, the decay of vortex sheet circulation due to turbulent diffusion is included. Finally, we modify the Gaussian wake model by incorporating the predicted shape and deflection of the curled wake, so that we can calculate the wake profiles behind yawed turbines. Model predictions are validated against large-eddy simulations and laboratory experiments for turbines with various operating conditions.

A quantitative understanding of the dominant mechanisms that govern the generation and decay of the counter-rotating vortex pair (CVP) produced by yawed wind turbines is needed to fully realize the potential of yawing for wind farm power maximization and regulation. Observations from large eddy simulations (LES) of yawed wind turbines in the turbulent atmospheric boundary layer and concepts from the aircraft trailing vortex literature inform a model for the shed vorticity and circulation. The model is formed through analytical integration of simplified forms of the vorticity transport equation. Based on an eddy viscosity approach, it uses the boundary-layer friction velocity as the velocity scale and the width of the vorticity distribution itself as the length scale. As with the widely used Jensen model for wake deficit evolution in wind farms, our analytical expressions do not require costly numerical integration of differential equations. The predicted downstream decay of maximum vorticity and total circulation agree well with LES results. We also show that the vorticity length scale grows linearly with downstream distance and find several power laws for the decay of maximum vorticity. These results support the notion that the decay of the CVP is dominated by gradual cancellation of the vorticity at the line of symmetry of the wake through cross-diffusion.

By analysing the Karman–Howarth equation for filtered-velocity fields in turbulent flows, we show that the two-point correlation between the filtered strain-rate and subfilter stress tensors plays a central role in the evolution of filtered-velocity correlation functions. Two-point correlation-based statistical a priori tests thus enable rigorous and physically meaningful studies of turbulence models. Using data from direct numerical simulations of isotropic and channel flow turbulence, we show that local eddy-viscosity models fail to exhibit the long tails observed in the real subfilter stress–strain-rate correlation functions. Stronger non-local correlations may be achieved by defining the eddy-viscosity model based on fractional gradients of order $0<\alpha <1$ (where $\alpha$ is the fractional gradient order) rather than the classical gradient corresponding to $\alpha =1$. Analyses of such correlation functions are presented for various orders of the fractional-gradient operators. It is found that in isotropic turbulence fractional derivative order $\alpha \sim 0.5$ yields best results, while for channel flow $\alpha \sim 0.2$ yields better results for the correlations in the streamwise direction, even well into the core channel region. In the spanwise direction, channel flow results show significantly more local interactions. The overall results confirm strong non-locality in the interactions between subfilter stresses and resolved-scale fluid deformation rates, but with non-trivial directional dependencies in non-isotropic flows. Hence, non-local operators thus exhibit interesting modelling capabilities and potential for large-eddy simulations although more developments are required, both on the theoretical and computational implementation fronts.

Based on a generalized local Kolmogorov–Hill equation expressing the evolution of kinetic energy integrated over spheres of size $\ell$ in the inertial range of fluid turbulence, we examine a possible definition of entropy and entropy generation for turbulence. Its measurement from direct numerical simulations in isotropic turbulence leads to confirmation of the validity of the fluctuation relation (FR) from non-equilibrium thermodynamics in the inertial range of turbulent flows. Specifically, the ratio of probability densities of forward and inverse cascade at scale $\ell$ is shown to follow exponential behaviour with the entropy generation rate if the latter is defined by including an appropriately defined notion of ‘temperature of turbulence’ proportional to the kinetic energy at scale $\ell$.

We study the dispersion characteristics of slightly buoyant droplets in a turbulent jet using large eddy simulations (LES). The droplet number density fields are represented using an Eulerian approach, with the dispersed phase modelled using the Fast-Eulerian method (Ferry & Balachandar, Intl J. Multiphase Flow, vol. 27, issue 7, pp. 1199–1226, 2001) that includes the droplet rise velocity. Radial concentration profiles and turbulent concentration fluxes for droplets of different sizes are analysed to quantify the ‘trajectory crossing effect’, when relative motions between particles and turbulent eddies tend to reduce turbulent diffusion. For finer LES grid resolutions, the model captures the differential, size-based dispersion characteristics of the droplets with the transverse dispersion of the larger droplet sizes suppressed, since trajectory crossing effects are explicitly resolved in LES. We examine a similarity solution model for the size-dependent radial concentration profiles based on a modified Schmidt number derived from the theory of turbulent diffusion of particles in the atmosphere proposed by Csanady (J. Atmos. Sci., vol. 20, pp. 201–208, 1963). The results are validated with the high resolution LES data and show good agreement. Then the size-dependent Schmidt number model is reformulated as a model for unresolved subgrid-scale trajectory crossing effects and used to calculate the subgrid concentration flux in a coarse LES of a turbulent jet, with slightly buoyant droplets injected at the centreline in the self-similar region of the jet. The results are compared to a simulation with higher grid resolution and a coarse simulation with a constant Schmidt number subgrid-scale model. We find that the subgrid model enhances the prediction accuracy of the concentration profiles and turbulent concentration flux for the coarse LES.

Drag for wall-bounded flows is directly related to the spatial flux of spanwise vorticity outward from the wall. In turbulent flows a key contribution to this wall-normal flux arises from nonlinear advection and stretching of vorticity, interpretable as a cascade. We study this process using numerical simulation data of turbulent channel flow at friction Reynolds number $Re_\tau =1000$. The net transfer from the wall of spanwise vorticity created by downstream pressure drop is due to two large opposing fluxes, one which is ‘down-gradient’ or outward from the wall, where most vorticity concentrates, and the other which is ‘up-gradient’ or toward the wall and acting against strong viscous diffusion in the near-wall region. We present evidence that the up-gradient/down-gradient transport occurs by a mechanism of correlated inflow/outflow and spanwise vortex stretching/contraction that was proposed by Lighthill. This mechanism is essentially Lagrangian, but we explicate its relation to the Eulerian anti-symmetric vorticity flux tensor. As evidence for the mechanism, we study (i) statistical correlations of the wall-normal velocity and of wall-normal flux of spanwise vorticity, (ii) vorticity flux cospectra identifying eddies involved in nonlinear vorticity transport in the two opposing directions and (iii) visualizations of coherent vortex structures which contribute to the transport. The ‘D-type’ vortices contributing to down-gradient transport in the log layer are found to be attached, hairpin-type vortices. However, the ‘U-type’ vortices contributing to up-gradient transport are detached, wall-parallel, pancake-shaped vortices with strong spanwise vorticity, as expected by Lighthill's mechanism. We discuss modifications to the attached eddy model and implications for turbulent drag reduction.

Aiming to study the rough-wall turbulent boundary layer structure over differently arranged roughness elements, an experimental study was conducted on flows with regular and random roughness. Varying planform densities of truncated cone roughness elements in a square staggered pattern were investigated. The same planform densities were also investigated in random arrangements. Velocity statistics were measured via two-component laser Doppler velocimetry and stereoscopic particle image velocimetry. Friction velocity, thickness, roughness length and zero-plane displacement, determined from spatially averaged flow statistics, showed only minor differences between the regular and random arrangements at the same density. Recent a priori morphometric and statistical drag prediction methods were evaluated against experimentally determined roughness length. Observed differences between regular and random surface flow parameters were due to the presence of secondary flows which manifest as high-momentum pathways and low-momentum pathways in the streamwise velocity. Contrary to expectation, these secondary flows were present over the random surfaces and not discernible over the regular surfaces. Previously identified streamwise-coherent spanwise roughness heterogeneity does not seem to be present, suggesting that such roughness heterogeneity is not necessary to sustain secondary flows. Evidence suggests that the observed secondary flows were initiated at the front edge of the roughness and sustained over irregular roughness. Due to the secondary flows, local turbulent boundary layer profiles do not scale with local wall shear stress but appear to scale with local turbulent shear stress above the roughness canopy. Additionally, quadrant analysis shows distinct changes in the populations of ejection and sweep events.

A large eddy simulation wall model is developed based on a formal interpretation of quasi-equilibrium that governs the momentum balance integrated in the wall-normal direction. The model substitutes the law-of-the-wall velocity profile for a smooth surface into the wall-normal integrated momentum balance, leading to a Lagrangian relaxation towards equilibrium (LaRTE) transport equation for the friction–velocity vector ${\boldsymbol u}_\tau (x,z,t)$. This partial differential equation includes a relaxation time scale governing the rate at which the wall stress can respond to imposed fluctuations due to the inertia of the fluid layer from the wall to the wall-model height. A priori tests based on channel flow direct numerical simulation (DNS) data show that the identified relaxation time scale ensures self-consistency with assumed quasi-equilibrium conditions. The new approach enables us to formally distinguish quasi-equilibrium from additional, non-equilibrium contributions to the wall stress. A particular model for non-equilibrium contributions is derived, motivated by laminar Stokes layer dynamics in the viscous sublayer when applying fast-varying pressure gradients. The new wall modelling approach is first tested in standard equilibrium channel flow in order to document various properties of the approach. The model is then applied in large eddy simulation of channel flow with a suddenly applied spanwise pressure gradient (SSPG). The resulting mean wall-stress evolution is compared with DNS with good agreement. At the onset of the SSPG, the laminar Stokes layer develops rapidly while the LaRTE portion of the stress has a delayed response due to its inherent relaxation dynamics. Results also highlight open challenges such as modelling the response of near-wall turbulence occurring above the viscous sublayer and at time scales faster than quasi-equilibrium conditions.

The spatio-temporal dynamics of separation bubbles induced to form in a fully developed turbulent boundary layer (with Reynolds number based on momentum thickness of the boundary layer of 490) over a flat plate is studied via direct numerical simulations. Two different separation bubbles are examined: one induced by a suction–blowing velocity profile on the top boundary and the other by a suction-only velocity profile. The latter condition allows reattachment to occur without an externally imposed favourable pressure gradient and leads to a separation bubble more representative of those occurring over airfoils and in diffusers. The suction-only separation bubble exhibits a range of clearly distinguishable modes, including a high-frequency mode and a low-frequency ‘breathing’ mode that has been observed in some previous experiments. The high-frequency mode is well characterized by classical frequency scalings for a plane mixing layer and is associated with the formation and shedding of spanwise-oriented vortex rollers. The topology associated with the low-frequency motion is revealed by applying dynamic mode decomposition to the data from the simulations and is shown to be dominated by highly elongated structures in the streamwise direction. The possibility of Görtler instability induced by the streamwise curvature on the upstream end of the separation bubble as the underlying mechanism for these structures and the associated low frequency is explored.

The recent Lagrangian relaxation towards equilibrium (LaRTE) approach (Fowler et al., J. Fluid Mech., vol. 934, 2022, A44) is a wall model for large-eddy simulations (LES) that isolates quasi-equilibrium wall-stress dynamics from non-equilibrium responses to time-varying LES inputs. Non-equilibrium physics can then be modelled separately, such as the laminar Stokes layers that form in the viscous region and generate rapid wall-stress responses to fast changes in the pressure gradient. To capture additional wall-stress contributions due to near-wall turbulent eddies, a model term motivated by the attached eddy hypothesis is proposed. The total modelled wall stress thus includes contributions from various processes operating at different time scales (i.e. the LaRTE quasi-equilibrium plus laminar and turbulent non-equilibrium wall stresses) and is called the multi-time-scale (MTS) wall model. It is applied in LES of turbulent channel flow subject to a wide range of unsteady conditions from quasi-equilibrium to non-equilibrium. Flows considered include pulsating and linearly accelerating channel flow for several forcing frequencies and acceleration rates, respectively. We also revisit the sudden spanwise pressure gradient flow (considered in Fowler et al., J. Fluid Mech., vol. 934, 2022, A44) to review how the newly introduced model features affect this flow. Results obtained with the MTS wall model show good agreement with direct numerical simulation data over a vast range of conditions in these various non-equilibrium channel flows. To further understand the MTS model, we also describe and test the instantaneous-equilibrium limit of the MTS wall model. In this limit, good wall-stress predictions are obtained with reduced model complexity but providing less complete information about the wall-stress physics.