Effect of an adverse pressure gradient on the streamwise Reynolds stress profile maxima in a turbulent boundary layer

It is widely accepted that in a turbulent boundary layer （TBL） with adverse pressure gradient （APG） an outer peak usually appears in the profile of streamwise Reynolds stress. However, the effect of APG on this outer peak is not clearly understood. In this paper, the effect of APG is analysed using the numerical and experimental results in the literature. Because the effect of upstream flow is inherent in the TBL, we first analyse this effect in TBLs with zero pressure gradient on flat plates. Under the individual effect of upstream flow, an outer peak already appears in the profile of streamwise Reynolds stress when the TBL continues developing in the streamwise direction. The APG accelerates the appearance of the outer peak, instead of being a trigger.     

A low Reynolds number variant of partially-averaged Navier-Stokes model for turbulence

A low Reynolds number (LRN) formulation based on the Partially Averaged Navier-Stokes (PANS) modelling method is presented, which incorporates improved asymptotic representation in near-wall turbulence modelling. The effect of near-wall viscous damping can thus be better accounted for in simulations of wall-bounded turbulent flows. The proposed LRN PANS model uses an LRN k-ε model as the base model and introduces directly its model functions into the PANS formulation. As a result, the inappropriate wall-limiting behavior inherent in the original PANS model is corrected. An interesting feature of the PANS model is that the turbulent Prandtl numbers in the k and ε equations are modified compared to the base model. It is found that this modification has a significant effect on the modelled turbulence. The proposed LRN PANS model is scrutinized in computations of decaying grid turbulence, turbulent channel flow and periodic hill flow, of which the latter has been computed at two different Reynolds numbers of Re = 10,600 and 37,000. In comparison with available DNS, LES or experimental data, the LRN PANS model produces improved predictions over the standard PANS model, particularly in the near-wall region and for resolved turbulence statistics. Furthermore, the LRN PANS model gives similar or better results - at a reduced CPU time - as compared to the Dynamic Smagorinsky model.     

Enhancement of finite Reynolds number effects: Inner-outer sublayer interaction in the turbulent boundary layer

Residual Reynolds number effects in the established data for the velocity profile in turbulent boundary layers (and in pipe or channel flows) are found to be remarkably large. We combine two eddy-viscosity models (with overlapping validity in the inertial sublayer) and show (both analytically and numerically) that this enhancement (which involves a viscous correlation length) arises from inner-outer sublayer interaction.     

Assessment of two‐equation turbulence modelling for high Reynolds number hydrofoil flows

This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two‐equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k–εmodel, the realizable k–εmodel, the standard k–ωmodel and the shear‐stress‐transport (SST) k–ωmodel. It has found that the realizable k–εturbulence model used with enhanced wall functions and near‐wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. Copyright © 2004 John Wiley & Sons, Ltd.     

The effect of the Reynolds number on the Reynolds stresses and vorticity in a turbulent far-wake

Using two orthogonal arrays of 16 X-wires, eight in the (x,y)-plane and eight in the (x,z)-plane, the effect of the Reynolds number in a turbulent plane far-wake has been investigated for two values of Re_θ (based on the free stream velocity and the momentum thickness), i.e. 1350 and 4600. It is observed that as the Reynolds number increases the magnitudes of the measured Reynolds stresses increase, as does the size of two-point vorticity correlation iso-contours. Discernible differences are also observed in probability density function, spectra and three-dimensional topologies. The Reynolds number dependence seems to vanish when Re_θ5000.     

低雷诺数下附面层组合抽吸方案对压气机特性影响的研究

为了分析附面层抽吸流动控制对低雷诺数下压气机特性的影响,本文采用数值方法模拟了低雷诺数下附面层组合抽吸方案对NASA Rotor 37跨音速压气机性能和稳定性的影响特点及作用机理。通过在该压气机转子叶片吸力面和机匣上分别设计附面层抽吸槽,探讨了组合抽吸方案对低雷诺数下(H=20km)压气机性能和稳定性的影响。结果表明:采用组合抽吸方案后,压气机峰值效率提高约1.3%;压气机最大增压比提高约2.5%;压气机转子的近失速点流量减小约14.6%。进一步分析作用机理发现,组合抽吸槽有效抑制了附面层径向涡向叶顶的运动和聚集,使叶顶附面层分离区减少约70%从而有效改善了压气机的流场特性。     

