ALBERT

Description: 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. © 2013 Cambridge University Press.

Description: The physics of the roughness sublayer are studied by direct numerical simulations (DNS) of an open-channel flow with sandgrain roughness. A double-averaging (DA) approach is used to separate the spatial variations of the time-averaged quantities and the turbulent fluctuations. The spatial inhomogeneity of velocity and Reynolds stresses results in an additional production term for the turbulent kinetic energy (TKE) - the 'wake production'; it is the excess wake kinetic energy (WKE), generated from the work of mean flow against the form drag, that is not directly dissipated into heat, but instead converted into turbulence. The wake production promotes wall-normal turbulent fluctuations and increases the pressure work, which ultimately leads to more homogeneous turbulence in the roughness sublayer, and to the increase of Reynolds shear stress and the drag on the rough wall. In the fully rough regime, roughness directly affects the generation of the wall-normal fluctuations, while in the transitionally rough regime, the region affected by roughness is separated from the region of intense generation of these fluctuations. The budget of the WKE and the connection between the wake and the turbulence suggest strong interactions between the roughness sublayer and the outer layer that are insensitive to the variation of the outer-layer conditions. Furthermore, the present results may have implications for the relationship between the roughness geometry and the flow dynamics in the region directly affected by roughness. © 2014 Cambridge University Press.

Description: We propose a new length scale as a basis for the modelling of subfilter motions in large-eddy simulations (LES) of turbulent flow. Rather than associating the model length scale with the computational grid, we put forward an approximation of the integral length scale to achieve a non-uniform flow coarsening through spatial filtering that reflects the local, instantaneous turbulence activity. Through the introduction of this grid-independent, solution-specific length scale it becomes possible to separate the problem of representing small-scale turbulent motions in a coarsened flow model from that of achieving an accurate numerical resolution of the primary flow scales. The formulation supports the notion of grid-independent LES, in which a prespecified reliability measure is used. We investigate a length-scale definition based on the resolved turbulent kinetic energy (TKE) and its dissipation. The proposed approach, which we call integral length-scale approximation (ILSA) model, is illustrated for turbulent channel flow at high Reynolds numbers and for homogeneous isotropic turbulence (HIT). We employ computational optimization of the model parameter based on various measures of subfilter activity, using the successive inverse polynomial interpolation (SIPI) and establish the efficiency of this route to subfilter modelling. © 2015 Cambridge University Press.

Description: Large-eddy simulations of temporally evolving turbulent mixing layers have been carried out. The effect of the initial conditions and the size of the computational box on the turbulent statistics and structures is examined in detail. A series of calculations was initialized using two different realizations of a spatially developing turbulent boundary-layer with their free streams moving in opposite directions. Computations initialized with mean flow plus random perturbations with prescribed moments were also conducted. In all cases, the initial transitional stage, from boundary-layer turbulence or random noise to mixing-layer turbulence, was followed by a self-similar period. The self-similar periods, however, differed considerably: the growth rates and turbulence intensities showed differences, and were affected both by the initial condition and by the computational domain size. In all simulations the presence of quasi-two-dimensional spanwise rollers was clear, together with ‘braid’ regions with quasi-streamwise vortices. The development of these structures, however, was different: if strong rollers were formed early (as in the cases initialized by random noise), a well-organized pattern persisted throughout the self-similar period. The presence of boundary layer turbulence, on the other hand, inhibited the growth of the inviscid instability, and delayed the formation of the roller–braid patterns. Increasing the domain size tended to make the flow more three-dimensional.