基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

低速压气机静叶非定常Clocking效应数值模拟

采用三维粘性非定常数值模拟方法研究了重复级低速压气机内部非定常流动,其中包括静叶Clocking(时序)效应以及尾迹-边界层、尾迹-尾迹和尾迹-泄漏流的相互作用.对静叶时序效应的结果分析显示,当两级静叶相对周向位置改变时最大效率变化只有0.1%,表明低速压气机中静叶时序效应不明显.另外,非定常结果流动图画显示,在非定常环境下上游周期性尾迹的通过对下游叶片边界层、尾迹以及叶尖泄漏流的发展都有一定影响.     

对转螺旋桨流场气动干扰数值模拟

基于结构网格动态面搭接技术,通过求解三维非定常雷诺平均纳维-斯托克斯(RANS)方程,对单独螺旋桨(SRP)滑流流场进行了数值模拟,结果与风洞试验数据吻合良好,在此基础上深入研究了对转螺旋桨（CRP）的非定常滑流效应及前后桨之间的气动干扰现象。计算结果表明:所采用的非定常RANS方法能够良好刻画对转螺旋桨滑流流场的发展变化特征,适用于分析前后桨气动干扰特性。前后桨桨尖涡的复杂相互作用是导致对转螺旋桨气动干扰的直接原因,后桨桨尖涡对于前桨桨尖涡具有一定耗散作用。在气动干扰影响下,对转螺旋桨气动力会产生周期性波动,波动周期与螺旋桨桨叶数相关, 6×6对转螺旋桨在1个旋转周期内产生了12次波动。气流穿过对转螺旋桨会发生两次加速,因此在相同工况下,对转螺旋桨的拉力系数和功率系数相比于单独螺旋桨分别增大约1倍,效率提升约1.5%。      

摆线桨非定常动态失速的三维效应

对摆线桨三维全尺寸模型进行了非定常数值模拟。验证了两种嵌套网格方法的可靠性与适用性,针对二维/三维情况下摆线桨非定常涡流动特性及诱导分离特性开展了数值分析,重点研究了展向气动力随方位角的分布与全局流场的诱导速度分布,并分析了尾迹捕捉精度与涡量耗散特征。结果表明:受摆线桨切向速度、展向诱导速度等三维效应影响,两种计算结果的瞬态气动力极值差值达到44%。三维动态失速涡的产生、脱落与再附明显弱于二维翼型情况,这对于摆线桨非定常气动特性有较大影响。      

柔性扑翼非定常流场的数值计算方法

提出一种将Delaunay图映射网格变形技术和非结构嵌套网格方法结合使用的策略,解决网格变形和嵌套网格单独用于柔性扑翼流场计算时需要网格再生的问题。该方法为嵌套网格中的每个嵌于背景网格的贴体非结构网格生成Delaunay背景图;每个时间步,根据扑翼的运动和变形规律移动背景图,再根据网格点和背景图的映射关系移动网格点,之后自动完成嵌套边界的定义和插值关系的建立。为方便嵌套关系的建立,嵌套网格进行分层管理。也研究了一种内存消耗少、效率较高的搜索算法,以及格心格式和格点格式统一的边界拓宽算法。非定常可压缩Navier-Stokes方程在非结构的动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代求解。几个扑翼算例的结果表明,该方法充分利用了Delaunay图映射网格变形方法的高效率,同时也发挥了嵌套网格处理大幅运动的优势;用于既有整体大幅扑动又有局部小变形的柔性扑翼流场计算,可取得令人满意的精度和效率。     

燃气二次喷射推力矢量控制系统启动过程非定常效应研究

采用数值模拟方法对燃气二次喷射推力矢量控制系统启动过程中非定常效应进行了研究.给出了详细的喷流瞬态干扰流场结构, 以及喷流瞬态干扰区随时间改变的流场细节特性,分析了流动参数和发动机内弹道性能的动态变化过程.研究表明,主/次流干扰流场内流动参数的分布及变化呈现复杂的多样性,启动过程中所呈现出的动态特性与系统稳定工作时有着明显的差异,侧向控制力、燃烧室平衡压强和发动机推力均有显著变化.在飞行器设计中,必须从安全控制的角度出发,考虑这一非定常效应的影响.     

非稳态流体网络方法在发动机空气冷却系统中的应用

将流体网络法与非稳态模拟方法相结合应用到发动机空气冷却系统沿程温度、压力、流量的非稳态特性预测中.本方法将空气系统简化为节点和元件组成的网络,节点和元件处建立守恒方程,针对传统解法的计算稳定性问题引入了压力修正和积分方法相结合的求解方法.程序应用到可参考对照的模型中得到验证.验证后的程序被应用到发动机轴部冷却通道空气系统非稳态特性预测中.      

上游尾迹与涡轮转子泄漏流相互作用数值模拟

叶轮机内部流动本质上是周期性非定常的,研究涡轮转子叶尖区域的非定常相互作用机理,对提高小展弦比高负荷涡轮性能具有重要意义.利用数值模拟方法研究了上游静子尾迹与涡轮转子叶尖泄漏流的非定常相互作用,分析了定常结果、时间平均结果以及瞬时时刻结果的流动图画.结果表明:上游静子尾迹与涡轮转子尖区二次流的相互作用能明显影响泄漏涡和机匣通道涡的时空演化规律,从而改变转子尖区的损失分布.上游尾迹在转子通道中传播时,诱导泄漏涡和通道涡区域出现周期性的扰动涡对,扰动涡对沿着泄漏涡和通道涡的轨迹向下游运动,使得转子尖区二次流结构呈现周期性变化.      

不同轴向间距下涡轮级内非定常流场的数值研究

针对不同轴向间距下涡轮级内非定常流动进行了数值研究.讨论了非定常条件下轴向间距变化对非定常性及流场分布的影响,探讨了转子和静子叶片之间的非定常相干的机理.结果表明,轴向间距L对转子内流动的非定常性有着强烈的影响.随着轴向间距的减小,静叶尾迹的强度及影响范围增加,静叶尾迹对动叶的周期性非定常影响还是主要出现在前缘附近.由于静叶尾迹的恢复效应与来流尾迹强度成正比,50%轴向间距下动叶出口流场的湍动能在动叶通道中央区域最小.      

湍流边界层厚度对三维空腔流动的影响

采用脱体涡模拟(DES)方法开展了不同湍流边界层厚度(TTBL)下的三维空腔非定常流动数值计算。空腔长、宽、深比例为5:1:1,来流马赫数为0.85,雷诺数为13.47×10⁶ m^-1,各工况湍流边界层厚度比值为1:2:4:8。研究结果表明,湍流边界层厚度对自由剪切层的发展、空腔底部静态压力分布、脉动压力及空腔流动类型均有重要影响,且随着边界层厚度的增大,下游剪切层覆盖的范围会增大,但是剪切层增长率降低;空腔前后静态压力压差减小、压力梯度下降;腔内局部测点的脉动压力声压级下降,各阶声压峰值频率向低频方向偏移;空腔流动类型往开式流动方向转换。