The objective of the present study is to provide a numerical database of thermal boundary layers and to contribute to the understanding of the dynamics of passive scalars at different Prandtl numbers. In this regard, a direct numerical simulation (DNS) of an incompressible zero-pressure-gradient turbulent boundary layer is performed with the Reynolds number based on momentum thickness ${\textit {Re}}_{\theta }$ ranging up to $1080$. Four passive scalars, characterized by the Prandtl numbers ${\textit {Pr}} = 1,2,4,6$ are simulated using the pseudo-spectral code SIMSON (Chevalier et al., SIMSON : a pseudo-spectral solver for incompressible boundary layer flows. Tech. Rep. TRITA-MEK 2007:07. KTH Mechanics, Stockholm, Sweden, 2007). To the best of our knowledge, the present DNS provides the thermal boundary layer with the highest Prandtl number available in the literature. It corresponds to that of water at ${\sim }24\,^{\circ }{\rm C}$, when the fluid temperature is considered as a passive scalar. Turbulence statistics for the flow and thermal fields are computed and compared with available numerical simulations at similar Reynolds numbers. The mean flow and scalar profiles, root-mean-squared velocity and scalar fluctuations, turbulent heat flux, turbulent Prandtl number and higher-order statistics agree well with the numerical data reported in the literature. Furthermore, the pre-multiplied two-dimensional spectra of the velocity and of the passive scalars are computed, providing a quantitative description of the energy distribution at the different length scales for various wall-normal locations. The energy distribution of the heat-flux fields at the wall is concentrated on longer temporal structures with increasing Prandtl number. This is due to the thinner thermal boundary layer as thermal diffusivity decreases and, thereby, the longer temporal structures exhibit a different footprint at the wall.

In the context of linear stability analysis, considering unsteady base flows is notoriously difficult. A generalisation of modal linear stability analysis, allowing for arbitrarily unsteady base flows over a finite time, is therefore required. The recently developed optimally time-dependent (OTD) modes form a projection basis for the tangent space. They capture the leading amplification directions in state space under the constraint that they form an orthonormal basis at all times. The present numerical study illustrates the possibility to describe a complex flow case using the leading OTD modes. The flow under investigation is an unsteady case of the Blasius boundary layer, featuring streamwise streaks of finite length and relevant to bypass transition. It corresponds to the state space trajectory initiated by the minimal seed; such a trajectory is unsteady, free from any spatial symmetry and shadows the laminar–turbulent separatrix for a finite time only. The finite-time instability of this unsteady base flow is investigated using the 8 leading OTD modes. The analysis includes the computation of finite-time Lyapunov exponents as well as instantaneous eigenvalues, and of the associated flow structures. The reconstructed instantaneous eigenmodes are all of outer type. They map unambiguously the spatial regions of largest instantaneous growth. Other flow structures, previously reported as secondary, are identified with this method as relevant to streak switching and to streamwise vortical ejections. The dynamics inside the tangent space features both modal and non-modal amplification. Non-normality within the reduced tangent subspace, quantified by the instantaneous numerical abscissa, emerges only as the unsteadiness of the base flow is reduced.

The presence of very-large-scale motions in wall-bounded turbulent flows is commonly associated with their footprint in the form of the superposition of the large scales at the wall and the additional amplitude modulation of small-scale near-wall turbulence. These two phenomena are currently understood to be interlinked, with the superposed large-scale velocity gradient causing the modulation of small-scale activity in the proximity of the wall. To challenge this idea, we devise a numerical strategy that selectively suppresses either superposition or amplitude modulation, in an effort to isolate and study the remaining phenomenon. Results from our direct numerical simulations indicate that a positive correlation between the amplitude of the small scales in the near-wall region and the large-scale signal in the outer flow persists even when near-wall large-scale motions are suppressed – i.e. in absence of superposition. Clearly, this kind of correlation cannot be caused by the near-wall large-scale velocity or its gradients, as both are absent. Conversely, when modulation is blocked, the near-wall footprints of the large scales seem to disappear. This study has been carried out on channel flows at friction Reynolds number $Re_\tau =1000$ in both standard simulation domains and minimal streamwise units (MSUs), where the streamwise fluctuation energy is enhanced. The consistency of the results obtained by the two approaches suggests that MSUs can capture correctly this kind of scale interaction at a much reduced cost.

