This is the accepted version of the paper.This version of the publication may differ from the final published version. Permanent repository link:

The statistical properties are presented for the direct numerical simulation (DNS) of a self-similar adverse pressure gradient (APG) turbulent boundary layer (TBL) at the verge of separation. The APG TBL has a momentum thickness based Reynolds number range from Re δ2 = 570 to 13800, with a self-similar region from Re δ2 = 10000 to 12300. Within this domain the average non-dimensional pressure gradient parameter β = 39, where for a unit density β = δ 1 P ′ e /τ w , with δ 1 the displacement thickness, τ w the mean shear stress at the wall, and P ′ e the farfield pressure gradient. This flow is compared to previous zero pressure gradient (ZPG) and mild APG TBL (β = 1) results of similar Reynolds number. All flows are generated via the DNS of a TBL on a flat surface with farfield boundary conditions tailored to apply the desired pressure gradient. The conditions for self-similarity, and the appropriate length and velocity scales are derived. The mean and Reynolds stress profiles are shown to collapse when non-dimensionalised on the basis of these length and velocity scales. As the pressure gradient increases, the extent of the wake region in the mean streamwise velocity profiles increases, whilst the extent of the log-layer and viscous sub-layer decreases. The Reynolds stress, production and dissipation profiles of the APG TBL cases exhibit a second outer peak, which becomes more pronounced and more spatially localised with increasing pressure gradient. This outer peak is located at the point of inflection of the mean velocity profiles, and is suggestive of the presence of a shear flow instability. The maximum streamwise velocity variance is located at a wall normal position of δ 1 of spanwise wavelength of 2δ 1 . In summary as the pressure gradient increases the flow has properties less like a ZPG TBL and more akin to a free shear layer.

Direct numerical simulation of turbulent flow in a straight square duct was performed in order to determine the minimal requirements for self-sustaining turbulence. It was found that turbulence can be maintained for values of the bulk Reynolds number above approximately 1100, corresponding to a friction-velocity-based Reynolds number of 80. The minimum value for the streamwise period of the computational domain measures around 190 wall units, roughly independently of the Reynolds number. Furthermore, we present a characterization of the flow state at marginal Reynolds numbers which substantially differs from the fully turbulent one. In particular, the marginal state exhibits a 4-vortex secondary flow structure alternating in time whereas the fully turbulent one presents the usual 8-vortex pattern. It is shown that in the regime of marginal Reynolds numbers buffer layer coherent structures play a crucial role in the appearance of secondary flow of Prandtl's second kind.

Statistically stationary and homogeneous shear turbulence (SS-HST) is investigated by means of a new direct numerical simulation code, spectral in the two horizontal directions and compact-finite-differences in the direction of the shear. No remeshing is used to impose the shear-periodic boundary condition. The influence of the geometry of the computational box is explored. Since HST has no characteristic outer length scale and tends to fill the computational domain, long-term simulations of HST are 'minimal' in the sense of containing on average only a few large-scale structures. It is found that the main limit is the spanwise box width, L z , which sets the length and velocity scales of the turbulence, and that the two other box dimensions should be sufficiently large (L x 2L z , L y L z ) to prevent other directions to be constrained as well. It is also found that very long boxes, L x 2L y , couple with the passing period of the shear-periodic boundary condition, and develop strong unphysical linearized bursts. Within those limits, the flow shows interesting similarities and differences with other shear flows, and in particular with the logarithmic layer of wall-bounded turbulence. They are explored in some detail. They include a self-sustaining process for large-scale streaks and quasi-periodic bursting. The bursting time scale is approximately universal, ∼ 20S −1 , and the availability of two different bursting systems allows the growth of the bursts to be related with some confidence to the shearing of initially isotropic turbulence. It is concluded that SS-HST, conducted within the proper computational parameters, is a very promising system to study shear turbulence in general.

The three-dimensional vortex clusters, and the structures based on the quadrant classification of the intense tangential Reynolds stress (Qs), are studied in direct numerical simulations of statistically stationary homogeneous shear turbulence (HST) at Taylor microscale Reynolds number Re λ ≈ 50-250, with emphasis on comparisons with turbulent channels (CHs). The Qs and vortex clusters in HST are found to be versions of the corresponding detached (in the sense of del Álamo et al. (J. Fluid Mech., vol. 561 (2006), pp. 329-358)) structures in CHs, although statistically symmetrised with respect to the substitution of sweeps by ejections and vice versa. In turn, these are more symmetric versions of the corresponding attached Qs and clusters. In both flows, only co-gradient sweeps and ejections larger than the local Corrsin scale are found to couple with the shear. They are oriented anisotropically, and are responsible for carrying most of the total Reynolds stress. Most large eddies in CHs are attached to the wall, but it is shown that this is probably a geometric consequence of their size, rather than the reason for their dynamical significance. Most small Q structures associated with different quadrants are far from each other in comparison to their size, but those that are close to each other tend to form quasi-streamwise trains of groups of a sweep and an ejection paired side by side in the spanwise direction, with a vortex cluster in between, generalising to three dimensions the corresponding arrangement of attached eddies in CHs. These pairs are organised around an inclined large-scale conditional vortex 'roller', and it is shown that the composite structure tends to be located at the interface between high-and low-velocity streaks, as well as in strong 'co-gradient' shear layers that separate streaks of either sign in which velocity is more uniform. It is further found that the conditional rollers are terminated by 'hooks' reminiscent of hairpins, both upright and inverted. The inverted hook weakens as the structures approach the wall, while the upright one changes little. At the same time, the inclination of the roller with respect to the mean velocity decreases from 45 • in HST to quasi-streamwise for † Email address for correspondence: jimenez@torroja.