湍流标度律的实验测量

实验测量了射流与圆柱尾流纵、横向的湍流速度脉动，计算了这两种方向的湍流速度结构函数的标度规律。结果表明，横向的湍流速度脉动比纵向有着更大的间歇性，这支持目前已有的一些关于均匀各向同性湍流的实验测量和数值计算的结果，而与另一些文献报道不同。此外，还对纵、横向湍流速度结构函数标度律的S－L模型进行了一些分析。     

壁湍流标度律的实验研究

对风洞中零压力梯度平板湍流边界层进行了实验研究，用热线风速仪测量了不同法向位置的脉动速度，研究了平板湍流边界层不同法向位置速度结构函数的标度律。     

用PIV实验研究边界层内结构函数标度律

湍流标度律的变化情况显示了湍流边界层内的流动状况。利用高分辨率、高帧率的粒子图像测速技术（PIV）对平板湍流边界层进行了实验研究。研究结果显示,当脉动速度增量结构函数的阶数为正整数且逐渐增大时,由实验所得的脉动速度增量结构函数的标度指数逐渐偏离K41和SL标度指数,出现奇异标度律;而当结构函数的阶数为正分数时,实验所得的标度指数和K41及SL标度指数较为接近;当结构函数的阶数由0趋向于-1时,实验的标度指数明显偏离K41和SL标度指数,出现奇异标度律。这也说明了用低阶结构函数的标度指数也能说明标度律的奇异性及边界层内存在间歇性。对于法向速度增量和涡量分量来说,不管结构函数的高阶或低阶的标度指数都不随壁面位置的不同而改变,具有普适性。     

壁面加热湍流标度律的实验研究

用热线风速仪测量了风洞中壁面加热和不加热平板湍流边界层不同法向位置的流向速度和温度的时间序列信号，研究了壁面加热的边条件和不同法向位置的剪切作用对速度和温度增量的p阶结构函数的非线性反常标度律的影响，实验结果表明，壁面加热的边界条件不同法向位置的剪切作用对层次结构模型中间歇参数β和最奇异指数γ的影响明显存在，并初步分析了在具有强剪切作用的壁湍流中这种影响产生的原因。     

用子波系数概率密度函数研究湍流多尺度结构的间歇性

用热线风速仪在充分发展平板湍流边界层中测量流向速度信号的时间序列并进行湍流多尺度结构统计特性的实验研究。用热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量风洞中平板湍流边界层不同法向位置的流向速度分量的时间序列信号,用子波分析对湍流脉动速度信号进行多尺度分解,将推广的自相似性(ESS)应用于子波分析。用子波系数的概率密度函数验证了多尺度相干结构对间歇性的影响,并证明与距壁面距离有关。对比研究了提取湍流边界层多尺度相干结构前后湍流多尺度结构函数间歇性的变化规律及概率密度函数的变化规律,验证了多尺度相干结构是湍流边界层产生奇异标度律的原因。当这些相干结构被分离出后,由子波系数计算的ESS标度律就回到K olm ogorov线性标度律,这一现象在整个湍流边界层都是很明显的。      

高温化学非平衡湍流边界层直接数值模拟

高超声速飞行器在特定工况下飞行时，其表面会存在湍流与化学非平衡耦合作用，但考虑化学反应的高温湍流边界层的研究工作还十分有限。选取高超声速楔形体头部斜激波后的流动状态，分别采用热完全气体模型、空气化学反应模型开展直接数值模拟(DNS)研究，对比分析了化学非平衡效应对湍流统计量、湍流脉动量的影响趋势，并研究了高温湍流边界层标度律。计算结果表明，以吸热离解反应为主的化学非平衡效应会使边界层平均温度和脉动温度显著降低，但对平均速度剖面、速度脉动、雷诺应力的影响较小；在高温化学非平衡湍流边界层的对数区，各流动参数的拓展自相似理论的相对标度指数基本符合标度律。     

