离心压气机中进口导叶尾迹的形态、输运及与叶轮流场的相互影响研究

对跨声速离心压气机中进口导叶与叶轮的相互影响作用进行了数值模拟研究,研究模型包括三种几何间距模型及对同一间距几何模型使用了两种不同的转静子交界面位置设定。计算结果表明：转静子交界面位置的设定对于多排叶片混合平面法定常计算性能会产生较大影响,5%叶轮轴向长度的变化就可造成1%的效率、2%的压比差别;进口导叶与叶轮的叶排间距大小对离心压气机级气动性能影响很小,但间距较小时会造成流场内较大的压力波动。在近叶排间距时,导叶切割叶轮激波会在导叶压力面侧产生一个高损失区,并沿导叶表面向上游传播;该高损失区由激波压力波造成。此外,激波存在与否,影响到两排叶片流场的相互影响作用强度：激波使叶轮对导叶流场的影响几倍甚至十几倍大于导叶对叶轮流场的影响;激波的影响作用使导叶压力面、吸力面侧压力波动由不同的原因造成。最后,在近叶排间距时,初始进入叶轮通道内且靠近主叶片压力面的尾迹相对靠近吸力面的尾迹片会滞后;完全进入叶轮通道后,近压力面侧尾迹的滞后更显著。     

高超声速飞行器表面温度分布与气动热耦合数值研究

针对高超声速飞行器热防护设计中的高温气体非平衡效应问题和气动热环境精确预测问题,基于流场的非平衡Navier-Stokes方程、表面的能量守恒方程和内部的热传导方程,考虑流场的非平衡效应、表面的热辐射效应、催化效应和烧蚀效应以及热防护层内部的热传导效应,建立了初步的表面温度分布与气动热的耦合计算方法,完善了高超声速飞行器气动物理流场计算软件(AEROPH_Flow)。在表面材料为碳-碳(C-C)条件下,对飞行高度为65km和飞行速度为8,10km/s的半球以及飞行高度为50km和飞行速度为8km/s的球锥模型,开展了表面温度分布与气动热的耦合计算,验证了计算方法和计算软件,分析了表面温度分布对气动热环境的影响。研究结果表明:表面温度分布对气动热的计算结果有较大影响,在气动热环境的预测中,不仅要考虑热化学非平衡效应和表面催化效应的影响,还要考虑表面温度分布的影响,最好是采用表面温度分布与气动热耦合计算的方法,以减小表面温度分布对气动热计算结果的影响。为此,需要发展完善非平衡流场/表面催化和烧蚀/热传导温度场(气/表/固)的计算模型、耦合求解技术和计算软件,实现对高超声速飞行器的真实飞行条件下高温气体非平衡效应和气动热环境的精确模拟。     

高焓激波风洞喷管流场非平衡特性研究

高焓激波风洞是开展高超声速流动研究的重要地面模拟设备,但其产生的高焓气流在喷管中的膨胀过程是一种典型的热化学非平衡流动,试验段特征参数通过直接实验测量难以完全确定。本文通过求解耦合双温度模型的轴对称Navier-Stokes方程,研究了高焓激波风洞中典型状态下气流的热化学非平衡流动特性,分析了焓值对非平衡特性的影响规律。结果表明,喷管出口自由流均匀区域达到出口截面直径的75%以上,能够为实验提供足够的空间;喷管出口自由流处于热化学非平衡状态,在喷管喉道后约1/5喷管长度处气流即已处于冻结流状态,组分浓度和振动温度随气流流动基本不变;焓值在8.4MJ/kg~19.5MJ/kg之间变化时,非平衡程度随着焓值的增加而增强,但是低焓值时非平衡程度的增强更加剧烈。     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

机匣处理作用下高负荷轴流压气机转/静子匹配设计研究

