Extension of the Kida law in turbulenceExtension de la loi de Kida en turbulence

We extend the validity range of Kida's log-stable law of stability index α=1.65 and intermittency parameter μ=0.2 to a new range of Reynolds number. This law describes intermittencies in fully developed turbulent flows or more precisely the p.d.f. of turbulence dissipation. Former measurements of the hyper-flatness factors of order 4, 5, 6 of turbulent velocity increments, coming from both experimental works and numerical simulations are used. We show that the power-law variation of these hyper-flatness factors with Taylor scale based Reynolds numbers Re_λ can be fitted, for Re_λ ranging from 35 to 750, by a log-stable law of stability index α=1.65 and intermittency parameter μ=0.21. To cite this article: N. Rimbert, O. Séro-Guillaume, C. R. Mecanique 331 (2003).     

Swirling-strength based large eddy simulation of turbulent flow around single square cylinder at low Reynolds numbers

In view of the fact that large scale vortices play the substantial role of momentum transport in turbulent flows, large eddy simulation(LES) is considered as a better simulation model. However, the sub-grid scale(SGS) models reported so far have not ascertained under what flow conditions the LES can lapse into the direct numerical simulation. To overcome this discrepancy, this paper develops a swirling strength based the SGS model to properly model the turbulence intermittency, with the primary characteristics that when the local swirling strength is zero, the local sub-grid viscosity will be vanished. In this paper, the model is used to investigate the flow characteristics of zero-incident incompressible turbulent flows around a single square cylinder(SC)at a low Reynolds number range Re ∈ [103, 104]. The flow characteristics investigated include the Reynolds number dependence of lift and drag coefficients, the distributions of time-spanwise averaged variables such as the sub-grid viscosity and the logarithm of Kolmogorov micro-scale to the base of 10 at Re = 2 500 and 104, the contours of spanwise and streamwise vorticity components at t = 170. It is revealed that the peak value of sub-grid viscosity ratio and its root mean square(RMS) values grow with the Reynolds number. The dissipation rate of turbulent kinetic energy is larger near the SC solid walls.The instantaneous factor of swirling strength intermittency(FSI) exhibits some laminated structure involved with vortex shedding.     

Direct numerical simulation of turbulent flows through concentric annulus with circumferential oscillation of inner wall

Periodic wall oscillations in the spanwise or circumferential direction can greatly reduce the friction drag in turbulent channel and pipe flows. In a concentric annulus, the constant rotation of the inner cylinder can intensify turbulence fluctuations and enhance skin friction due to centrifugal instabilities. In the present study, the effects of the periodic oscillation of the inner wall on turbulent flows through concentric annulus are investigated by the direct numerical simulation (DNS). The radius ratio of the inner to the outer cylinders is 0.1, and the Reynolds number is 2 225 based on the bulk mean velocity U_m and the half annulus gap H. The influence of oscillation period is considered. It is found that for short-period oscillations, the Stokes layer formed by the circumferential wall movement can effectively inhibit the near-wall coherent motions and lead to skin friction reduction, while for long-period oscillations, the centrifugal instability has enough time to develop and generate new vortices, resulting in the enhancement of turbulence intensity and skin friction.     

Turbulence in a three‐dimensional wall‐bounded shear flow

A new turbulent flow with distinct three‐dimensional characteristics has been designed in order to study the impact of mean‐flow skewing on the turbulent coherent vortices and Reynolds‐averaged statistics. The skewing of a unidirectional plane Couette flow was achieved by means of a spanwise pressure gradient. Direct numerical simulations of the statistically steady Couette–Poiseuille flow enabled in‐depth explorations of the turbulence field in the skewed flow. The imposition of a modest spanwise gradient turned the mean flow about 8° away from the original Couette flow direction and this turning angle remained nearly the same over the entire cross section. Nevertheless, a substantial non‐alignment between the turbulent shear stress angle and the mean velocity gradient angle was observed. The structure parameter turned out to slightly exceed that in the pure Couette flow, contrary to the observations made in some other three‐dimensional shear flows. Coherent flow structures, which are known to be associated with the Reynolds shear stress in near‐wall regions, were identified by the λ₂‐criterion. Instantaneous and ensemble‐averaged vortices resembled those found in the unidirectional Couette flow. In the skewed flow, however, the vortex structures were turned to align with the local mean‐flow direction. The conventional symmetry between Case 1 and Case 2 vortices was broken due to the mean‐flow three‐dimensionality. The turning of the coherent vortices and the accompanying symmetry‐breaking gave rise to secondary and tertiary turbulent shear stress components. By averaging the already ensemble‐averaged shear stresses associated with Case 1 and Case 2 vortices in the homogeneous directions, a direct link between the educed near‐wall structures and the Reynolds‐averaged turbulent stresses was established. These observations provide evidence in support of the hypothesis that the structural model proposed for two‐dimensional turbulent boundary layers remains valid also in flows with moderate mean three‐dimensionality. Copyright © 2009 John Wiley & Sons, Ltd.     

