壁湍流多尺度相干结构复涡黏模型的实验研究

在湍流相干结构动力学方程中,非相干结构成分对相干结构贡献的雷诺应力的模型为涡黏性模型,即涡黏性系数乘以相干结构平均速度变形率的形式.基于非相干结构成分对相干结构贡献的雷诺应力与相干结构速度变形率之间存在相位差的事实,在理论上提出了非相干结构成分对相干结构贡献的雷诺应力复涡黏性模型的假设.应用热线测速技术,在低速风洞中对湍流边界层非相干结构成分对相干结构贡献的雷诺应力与相干结构法向速度变形率之间的相位关系进行了实验测量.通过分析湍流相干结构猝发过程中非相干结构成分对相干结构贡献的雷诺应力与相干结构速度变形率之间的相位关系,研究了相干结构雷诺应力分量与流向速度法向梯度之间的相位差沿湍流边界层法向的变化规律,肯定了湍流相干结构复涡黏性系数模型的合理性.     

Experimental research on complex eddy viscosity modeling of multi-scale coherent structures in wall turbulence

在湍流相干结构动力学方程中,非相干结构成分对相干结构贡献的雷诺应力的模型为涡黏性 模型,即涡黏性系数乘以相干结构平均速度变形率的形式. 基于非相干结构成分对相干结构贡 献的雷诺应力与相干结构速度变形率之间存在相位差的事实,在理论上提出了非相干结构成 分对相干结构贡献的雷诺应力复涡黏性模型的假设. 应用热线测速技术,在低速风洞中对湍 流边界层非相干结构成分对相干结构贡献的雷诺应力与相干结构法向速度变形率之间的相位 关系进行了实验测量. 通过分析湍流相干结构猝发过程中非相干结构成分对相干结构贡献的 雷诺应力与相干结构速度变形率之间的相位关系,研究了相干结构雷诺应力分量与流向速度 法向梯度之间的相位差沿湍流边界层法向的变化规律,肯定了湍流相干结构复涡黏性系数模 型的合理性.     

湍流边界层复涡黏模式的实验研究

在开口式循环水槽底部湍流边界层外区引入周期性扰动。利用X型热膜探针在扰动下游进行测量。用实验的方法研究了周期性大尺度结构下壁湍流涡黏模式中涡黏系数的形式，结果发现和周期性扰动对应的变形率及与之对应的雷诺应力间存在着相位差。这是目前许多最终导致涡黏系数的湍流模式理论都没有考虑到的一个重要因素。     

圆柱尾流雷诺应力与平均运动变形率的空间弛豫效应

针对圆柱尾流中沿流向存在的Karman涡街周期性涡旋结构,对湍流雷诺应力与平均运动变形率之间的空间弛豫效应进行了实验研究.在回流式水槽中,放入不同直径的圆柱模型,获得不同雷诺数下的圆柱尾流,利用二维高时间分辨率粒子图像测速(TRPIV)技术测量圆柱尾流二维瞬时速度空间分布图像的时间序列.经过数字图像处理,获得二维雷诺应...     

含有周期涡结构的湍流雷诺应力与平均速度梯度之间的迟滞效应

针对存在时空周期涡旋结构的湍流场，对湍流雷诺应力与平均速度梯度之间的空间迟滞效应进行了实验研究，发现二者流向相位差呈周期变化．在低速回流式水槽中，利用二维高时间分辨率粒子图像测速(TimeResolved Particle Image Velocimetry, TRPIV)技术，对Re=324的圆柱尾流流场进行测量．经过49个周期166个相位的8 215个PIV瞬时流场，经过周期相位平均，得到一个周期内不同相位典型的湍流雷诺应力和平均速度梯度的空间分布．利用湍流雷诺应力相位平均图像和平均速度梯度图像在不同时间相位下的空间互相关函数最大值对应的流向空间距离，得到湍流雷诺应力与平均速度梯度之间沿流向的空间相位差，并绘制流向相位差随周期相位的演变过程．本文验证了湍流复涡黏模型的合理性，为建立符合周期涡旋结构物理机理的湍流模型研究提供了实验依据．     

壁湍流猝发过程中速度分量的相位差对雷诺应力影响的实验研究