Well-resolved direct numerical simulations (DNS) have been performed of the flow in a smooth circular pipe of radius $R$ and axial length $10{\rm \pi} R$ at friction Reynolds numbers up to $Re_\tau =5200$ using the pseudo-spectral code OPENPIPEFLOW. Various turbulence statistics are documented and compared with other DNS and experimental data in pipes as well as channels. Small but distinct differences between various datasets are identified. The friction factor $\lambda$ overshoots by $2\,\%$ and undershoots by $0.6\,\%$ the Prandtl friction law at low and high $Re$ ranges, respectively. In addition, $\lambda$ in our results is slightly higher than in Pirozzoli et al. (J. Fluid Mech., vol. 926, 2021, A28), but matches well the experiments in Furuichi et al. (Phys. Fluids, vol. 27, issue 9, 2015, 095108). The log-law indicator function, which is nearly indistinguishable between pipe and channel up to $y^+=250$, has not yet developed a plateau farther away from the wall in the pipes even for the $Re_\tau =5200$ cases. The wall shear stress fluctuations and the inner peak of the axial turbulence intensity – which grow monotonically with $Re_\tau$ – are lower in the pipe than in the channel, but the difference decreases with increasing $Re_\tau$. While the wall value is slightly lower in the channel than in the pipe at the same $Re_\tau$, the inner peak of the pressure fluctuation shows negligible differences between them. The Reynolds number scaling of all these quantities agrees with both the logarithmic and defect-power laws if the coefficients are properly chosen. The one-dimensional spectrum of the axial velocity fluctuation exhibits a $k^{-1}$ dependence at an intermediate distance from the wall – also seen in the channel. In summary, these high-fidelity data enable us to provide better insights into the flow physics in the pipes as well as the similarity/difference among different types of wall turbulence.

In this study, a new well-resolved large-eddy simulation of an incompressible near-equilibrium adverse-pressure-gradient (APG) turbulent boundary layer (TBL) over a flat plate is presented. In this simulation, we have established a near-equilibrium APG over a wide Reynolds-number range. In this so-called region of interest, the Rotta–Clauser pressure-gradient parameter $\beta$ exhibits an approximately constant value of around 1.4, and the Reynolds number based on momentum thickness reaches ${\textit {Re}}_{\theta }=8700$. To the best of the authors’ knowledge, this is to date the highest ${\textit {Re}}_{\theta }$ achieved for a near-equilibrium APG TBL under an approximately constant moderate APG. We evaluated the self-similarity of the outer region using two scalings, namely the Zagarola–Smits and an alternative scaling based on edge velocity and displacement thickness. Our results reveal that outer-layer similarity is achieved, and the viscous scaling collapses the near-wall region of the mean flow in agreement with classical theory. Spectral analysis reveals that the APG displaces some small-scale energy from the near-wall to the outer region, an effect observed for all the components of the Reynolds-stress tensor, which becomes more evident at higher Reynolds numbers. In general, the effects of the APG are more noticeable at lower Reynolds numbers. For instance, the outer peak of turbulent-kinetic-energy (TKE) production is less prominent at higher $Re$. Although the scale separation increases with ${\textit {Re}}$ in zero-pressure-gradient TBLs, this effect becomes accentuated by the APG. Despite the reduction of the outer TKE production at higher Reynolds numbers, the mechanisms of energisation of large scales are still present.