可压缩湍流的多尺度分析

作者的研究团队近几年在可压缩湍流的多尺度性质方面开展了系统的研究工作.通过多过程分解,研究了可压缩湍流中的速度和热力学量的剪切部分、胀压部分、伪声模态、声模态、熵模态的多尺度性质,并总结了各类可压缩条件下的速度和热力学量的谱的标度律.通过滤波方法,研究了动能和热力学量的多尺度传输现象,并重点分析了可压缩性对动能胀压部分...     

磁流体湍流中的结构和能量传输

磁流体湍流中存在多种相干结构,包括电流和涡旋结构。文章主要对近几年有关磁流体湍流相干结构及其相关的能量传输的一些工作进行了介绍。不同于中性流体湍流中的管状结构,磁流体中的相干结构多呈现为片状,强电流片附近也存在强涡片结构。磁流体湍流中的能量级串使能量跨尺度从大尺度传递至小尺度;与中性流体湍流能量传输主要集中在相近尺度上不同,磁流体湍流除了动能传输,还有磁能传输,以及动能与磁能之间的跨尺度互相转换。研究表明,磁流体湍流中的能量传输及耗散具有强间歇性,集中在较小区域范围,且与相干结构存在关联。在存在背景磁场的情况下,磁流体湍流中的相干结构在背景磁场方向上被拉长,而同样方向上的能量传输受到抑制。在高马赫数的情况下,激波的产生会使能量传输明显增强,同时带来很强的间歇性,结构函数标度律远远偏离线性标度关系,并且随着阶数的增加,标度指数会达到饱和值。     

湍流边界层相干结构空间拓扑形态的 层析TRPIV测量

用层析TRPIV(time-resolved paticle image velocimetry)技术精细测量了水洞中平板湍流边界层三分量速度的时空序列信号,提出了空间局部平均多尺度速度结构函数的新概念描述湍流多尺度涡结构的空间拉伸、压缩、剪切变形和旋转.用空间局部平均多尺度速度结构函数对湍流脉动速度进行了空间多尺度分解.用空间局部平均多尺度速度结构函数的新概念,根据湍流多尺度涡结构在空间流向的拉伸和压缩特征,提出了新的湍流相干结构条件采样方法,检测并提取了层析TRPIV数据中相干结构的“喷射”和“扫掠”事件中的速度、涡量等物理量的空间拓扑形态.发现在喷射和扫掠事件中均存在一对反向旋转的准流向“马蹄形”涡结构.      

壁湍流调制雷诺应力和拟序结构应变率空间形态的层析PIV实验研究

采用高时间分辨率的层析PIV技术,测量了水洞壁湍流三维速度分量空间分布的时间序列,应用局部平均速度结构函数概念对流向速度信号进行多尺度分解,以流向速度分量的局部平均速度结构函数过零点作为特征量检测壁湍流中拟序结构猝发的喷射和扫掠过程,应用空间相位平均技术提取拟序结构猝发的喷射和扫掠过程各速度分量、涡量分量、拟序结构速度应变率分量以及调制雷诺应力分量的空间相位平均拓扑形态。为了引入平衡态涡粘模型假设模拟调制雷诺应力,研究了拟序结构猝发过程中调制雷诺应力分量和速度变形率分量的空间分布形态,发现两者之间的空间相位分布不一致。由于存在时空相位不同步性,说明需要考虑大尺度拟序涡结构引起动量传递的时空弛豫效应。应用经典的线性平衡态下的Boussinesq涡粘模型不能准确地描述壁湍流拟序结构动量传递非平衡现象的物理机理。对于壁湍流拟序结构动力学方程中调制雷诺应力的模拟,应采用包含时空相位信息的复涡粘张量模型。由于雷诺应力与速度变形率的时空不同步性,对非平衡非局部湍流场的数值模拟提出了一个挑战性的问题,建议采用包含相位信息的复涡粘张量模型来模拟雷诺应力张量,从而更加符合雷诺应力演化的物理机理,这一模型有可能成为一个很有发展前景的封闭模型,从而更加准确地预测工业领域中广泛存在的非平衡湍流。     

