Conditionally-averaged structures in wall-bounded turbulent flows

The quadrant-splitting and the wall-shear detection techniques were used to obtain ensemble-averaged wall layer structures. The two techniques give similar results for Q4 events, but the wall-shear method leads to smearing of Q2 events. Events were found to maintain their identity for very long times. The ensemble-averaged structures scale with outer variables. Turbulence producing events were associated with one dominant vortical structure rather than a pair of counter-rotating structures. An asymmetry-preserving averaging scheme was devised that allowed a picture of the average structure which more closely resembles the instantaneous one, to be obtained

On the large-eddy simulation of transitional wall-bounded flows

The structure of the subgrid scale fields in plane channel flow has been studied at various stages of the transition process to turbulence. The residual stress and subgrid scale dissipation calculated using velocity fields generated by direct numerical simulations of the Navier-Stokes equations are significantly different from their counterparts in turbulent flows. The subgrid scale dissipation changes sign over extended areas of the channel, indicating energy flow from the small scales to the large scales. This reversed energy cascade becomes less pronounced at the later stages of transition. Standard residual stress models of the Smagorinsky type are excessively dissipative. Rescaling the model constant improves the prediction of the total (integrated) subgrid scale dissipation, but not that of the local one. Despite the somewhat excessive dissipation of the rescaled Smagorinsky model, the results of a large eddy simulation of transition on a flat-plate boundary layer compare quite well with those of a direct simulation, and require only a small fraction of the computational effort. The inclusion of non-dissipative models, which could lead to further improvements, is proposed

Numerical simulation of roughness effects on the flow past a circular cylinder