用IFA300恒温热线风速仪和×形二分量热线探针,以采样间隔小于最小湍流时间尺度的分辨率,精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、展向速度分量的时间序列信号.用子波分析辨识壁湍流相干结构猝发事件的能量最大准则,确定壁湍流相干结构猝发事件的时间尺度;用条件相位平均技术提取了相干结构猝发过程中流向、展向脉动速度分量条件相位平均波形,用互相关方法研究了相干结构猝发过程中流向、展向脉动速度分量条件相位平均波形的相位差关系及其对雷诺应力的影响,发现在缓冲层和对数律区,展向脉动速度与流向脉动速度的条件相位平均波形具有不同的相位;当两者相位基本一致时,雷诺应力达到正的最大值,此时湍流相干结构的产生非常活跃;当两者相位差分别集中在90°和270°附近时,雷诺应力的幅值减小并接近于零,此时湍流相干结构的产生和猝发都得到了抑制.     

雷诺应力各向异性涡黏模型的层析TRPIV测量

利用层析TRPIV测量水洞中平板湍流边界层3D-3C速度场的高分辨率时间序列数据库. 提出了空间局部平均多尺度速度结构函数的新概念, 描述湍流多尺度涡结构的空间拉伸、压缩、剪切变形和旋转. 用空间局部平均多尺度速度结构函数对湍流脉动速度进行了空间多尺度分解. 用空间流向局部平均多尺度速度结构函数, 根据湍流多尺度涡结构在流向的拉伸和压缩物理特征, 提出了新的湍流相干结构条件采样方法, 检测并提取了层析TRPIV数据中相干结构“喷射”和“扫掠”事件中的脉动速度、平均速度变形率、雷诺应力等物理量的空间拓扑形态. 通过研究平均速度变形率各分量与雷诺应力各分量之间的空间相位差异,肯定了壁湍流相干结构雷诺应力各向异性复涡黏模型的合理性.      

可压缩均匀剪切湍流的二阶矩模拟

基于对可压缩湍流中脉动压力场和脉动速度场特征的理论分析以及ＤＮＳ结果,建立了可均匀剪切湍流中压力－变形率关联的压缩性修正模式,应用这个模式,加上Ｓａｒｋａｒ等建立的脉动体胀率项（ｄｉｌａｔａｔｉｏｎａｌ ｔｅｒｍｓ）的模式,预测可压缩均匀剪切湍流随时间的发展,所得雷诺应力各是性张量的平衡值与Ｂｌａｉｓｄｅｌｌ等的ＤＮＳ数据非常一致。这个模式准确地预测出均匀剪切湍流中压缩性导致的雷诺应力结构的“流向     

湍流边界层相干结构剩余脉动雷诺应力项的实验测量

用IFA300恒温热线风速仪和X形二分量热线探针以采样间隔小于湍流耗散时间尺度的分辨率精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、法向速度分量的时间序列信号。用条件采样和相位平均技术提取了相干结构猝发过程中相干结构剩余脉动雷诺应力和随机脉动对相干结构贡献的雷诺应力的条件相位平均波形。基于理论上对湍流相干结构动量方程中随机脉动对相干结构贡献的雷诺应力和相干结构剩余脉动雷诺应力项的分析,对两种雷诺应力项进行了对比研究。研究发现,相干结构剩余脉动雷诺应力项在数值上具有和随机脉动对相干结构贡献的雷诺应力相同的数量级,表明在相干结构动力学模型方程中,相干结构剩余脉动雷诺应力项并不像以前估计的那样可以忽略不计。     

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

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

Lower-Dimensional Manifold (Algebraic) Representation of Reynolds Stress Closure Equations