The application of drag-control strategies on canonical wall-bounded turbulence, such as periodic channel and zero- or adverse-pressure-gradient boundary layers, raises the question on how to distinguish consistently the origin of control effects under different reference conditions. We employ the RD identity (Renard & Deck, J. Fluid Mech., vol. 790, 2016, pp. 339–367) to decompose the mean friction drag and investigate the control effects of uniform blowing and suction applied to an NACA4412 airfoil at chord Reynolds numbers $Re_c=200\,000$ and $400\,000$. The connection of the drag reduction/increase by using blowing/suction with the turbulence statistics (including viscous dissipation, turbulence kinetic energy production and spatial growth of the flow) across the boundary layer, subjected to adverse or favourable pressure gradients, is examined. We found that the inner and outer peaks of the contributions associated with the friction-drag generation show good scaling with either inner or outer units, respectively. They are also independent of the Reynolds number, control scheme and intensity of the blowing/suction. The small- and large-scale structures are separated with an adaptive scale-decomposition method, namely the empirical mode decomposition (EMD), which aims to analyse the scale-specific contribution of turbulent motions to friction-drag generation. Results unveil that blowing on the suction side of the airfoil is able to enhance the contribution of large-scale motions and to suppress that of small scales; however, suction behaves contrarily. The contributions related to cross-scale interactions remain almost unchanged with different control strategies.

Two models based on convolutional neural networks are trained to predict the two-dimensional instantaneous velocity-fluctuation fields at different wall-normal locations in a turbulent open-channel flow, using the wall-shear-stress components and the wall pressure as inputs. The first model is a fully convolutional neural network (FCN) which directly predicts the fluctuations, while the second one reconstructs the flow fields using a linear combination of orthonormal basis functions, obtained through proper orthogonal decomposition (POD), and is hence named FCN-POD. Both models are trained using data from direct numerical simulations at friction Reynolds numbers $Re_{\tau } = 180$ and 550. Being able to predict the nonlinear interactions in the flow, both models show better predictions than the extended proper orthogonal decomposition (EPOD), which establishes a linear relation between the input and output fields. The performance of the models is compared based on predictions of the instantaneous fluctuation fields, turbulence statistics and power-spectral densities. FCN exhibits the best predictions closer to the wall, whereas FCN-POD provides better predictions at larger wall-normal distances. We also assessed the feasibility of transfer learning for the FCN model, using the model parameters learned from the $Re_{\tau }=180$ dataset to initialize those of the model that is trained on the $Re_{\tau }=550$ dataset. After training the initialized model at the new $Re_{\tau }$, our results indicate the possibility of matching the reference-model performance up to $y^{+}=50$, with $50\,\%$ and $25\,\%$ of the original training data. We expect that these non-intrusive sensing models will play an important role in applications related to closed-loop control of wall-bounded turbulence.

We investigate spanwise-coherent structures in the turbulent flow around airfoils, motivated by their connection with trailing-edge noise. We analyse well-resolved large-eddy simulations (LES) of the flow around NACA 0012 and NACA 4412 airfoils, both at a Reynolds number of 400 000 based on the chord length. Spectral proper orthogonal decomposition performed on the data reveals that the most energetic coherent structures are hydrodynamic waves, extending over the turbulent boundary layers around the airfoils with significant amplitudes near the trailing edge. Resolvent analysis was used to model such structures, using the mean field as a base flow. We then focus on evaluating the dependence of such structures on the domain size, to ensure that they are not an artefact of periodic boundary conditions in small computational boxes. To this end, we performed incompressible LES of a zero-pressure-gradient turbulent boundary layer, for three different spanwise sizes, with the momentum-thickness Reynolds number matching those near the airfoils trailing edge. The same coherent hydrodynamic waves were observed for the three domains. Such waves are accurately modelled as the most amplified flow response from resolvent analysis. The signature of such wide structures is seen in non-premultiplied spanwise wavenumber spectra, which collapse for the three computational domains. These results suggest that the spanwise-elongated structures are not domain-size dependent for the studied simulations, indicating thus the presence of very wide structures in wall-bounded turbulent flows.