基于互相关的光流测速法在湍流热对流中的应用

同时准确地测量大尺度流动结构和小尺度湍流统计特性对于理解湍流的能量传递规律有着重要意义。尽管粒子成像测速（PIV）技术得到了广泛的应用，但由于算法原理的要求以及商用相机分辨率的限制，PIV技术在同时获取大、小尺度流动特性方面并不实用。另一方面，近年来逐渐引起关注的光流测速法虽然可以实现单像素分辨率的大区域速度场测量，但其准确性通常受限于应用场景而难以保证。本文采用一种基于互相关的光流测速法来分析RayleighBénard对流的粒子图像，进而同时获得全局速度场和小尺度湍流特性。结果表明，这种方法不仅可以获得与PIV技术一致的全局场和流动强度，还能有效地解析出湍流耗散区的速度结构函数。进一步的，分别基于速度结构函数和定义直接计算而获得的能量耗散率不仅在数值大小上互相吻合，还与直接数值模拟的结果以及理论预测相一致。这些结果表明基于互相关的光流测速法在同时测量大尺度流动结构和小尺度统计特性上有很好的应用潜力。     

均匀各向同性大气湍流对声爆传播特性的影响

基于离散Fourier模态有限和生成的随机大气湍流场，采用修正波形参数方法，开展了均匀各向同性大气湍流对典型超声速客机声爆传播特性的影响分析。计算采用的超声速客机模型为自行生成的简化超声速公务机模型。首先，应用航空工业空气动力研究院自主研制的CFD软件ARI_Overset在三维空间求解Navier-Stokes方程，得到作为声爆远场传播初始值的近场空间压力分布；其次，基于Von Karman能量谱，采用离散Fourier模态有限和形式生成随机均匀各向同性大气湍流场；最后，采用修正波形参数方法模拟了声爆信号在随机速度湍流场中的传播过程。数值结果表明：各向同性大气湍流对于地面声爆特征有重要影响。与无湍流状态相比，均匀各向同性大气湍流使得地面声爆特征增强的概率约为55%，使得地面声爆特征减弱的概率约为45%，故总体而言大气湍流效应更倾向于增强地面声爆特征；均匀各向同性大气湍流对于声爆传播路径影响相对较小，但是这一变化仍会导致地面信号接收点的不确定性。     

湍流射流扩散火焰的层流火焰面模拟

采用双尺度κ-ε湍流模型，标量联合的概率密度函数（PDF）输运方程和层流火争面模型相结合，模拟湍流射流扩散火焰，比较了几种不同κ-ε湍流模型对轴对称湍流射流速度场模拟的结果，建立了混合物分数和湍流频率的联合PDF输运方程，并假定湍流频率和标量 线性关系，得到标量耗散率的概率分布，化学组分的浓度和温度由层流火焰面系综的平均得到，计算结果与直接用化学反应动力学计算的结果及实验数据做了比较。     

吹吸扰动控制猝发过程中雷诺应力的实验

采用热线测速技术和双丝边界层热线探针,在风洞中以高于对应最小湍流时间尺度的分辨率,精细测量了平板湍流边界层施加不同频率周期性质量引射前后,不同流向和法向位置的瞬时流向、法向和展向速度分量的时间序列信号.采用流向脉动速度局部平均结构函数模极大值法检测了壁湍流猝发过程中拟序结构不同速度分量和雷诺应力分量的条件相位平均波形.发现在扰动源下游一定范围内,摩擦因数和雷诺应力幅值有所减小,扰动改变了法向或展向脉动速度与流向脉动速度的条件相位平均波形的相位差,使其保持在π/2的奇数倍左右,此时雷诺应力几乎为0,改变了湍流边界层原有雷诺应力的产生和扩散输运机制,抑制了湍流雷诺应力的产生和输运.      