For complex turbulent flows, Reynolds stress closure modeling (RSCM) is the lowest level at which models can be developed with some fidelity to the governing Navier–Stokes equations. Citing computational burden, researchers have long sought to reduce the seven-equation RSCM to the so-called algebraic Reynolds stress model which involves solving only two evolution equations for turbulent kinetic energy and dissipation. In the past, reduction has been accomplished successfully in the weak-equilibrium limit of turbulence. In non-equilibrium turbulence, attempts at reduction have lacked mathematical rigor and have been based on ad hoc hypotheses resulting in less than adequate models.?In this work we undertake a formal (numerical) examination of the dynamical system of equations that constitute the Reynolds stress closure model to investigate the following questions. (i) When does the RSCM equation system formally permit reduced representation? (ii) What is the dimensionality (number of independent variables) of the permitted reduced system? (iii) How can one derive the reduced system (algebraic Reynolds stress model) from the full RSCM equations? Our analysis reveals that a lower-dimensional representation of the RSCM equations is possible not only in the equilibrium limit, but also in the slow-manifold stage of non-equilibrium turbulence. The degree of reduction depends on the type of mean-flow deformation and state of turbulence. We further develop two novel methods for deriving algebraic Reynolds stress models from RSCM equations in non-equilibrium turbulence. The present work is expected to play an important role in bringing much of the sophistication of the RSCM into the realm of two-equation algebraic Reynolds stress models. Another objective of this work is to place the other algebraic stress modeling efforts in the lower-dimensional modeling context. Received 19 November 1999 and accepted 3 August 2000     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

Rotation invariant constitutive relation for Reynolds stress structure parameter

A new Reynolds stress constitutive formula is constructed using the firstorder statistics of turbulent fluctuations instead of the mean strain rate. It includes zero empirical coefficients. The formula is validated with the direct numerical simulation(DNS) data of turbulent channel flow at Reτ =180. The Reynolds stresses given by the proposed formula agree very well with the DNS results. The good agreement persists even after the multi-angle rotation of the coordinate system, indicating the rotation invariance of the formula. The autocorrelation of the fluctuating velocity rather than the mean strain rate is close to the essence of the Reynolds stress.     

A dynamic global-coefficient mixed subgrid-scale model for large-eddy simulation of turbulent flows

A dynamic global-coefficient mixed subgrid-scale eddy-viscosity model for large-eddy simulation of turbulent flows in complex geometries is developed. In the present model, the subgrid-scale stress is decomposed into the modified Leonard stress, cross stress, and subgrid-scale Reynolds stress. The modified Leonard stress is explicitly computed assuming a scale similarity, while the cross stress and the subgrid-scale Reynolds stress are modeled using the global-coefficient eddy-viscosity model. The model coefficient is determined by a dynamic procedure based on the global-equilibrium between the subgrid-scale dissipation and the viscous dissipation. The new model relieves some of the difficulties associated with an eddy-viscosity closure, such as the nonalignment of the principal axes of the subgrid-scale stress tensor and the strain rate tensor and the anisotropy of turbulent flow fields, while, like other dynamic global-coefficient models, it does not require averaging or clipping of the model coefficient for numerical stabilization. The combination of the global-coefficient eddy-viscosity model and a scale-similarity model is demonstrated to produce improved predictions in a number of turbulent flow simulations.     

Anisotropic Organised Eddy Simulation for the prediction of non-equilibrium turbulent flows around bodies

The unsteady turbulent flow around bodies at high Reynolds number is predicted by an anisotropic eddy-viscosity model in the context of the Organised Eddy Simulation (OES). A tensorial eddy-viscosity concept is developed to reinforce turbulent stress anisotropy, that is a crucial characteristic of non-equilibrium turbulence in the near-region. The theoretical aspects of the modelling are investigated by means of a phase-averaged PIV in the flow around a circular cylinder at Reynolds number 1.4×10⁵. A pronounced stress–strain misalignment is quantified in the near-wake region of the detached flow, that is well captured by a tensorial eddy-viscosity concept. This is achieved by modelling the turbulence stress anisotropy tensor by its projection onto the principal matrices of the strain-rate tensor. Additional transport equations for the projection coefficients are derived from a second-order moment closure scheme. The modification of the turbulence length scale yielded by OES is used in the Detached Eddy Simulation hybrid approach. The detached turbulent flows around a NACA0012 airfoil (2-D) and a circular cylinder (3-D) are studied at Reynolds numbers 10⁵ and 1.4×10⁵, respectively. The results compared to experimental ones emphasise the predictive capabilities of the OES approach concerning the flow physics capture for turbulent unsteady flows around bodies at high Reynolds numbers.