高阶精度方法下的湍流生成项对低速流动数值模拟的影响研究

基于雷诺平均的Navier-Stokes方程和结构网格技术,采用五阶空间离散精度的WCNS格式,开展了SST两方程模型不同湍流生成项对低速流动数值模拟的计算分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR7301两段翼型和Trap Wing高升力构型,研究内容主要包括不同湍流生成项对收敛历程、边界层湍流粘性系数分布、边界层速度分布、压力系数分布、气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于低速二维流动,不同湍流生成项对收敛历程有比较明显的影响,对附面层湍流粘性系数分布和速度型影响不明显,不同湍流生成项主要影响主翼前缘的吸力峰值,进而影响升力系数和压差阻力系统;对于低速三维流动,不同湍流生成项对低速流动的收敛特性影响不明显,对翼梢涡的模拟精度有比较明显的影响,进而影响翼梢站位的压力分布和总体气动特性。     

无源腔结构对大膨胀比单膨胀斜面喷管的影响

开展了基于无源腔的单膨胀斜面喷管(SERN)流动控制方法研究,采用标准k-ε两方程湍流模型,通过求解Reynolds平均Navier-Stokes方程,数值研究了过膨胀状态下无源腔结构对SERN内部流场和性能的影响.计算结果表明:与原型SERN相比,落压比为8时,无源腔结构的存在促进SERN上膨胀斜面流动的分离,改变SERN内部流场结构,影响SERN上膨胀斜面压力分布,从而使SERN的轴向推力系数提高1.75%,同时不会带来低头力矩的恶化.在落压比为3~13范围内,无源腔结构的存在促进了无源腔附近附面层流动分离,减少了SERN内激波串数目,并使SERN上膨胀斜面压力峰值增加,提高SERN轴向推力系数.      

低雷诺数下压气机的二次流旋涡结构

为了探究高空低雷诺数条件下跨声速压气机的流动规律,对NASA Rotor37进行单通道数值模拟,探索其在低雷诺数进气条件下二次流的旋涡结构.研究发现:马蹄涡压力面分支诱发压力面角区诱导涡,壁角涡形成了顺流和逆流的两段式结构,脱落涡由叶根角区发展起来后不断从尾缘脱落,泄漏涡近失速点仅局部破裂不是失稳触发的主要原因.通道中的激波系诱发了吸力面和压力面的两个径向涡,压力面径向涡构成闭合的气泡式分离,吸力面径向涡在叶顶的破碎诱导产生分离涡,触发了低雷诺数下压气机的失稳.流场旋涡结构由马蹄涡、壁角涡、径向涡、泄漏涡、分离涡、脱落涡6个大尺度旋涡以及其他小尺度旋涡组成.      

航空发动机主燃烧室燃油总管流动换热特性研究

为优化发动机燃油管路设计并预防燃油结焦,对发动机典型工况下燃油总管和分管内部流动进行了数值研究,给出了速度和压力分布情况。对由于热气流的加热作用引起的燃油温升效应进行了传热计算。对总管及分管进行了热防护设计,并比较了不同热防护措施的隔热效果。结果表明:在计算工况下,油管出口质量流量满足设计要求;内外混合隔热方式可有效降低燃油温升,减小燃油过热引起结焦的可能性。提出了改善管内流动与换热特性的措施,例如增加空气隔热层,为恶劣热负荷条件下燃油管路热防护设计提供理论基础。     

涡轮叶片叶顶间隙变化减敏研究

采用数值方法联合标准k-ω两方程湍流模型求解雷诺平均Navier-Stokes方程组,研究了平顶和分别采用被动、主动以及复合间隙控制方法的凹槽顶、平顶喷气、凹槽顶喷气4种不同叶顶结构下涡轮性能对叶顶间隙变化的敏感性,分析了间隙变化对间隙流场结构及损失的影响,研究了涡轮性能随叶顶间隙变化的规律,并对采用不同间隙控制方法的叶顶结构所带来的效率收益对攻角变化敏感性进行了评估.结果表明:在不同间隙控制方法下叶顶泄漏涡强度随间隙的增加速度变缓,尤其在复合间隙控制方法下更为明显,从而使得在复合间隙控制方法下涡轮性能对间隙变化的敏感度最低;仅仅涡轮叶片前缘区域泄漏流动对攻角变化比较敏感,使得凹槽顶结构具有最佳的攻角变化适应性.      

纳秒脉冲等离子体分离流控制频率优化及涡运动过程分析

将纳秒脉冲驱动的介质阻挡放电等离子体激励器应用到NASA SC（2）-0712翼型上,在迎角分别为15°和20°时,开展了在不同雷诺数下的分离流动控制研究。通过模型表面静压测量,得到了不同激励频率下的分离流动控制效果。对翼型表面压力进行分布积分,得到了在不同雷诺数和激励频率下的升力系数,表明分离流的控制效果有一个较宽的激励频率范围,只要激励频率落在相应的频带范围内,均能实现有效的分离抑制。流动显示结果表明,分离流的控制在瞬时表现为放电后可形成大尺度旋涡拟序结构。旋涡的周期性产生、运动和演化造成了分离剪切流动的动态变化过程,从而促进了高/低速气流的动态掺混。     