用二阶矩两相湍流模型模拟鼓泡床内气泡—液体湍流两相流动

建立了二阶矩气液两相湍流模型，模拟了不同工况下二维矩形断面气液鼓泡床中气泡－液体湍流两相流动，给出了气泡和液体的速度场、气泡体积分数和两相雷诺应力分布，基本工况的模拟结果与美国俄亥俄州立大学PIV测量结果符合较好。文中研究了不同气体表观速度对两相流场的影响。模拟结果显示了鼓泡床内液体的回流流动和气泡的上升运动、各向异性的上湍流，气泡湍流脉动比液体的强，以及气泡体积分数和两相湍流强度随着气体表观速度的增大而增大等规律。     

湍流参数对复杂形状涡轮盘腔流场的数值研究

以罗罗公司未来某型双轴发动机的高压涡轮旋转盘腔为计算模型，目的是要对复杂的内部旋转空腔流场有一个详细的了解，通过CFD(FLUENT6．0．12)的研究为即将进行的放大比例的旋转涡轮盘腔试验提供指导。计算中采用了Spalart-Allmaras湍流模型，结果表明湍流参数相同的流动，其流体结构是非常相似的；流场主要由湍流参数λT决定涡旋的大小和数量；对于湍流参数小于0．199的旋转流动，盘腔内存在着自由涡流。     

基于人工神经网络的可压缩湍流大涡模拟模型

在国家数值风洞（NNW）工程项目的指导下，空间人工神经网络（SANN）模型被用于强可压缩湍流大涡模拟（LES）研究，其中流场的湍流马赫数分别为0.6、0.8、1.0。基于湍流的多尺度空间结构特性和人工神经网络方法发展的高精度空间神经网络（SANN）模型适用于不可压缩湍流和弱可压缩湍流。对于强可压缩湍流，流场中会出现激波结构，给大涡模拟带来了挑战。本文的研究结果表明：SANN模型适用于强可压缩湍流的大涡模拟。在先验分析中，SANN模型预测的亚格子应力和亚格子热流的相关系数超过0.995，远远高于梯度模型和近似反卷积模型等传统模型；传统模型的相对误差大于30%，而SANN模型在这方面有很大的改进，相对误差低于11%。在后验分析中，与隐式大涡模拟（ILES）、动态Smagorinsky模型（DSM）、动态混合模型（DMM）相比，SANN模型能更精确地预测能谱、各类湍流统计特性以及瞬态流动结构。因此，基于湍流多尺度空间结构特性的人工神经网络模型加深了人们对强可压缩湍流亚格子建模的认识，同时可以服务于NNW工程的流体力学模型构造。     

加力燃烧室单相燃烧模型及其应用

提出了描述航空发动机加力燃烧室内气相湍流燃烧的模型方程，并考虑了加力喷油效应。通过结合双时间尺度ｋ－ε湍流模型、热通量高温辐射模型及ＳＩＭＰＬＥ方法，对实际运行的加力燃烧室内的湍流燃烧进行了数值模拟，得到了数值流场图及温度、组分浓度的等值线分布，给出了燃烧效率和总压恢复系数，数值解证明了湍流燃烧模型的可解性与合理性。     

航空发动机燃烧室流动数值计算中湍流模型的比较

分别采用标准k-ε湍流模型、RNGk-ε湍流模型、Realizablek-ε湍流模型以及雷诺应力模型,对某型航空发动机燃烧室流动进行了数值计算.近壁处理采用标准壁面函数法,计算得到速度矢量分布以及质量流量、湍流粘度比和湍流强度等参数.四种湍流模型计算的总体流动差别较小,但流场的细节有较明显的不同.标准k-ε模型、Realizable k-ε模型和雷诺应力模型的湍流粘性计算结果较为接近,而RNGk-ε模型计算的湍流粘性较小.     

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。