The Relaminarization Mechanisms of Turbulent Channel Flow at Low Reynolds Numbers

The mechanisms of laminarization in wall-bounded flows have been investigated by performing direct numerical simulations (DNS) of turbulent channel flows. By decreasing Reynolds numbers systematically, the effects of the low Reynolds number are studied in connection with the near-wall turbulent structure and turbulent statistics. At approximately the critical Reynolds number, the turbulent skin friction is reduced, and the turbulent structure changes qualitatively in the very near-wall region. Instantaneous turbulent structures reveal that streamwise vortices, the cores of which are at y+ 10, disappear, although low speed streaks and Reynolds shear stress are still produced by larger streamwise vortices located in the buffer region y+ > 10. Sweep motions induced by these vortical structures are shifted toward the center of a channel and also significantly deterred, which may heighten the effects of the viscous sublayer over most of the channel section and suppress the regeneration mechanisms of new streamwise vortices in the very near-wall region. To investigate the details of how large-scale coherent vortices affect the viscous sublayer and the relevant small-scale streamwise vortices, a body force is virtually imposed in the wall-normal direction to enhance the large streamwise vortices. As a result, it is found that when they are sufficiently enhanced, the small-scale vortices reappear, and the sweep events are again dominant in the viscous sublayer.     

DNS of turbulent heat transfer in a channel flow with a high spatial resolution

Direct numerical simulations of turbulent heat transfer in a channel flow are performed to investigate the effects of Reynolds and Prandtl numbers on higher-order turbulence statistics such as a turbulent Prandtl number and the budget for the dissipation rate of the temperature variance. The Reynolds numbers based on the friction velocity and the channel half width are 180 and 395, and the molecular Prandtl numbers Pr’s 0.71–10.0. Careful attention is paid to ensure accuracy of the higher-order statistics through the use of a high spatial resolution comparable to Batchelor length scale. The wall-asymptotic value of the turbulent Prandtl number is mostly independent of Reynolds number for the current range of Pr’s. The budget for the dissipation rate of the temperature variance has been computed, and the negligible effect of a Reynolds number on the sum of all source and sink terms in near-wall region in the current computational range is found. This result is quite similar to the one in the budget for the dissipation rate of turbulent energy. In addition, a priori test for existing models is also performed to assess the Pr dependence on the individual terms and their summations in the budget.     

Large eddy simulation of temporally developing juncture flows

Large eddy simulation (LES) results are reported for temporally developing solid–solid and solid–rigid-lid juncture flows. A MacCormack-type scheme that is second-order in time, and fourth-order in space for the convective terms and second-order in space for the viscous terms, is used. The simulations are obtained for a low subsonic Mach number. The subgrid-scale stresses (SGS) are modeled using the dynamic modeling procedure. The turbulent flow field generated on a flat-plate boundary layer is used to initialize the juncture flow simulations. The results of the flat-plate boundary layer simulations are validated with experimental and direct numerical simulations (DNS) data. In juncture flow simulations, the presence of an adjacent solid-wall/rigid-lid boundary altered the mean and the turbulent field, setting up gradients in the anisotropy of normal Reynolds stresses resulting in the formation of turbulence-induced secondary vortices. The relative size of these secondary vortices and the distribution of mean and turbulent quantities are in qualitative agreement with the experimental observations for the solid–solid juncture. The overall distribution of the mean and turbulence quantities showed close resemblance between the solid–solid and the solid–rigid-lid junctures; except for the absence of a second vortical region near the rigid-lid boundary. In agreement with the experimental observations, it was found that the normalized anisotropy term exhibited similarity when plotted against the distance from the boundary, regardless of the type of boundary, i.e. solid-wall or rigid-lid. The turbulent kinetic energy increased near the rigid-lid boundary. While the surface normal velocity fluctuations decreased to zero at the rigid-lid boundary, the other two velocity components showed an increase in their energy, which is also consistent with the experimental observations. © 1998 John Wiley & Sons, Ltd.     