为了探索机匣处理作用下转/静子的轴向匹配方法以进一步提高压气机级的失速裕度,研究了静子的叶型安装角及"弯"、"掠"规律对压气机性能的影响,针对机匣处理与优化静子的组合结构进行了非定常数值模拟,阐述了该结构的扩稳机理以及压气机新的失速机制。研究结果表明,在机匣处理作用下,静子成为压气机失速的触发因素,通过对静子叶型安装角及"弯"、"掠"规律的优化均可进一步提高压气机级的失速裕度,其中改变静子"弯"型对压气机级失速裕度的改善最大。组合应用机匣处理与尖部反弯根部正弯静子后,压气机效率基本不变,失速裕度提升了80.2%,较单独使用机匣处理提升30.9%。在该组合结构作用下,压气机的失速由静子触发,静子叶根吸力面在激波作用下发生附面层分离,且与轮毂表面附面层相互作用形成角区涡,接近失速边界时,静子叶根形成"前缘溢流,尾缘反流"现象,造成静子通道的大范围堵塞,诱发压气机失速。压气机级的扩稳应充分考虑机匣处理的影响,对静子进行优化设计。     

预旋系统稳态和非稳态计算对比研究

针对预旋系统由于旋转引起流场和温度场参数周期性瞬态变化的问题，分别采用滑移网格的瞬态法和固定转子相位的稳态法进行数值求解，并在稳态计算中通过改变转子相位来近似模拟非稳态问题的时空变化特性。通过稳态空间平均结果与非稳态时均结果的对比，以期为预旋系统非稳态问题的低成本求解提供方法依据。结果表明：稳态计算结果与非稳态计算结果相比，周期性波动频率一致，接受孔进口处，稳态计算的压力波动幅度小39%左右，温度波动幅度小15%左右；多个相位的稳态空间平均结果与非稳态时均结果相比，压力高0.2%，温度低0.1%；采用多个计算周期数与单个周期的非稳态计算结果差异微小。当采用喷嘴出口气流中心正对接受孔迎风面前缘的转子相位时，稳态计算结果与非稳态计算时均结果最接近。     

核热推进地面试验技术研究

核热推进具有比冲高、推力大等特点,在载人深空探测和星际货运任务上具有广阔应用前景。核热推进技术的研发需要进行大量地面试验。首先回顾了美国与俄罗斯的核热推进地面试验技术的发展,对地面试验技术进行分类总结。然后基于一种小型核火箭方案,研究了燃料元件非核试验、燃料元件辐照考验试验和带核整机地面试验等关键地面试验技术,并提出了初步试验方案。最后对我国核热推进地面试验的发展提出了一些建议。     

基于自适应非结构嵌套网格的旋翼流场模拟

针对直升机旋翼CFD仿真的复杂性,提出了改进的适合于格心格式求解器的非结构嵌套网格算法。采用自适应网格技术在旋翼流场仿真的整个过程中进行网格的自适应加密和疏化操作,以更好地捕捉桨尖涡等流动细节。对于频繁的自适应过程中产生的大量重复点和无用点,采用了高效的交替数字树算法(Alternating digital tree,ADT)和标记-删除-移动算法(Mark,delete,move,MDM)进行删除,节约了不必要的存储。针对格心格式的求解器,采用了基于梯度的网格间插值方式,简化了网格间数值传递的复杂性,同时不降低求解器的精度。对CaradonnaTung旋翼悬停算例和HLISHAPE 7A旋翼悬停算例进行了模拟验证,计算值与实验值吻合,表明本文建立的方法具有良好的鲁棒性和有效性。最后,与未采用自适应时求解器对桨尖涡的捕捉效果进行了对比,结果表明本文所采用的方法可以明显地提高求解器对桨尖涡的捕捉。     

考虑非理想噪声干扰的直升机飞行动力学模型辨识方法

利用MATLAB/SIMULINK平台搭建了直升机仿真实验系统，并以UH-60直升机为对象实现了考虑非理想噪声干扰的仿真飞行实验。建立了非理想噪声数学模型并基于增广最小二乘法设计 了一套同时考虑模型参数和噪声参数的综合辨识算法，在此基础上，利用仿真实验数据实现了UH-60直升机纵向飞行动力学模型的辨识与验证。最后，通过与普通最小二乘法的对比验证了本文方法的优越性。 