湍流模型对叶片表面换热计算的影响

为了比较湍流模型对涡轮叶片外换热计算结果的影响,采用五种湍流模型及两种壁面函数计算了NASA-MarkII导向叶片表面换热并与实验数据进行了对比,结果表明:在相同边界条件下不同湍流模型的计算结果有很明显的差别,即使是同一种湍流模型,如果采用不同的壁面处理函数其计算结果也是有很大差别的;某些湍流模型的计算值只是在某个区域较为理想,还不能找到在整个叶片表面计算结果与实验数据较为接近的湍流模型.在尚无普适性较好的湍流模型的情况下,研究在不同计算域采用不同湍流模型和壁面处理函数的计算技术,是一种较好的可行的方案.      

温度边界层对壁湍流多尺度相干结构的影响

用热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量了风洞中壁面加热平板湍流边界层不同法向位置流向速度分量的时间序列信号, 与常温情况相比, 热对流加速了壁湍流中流体的动量、质量和能量的交换, 加快了湍流边界层的发展;用子波分析方法研究了壁面加热对壁湍流多尺度相干结构喷射和扫掠条件相位平均波形的影响, 显示热对流明显增大了缓冲层相干结构扫掠过程和喷射过程的强度, 是增强壁湍流中流体的动量、质量和能量交换的根本原因.      

高超声速飞行器-进气道一体化热流数值计算

采用CFD(计算流体动力学)技术, 开展了飞行器前体/发动机一体化气动热环境分析.对层流区、转捩区和湍流区分别采用计算模型, 在湍流区利用压缩性修正的SSGZ-Jk-ε湍流模型, 在转捩区引入代数型转捩因子模型描述边界层由层流逐渐过渡为完全湍流的流动过程.计算了前体和内通道的表面热流, 并与实验结果进行了对比.结果表明所采用的计算方法可以较好地预测前体及发动机内通道热流率, 流动状态、几何结构及激波入射对热流值影响较大.      

PSE方法在后部加载叶栅绕流稳定性分析及转捩预测中的应用

采用数值求解与实验测量相结合的方法, 分析后加载叶栅绕流的稳定性并预测转捩.推导出正交曲线坐标系下三维扰动波的抛物化稳定性方程, 在风洞实验中测量叶栅叶片表面的静压分布, 并结合测量结果对设计叶栅的流动稳定性进行特征值分析.分析结果表明:(1)后加载叶型能有效推迟转捩, 提高叶栅的变冲角特性;(2)当冲角变化时叶顶流动稳定性的变化最显著, 负冲角流动稳定, 正冲角流动最不稳定.      

高超声速钝楔边界层转捩大涡模拟

 以五阶迎风和八阶对称格式混合差分格式求解三维可压缩滤波Navier-Stokes方程,对来流Mach数为6.0、半锥角5°的高超声速空间发展钝楔边界层转捩至完全湍流进行了大涡模拟。时间推进采用紧致存储三阶Runge-Kutta方法,亚格子尺度模型为Driest因子修正的Smagorinsky涡黏性模型。通过定常流场入口边界附近吹／吸引入不稳定扰动斜波的方法数值模拟得到了层流失稳转捩直至完全湍流的空间发展全过程。对扰动的线性、非线性增长以及湍流斑的形成和发展进行了分析,给出了转捩及完全湍流下的速度相关量统计并与实验、DNS结果进行了对比分析,计算结果与理论及实验吻合。     

A Nonlinear Sub-grid Scale Model for Compressible Turbulent Flow

本文采用直接过滤的Navier-Stokes方程组作为可压缩湍流大涡模拟控制方程组。方程组中因过滤产生的高阶相关项用Taylor级数展开近似,但仅保留级数的一阶导数项。这样产生的误差相当于丢失了模型的部分耗散作用,本文用一种动力学非线性亚格子模型来补偿丢失的耗散影响。本文根据Caylay-Hamilton定理导出了一个非线性亚格子模型,模型中包含的系数由动力学模式确定。与传统的动力学Smagorinsky模式相比,这种动力学非线性模型(DNM)稳定性更好,且不需要在统计均匀方向进行统计平均来计算模型系数,因此减少了计算开销。本文用这种非线性亚格子模型对绕双椭球高超声速湍流进行了模拟,并将所得结果与实验值、计算值及理论值进行了对比。结果表明,本文模型可以有效地模拟可压缩湍流流场。     

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