Direct numerical simulations, performed with a high-order spectral-element method, are used to study coherent structures in turbulent pipe flow at friction Reynolds numbers $Re_{\tau } = 180$ and $550$. The database was analysed using spectral proper orthogonal decomposition (SPOD) to identify energetically dominant coherent structures, most of which turn out to be streaks and quasi-streamwise vortices. To understand how such structures can be modelled, the linear flow responses to harmonic forcing were computed using the singular value decomposition of the resolvent operator, using the mean field as a base flow. The SPOD and resolvent analysis were calculated for several combinations of frequencies and wavenumbers, allowing the mapping out of similarities between SPOD modes and optimal responses for a wide range of relevant scales in turbulent pipe flows. In order to explore physical reasons behind the agreement between both methods, an indicator of lift-up mechanism in the resolvent analysis was introduced, activated when optimal forcing is dominated by the wall-normal and azimuthal components, and associated response corresponds to streaks of streamwise velocity. Good agreement between leading SPOD and resolvent modes is observed in a large region of parameter space. In this region, a significant gain separation is found in resolvent analysis, which may be attributed to the strong amplification associated with the lift-up mechanism, here understood as nonlinear forcing terms leading to the appearance of streamwise vortices, which in turn form high-amplitude streaks. For both Reynolds numbers, the observed concordances were generally for structures with large energy in the buffer layer. The results highlight resolvent analysis as a pertinent reduced-order model for coherent structures in wall-bounded turbulence, particularly for streamwise elongated structures corresponding to near-wall streamwise vortices and streaks.

Recent progress in understanding subcritical transition to turbulence is based on the concept of the edge, the manifold separating the basins of attraction of the laminar and the turbulent state. Originally developed in numerical studies of parallel shear flows with a linearly stable base flow, this concept is adapted here to the case of a spatially developing Blasius boundary layer. Longer time horizons fundamentally change the nature of the problem due to the loss of stability of the base flow due to Tollmien–Schlichting (TS) waves. We demonstrate, using a moving box technique, that efficient long-time tracking of edge trajectories is possible for the parameter range relevant to bypass transition, even if the asymptotic state itself remains out of reach. The flow along the edge trajectory features streak switching observed for the first time in the Blasius boundary layer. At long enough times, TS waves co-exist with the coherent structure characteristic of edge trajectories. In this situation we suggest a reinterpretation of the edge as a manifold dividing the state space between the two main types of boundary layer transition, i.e. bypass transition and classical transition.

Isolated patches of turbulence in transitional straight pipes are sustained by a strong instability at their upstream front, where the production of turbulent kinetic energy (TKE) is up to five times higher than in the core. Direct numerical simulations presented in this paper show no evidence of such strong fronts if the pipe is bent. We examine the temporal and spatial evolution of puffs and slugs in a toroidal pipe with pipe-to-torus diameter ratio $\unicode[STIX]{x1D6FF}=D/d=0.01$ at several subcritical Reynolds numbers. Results show that the upstream overshoot of TKE production is at most one-and-a-half times the value in the core and that the average cross-flow fluctuations at the front are up to three times lower if compared to a straight pipe, while attaining similar values in the core. Localised turbulence can be sustained at smaller energies through a redistribution of turbulent fluctuations and vortical structures by the in-plane Dean motion of the mean flow. This asymmetry determines a strong localisation of TKE production near the outer bend, where linear and nonlinear mechanisms optimally amplify perturbations. We further observe a substantial reduction of the range of Reynolds numbers for long-lived intermittent turbulence, in agreement with experimental data from the literature. Moreover, no occurrence of nucleation of spots through splitting could be detected in the range of parameters considered. Based on the present results, we argue that this mechanism gradually becomes marginal as the curvature of the pipe increases and the transition scenario approaches a dynamical switch from subcritical to supercritical.