Direct Numerical Simulation of Self-Similar Turbulent Boundary Layers in Adverse Pressure Gradients

Direct numerical simulations of the Navier–Stokes equations have been carried out with the objective of studying turbulent boundary layers in adverse pressure gradients. The boundary layer flows concerned are of the equilibrium type which makes the analysis simpler and the results can be compared with earlier experiments and simulations. This type of turbulent boundary layers also permits an analysis of the equation of motion to predict separation. The linear analysis based on the assumption of asymptotically high Reynolds number gives results that are not applicable to finite Reynolds number flows. A different non-linear approach is presented to obtain a useful relation between the freestream variation and other mean flow parameters. Comparison of turbulent statistics from the zero pressure gradient case and two adverse pressure gradient cases shows the development of an outer peak in the turbulent energy in agreement with experiment. The turbulent flows have also been investigated using a differential Reynolds stress model. Profiles for velocity and turbulence quantities obtained from the direct numerical simulations were used as initial data. The initial transients in the model predictions vanished rapidly. The model predictions are compared with the direct simulations and low Reynolds number effects are investigated.     

Hairpin vortices in the largest scale of turbulent boundary layers

We have conducted direct numerical simulations of a turbulent boundary layer for the momentum-thickness-based Reynolds number ${\mathit{Re}}_{\theta }$ = 180–4600. To extract the largest-scale vortices, we coarse-grain the fluctuating velocity fields by using a Gaussian filter with the filter width comparable to the boundary layer thickness. Most of the largest-scale vortices identified by isosurfaces of the second invariant of the coarse-grained velocity gradient tensor are similar to coherent vortices observed in low-Reynolds-number regions, that is, hairpin vortices or quasi-streamwise vortices inclined to the wall. We also develop a percolation analysis to investigate the threshold-dependence of the isosurfaces and objectively identify the largest-scale hairpin vortices in terms of the coarse-grained vorticity, which leads to the quantitative evidence that they never disappear even in fully developed turbulent regions. Hence, we conclude that hairpin vortices exist in the largest-scale structures irrespective of the Reynolds number.     

Numerical prediction of turbulent thermocapillary convection in superposed fluid layers with a free interface

The present paper introduces a new numerical method for predicting the characteristics of thermocapillary turbulent convection in a differentially-heated rectangular cavity with two superposed and immiscible fluid layers. The unsteady Reynolds form of the Navier–Stokes equations and energy equation are solved by using the control volume approach on a staggered grid system using SIMPLE algorithm. The turbulence quantities are predicted by applying the standard k–ε turbulence model. The level set formulation is applied for predicting the topological changes of the interface separating the two fluid layers and to provide an accurate and robust modeling of the interfacial normal and tangential stresses. The computational results obtained showed good agreement when compared with the previous experimental, numerical and analytical benchmark data for different validation cases in both laminar and turbulent regimes. The present numerical method is then applied to predict the velocity and temperature distribution in two immiscible liquid layers with undeformable interface for a wide range of Marangoni numbers. The laminar-turbulent transition is demonstrated by obtaining the turbulence features at high interfacial temperature gradient which is characterized by high Marangoni number. The effect of increasing Marangoni number on the interface dynamics in turbulent regime is also investigated.     

A Bypass Transition Model Based on the Intermittency Function

An intermittency model that is formulated in local variables is proposed for representing bypass transition in Reynolds-Averaged Navier-Stokes (RANS) computations. No external data correlation is used to fix transition. Transition is initiated by diffusion, and a source term carries it to completion. A sink term is created to predict the laminar region before transition, then it vanishes in the turbulent region. Both the source and sink are functions of a wall-distance Reynolds number and turbulence scale. A modification is introduced to predict transition in separated boundary layers. The transition model is incorporated with the k?ω RANS model. The present model is implemented into a general purpose, computational fluid dynamics (CFD) code. The model is validated with several test cases. Decent agreement with the available data is observed in a range of flows.     

