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

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

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

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

大扩张角涡轮过渡段非定常流动数值仿真

为进一步减小涡轮过渡段流动损失,深入了解涡轮过渡段中的非定常损失机理,开展了大扩张角过渡段研究。在过渡段的非定常流动机理研究中,过渡段进口流场的最显著特点是:转子泄漏涡、通道涡和尾迹。采用数值方法对大扩张角涡轮过渡段进行3维非定常数值仿真。结果表明:支板尾缘部分的静压波动小于支板前缘部分的;高压涡轮静子尾迹被转子切割后进入转子通道中向下游传播并在过渡段内形成尾迹通道,尾迹在过渡段内的时空演化是过渡段内损失的主要来源;过渡段支板表面负荷分布发生明显的周期性变化,支板表面承受较强的非定常力,在过渡段设计中必须考虑。     

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

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

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

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

高速进气道出口快速撤锥再起动特性研究

为了探寻出口快速撤锥过程对高速进气道起动性能的影响,对二元高速进气道出口快速撤锥过程开展了非定常数值仿真研究。对比了进气道在一个喘振周期中从不同状态点处进行快速撤锥所表现出的起动性能的差异。研究表明,当来流马赫数低于进气道的加速自起动马赫数时,对于高速进气道因出口堵锥节流所引起的不起动流态(喘振),从不同喘振状态开始撤锥,进气道表现出不同的起动特性,而进气道最终能否顺利建立起动流态与快速撤锥时机的选取有关。     

关于变形网格“几何守恒律”概念的讨论

在国内外文献中应用包含有变形的动网格进行流动数值模拟时,要满足所谓的几何守恒律条件;通过对这一概念最早提出推导过程和近年来应用发展情况分析以后发现,逻辑上存在不适当之处.从流体力学基本理论出发,对包含有变形的动网格技术的理论基础进行讨论,证实几何守恒律是流体力学控制方程的伴随方程或退化方程;通过简单模型分析有限体积方法离散过程,发现计算过程中不满足几何守恒律所引起的非物理现象本质上是目前有限体积离散处理中微元体积计算方法不符合物理定律所致;因此在变形动网格计算方法中不存在必须要满足的几何守恒律.最后根据理论提出物理背景清晰的微元体积计算方法,数值验证可行.     

昆虫翅悬停飞行中的三维空间旋涡结构

昆虫(果蝇)悬停飞行中,翅膀按照特定的拍动方式往复运动,产生非定常高升力维持身体的平衡.研究昆虫高升力机理,需要探索拍动翼运动引发的三维空间非定常流场的特性,尤其是三维空间非定常涡的发展变化过程.本文将氢气泡流动显示技术应用于动态模型实验,定性的观察拍动翼前缘涡(LEV)的发展破裂过程.并利用数字体视粒子图像测速,DSPIV(Digital Steroscopic Particle Image Velocimetry)技术,测得了拍动翼运动瞬时相位和相位平均的三维空间流场信息(速度向量场、截面涡量场、空间涡量场,以及三维空间流线),揭示了拍动翼展向流动的存在,并结合定性和定量方法多角度说明了前缘涡沿展向发展到破裂的流动结构,并说明了侧缘涡与前缘涡的相互影响.测量结果表明:在雷诺数960的情况下,拍动翼运动至相位时,翼面上前缘涡在距翼根约60%展长的位置发生破裂;翼根至破裂点之间,展向流动稳定,指向翼梢;破裂点以后,展向流改变方向,指向翼根.     

SDM标模大迎角动态特性试验研究

介绍了CARDC高速所进行的SDM模型大迎角的动导数和大迎角大振幅动态试验情况与试验结果.SDM标模的试验结果表明:试验结果与国外文献数据具有较好的一致性.试验结果随马赫数和迎角的变化表现出明显的非线性,而减缩频率的影响并不显著.天平与测试系统的重复性精度较高.在试验条件下,减缩频率、振幅、平均迎角及来流肘数均对模型的动态气动特性有影响,其影响量值的大小与具体的试验条件有关.在迎角约为15°～75°范围内,试验曲线存在着迟滞环.     

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

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