The turbulent/non-turbulent interface in an adverse pressure gradient turbulent boundary layer

The turbulent/non-turbulent interface (TNTI) in an adverse pressure gradient (APG, β = 1.45) turbulent boundary layer (TBL) is explored here by using direct numerical simulation (DNS) data; β is the Clauser pressure gradient parameter. For comparison, the DNS data for a zero pressure gradient (ZPG) TBL is included. The interface is extracted with an approach based on enstrophy criteria. Depending on the enstrophy, the outer boundary layer flow can be classified into the free stream, boundary layer wake, and intermittent flow regimes. The fractal dimension of the interface is obtained by using the box-counting algorithm, and was found to be constant over a long range of box sizes. The TNTI shows a monofractal behavior. The geometric complexity of a TNTI can be determined in terms of the genus, which is defined as the number of handles in a geometric object. We examine the volume and projection area of the genus of the TNTI to analyze the entrainment process. The geometric complexity of the APG TBL interface and the local entrainment are greater than those of the ZPG TBL, as is evident in the increases in the genus near the interface. The local entrainment velocity is dominantly affected by the viscous diffusion at the interface.     

Numerical prediction of locally forced turbulent boundary layer

An unsteady numerical simulation was performed to analyze flow structure behind a local suction/blowing in a flat-plate turbulent boundary layer. The local forcing was given to the boundary layer flow by means of a sinusoidally oscillating jet. A version of the unsteady k––f_μ model [Fluid Dyn. Res. 26 (6) (2000) 421] was employed. The Reynolds number based on the momentum thickness was about Re_θ=1700. The forcing frequency was varied in the range 0.011f⁺0.044 with a fixed forcing amplitude A_o=0.4. The predicted results were compared and validated with the experimental data. It was shown that the unsteady locally forced boundary layer flow is predicted well by the k––f_μ model. The time-dependent numerical flow visualizations were demonstrated during one period of the local forcing. The effect of the pitch angle of local forcing on the reduction of skin friction was examined.     

Hybrid LES-RANS: An approach to make LES applicable at high Reynolds number

The main bottle neck for using large eddy simulations (LES) at high Reynolds number is the requirement of very fine meshes near walls. Hybrid LES-Reynolds-averaged Navier-Stokes (RANS) was invented to get rid of this limitation. In this method, unsteady RANS (URANS) is used near walls and away from walls LES is used. The matching between URANS and LES takes place in the inner log-region. In the present paper, a method to improve standard LES-RANS is evaluated. The improvement consists of adding instantaneous turbulent fluctuations (forcing conditions) at the matching plane in order to provide the equations in the LES region with relevant turbulent structures. The fluctuations are taken from a DNS of a generic boundary layer. Simulations of fully developed channel flow and plane asymmetric diffuser flow are presented. Hybrid LES-RANS is used both with and without forcing conditions.     

DNS of a spatially developing turbulent boundary layer with passive scalar transport

A direct numerical simulation (DNS) of a spatially developing turbulent boundary layer over a flat plate under zero pressure gradient (ZPG) has been carried out. The evolution of several passive scalars with both isoscalar and isoflux wall boundary condition are computed during the simulation. The Navier–Stokes equations as well as the scalar transport equation are solved using a fully spectral method. The highest Reynolds number based on the free-stream velocity U_∞ and momentum thickness θ is Re_θ=830, and the molecular Prandtl numbers are 0.2, 0.71 and 2. To the authors’ knowledge, this Reynolds number is to date the highest with such a variety of scalars. A large number of turbulence statistics for both flow and scalar fields are obtained and compared when possible to existing experimental and numerical simulations at comparable Reynolds number. The main focus of the present paper is on the statistical behaviour of the scalars in the outer region of the boundary layer, distinctly different from the channel-flow simulations. Agreements as well as discrepancies are discussed while the influence of the molecular Prandtl number and wall boundary conditions is also highlighted. A Pr scaling for various quantities is proposed in outer scalings. In addition, spanwise two-point correlation and instantaneous fields are employed to investigate the near-wall streak spacing and the coherence between the velocity and the scalar fields. Probability density functions (PDF) and joint probability density functions (JPDF) are shown to identify the intermittency both near the wall and in the outer region of the boundary layer. The present simulation data will be available online for the research community.     

