Eddy Structures in a Simulated Low Reynolds Number Turbulent Boundary Layer

Coherent structures in a tripped turbulent boundary layer havebeen analysed by applying a new conditional sampling algorithm tolarge-eddy simulation (LES) data. The space-time development of theevents and structural characteristics were examined. The new conditionalsampling scheme is shown to be very effective in the eduction of thecoherent eddy structures, allowing the simultaneous detection, trackingand averaging of several three-dimensional (3-D) structures. Alignmentof the triggering events in all spatial directions and flip-averaging(spatial reflection of the samples) enhance significantly the extractionof detailed features of the structures. The detection method minimisesthe smearing of the spatial details and avoids imposing artificialsymmetry, which is often intrinsic to other conditional samplingschemes.The results show the existence of cigar-shapedstreamwise vortices which are directly associated with negative pressurefluctuation peaks (positive source term of the pressure Poissonequation). They are inclined at 12^° to the wall andtilted laterally at an angle of 7^° to the streamwisedirection. The streamwise vortices induce ejections and sweeps throughan advection mechanism due to the tilting or inclination of thevortices. There is no evidence of hairpin vortices in either theconditional averages or instantaneous flow fields. Near-wall shearlayers are found to be related to the positive pressure fluctuationpeaks as a result of complex interactions between ejection and sweepevents. The shear layer structure has an inclination of10^°, being located between two tilted, streamwisevortices of opposite direction of rotation.The presentresults are very close to other Direct Numerical Simulation studiesusing different conditional sampling schemes. Conceptual models for thestreamwise vortices and shear layer structures are proposed to accountfor the results.     

Reynolds stress and vorticity in turbulent wall flows

Experimental measurements in a boundary layer and a large-eddy simulation of plane channel flow have been used to study the dynamics of vorticity and mass transport in the nearwall region. It was found that Reynolds stress generation occurs in the vicinity of quasi-streamwise vortices, and that smoke particles tend to be ejected from the wall near these vortical structures.     

雷诺数对三角翼绕流的影响

应用染色液流动显示技术研究了雷诺数对６０°尖前缘三角翼前缘涡破裂位置、背风面流动结构等的影响，并详细分析了背风面流动随攻角的变化．     

Reynolds Number Dependence of Velocity Structure Functions in a Turbulent Pipe Flow

Low-order moments of the increments δu andδv where u and v are the axial and radial velocity fluctuations respectively, have been obtained using single and X-hot wires mainly on the axis of a fully developed pipe flow for different values of the Taylor microscale Reynolds numberR _λ. The mean energy dissipation rate〉ε〈 was inferred from the uspectrum after the latter was corrected for the spatial resolution of the hot-wire probes. The corrected Kolmogorov-normalized second-order structure functions show a continuous evolution withR _λ. In particular, the scaling exponentζ _v , corresponding to the v structure function, continues to increase with R _λ in contrast to the nearly unchanged value of ζ _u . The Kolmogorov constant for δu shows a smaller rate of increase with R _λ than that forδv. The level of agreement with local isotropy is examined in the context of the competing influences ofR _λ and the mean shear. There is close but not perfect agreement between the present results on the pipe axis and those on the centreline of a fully developed channel flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

Statistical Structure at the Wall of the High Reynolds Number Turbulent Boundary Layer

