逆压梯度边界层未分离情况下的湍流模式适用性研究

边界层逆压梯度作用下的流动是许多工程中的一个基础问题,由于逆压梯度作用,流动形态复杂,使得数值模拟有很大的难度。基于雷诺平均纳维-斯托克斯RANS(Reynolds Averaged Navier-Stokes)方程对二维平板逆压梯度边界层作数值计算研究,选取6种代表性的湍流模式,得到局部摩擦系数的数值解,与实验值比较,发现k-ω模式具有很好的精度。基于该湍流模式,给出了湍动能分布,该结果有助于认识逆压梯度边界层流动的复杂特征。     

抽吸和压力梯度在层流边界层转捩过程中的作用

用空间模式的二次稳定性理论研究了抽吸和压力梯度对边界层三维亚谐扰动流动稳定性的影响．数值结果表明，固体边界上的抽吸有明显的层流控制作用，逆压梯度则有较强的不稳定作用．     

抽吸和压力梯度在层流边界层转捩过程中的作用

用空间模式的二次稳定性理论研究了抽吸和压力梯度对边界层三维亚谐扰动流动稳定性的影响．数值结果表明，固体边界上的抽吸有明显的层流控制作用，逆压梯度则有较强的不稳定作用．     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

k-ω SST湍流模式三维激波分离流动修正

“激波?边界层分离”是航空气动领域的典型湍流非平衡流动问题, 准确模拟激波分离对于跨声速飞行器气动性能评估和优化设计具有重要意义. 然而传统涡黏性湍流模式中涡黏性系数的定义方式并不适用于非平衡流动, k-ω SST湍流模式为此引入的Bradshaw假设在应用于三维强逆压梯度和较大分离流动时反而限制了雷诺应力的生成, 导致包括k-ω SST在内的常用涡黏性湍流模式均无法对此类流动进行准确模拟. 同时, 现有的非线性雷诺应力本构关系也并不能有效提高模拟精度. 为此, 针对k-ω SST模式分别提出了基于Bradshaw假设和基于长度尺度的两种激波分离流动修正方法. 前者通过提高Bradshaw常数的方式放宽了对雷诺应力生成的限制, 后者则从湍流长度尺度概念出发, 利用混合长度理论、湍动能生成/耗散之比和一种新定义的长度尺度之比构造了ω方程耗散项修正函数, 提高了模式在三维激波分离流动中的建模长度尺度. 两种方法对ONERA M6机翼跨声速大攻角流动均能得到较雷诺应力模式更好的模拟结果. 进一步的雷诺应力分析表明, 三维激波分离流动中“主雷诺应力分量”的概念不再成立, 各雷诺应力分量大小接近. 网格收敛性分析、对其他攻角状态的验证以及湍流平板边界层壁面律验证进一步确认了所提出的两种修正方法的合理性、有效性和通用性.      

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

压力梯度对边界层两相流动稳定性的影响

分析了有压力梯度的边界层两相流动稳定性,推导出类似于Saffman理论的修正的稳定性方程,数值计算采用高精度的谱方法。结果说明,压力梯度对边界层两相流动稳定性有显著的影响,顺压梯度增强流动稳定性,而逆压梯度则促进流动失稳。在不同的压力梯度和浓度下,Stokes数对流动稳定性的影响是一致的,存在一个临界Stokes数,小Stokes数促进流动失稳,而大Stokes数则提高临界雷诺数,抑制流动失稳的最佳Stokes数为10的量级。     

应用于激波/边界层相互作用的非线性湍流模式

选择8个近年来有代表性的非线性湍流模式，研究2个跨声速激波／边界层相互作用问题．采用的非线性湍流模式包括4个二阶模式和4个三阶模式．2个跨声速激波／边界层相互作用的流动是轴对称圆弧突起绕流和二维管道突起流动．通过数值计算结果和实验结果的比较，对有关的非线性湍流模式进行评估和分析．计算结果表明，非线性模式的模化系数与平均流动应变不变量以及涡量不变量有关，反映了湍流的各向异性，比线性模式优越得多．     

对逆压区湍流边界层拟序结构的观测

用氢气泡法观测了水槽中二元扩散段逆压区湍流边界层的拟序结构。正面观测所得的快斑与慢斑和侧面观测所得的涡状运动等图象都显著地与于无压力梯度时的不同。慢斑较宽,慢斑内部流动平静,慢斑边缘产生紊乱流动。涡较大,可看到涡的诱导作用。     

管道凸起激波/湍流边界层干扰的数值模拟

采用基于部分平均的湍流方程组,模拟了管道凸起流动中两种不同强度的激波-湍流边界层干扰现象,对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布、边界层厚度和速度型与实验值比较吻合很好.结果表明基于部分平均的湍流方程组既能够准确模拟小分离的激波/湍流边界层干扰流动又能够准确模拟大分离的激波/湍流边界层干扰流动,较好地预测了典型的λ激波结构.     

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.     