On near-wall turbulence-generating events in a turbulent boundary layer on a riblet surface

The boundary layer over a drag reducing riblet surface is investigated using hot-wire anemometry and flow visualisation. The concept of a riblet sublayer is introduced, and a definition is proposed in terms of a region of reduced turbulence energy production formed near the wall by the addition of riblets. The hot wire records are examined using a modified form of quadrant analysis, and results obtained over plain and riblet surfaces are compared. Close to the wall, the addition of riblets produces a marked reduction in the occurrence of ejection (2nd quadrant) events. A corresponding increase in the incidence of sweep (4th quadrant) events is accompanied by the development of a strong tendency toward a preferred event duration, and a preferred interval between events. These changes diminish rapidly with distance from the surface, becoming almost undetectable beyondy ⁺=40. They are discussed in the light of flow visualisation results, and interpreted in terms of mechanisms associated with the interaction between the riblets and the inner boundary layer flow structures. A conceptual model of the flow mechanisms in the riblet sublayer is proposed.     

LES of the adverse-pressure gradient turbulent boundary layer

We describe large-eddy simulations (LES) of the flat-plate turbulent boundary layer in the presence of an adverse pressure gradient. The stretched-vortex subgrid-scale model is used in the domain of the flow coupled to a wall model that explicitly accounts for the presence of a finite pressure gradient. The LES are designed to match recent experiments conducted at the University of Melbourne wind tunnel where a plate section with zero pressure gradient is followed by section with constant adverse pressure gradient. First, LES are described at Reynolds numbers based on the local free-stream velocity and the local momentum thickness in the range 6560–13,900 chosen to match the experimental conditions. This is followed by a discussion of further LES at Reynolds numbers at approximately 10 times and 100 times these values, which are well out of range of present day direct numerical simulation and wall-resolved LES. For the lower Reynolds number runs, mean velocity profiles, one-point turbulent statistics of the velocity fluctuations, skin friction and the Clauser and acceleration parameters along the streamwise, adverse pressure-gradient domain are compared to the experimental measurements. For the full range of LES, the relationship of the skin-friction coefficient, in the form of the ratio of the local free-stream velocity to the local friction velocity, to both Reynolds number and the Clauser parameter is explored. At large Reynolds numbers, a region of collapse is found that is well described by a simple log-like empirical relationship over two orders of magnitude. This is expected to be useful for constant adverse-pressure gradient flows. It is concluded that the present adverse pressure gradient boundary layers are far from an equilibrium state.     

随机湍流工况低雷诺数风力机翼型优化研究

为改善小型风力机随机湍流工况适应性,以NACA0012翼型为研究对象,采用非嵌入式概率配置点法,获得随机湍流工况下小型风力机叶片翼型运行攻角分布规律;在气动优化中耦合层流分离预测,基于Transition SST模型、拉丁超立方试验设计、Kriging模型和带精英策略非支配排序遗传算法NSGA-II进行高湍流低雷诺数风力机翼型气动优化。结果表明,优化翼型叶片平均风能捕获效率分别提高3.01%和4.76%,标准差分别降低4.76%和14.93%,优化翼型湍流适应性增强。该方法将翼型设计与湍流风况相匹配,为湍流工况低雷诺数翼型及小型风力机设计提供参考。     

Investigation of a hypersonic crossing shock wave/turbulent boundary layer interaction

A combined theoretical and experimental study is presented for the interaction between crossing shock waves generated by (10°, 10°) sharp fins and a flat plate turbulent boundary layer at Mach 8.3. The theoretical model is the full 3-D mean compressible Reynolds-averaged Navier-Stokes RANS) equations incorporating the algebraic turbulent eddy viscosity model of Baldwin and Lomax. A grid refinement study indicated that adequate resolution of the flowfield has been achieved. Computed results agree well with experiment for surface pressure and surface flow patterns and for pitot pressure and yaw angle profiles in the flowfield. The computations, however, significantly overpredict surface heat transfer. Analysis of the computed flowfield results indicates the formation of complex streamline and wave structures within the interaction region.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

The effect of particles on fluid turbulence in a turbulent boundary layer over a cylinder