Three methods are evaluated to estimate the streamwise velocity fluctuations of a zero-pressure-gradient turbulent boundary layer of momentum-thickness-based Reynolds number up to $Re_{\unicode[STIX]{x1D703}}\simeq 8200$, using as input velocity fluctuations at different wall-normal positions. A system identification approach is considered where large-eddy simulation data are used to build single and multiple-input linear and nonlinear transfer functions. Such transfer functions are then treated as convolution kernels and may be used as models for the prediction of the fluctuations. Good agreement between predicted and reference data is observed when the streamwise velocity in the near-wall region is estimated from fluctuations in the outer region. Both the unsteady behaviour of the fluctuations and the spectral content of the data are properly predicted. It is shown that approximately 45 % of the energy in the near-wall peak is linearly correlated with the outer-layer structures, for the reference case $Re_{\unicode[STIX]{x1D703}}=4430$. These identified transfer functions allow insight into the causality between the different wall-normal locations in a turbulent boundary layer along with an estimation of the tilting angle of the large-scale structures. Differences in accuracy of the methods (single- and multiple-input linear and nonlinear) are assessed by evaluating the coherence of the structures between wall-normally separated positions. It is shown that the large-scale fluctuations are coherent between the outer and inner layers, by means of an interactions which strengthens with increasing Reynolds number, whereas the finer-scale fluctuations are only coherent within the near-wall region. This enables the possibility of considering the wall-shear stress as an input measurement, which would more easily allow the implementation of these methods in experimental applications. A parametric study was also performed by evaluating the effect of the Reynolds number, wall-normal positions and input quantities considered in the model. Since the methods vary in terms of their complexity for implementation, computational expense and accuracy, the technique of choice will depend on the application under consideration. We also assessed the possibility of designing and testing the models at different Reynolds numbers, where it is shown that the prediction of the near-wall peak from wall-shear-stress measurements is practically unaffected even for a one order of magnitude change in the corresponding Reynolds number of the design and test, indicating that the interaction between the near-wall peak fluctuations and the wall is approximately Reynolds-number independent. Furthermore, given the performance of such methods in the prediction of flow features in turbulent boundary layers, they have a good potential for implementation in experiments and realistic flow control applications, where the prediction of the near-wall peak led to correlations above 0.80 when wall-shear stress was used in a multiple-input or nonlinear scheme. Errors of the order of 20 % were also observed in the determination of the near-wall spectral peak, depending on the employed method.

Motivated by the relevance of edge state solutions as mediators of transition, we use direct numerical simulations to study the effect of spatially non-uniform viscosity on their energy and stability in minimal channel flows. What we seek is a theoretical support rooted in a fully nonlinear framework that explains the modified threshold for transition to turbulence in flows with temperature-dependent viscosity. Consistently over a range of subcritical Reynolds numbers, we find that decreasing viscosity away from the walls weakens the streamwise streaks and the vortical structures responsible for their regeneration. The entire self-sustained cycle of the edge state is maintained on a lower kinetic energy level with a smaller driving force, compared to a flow with constant viscosity. Increasing viscosity away from the walls has the opposite effect. In both cases, the effect is proportional to the strength of the viscosity gradient. The results presented highlight a local shift in the state space of the position of the edge state relative to the laminar attractor with the consequent modulation of its basin of attraction in the proximity of the edge state and of the surrounding manifold. The implication is that the threshold for transition is reduced for perturbations evolving in the neighbourhood of the edge state in the case that viscosity decreases away from the walls, and vice versa.

Swirl-switching is a low-frequency oscillatory phenomenon which affects the Dean vortices in bent pipes and may cause fatigue in piping systems. Despite thirty years worth of research, the mechanism that causes these oscillations and the frequencies that characterise them remain unclear. Here we show that a three-dimensional wave-like structure is responsible for the low-frequency switching of the dominant Dean vortex. The present study, performed via direct numerical simulation, focuses on the turbulent flow through a $90^{\circ }$ pipe bend preceded and followed by straight pipe segments. A pipe with curvature 0.3 (defined as ratio between pipe radius and bend radius) is studied for a bulk Reynolds number $Re=11\,700$, corresponding to a friction Reynolds number $Re_{\unicode[STIX]{x1D70F}}\approx 360$. Synthetic turbulence is generated at the inflow section and used instead of the classical recycling method in order to avoid the interference between recycling and swirl-switching frequencies. The flow field is analysed by three-dimensional proper orthogonal decomposition (POD) which for the first time allows the identification of the source of swirl-switching: a wave-like structure that originates in the pipe bend. Contrary to some previous studies, the flow in the upstream pipe does not show any direct influence on the swirl-switching modes. Our analysis further shows that a three-dimensional characterisation of the modes is crucial to understand the mechanism, and that reconstructions based on two-dimensional POD modes are incomplete.