The Effect of Pressure Gradient on Boundary Layer Receptivity

This paper presents numerical results for the receptivity of three laminar boundary layers with zero (ZPG), adverse (APG) and favourable (FPG) pressure gradients. Each boundary layer is subjected to a series of simple freestream waveforms which can be considered as constituent parts of either an isotropic or a non-isotropic turbulent freestream. Each freestream waveform has a single frequency in each spatial direction and is divided into two mutually perpendicular components. The first component has a zero spanwise velocity and hence lies in the streamwise normal plane whereas the second component lies in a plane which is perpendicular both to this plane and the spatial frequency vector. High boundary layer receptivities are only obtained for a minority of these waveforms and so only the resulting flow structures for these waveforms are considered in detail. The dominant flow structures are identified as either Tollmien Schlichting (T-S) waves or streaky structures. The streaky structures can be induced by both freestream components, but the response to the second component, which results in streamwise vortices in the freestream, is considerably stronger and occurs over a much larger streamwise frequency range. The boundary layer is only receptive to a relatively narrow band of spanwise wavelengths ranging from approximately one to four times the local boundary layer thickness. The APG leads to receptivities which are more than double those for the FPG case. The ratio of the freestream fluctuation streamwise wavelength to the distance from the plate leading edge is identified as an important influential parameter for receptivity leading to streaks. Significant T-S activity is only observed for APG, but is also detected for ZPG.     

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 study of adverse pressure gradient generation over a flat plate using a rotating cylinder

Generating an adverse pressure gradient (APG), using a rotating cylinder in the proximity of a plane wall under a laminar freestream flow, is studied numerically in this work. The magnitude of the generated APG is a function of the gap, G, between the cylinder and the wall, and the rotational speed of the cylinder, Ω. The flow in such a configuration is characterized by periodic transient vortex shedding at high Reynolds number. A numerical model for the computation of the transient flow for this configuration is developed using the ANSYS CFD simulation tool. The model is validated against published experimental and numerical data for similar flow configurations and excellent agreement is observed. A parametric study is carried out for different combinations of G and Ω for two different Reynolds numbers of 200 and 1000 to examine the development of the resulting separation bubble due to the generated APG. The mechanism of the boundary layer separation over the plane wall and the corresponding wake dynamics is investigated. Results are presented in terms of the distribution of the pressure coefficient as well as skin friction coefficient along the wall and flow patterns around and downstream of the cylinder in the proximity of the wall. The results of these computations confirm that using a rotating cylinder over a plane wall in a freestream flow is an effective technique to generate a controlled range of adverse pressure gradients.     

Revisiting History Effects in Adverse-Pressure-Gradient Turbulent Boundary Layers

The goal of this study is to present a first step towards establishing criteria aimed at assessing whether a particular adverse-pressure-gradient (APG) turbulent boundary layer (TBL) can be considered well-behaved, i.e., whether it is independent of the inflow conditions and is exempt of numerical or experimental artifacts. To this end, we analyzed several high-quality datasets, including in-house numerical databases of APG TBLs developing over flat-plates and the suction side of a wing section, and five studies available in the literature. Due to the impact of the flow history on the particular state of the boundary layer, we developed three criteria of convergence to well-behaved conditions, to be used depending on the particular case under study. (i) In the first criterion, we develop empirical correlations defining the R e _?? -evolution of the skin-friction coefficient and the shape factor in APG TBLs with constant values of the Clauser pressure-gradient parameter β = 1 and 2 (note that β = δ ^?/τ _w dP _e /dx, where δ ^? is the displacement thickness, τ _w the wall-shear stress and dP _e /dx the streamwise pressure gradient). (ii) In the second one, we propose a predictive method to obtain the skin-friction curve corresponding to an APG TBL subjected to any streamwise evolution of β, based only on data from zero-pressure-gradient TBLs. (iii) The third method relies on the diagnostic-plot concept modified with the shape factor, which scales APG TBLs subjected to a wide range of pressure-gradient conditions. These three criteria allow to ensure the correct flow development of a particular TBL, and thus to separate history and pressure-gradient effects in the analysis.     

Higher-order and length-scale statistics from DNS of a decelerated planar wall-bounded turbulent flow

A DNS database is presented to document third- and fourth-order moments and their budgets for fully developed plane-channel flow and for strained plane-channel flow. The effect of straining has a similar effect on statistics, such as the skewness and flatness of velocities, as in an adverse-pressure-gradient (APG) boundary layer. In addition to higher-order statistics, some modeling implications are also described, including issues related to decomposition of the velocity–pressure gradient correlations, the assumption of dissipation isotropy and a fourth-order turbulence time scale. An analysis of two-point correlations along the inhomogeneous direction is made to include the effects of APG straining. This reveals the advantage of a one-point length-scale based on wall-normal velocity fluctuations as an alternative to the traditionally used length-scale based on turbulence kinetic energy. The present study should prove to be useful for turbulence modelers in need of data to develop and assess higher-than-second-order Reynolds-averaged Navier–Stokes closure models.     