In the present work large eddy simulations of the flow past a rough cylinder are performed at a Reynolds number of Re = 4.2 × 105 and an equivalent sand-grain surface roughness height ks = 0.02D. In order to determine the effects of the surface roughness on the boundary layer transition and as a consequence on the wake topology, results are compared to those of the smooth cylinder. It is shown that surface roughness triggers the transition to turbulence in the boundary layer, thus leading to an early separation caused by the increased drag and momentum deficit. Thus, the drag coefficient increases up to CD 1.122 (if compared to the smooth cylinder it should be about CD 0.3 - 0.5). The wake topology also changes and resembles more the subcritical wake observed for the smooth cylinder at lower Reynolds numbers than the expected critical wake at this Reynolds number.Peer ReviewedPostprint (author's final draft

Numerical simulation of pulsating turbulent channel flow

Direct and large-eddy simulations of the Navier–Stokes equations are used to study the pulsating flow in a channel. The cases examined span a wide range of frequencies of the driving pressure gradient, and encompass different physical behaviors, from the quasi-Stokes flow observed at high frequencies, to a quasisteady behavior at the lowest ones. The validity of the dynamic Smagorinsky model to study this kind of unsteady flow is established by a posteriori comparison with direct simulations and experimental data. It is shown that the fluctuations generated in the near-wall region by the unsteady pressure gradient do not propagate beyond a certain distance l t from the wall, which can be estimated quite accurately by a simple eddyviscosity argument. No substantial departure from the Stokes regime at very high frequency (ω + as high as 0.1) is observed. The time-dependent characteristics of the flow are examined in detail, as well as the topology of the coherent structures

Large-eddy simulation of three-dimensional dunes in a steady, unidirectional flow. Part 2. Flow structures

We performed large-eddy simulations of the flow over a series of three-dimensional (3D) dunes at laboratory scale. The bedform three-dimensionality was imposed by shifting a standard two-dimensional (2D) dune shape in the streamwise direction according to a sine wave. The turbulence statistics were discussed in Part 1 of this article (Omidyeganeh & Piomelli, J. Fluid Mech., vol. 721, 2013, pp. 454–483). Coherent flow structures and their statistics are discussed concentrating on two cases with the same crestline amplitudes and wavelengths but different crestline alignments: in-phase and staggered. The present paper shows that the induced large-scale mean streamwise vortices are the primary factor that alters the features of the instantaneous flow structures. Wall turbulence is insensitive to the crestline alignment; alternating high- and low-speed streaks appear in the internal boundary layer developing on the stoss side, whereas over the node plane (the plane normal to the spanwise direction at the node of the crestline), they are inclined towards the lobe plane (the plane normal to the spanwise direction at the most downstream point of the crestline) due to the mean spanwise pressure gradient. Spanwise vortices (rollers) generated by Kelvin–Helmholtz instability in the separated shear layer appear regularly over the lobe with much larger length scale than those over the saddle (the plane normal to the spanwise direction at the most upstream point of the crestline). Rollers over the lobe may extend to the saddle plane and affect the reattachment features; their shedding is more frequent than in 2D geometries. Vortices shed from the separated shear layer in the lobe plane undergo a three-dimensional instability while being advected downstream, and rise toward the free surface. They develop into a horseshoe shape (similar to the 2D case) and affect the whole channel depth, whereas those generated near the saddle are advected downstream and toward the bed. When the tip of such a horseshoe reaches the free surface, the ejection of flow at the surface causes ‘boils’ (upwelling events on the surface). Strong boil events are observed on the surface of the lobe planes of 3D dunes more frequently than in the saddle planes. They also appear more frequently than in the corresponding 2D geometry. The crestline alignment of the dune alters the dynamics of the flow structures, in that they appear in the lobe plane and are advected towards the saddle plane of the next dune, where they are dissipated. Boil events occur at a higher frequency in the staggered alignment, but with less intensity than in the in-phase alignment

Large-eddy simulation of three-dimensional dunes in a steady, unidirectional flow. Part 1. Turbulence statistics

We performed large-eddy simulations of flow over a series of three-dimensional dunes at laboratory scale (Reynolds number based on the average channel depth and streamwise velocity was 18 900) using the Lagrangian dynamic eddy-viscosity subgrid-scale model. The bedform three-dimensionality was imposed by shifting a standard two-dimensional dune shape in the streamwise direction according to a sine wave. The statistics of the flow are discussed in 10 cases with in-phase and staggered crestlines, different deformation amplitudes and wavelengths. The results are validated qualitatively against experiments. The three-dimensional separation of flow at the crestline alters the distribution of wall pressure, which in turn may cause secondary flow across the stream, which directs low-momentum fluid, near the bed, toward the lobe (the most downstream point on the crestline) and high-momentum fluid, near the top surface, toward the saddle (the most upstream point on the crestline). The mean flow is characterized by a pair of counter-rotating streamwise vortices, with core radius of the order of the flow depth. However, for wavelengths smaller than the flow depth, the secondary flow exists only near the bed and the mean flow away from the bed resembles the two-dimensional case. Staggering the crestlines alters the secondary motion; the fastest flow occurs between the lobe and the saddle planes, and two pairs of streamwise vortices appear (a strong one, centred about the lobe, and a weaker one, coming from the previous dune, centred around the saddle). The distribution of the wall stress and the focal points of separation and attachment on the bed are discussed. The sensitivity of the average reattachment length, depends on the induced secondary flow, the streamwise and spanwise components of the channel resistance (the skin friction and the form drag), and the contribution of the form drag to the total resistance are also studied. Three-dimensionality of the bed increases the drag in the channel; the form drag contributes more than in the two-dimensional case to the resistance, except for the staggered-crest case. Turbulent-kinetic energy is increased in the separated shear layer by the introduction of three-dimensionality, but its value normalized by the plane-averaged wall stress is lower than in the corresponding two-dimensional dunes. The upward flow on the stoss side and higher deceleration of flow on the lee side over the lobe plane lift and broaden the separated shear layer, respectively, affecting the turbulent kinetic energy

A dynamic subfilter-scale stress model for large eddy simulations based on physical flow scales

We propose a new definition of the length scale in an eddy-viscosity model for large-eddy simulations (LES). This formulation extends and generalizes a previous proposal [Piomelli, Rouhi and Geurts, Proc. ETMM10, 2014], in which the LES length scale was expressed in terms of the integral length-scale of turbulence determined by the flow characteristics and explicitly decoupled from the simulation grid; this approach was named Integral Length-Scale Approximation (ILSA). As in the original ILSA, the model coefficient was determined by the user, and required to maintain a desired contribution of the unresolved, subfilter scales (SFS) to the global transport. We propose a local formulation (local ILSA) in which the model coefficient is local in space, allowing a precise control over SFS activity as a function of location. This new formulation preserves the properties of the global model; application to channel flow and backward-facing step verifies its features and accuracy.</p

Large-eddy simulations of the flow on an aerofoil with leading-edge imperfections

We performed large-eddy simulations of the flow over an aerofoil to understand the effects of leading-edge roughness designed to mimic ice accretion. The roughness elements protrude outside the boundary layer, which, near the leading edge, is very thin; thus, the configuration does not represent a classical rough-wall boundary layer, but rather the flow over macroscopic obstacles. A grid convergence study is conducted and results are validated by comparison to numerical and experimental studies in the literature. The main effect of the obstacles is to accelerate transition to turbulence. Significant variations in structure generation are observed for different roughness shapes. The three-dimensionality of the irregularities has a strong impact on the flow: it creates alternating regions of high-speed (‘peaks’) and low-speed (‘valleys’) regions, a phenomenon termed ‘channelling’. The valley regions resemble a decelerating boundary layer: they exhibit considerable wake and higher levels of Reynolds stresses. The peak regions, on the other hand, are more similar to an accelerating one. Implications of the channelling phenomenon on turbulence modelling are discussed.VK acknowledges the financial support by Mitacs, Bombardier Aerospace and CARIC/CRIAQ. UP acknowledges the support from the Natural Science and Engineering Research Council of Canada (NSERC) under the Discovery Grant program, and the Canada Research Chair program. This research was enabled in part by computational support provided by Compute Ontario (computeontario.ca) and Southern Ontario Smart Computing Innovation Platform (SOSCIP) (www.soscip.org).Peer ReviewedPostprint (author's final draft