With the motivation of determining the critical roughness size, a global stability and sensitivity analysis of a three-dimensional Falkner–Skan–Cooke (FSC) boundary layer with a cylindrical surface roughness is performed. The roughness size is chosen such that breakdown to turbulence is initiated by a global version of traditional secondary instabilities of the cross-flow (CF) vortices instead of an immediate flow tripping at the roughness. The resulting global eigenvalue spectra of the systems are found to be very sensitive to numerical parameters and domain size. This sensitivity to numerical parameters is quantified using the $\unicode[STIX]{x1D700}$-pseudospectrum, and the dependency on the domain is analysed through an impulse response, structural sensitivity analysis and an energy budget. It is shown that while the frequencies remain relatively unchanged, the growth rates increase with domain size, which originates from the inclusion of stronger CF vortices in the baseflow. This is reflected in a change in the rate of advective energy transport by the baseflow. It is concluded that the onset of global instability in a FSC boundary layer as the roughness height is increased does not correspond to an immediate flow tripping behind the roughness, but occurs for lower roughness heights if sufficiently long domains are considered. However, the great sensitivity results in an inability to accurately pinpoint the exact parameter values for the bifurcation, and the large spatial growth of the disturbances in the long domains eventually becomes larger than can be resolved using finite-precision arithmetic.

Direct numerical simulations (DNS) of particle-laden turbulent flow in straight, mildly curved and strongly bent pipes are performed in which the solid phase is modelled as small heavy spherical particles. A total of seven populations of dilute particles with different Stokes numbers, one-way coupled with their carrier phase, are simulated. The objective is to examine the effect of the curvature on micro-particle transport and accumulation. It is shown that even a slight non-zero curvature in the flow configuration strongly impact the particle concentration map such that the concentration of inertial particles with bulk Stokes number $0.45$ (based on bulk velocity and pipe radius) at the inner bend wall of mildly curved pipe becomes $12.8$ times larger than that in the viscous sublayer of the straight pipe. Near-wall helicoidal particle streaks are observed in the curved configurations with their inclination varying with the strength of the secondary motion of the carrier phase. A reflection layer, as previously observed in particle laden turbulent S-shaped channels, is also apparent in the strongly curved pipe with heavy particles. In addition, depending on the curvature, the central regions of the mean Dean vortices appear to be completely depleted of particles, as observed also in the partially relaminarised region at the inner bend. The turbophoretic drift of the particles is shown to be affected by weak and strong secondary motions of the carrier phase and geometry-induced centrifugal forces. The first- and second-order moments of the velocity and acceleration of the particulate phase in the same configurations are addressed in a companion paper by the same authors. The current data set will be useful for modelling particles advected in wall-bounded turbulent flows where the effects of the curvature are not negligible.

The modal instability encountered by the incompressible flow inside a toroidal pipe is studied, for the first time, by means of linear stability analysis and direct numerical simulation (DNS). In addition to the unquestionable aesthetic appeal, the torus represents the smallest departure from the canonical straight pipe flow, at least for low curvatures. The flow is governed by only two parameters: the Reynolds number $\mathit{Re}$ and the curvature of the torus ${\it\delta}$, i.e. the ratio between pipe radius and torus radius. The absence of additional features, such as torsion in the case of a helical pipe, allows us to isolate the effect that the curvature has on the onset of the instability. Results show that the flow is linearly unstable for all curvatures investigated between 0.002 and unity, and undergoes a Hopf bifurcation at $\mathit{Re}$ of about 4000. The bifurcation is followed by the onset of a periodic regime, characterised by travelling waves with wavelength $\mathit{O}(1)$ pipe diameters. The neutral curve associated with the instability is traced in parameter space by means of a novel continuation algorithm. Tracking the bifurcation provides a complete description of the modal onset of instability as a function of the two governing parameters, and allows a precise calculation of the critical values of $\mathit{Re}$ and ${\it\delta}$. Several different modes are found, with differing properties and eigenfunction shapes. Some eigenmodes are observed to belong to groups with a set of common characteristics, deemed ‘families’, while others appear as ‘isolated’. Comparison with nonlinear DNS shows excellent agreement, confirming every aspect of the linear analysis, its accuracy, and proving its significance for the nonlinear flow. Experimental data from the literature are also shown to be in considerable agreement with the present results.