Adverse pressure gradient turbulent flows over rough walls

An experimental study of a fully developed turbulent channel flow and an adverse pressure gradient (APG) turbulent channel flow over smooth and rough walls has been performed using a particle image velocimetry (PIV) technique. The rough walls comprised two-dimensional square ribs of nominal height, k = 3 mm and pitch, p = 2k, 4k and 8k. It was observed that rib roughness enhanced the drag characteristics, and the degree of enhancement increased with increasing pitch. Similarly, rib roughness significantly increased the level of turbulence production, Reynolds stresses and wall-normal transport of turbulence kinetic energy and Reynolds shear stress well beyond the roughness sublayer. On the contrary, the distributions of the eddy viscosity, mixing length and streamwise transport of turbulence kinetic energy and Reynolds shear stress were reduced by wall roughness, especially in the outer layer. Adverse pressure gradient produced a further reduction in the mean velocity (in comparison to the results obtained in the parallel section) but increased the wall-normal extent across which the mean flow above the ribs is spatially inhomogeneous in the streamwise direction. APG also reinforced wall roughness in augmenting the equivalent sand grain roughness height. The combination of wall roughness and APG significantly increased turbulence production and Reynolds stresses except in the immediate vicinity of the rough walls. The transport velocities of the turbulence kinetic energy and Reynolds shear stress were also augmented by APG across most part of the rough-wall boundary layer. Further, APG enhanced the distributions of the eddy viscosity across most of the boundary layer but reduced the mixing length outside the roughness sublayer.     

Higher-order moments and spectra of velocity fluctuations in adverse-pressure-gradient turbulent boundary layer

The characteristics of turbulent boundary layer flows with adverse pressure gradients (APGs) differ significantly from those of canonical boundary layers. We have investigated the effects of an APG on the higher-order moments and spectra of velocity fluctuations. The local wavelet spectra reveal a large difference in energy-containing frequencies of streamwise and wall-normal components of turbulent velocities, which results in smaller Reynolds shear stress production. Moreover, an analysis of bispectra in the Fourier space has revealed that non-local interactions, consisting of streamwise fluctuating velocity with low frequency and wall-normal velocity with high frequency, occur in the APG flow. However, the small-scale motions are not affected by imposing an APG.     

Development of a generalized multi‐layer model for 3‐D simulation of free surface flows

A mathematical model was developed for three‐dimensional (3‐D) simulation of free surface flows. In this model, the flow depth is divided into a number of layers and shallow water equations are integrated in each layer to derive the hydrodynamic equations. To give a general form to this model, each layer is assumed to be non‐horizontal with varying thickness in the flow domain. A non‐orthogonal curvilinear coordinate system is employed in the model, to allow for flexibility in dealing with the irregular geometry of natural watercourses. Due to the similarity in governing equations, two‐dimensional (2‐D) depth averaged programs can be developed into a multi‐layer model. The development for a depth averaged program and its numerical scheme is described in this paper. Experimental data and semi‐analytical solutions are used to evaluate the performance of the model. Three different cases of open channel flow are tested: 1‐flow in a straight open channel, 2‐the flow development region in a channel, and 3‐flow in a meandering channel. It is shown that the model has the capability to predict velocity distribution and secondary flows in complex 3‐D flow conditions. Copyright © 2004 John Wiley & Sons, Ltd.     

Characterizing developing adverse pressure gradient flows subject to surface roughness

An experimental study was conducted to examine the effects of surface roughness and adverse pressure gradient (APG) on the development of a turbulent boundary layer. Hot-wire anemometry measurements were carried out using single and X-wire probes in all regions of a developing APG flow in an open return wind tunnel test section. The same experimental conditions (i.e., T _∞, U _ref, and C _p) were maintained for smooth, k ⁺ = 0, and rough, k ⁺ = 41–60, surfaces with Reynolds number based on momentum thickness, 3,000 < Re _θ < 40,000. The experiment was carefully designed such that the x-dependence in the flow field was known. Despite this fact, only a very small region of the boundary layer showed a balance of the various terms in the integrated boundary layer equation. The skin friction computed from this technique showed up to a 58% increase due to the surface roughness. Various equilibrium parameters were studied and the effect of roughness was investigated. The generated flow was not in equilibrium according to the Clauser (J Aero Sci 21:91–108, 1954) definition due to its developing nature. After a development region, the flow reached the equilibrium condition as defined by Castillo and George (2001), where Λ = const, is the pressure gradient parameter. Moreover, it was found that this equilibrium condition can be used to classify developing APG flows. Furthermore, the Zagarola and Smits (J Fluid Mech 373:33–79, 1998a) scaling of the mean velocity deficit, U _∞δ*/δ, can also be used as a criteria to classify developing APG flows which supports the equilibrium condition of Castillo and George (2001). With this information a ‘full APG region’ was defined.