壁湍流相干结构和减阻控制机理

剪切湍流中相干结构的发现是上世纪湍流研究的重大进展之一,这些大尺度的相干运动在湍流的动力学过程中起重要作用,也为湍流的控制指出了新的方向.壁湍流高摩擦阻力的产生与近壁区流动结构密切相关,基于近壁区湍流动力学过程的减阻控制方案可以有效降低湍流的摩擦阻力,但是随着雷诺数的升高, 这些控制方案的有效性逐渐降低.近年来研究发现, 在高雷诺数情况下外区存在大尺度的相干运动,这种大尺度运动对近壁区湍流和壁面摩擦阻力的产生有重要影响,为高雷诺数湍流减阻控制策略的设计提出了新的挑战.该文将对壁湍流相干结构的研究历史加以简单的回顾,重点介绍近壁区相干结构及其控制机理、近年来高雷诺数外区大尺度运动的研究进展,在此基础上提出高雷诺数减阻控制研究的关键科学问题.      

A modified STRUCT model for efficient engineering computations of turbulent flows in hydro-energy machinery

A modified STRUCT (MST) turbulence model for efficient engineering computations of turbulent flows in hydro-energy machinery is proposed in this paper. The MST model switches between URANS and LES-like modes using a new damping function to adjust the turbulent viscosity. Compared with the original STRUCT method, the modifications are as follows: (1) the BSL k-ω model with the Spalart-Shur correction is chosen as the new baseline to improve the sensitivity to rotation and curvature; (2) a new adaptive time-scale ratio is proposed to avoid the arbitrariness of geometric averaging operation in the original method; (3) the normalized helicity is introduced into the new damping function to detect the energy backscatter phenomenon. Five classical high Reynolds number flow cases are tested. The results show that the turbulent viscosity of the MST model is reasonably reduced in the massively separated regions and LES-like mode is activated, which captures more turbulent vortices and fluctuations on the URANS grids. With high efficiency and robustness, the MST model inherits the advantages of the original STRUCT method and improves the prediction accuracy of the turbulence with rotation and curvature, which enables efficient engineering computations of turbulent flows in hydro-energy machinery.     

Compound wall treatment with low‐Re turbulence model

A generalized treatment for the wall boundary conditions relating to turbulent flows is developed that blends the integration to a solid wall with wall functions. The blending function ensures a smooth transition between the viscous and turbulent regions. An improved low Reynolds number k?ε model is coupled with the proposed compound wall treatment to determine the turbulence field. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. Computations with fine and coarse meshes of a few flow cases yield appreciably good agreement with the direct numerical simulation and experimental data. The method is recommended for computing the complex flows where computational grids cannot satisfy a priori the prerequisites of viscous/turbulence regions. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation of an axisymmetric separated and reattached flow over a longitudinal blunt circular cylinder

This article presents a numerical investigation of turbulent flow in an axisymmetric separated and reattached flow over a longitudinal blunt circular cylinder. The governing equations were discretized by the finite-volume method and SIMPLER method was applied to solve the equations on a staggered grid. The turbulent flow was numerically simulated using the standard k–ε, Abe–Kondoh–Nagano (AKN) and Shear Stress Transport (SST) turbulence models. The comparisons made between numerical results and experimental measurements showed that the SST model is superior to other models in the present calculation.Computations were performed for three different Reynolds numbers of 6000, 10 000 and 20 000 based on the cylinder diameter. To our knowledge, this study represents the first numerical investigation of the present flow configuration. The computational results were validated with the available experimental data of reattachment length, mean velocity distribution and wall static pressure coefficient in the turbulent blunt circular cylinder flows. Further, other characteristics of the flow, such as turbulent kinetic energy, pressure, streamlines, and the velocity vectors are discussed.The results show that the main characteristics of the turbulence flow in the separation region, such as reattachment length or velocity profiles, are nearly independent of the Reynolds number. The obtained results showed that a secondary separation bubble may appear in the main separation bubble near the leading edge. Furthermore, it was found that the turbulent kinetic energy has a large effect on the formation of the secondary bubble.     

LES-based filter-matrix lattice Boltzmann model for simulating fully developed turbulent channel flow