非结构变形网格和离散几何守恒律

数值模拟流固耦合问题或多体分离问题的非定常流动时,常采用基于任意拉格朗日-欧拉（ALE）方程的有限体积法,涉及到变形网格和离散几何守恒律。在对变形网格算法进行综述时,按照构造思想分为物理比拟、椭圆光顺、插值、运动子网格（MSA）及其混合法共5类,分别介绍了基本原理、研究现状和适用范围,通过算例比较表明,径向基函数（RBFs）和运动子网格相结合的混合方法既有很好的变形能力,也有较高的计算效率,值得进一步发展和推广。在介绍了离散几何守恒律（DGCL）概念之后,采用二维几何模型进行分析,指出其机理是离散过程中体积增量与网格面元扫过的体积不相等造成的,把目前国内外应用的算法分为面积修正法、给定速度的面积修正法、速度修正法和体积修正法共4类,对其应用范围和存在的问题进行讨论,认为提出的体积修正算法既可以保证流固界面条件,也可以用于时间多层格式。     

基于POD和DMD方法的跨声速抖振模态分析

跨声速抖振现象是由于非定常跨声速流动中激波的自激振荡而引起的结构强迫振荡,这种现象在跨声速飞行器中普遍存在,对飞机的结构强度和疲劳寿命有不利影响。基于模态分解的分析方法是进一步发展抖振控制手段的有效工具。本文通过两类典型模态分析方法（本征正交分解（POD）和动态模态分解（DMD））对OAT15A翼型的跨声速抖振现象进行分析,通过对模态频率、翼面压力分布、流场重构误差等方面的研究,将两种模态分解方法进行对比。发现基于频率特征的DMD方法能够准确捕捉抖振的临界稳定特征和抖振主频的典型模态,同时能够更准确反映流场变量在激波间断附近随时间的变化过程;而POD方法尽管在流场重构时具有较小的总体误差,但对激波附近压强随时间的变化历程拟合较差。     

存在烧蚀的高超声速压缩拐角流动数值研究

耦合求解热化学非平衡流控制方程和烧蚀壁面边界条件,进行存在石墨烧蚀的压缩拐角流场数值模拟。流场化学反应采用16组元(N2,O 2,NO,N,O,NO +,N+2,O +,N+,CO,CO 2,C,C2,C3,CN,e-)29个反应的非平衡模型,热力非平衡的双温度模型下,不同反应采用不同控制温度。石墨材料表面反应包括碳的氧化反应、碳催化的O 原子复合反应和碳的升华反应。对15°、18°、24°压缩拐角模型,在自由流 Ma =10~30,总焓值6~55 MJ/kg 范围,分别进行无烧蚀的壁面催化与非催化条件和石墨烧蚀条件下的流场计算,分析各类条件下的流场结构、流动分离特性以及流场热化学参数分布特点,研究壁面条件对流动特性的影响。结果表明：流动分离可能性和分离区范围随着压缩拐角斜面倾角增大而增大,随来流马赫数增大而减小;相对于低壁温条件,无烧蚀的辐射平衡壁温和壁面烧蚀条件下流动分离区增大,斜面上压力、摩阻和热流峰值点也有所后移。     

激励器结构对三电极等离子体高能合成射流流场及其冲量特性的影响

等离子体激励器以其结构简单、响应速度快、环境适应性强等优势,已成为主动流动控制技术和流体力学研究的前沿与热点。相比于传统两电极激励器,三电极等离子体高能合成射流激励器具有更高的能量效率,形成射流冲量更大,有望成为新型快响应直接力产生装置。为揭示激励器结构对射流流场和冲量特性的影响规律,进而优化激励器结构参数,利用电参数测量装置、高速阴影系统及自主设计的单丝扭摆式微冲量测量系统对不同射流孔径、腔体体积和电极间距的三电极激励器放电特性、射流流场及其冲量进行了实验研究。为对比激励器在不同工况条件下的工作特性,定义无量纲能量沉积ε和无量纲射流冲量 I *,并分析了激励器结构参数对ε和 I *的影响。结果表明对于给定无量纲能量沉积ε,激励器存在最优射流孔径;激励器无量纲能量沉积ε和无量纲射流冲量I *随腔体体积增加而减小,随激励器电极间距增加而增加;射流强度及其流场影响区域随腔体体积增加而减小,随激励器电极间距增加而增加。对比不同腔体体积和电极间距工况条件下 I *随ε的变化可知,为设计具有较好射流冲量水平的激励器,在相同无量纲能量沉积ε条件下,应尽量增大激励器无量纲射流冲量 I *。当设计激励器无量纲能量沉积ε小于初始工况时,应增大初始工况激励器腔体体积使无量纲能量沉积ε降低至设计值;当设计激励器无量纲能量沉积ε大于初始工况时,应增大初始工况激励器电极间距使无量纲能量沉积ε增加至设计值,使设计激励器具有较好的射流冲量水平。