Recent experimental studies have shown rich transition behaviour in rotating plane Couette flow (RPCF). In this paper we study the transition in supercritical RPCF theoretically by determination of equilibrium and periodic orbit tertiary states via Floquet analysis on secondary Taylor vortex solutions. Two new tertiary states are discovered which we name oscillatory wavy vortex flow (oWVF) and skewed vortex flow (SVF). We present the bifurcation routes and stability properties of these new tertiary states and, in addition, we describe a bifurcation procedure whereby a set of defected wavy twist vortices is approached. Further to this, transition scenarios at flow parameters relevant to experimental works are investigated by computation of the set of stable attractors which exist on a large domain. The physically observed flow states are shown to share features with states in our set of attractors.

Two modal decomposition techniques are employed to analyse the stability of wind turbine wakes. A numerical study on a single wind turbine wake is carried out focusing on the instability onset of the trailing tip vortices shed from the turbine blades. The numerical model is based on large-eddy simulations (LES) of the Navier–Stokes equations using the actuator line (ACL) method to simulate the wake behind the Tjæreborg wind turbine. The wake is perturbed by low-amplitude excitation sources located in the neighbourhood of the tip spirals. The amplification of the waves travelling along the spiral triggers instabilities, leading to breakdown of the wake. Based on the grid configurations and the type of excitations, two basic flow cases, symmetric and asymmetric, are identified. In the symmetric setup, we impose a 120° symmetry condition in the dynamics of the flow and in the asymmetric setup we calculate the full 360° wake. Different cases are subsequently analysed using dynamic mode decomposition (DMD) and proper orthogonal decomposition (POD). The results reveal that the main instability mechanism is dispersive and that the modal growth in the symmetric setup arises only for some specific frequencies and spatial structures, e.g. two dominant groups of modes with positive growth (spatial structures) are identified, while breaking the symmetry reveals that almost all the modes have positive growth rate. In both setups, the most unstable modes have a non-dimensional spatial growth rate close to $\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}}\pi /2$ and they are characterized by an out-of-phase displacement of successive helix turns leading to local vortex pairing. The present results indicate that the asymmetric case is crucial to study, as the stability characteristics of the flow change significantly compared to the symmetric configurations. Based on the constant non-dimensional growth rate of disturbances, we derive a new analytical relationship between the length of the wake up to the turbulent breakdown and the operating conditions of a wind turbine.

A recent assessment of available direct numerical simulation (DNS) data from turbulent boundary layer flows (Schlatter & Örlü, J. Fluid Mech., vol. 659, 2010, pp. 116–126) showed surprisingly large differences not only in the skin friction coefficient or shape factor, but also in their predictions of mean and fluctuation profiles far into the sublayer. While such differences are expected at very low Reynolds numbers and/or the immediate vicinity of the inflow or tripping region, it remains unclear whether inflow and tripping effects explain the differences observed even at moderate Reynolds numbers. This question is systematically addressed by re-simulating the DNS of a zero-pressure-gradient turbulent boundary layer flow by Schlatter et al. (Phys. Fluids, vol. 21, 2009, art. 051702). The previous DNS serves as the baseline simulation, and the new DNS with a range of physically different inflow conditions and tripping effects are carefully compared. The downstream evolution of integral quantities as well as mean and fluctuation profiles is analysed, and the results show that different inflow conditions and tripping effects do indeed explain most of the differences observed when comparing available DNS at low Reynolds number. It is further found that, if transition is initiated inside the boundary layer at a low enough Reynolds number (based on the momentum-loss thickness) , all quantities agree well for both inner and outer layer for . This result gives a lower limit for meaningful comparisons between numerical and/or wind tunnel experiments, assuming that the flow was not severely over- or understimulated. It is further shown that even profiles of the wall-normal velocity fluctuations and Reynolds shear stress collapse for higher irrespective of the upstream conditions. In addition, the overshoot in the total shear stress within the sublayer observed in the DNS of Wu & Moin (Phys. Fluids, vol. 22, 2010, art. 085105) has been identified as a feature of transitional boundary layers.