In this paper, a three-dimensional filter-matrix lattice Boltzmann (FMLB) model based on large eddy simulation (LES) was verified for simulating wall-bounded turbulent flows. The Vreman subgrid-scale model was employed in the present FMLB–LES framework, which had been proved to be capable of predicting turbulent near-wall region accurately. The fully developed turbulent channel flows were performed at a friction Reynolds number Re_τ of 180. The turbulence statistics computed from the present FMLB–LES simulations, including mean stream velocity profile, Reynolds stress profile and root-mean-square velocity fluctuations greed well with the LES results of multiple-relaxation-time (MRT) LB model, and some discrepancies in comparison with those direct numerical simulation (DNS) data of Kim et al. was also observed due to the relatively low grid resolution. Moreover, to investigate the influence of grid resolution on the present LES simulation, a DNS simulation on a finer gird was also implemented by present FMLB–D3Q19 model. Comparisons of detailed computed various turbulence statistics with available benchmark data of DNS showed quite well agreement.     

DNS of viscoelastic turbulent channel flow with rectangular orifice at low Reynolds number

Direct numerical simulations of turbulent viscoelastic-fluid flow in a channel with a rectangular orifice were performed to investigate the influence of viscoelasticity on turbulence statistics and turbulent structures downstream of the orifice. The geometry considered is periodic rectangular orifices with 1:2 expansion. The constitutive equation follows the Giesekus model, valid for polymer (or surfactant) solutions, which are generally capable of reducing the turbulent frictional drag in a smooth channel. The friction Reynolds number and the Weissenberg number were set to 100 and 20-30, respectively. A drag reduction of about 20% was achieved in the viscoelastic flows. The onset Reynolds number for the transition from a symmetric to an asymmetric state was found to be shifted to higher values than that for the Newtonian flow. In the viscoelastic flow, the turbulent kinetic energy was decreased and fewer turbulent eddies were observed, as the Kelvin-Helmholtz vortices were quickly damped. Away from the orifice, quasi-streamwise vortices in the viscoelastic flow were sustained for a longer period, accompanied by energy exchange from elastic energy of the viscoelastic fluid to kinetic energy.     

Coherent Structures and their Frequency Signature in the Separated Shear Layer on the Sides of a Square Cylinder

The purpose of the present work is to study the various specific time scales of the turbulent separating flow around a square cylinder, in order to determine the Reynolds number effect on the separating shear layer, where occurs a transition to turbulence. Unsteady analysis based on large eddy simulation (LES) at intermediate Reynolds numbers and laser doppler velocimetry (LDV) measurements at high Reynolds numbers are carried out. The Reynolds number, based on the cylinder diameter D and the inflow velocity U _o , is ranging from Re?=?50 to Re?=?300,000. A special focus is performed on the coherent structures developing on the sides and in the wake of a square cylinder. For a large Reynolds number range above Re?≈?1,000, both signatures of Von Karman (VK) and Kelvin–Helmholtz (KH) type vortical structures are found on velocity time samples. The combination of their frequency signature is studied based on Fourier and wavelet analysis. In the present study, We observe the occurrence of KH pairings in the separating shear layer on the side of the cylinder, and confirm the intermittency nature of such a shear flow. These issues concerning the structure of the near wake shear layer which were addressed for the round cylinder case in a recent experimental publication (Rajagopalan and Antonia, Exp Fluids 38:393–402, 2005) are of interest in the present flow configuration as well.     

Direct numerical simulation of turbulent flow in a spanwise rotating square duct at high rotation numbers

In this paper, direct numerical simulations have been performed to study the effects of Coriolis force on the turbulent flow field confined within a square duct subjected to spanwise system rotations at high rotation numbers. In response to the system rotation, secondary flows appear as large streamwise counter-rotating vortices, which interact intensely with the four boundary layers and have a significant impact on flow statistics, velocity spectra and coherent structures. It is observed that at sufficiently high rotation numbers, a Taylor–Proudman region appears and complete laminarization is almost reached near the top and side walls. The influence of large organized secondary flows on the production rate and re-distribution of turbulent kinetic energy has been investigated through a spectral analysis. It is observed that the Coriolis force dominates the transport of Reynolds stresses and turbulent kinetic energy, and forces the spectra of streamwise and vertical velocities to synchronize within a wide range of scales.