基于MSFLE与Liutex的涡合并实验研究

应用具有拉格朗日性质的运动单帧长曝光(MSFLE)图像测量方法,对矩形管内湍流边界层涡结构合并现象开展了实验研究。研究了雷诺数Reθ在97到194之间管道湍流边界层流向-法向平面内涡合并现象。在实验中,测量系统以与涡运动速度相近的速度匀速移动,采用连续长曝光记录示踪粒子运动轨迹,以捕捉涡随时空的演变过程,骨架提取图像处理算法获得示踪粒子速度,并应用Liutex物理量表征涡的旋转强度。研究表明:MSFLE方法测量装置简单,对实验条件要求低,可从拉格朗日视角直观测量湍流边界层涡结构及周围流场的时空演变过程,MSFLE图像测量方法与Liutex涡识别算法结合可以很好地应用于湍流边界层涡结构的可视化与量化。矩形管内湍流边界层涡合并的条件是两个涡相邻、强度和尺寸基本相同且为同向旋转涡,合并中两个涡的强度呈反向变化,合并生成的新涡的强度和尺寸基本为初始合并时两个涡之和,且旋转方向与两个涡同向。     

基于释热率分析的钝体贫油熄火过程分析

针对锥形钝体稳定的甲烷-空气预混湍流火焰复杂的熄火过程,采用大涡模拟(LES)与输运方程概率密度函数(TPDF)湍流燃烧模型相结合的模拟方法,研究远离熄火、近熄火及熄火点的火焰和释热率(HRR)数值变化情况,定量分析熄火判据。结果表明:冷态速度模拟结果和实验的相对均方根误差在10%以内,热态误差在20%以内;释热率是否出现在OH和CH₂O重叠的区域,是判断熄火的一个重要参数;远离熄火时,释热率高的区域主要沿内侧剪切层出现;近熄火工况下,释热率在流向轴上闭合,回流区下游也出现较大的释热率;熄火点工况下,较大释热率的区域在回流区下游和上游均有出现;模拟预测的熄火情况和实验PLIF结果一致;平均释热率可作为判断熄火的定量依据,即当钝体后方0.2d处内侧剪切层平均释热率与回流区平均释热率的比值小于4时,发生熄火。      

航空弧齿锥齿轮副风阻功率损失分析与优化

基于CFD理论,利用Fluent求解软件,借助超级计算机强大的并行运算能力对航空弧齿锥齿轮副风阻功率损失进行仿真计算。采用局部综合法建立弧齿锥齿轮副三维模型,选用RNG k-ε湍流模型,考虑平均流动中的旋流流动情况,与标准k-ε模型相比,RNG通过修正湍流黏度并很好地处理了高应变率以及流线弯曲程度较大的流动。齿轮边界运动通过UDF(user-defined functions)函数驱动,同时采用动网格模拟流场形状由于边界运动而随时间改变问题。最后得出无挡风罩和不同挡风罩配置下的齿轮副风阻功率损失,证实了合理安装挡风罩能够有效降低齿轮风阻损失,并分析多组仿真实验间的减速器内流场压力、速度、湍流动能云图变化,得出了最优化的挡风罩配置,以求最小化风阻功率损失,文中减阻效果最好的挡风罩能降低55.3%的齿轮风阻损失,此时挡风罩间隙为1 mm,为工程实际应用挡风罩的设计提供了参考。     

应用三维粒子动态分析仪测量三分支联接流场

应用三维粒子动态分析仪 (3D PDA)测得了矩形截面三分支联接的三维速度分布。在测量中 ,经过理论分析与对比试验 ,选择了蚊香烟雾作为散射粒子。实验结果表明 ,在只有支管进气的情况下 ,在总管的封闭端和支管下游总管壁面处产生了回流 ;以及气流从支管进入总管后由于截面扩张在总管横截面上所诱发的二次流 ,并且二次流现象只发生在支管下游附近的总管横截面上。     

SAS与PANS模型在圆柱绕流中的应用比较

分别运用SST(Shear Stress Transport)、SST-SAS(Scale-Adaptive Simulation)、两种变fk(模化湍动能的比例)函数的SST-PANS(Partially Averaged Navier-Stokes)湍流模型对Re=3900的圆柱绕流进行了数值研究,重点从湍流结构捕捉、气动力计算、涡黏性控制等方面,比较了SAS与PANS两类RANS/LES混合模型的计算能力,并通过不同网格计算分析了模型的网格敏感性。数值结果表明:SAS及两种变fk方法的PANS模型均具有求解小尺度涡运动的能力,并能较好地反映出绕流尾迹的三维非定常特性,同时PANS模型能捕捉到更多的非定常结构;SAS模型中自适应尺度Lvk立足于当地流动,对网格依赖较小,计算的湍动粘度分布更合理,能够更好地计算剪切层及回流区;两种PANS模型网格独立性较差,出现了雷诺应力不足的现象;类DES可变fk函数构造相对简单,所得fk分布更准确,使用tanh函数计算的尾迹区fk值偏低,对流场调控能力稍差。     

分区和有限容积法计算长尾喷管中的湍流流动

为了给长尾喷管热防护结构设计提供热边界条件，从轴对称二维雷诺平均N－S方程出发，采用分区算法和有限容积法计算了长尾喷管中的湍流流动，分析了流场的结构特性，提出了长尾喷管中燃气温度，压强和马赫数分布的某些规律以及燃气温度沿壁面的分布，计算和分析结果提出了长尾喷管热防护结构设计应注意的事项。     

环型超声速空气引射器零二次流流场数值研究

采用二维轴对称雷诺平均方程和标准κ-ε双方程湍流模型，数值研究了环型超声速引散器零二次流的流场结构及盲腔压强的变化，空间上采用二阶迎风格式进行耦合求解，时间上采用显示的Runge-Kutta方法进行迭代推进。结果表明，引射器几何参数不变的情况下，启动后的盲腔压强与引射气流总压之比为一常数；喷管马赫数不变情况下，喷管出口面积与混合室入口面积比越小，盲腔压强越低，扩压器性能越好，启动要求的总压越低，对超声速空气引射器的设计具有指导意义。     

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

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

旋涡运动与湍流模型

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

风对建筑绕流流动的数值方法

本文利用k-ε湍流模型及SIMPLE方法,在二维条件下,就风对建筑绕流流动进行了数值模拟,建筑及地面的边界条件采用以低Re数k-ε湍流模型中导出的壁函数处理近壁湍流。自由边界采用外推格式。对于计算域中的建筑物则利用控制粘性系数的方法处理。并对计算结果和风洞实验进行了比较。 为了研究建筑物之间的气流及压力分布,本文还对并列两栋建筑的不同情况进行了数值计算和分析。