The one and two-point statistical structure of very high Reynolds number turbulence in the surface layer near a rigid `wall' is analysed. The essential mechanisms for turbulent eddies impinging on the wall are studied using linearised rapid distortion theory, which show how the mean shear and blocking actions of the surface act first independently and then, over the life time of the eddy, interactively. Previous analytical results are reinterpreted and some new results are derived to show how the integral length scales, cross correlations and spectra of the different components of the turbulence are distorted depending on the form of the spectra of eddies above the surface layer and how they are related to motions of characteristic eddy structures near the surface. These results are applied to derive some quantitative and qualitative predictions in the surface layers (SL), where the eddies are affected by local shear dynamics, and in the `eddy surface layer' (ESL) where quasi independents loping elongated eddies interact directly with the wall, and where there is a large range of wave number within which the spectra of the horizontal velocity components are proportional to k ⁻¹. The longest eddies in the boundary layer occur near the wall. Field experiments agree with the theoretical model predictions for the quite different forms for the spectra, cospectra and cross correlations for the vertical and horizontal components of the velocity field. By showing that in SL the energy exchange between the large and small scale eddies is local(`staircase') energy cascade, whereas in ESL there is a direct nonlocal (`elevator-like')energy transfer to the small scales, it is shown why the thickness of the ESL increases over rougher surfaces and as the Reynolds number decreases. This revised version was published online in August 2006 with corrections to the Cover Date.     

Large-Eddy Simulation of the Flow Over a Circular Cylinder at Reynolds Number 3900 Using the OpenFOAM Toolbox

The flow over a circular cylinder at Reynolds number 3900 and Mach number 0.2 was predicted numerically using the technique of large-eddy simulation. The computations were carried out with an O-type curvilinear grid of size of 300 × 300 × 64. The numerical simulations were performed using a second-order finite-volume method with central-difference schemes for the approximation of convective terms. A conventional Smagorinsky and a dynamic k-equation eddy viscosity sub-grid scale models were applied. The integration time interval for data sampling was extended up to 150 vortex shedding periods for the purpose of obtaining a fully converged mean flow field. The present numerical results were found to be in good agreement with existing experimental data and previously obtained large-eddy simulation results. This gives an indication on the adequacy and accuracy of the selected large-eddy simulation technique implemented in the OpenFOAM toolbox.     

Numerical study of turbulent plane Couette flow at low Reynolds number

A direct numerical simulation of fully developed turbulent plane Couette flow has been performed. Unsteady large-scale structures, which contributed to the instantaneous energy level, were observed. These evolving and drifting vortical structures vanished after time-averaging, and the resulting mean velocity and streamwise turbulence intensity compared favourably with recent laboratory data.     

方柱绕流的数值模拟

采用有限差分法,对雷诺数为2.2×10~4的方柱绕流进行了大涡模拟(简称LES)。运用时间分裂控制(Split-Operator)法,将N-S方程分为对流步、扩散步和传播步。对Smagorinsky假设在近壁区的发散问题用两层模型进行处理。对流项用迎风—中心差分格式模拟,压力方程用SOR法迭代求解。计算得到的沿对称线的时均顺流向速度与文献上的实验结果进行了比较,结果吻合较好,同时还对绕方柱流的流场结构进行了分析研究。     

The Use of DNS to Define Stress Producing Events for Turbulent Flow over a Smooth Wall

Over the past 15 years direct numerical simulations (DNS) of turbulent flow and particle image velocimetry (PIV) have provided the opportunity to obtain information about a turbulent velocity field simultaneously at a large number of locations. This paper gives a personal viewpoint of how these techniques are providing new insights about the Reynolds stress producing structures in turbulence generated by flow over a smooth boundary. This revised version was published online in August 2006 with corrections to the Cover Date.     

Simulation and Modelling of Turbulent Trailing-Edge Flow

Computations of turbulent trailing-edge flow have been carried out at a Reynolds number of 1000 (based on the free-stream quantities and the trailing-edge thickness) using an unsteady 3D Reynolds-Averaged Navier–Stokes (URANS) code, in which two-equation (k–ε) turbulence models with various low-Re near wall treatments were implemented. Results from a direct numerical simulation (DNS) of the same flow are available for comparison and assessment of the turbulence models used in the URANS code. Two-dimensional URANS calculations are carried out with turbulence mean properties from the DNS used at the inlet; the inflow boundary-layer thickness is 6.42 times the trailing-edge thickness, close to typical turbine blade flow applications. Many of the key flow features observed in DNS are also predicted by the modelling; the flow oscillates in a similar way to that found in bluff-body flow with a von Kármán vortex street produced downstream. The recirculation bubble predicted by unsteady RANS has a similar shape to DNS, but with a length only half that of the DNS. It is found that the unsteadiness plays an important role in the near wake, comparable to the modelled turbulence, but that far downstream the modelled turbulence dominates. A spectral analysis applied to the force coefficient in the wall normal direction shows that a Strouhal number based on the trailing-edge thickness is 0.23, approximately twice that observed in DNS. To assess the modelling approximations, an a priori analysis has been applied using DNS data for the key individual terms in the turbulence model equations. A possible refinement to account for pressure transport is discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Effects of the Similarity Model in Finite-Difference LES of Isotropic Turbulence Using a Lagrangian Dynamic Mixed Model

A Lagrangian dynamic formulation of the mixed similarity subgrid (SGS) model for large-eddy simulation (LES) of turbulence is proposed. In this model, averaging is performed over fluid trajectories, which makes the model applicable to complex flows without directions of statistical homogeneity. An alternative version based on a Taylor series expansion (nonlinear mixed model) is also examined. The Lagrangian models are implemented in a finite difference code and tested in forced and decaying isotropic turbulence. As comparison, the dynamic Smagorinsky model and volume-averaged formulations of the mixed models are also tested. Good results are obtained, except in the case of low-resolution LES (32³) of decaying turbulence, where the similarity coefficient becomes negative due to the fact that the test-filter scale exceeds the integral scale of turbulence. At a higher resolution (64³), the dynamic similarity coefficient is positive and good agreement is found between predicted and measured kinetic energy evolution. Compared to the eddy viscosity term, the similarity or the nonlinear terms contribute significantly to both SGS dissipation of kinetic energy and SGS force. In order to dynamically test the accuracy of the modeling, the error incurred in satisfying the Germano identity is evaluated. It is found that the dynamic Smagorinksy model generates a very large error, only 3% lower than the worst-case scenario without model. Addition of the similarity or nonlinear term decreases the error by up to about 50%, confirming that it represents a more realistic parameterization than the Smagorinsky model alone.     

Large-Eddy Simulation of a Turbulent Hydrogen Diffusion Flame

In this work a large-eddy simulation (LES) of a turbulent hydrogen jet diffusion flame is presented. The numerical method handles fluctuations of density in space and in time, but assumes density to be independent of pressure (incompressibility). The chemical composition of the fluid is described by solving the filtered transport equation for mixture fraction f. Density, viscosity and temperature are evaluated assuming chemical equilibrium. To account for sub-grid fluctuations of f, its sub-grid distribution is presumed to have the shape of a β-function. The results of the simulation are discussed extensively. The influence of inlet boundary conditions is addressed and radial profiles at different axial positions are shown for a complete set of one-point statistical data. Agreement of numerical results and experimental data is very good. Furthermore, a comparison of Reynolds- and Favre-averages is done and energy spectra at different locations in the flame are discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

高雷诺数流动的加罚有限元数值分析

在传统的伽辽金变分为有限元数值逼近思想的基础上，本文采用改进的加罚有限元法对粘性不可压缩紊流动进行数值计算。在高雷诺数流动时，为避免对流效应过强产生的数值计算的振荡，对标准权函数引入迎风修正项，同时采用雷诺数加载法，保证数值解的收敛性。在有限元方法离散过程中，采用有效的隐式压力－显式速度方式，以准确的速度场确保获得压力的稳定性。速度压力项选用不同阶次的插值函数。实践证明：当控制方法各对流项假扩散被     

From Two-Point Closures of Isotropic Turbulence to LES of Shear Flows

We first recall the EDQNM two-point closure approach of three-dimensional isotropic turbulence. It allows in particular prediction of the infrared kinetic-energy dynamics (with ak ₄ backscatter) and the associated time-decay law of kinetic-energy, useful in particular for one-point closure modelling. Afterwards, we show how the spectral eddy viscosity concept may be used for large-eddy simulations: we introduce the plateau-peak model and the spectral-dynamic models. They are applied to decaying isotropic turbulence, and allow recovery of the EDQNM infrared energy dynamics. Anew infrared k ₂ law for the pressure spectrum, predicted by the closure, is also well verified. Assuming that subgrid scales are not too far from isotropy, the spectral-dynamic model is applied to the channel flow at h ₊= 390, with statistics in very good agreement with DNS, while reducing considerably the computational time. We study with the aid of DNS and LES the case of the channel rotating about an axis of spanwise direction. The calculations allow to recover the universal linear behaviour of the mean velocity profile, with a local Rossby number equal to −1. We present also LES (using the Grenoble Filtered Structure-Function Model), of a turbulent boundary layer passing over a cavity. Finally, we make some remarks on the future of LES for industrial applications. This revised version was published online in July 2006 with corrections to the Cover Date.     

Three‐dimensional numerical simulation of compound meandering open channel flow by the Reynolds stress model

Turbulent flow in a compound meandering open channel with seminatural cross sections is one of the most complicated turbulent flows as the flow pattern is influenced by the combined action of various forces, such as centrifugal force, pressure, and shear stresses. In this paper, a three‐dimensional (3D) Reynolds stress model (RSM) is adopted to simulate the compound meandering channel flows. Governing equations of the flow are solved numerically with finite‐volume method. The velocity fields, wall shear stresses, and Reynolds stresses are calculated for a range of input conditions. Good agreement between the simulated results and measurements indicates that RSM can successfully predict the complicated flow phenomenon. Copyright © 2008 John Wiley & Sons, Ltd.     

Investigation of Anisotropy-Resolving Turbulence Models by Reference to Highly-Resolved LES Data for Separated Flow

The predictive properties of several non-linear eddy-viscosity models are investigated by reference to highly-resolved LES data obtained by the authors for an internal flow featuring massive separation from a curved surface. The test geometry is a periodic segment of a channel constricted by two-dimensional (2D) `hills' on the lower wall. The mean-flow Reynolds number is 21560. Periodic boundary conditions are applied in the streamwise and spanwise directions. This makes the statistical properties of the simulated flow genuinely 2D and independent from boundary conditions, except at the walls. The simulation was performed on a high-quality, 5M-node grid. The focus of the study is on the exploitation of the LES data for the mean-flow, Reynolds stresses and macro-length-scale. Model solutions are first compared with the LES data, and selected models are then subjected to a-priori studies designed to elucidate the role of specific model fragments in the non-linear stress-strain/vorticity relation and their contribution to observed defects in the mean-flow and turbulence fields. The role of the equation governing the length-scale, via different surrogate variables, is also investigated. It is shown that, while most non-linear models overestimate the separation region, due mainly to model defects that result in insufficient shear stress in the separated shear layer, model forms can be derived which provide a satisfactory representation of the flow. One such model is identified. This combines a particular quadratic constitutive relation with a wall-anisotropy term, a high-normal-strain correction and a new form of the equation for the specific dissipation ω = ∈/k. This revised version was published online in July 2006 with corrections to the Cover Date.     

Direct Numerical Simulation of Transitional Turbulent Flow in a Closed Rotor-Stator Cavity

The transitional turbulent regime in confined flow between a rotating and a stationary disc is studied using direct numerical simulation. Besides its fundamental importance as a three-dimensional prototype flow, such flows frequently arise in many industrial devices, especially in turbomachinary applications. The present contribution extends the DNS simulation into the turbulent flow regime, to a rotational Reynolds number Re =3 × 10⁵. An annular rotor-stator cavity of radial extension ΔR and height H, is considered with L = 4.72(L = ΔR/H) and Rm = 2.33 (Rm = (R ₁+ R ₀)/ΔR). The direct numerical simulation is performed by integrating the time-dependent Navier–Stokes equations until a statistically steady state is reached. A three-dimensional spectral method is used with the aim of providing both very accurate instantaneous fields and reliable statistical data. The instantaneous quantities are analysed in order to enhance our knowledge of the physics of turbulent rotating flows. Also, the results have been averaged so as to provide target turbulence data for any subsequent modelling attempts at reproducing the flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

静电旋风分离器气相流场的数值模拟

对一内部安装电晕极的切向进口旋风分离器,以三维贴体坐标为基础,应用Bradshaw的修正k-ε湍流模型,用非交错的SIMPLE算法对静电旋风分离器的气相流场进行求解,计算结果与实验结果进行了对比。分析了电晕极的不同安装位置对旋风分离器流场的影响。从流场的角度来看,电晕极安装在筒体与排气管之间并靠近排气管的位置有利于提高静电旋风分离器的分离效率。