The turbulent fluid and particle interaction in the turbulent boundary layer for cross flow over a cylinder has been experimentally studied. A phase-Doppler anemometer was used to measure the mean and fluctuating velocities of both phases. Two size ranges of particles (30μm–60μm and 80μm–150μm) at certain concentrations were used for considering the effects of particle sizes on the mean velocity profiles and on the turbulent intensity levels. The measurements clearly demonstrated that the larger particles damped fluid turbulence. For the smaller particles, this damping effect was less noticeable. The measurements further showed a delay in the separation point for two phase turbulent cross flow over a cylinder. The project supported by the National Natural Science Foundation of China     

Fluid forces on a very low Reynolds number airfoil and their prediction

This paper presents the measurements of mean and fluctuating forces on an NACA0012 airfoil over a large range of angle (α) of attack (0-90°) and low to small chord Reynolds numbers (Re_c), 5.3 × 10³-5.1 × 10⁴, which is of both fundamental and practical importance. The forces, measured using a load cell, display good agreement with the estimate from the LDA-measured cross-flow distributions of velocities in the wake based on the momentum conservation. The dependence of the forces on both α and Re_c is determined and discussed in detail. It has been found that the stall of an airfoil, characterized by a drop in the lift force and a jump in the drag force, occurs at Re_c ? 1.05 × 10⁴ but is absent at Re_c = 5.3 × 10³. A theoretical analysis is developed to predict and explain the observed dependence of the mean lift and drag on α.     

Large-scale spanwise periodicity in a turbulent boundary layer induced by highly ordered and directional surface roughness

The effect of converging–diverging riblet-type surface roughness (riblets arranged in a ‘herringbone’ pattern) are investigated experimentally in a zero pressure gradient turbulent boundary layer. For this initial parametric investigation three different parameters of the surface roughness are analysed in detail; the converging–diverging riblet yaw angle α, the streamwise fetch or development length over the rough surface F_x and the viscous-scaled riblet height h⁺. It is observed that this highly directional surface roughness pattern induces a large-scale spanwise periodicity onto the boundary layer, resulting in a pronounced spanwise modification of the boundary layer thickness. Hot-wire measurements reveal that above the diverging region, the local mean velocity increases while the turbulent intensity decreases, resulting in a thinner overall boundary layer thickness in these locations. The opposite situation occurs over the converging region, where the local mean velocity is decreased and the turbulent intensity increases, producing a locally thicker boundary layer. Increasing the converging–diverging angle or the viscous-scaled riblet height results in stronger spanwise perturbations. For the strongest convergent–divergent angle, the spanwise variation of the boundary layer thickness between the diverging and converging region is almost a factor of two. Such a large variation is remarkable considering that the riblet height is only 1% of the unperturbed boundary layer thickness. Increasing the fetch seems to cause the perturbations to grow further from the surface, while the overall strength of the induced high and low speed regions remain relatively unaltered. Further analysis of the pre-multiplied energy spectra suggests that the surface roughness has modified or redistributed the largest scale energetic structures.     

自然对流边界层中湍流的发生

自然对流边界层中从层流到湍流的转捩经历了浮力振型、无摩擦振型和黏性振型的三重流动不稳定性相继产生的前转捩过程,以及近壁迅速出现强湍流源,随之平缓地向自模拟的湍流边界层过渡的热转捩过程.浮力振型在修正Grashof数G>40时开始失稳并成为主要振型,在振幅分布中3种振型的临界层位置处出现3个峰值;在G>100时浮力振型消失,无摩擦振型失稳并成为主要振型,振幅分布中在近壁区还出现黏性振型的峰值;在G>170时无摩擦振型经非线性演化在外层形成较弱的湍流,但内层黏性应力仍远高于湍流应力,振幅分布中仅有与黏性振型相应的峰值,在频谱中黏性振型的基频、第一、第二、第三阶亚谐频随G的增加相继出现,此时黏性不稳定波的高频成分已转化为湍流,但低频成分仍按线性规律增长,直至湍流惯性子区开始形成;至G>800时黏性振型消失,并在G=850附近时近壁区出现强湍流源,湍流应力、湍能产生项和近壁湍流热流率剧增.在热转捩后期,湍流应力和湍能产生项明显下降,流动在内外层趋于平衡.