计算气动光学研究进展

计算气动光学是采用数值计算的方法研究高速空气动力流场对光波传输和光学成像影响的一门交叉学科,在红外成像制导设计等工程应用中发挥着重要的作用。总结了计算气动光学研究的三种基本计算方法,包括基于CFD简化方法的光学统计估算、基于RANS流场的计算和基于LES/DNS瞬态流场的计算,对这三种计算方法的基本理论、计算过程进行了介绍,综述了新近的一些计算气动光学研究进展,最后指出了未来的研究重点和方向。     

超高速干气密封扰流效应及抑扰机制

干气密封在高速时优异的动压性能使其应用范围从传统的压缩机、离心机等中高速设备逐渐扩大到航空发动机、（微型）燃气轮机等超高速设备中。基于实际超高速工况特点,对转速范围为10 000~120 000 r/min时的干气密封性能进行了系统性仿真计算,结果发现：在一定几何参数和工况参数下,类似于气浮轴承的微振动现象,干气密封会出现疑似受气体压力波动流影响的开启力、泄漏量与转速非正相关变化的扰流现象,尤其在高压、大膜厚、小槽深时的扰流效应愈加显著;在转速持续增大过程中,干气密封微尺度流场会出现二次拐点现象,且一次拐点发生转速与设计参数有关,而二次拐点发生转速基本约为90 000 r/min。同时结合导流织构的设计思路,进一步研究了超高速下干气密封槽底导流织构的驱动导流效应,结果表明：加设导流织构后,承载效果明显提高,拐点发生工况延后且压力波动区域被压缩。表明导流织构具有良好的抑制扰流、维持开启力与转速持续正相关的作用,在此基础上,进一步阐释了导流织构的抑扰机制,以期为突破干气密封在超高速工况下的应用壁垒提供新思路。     

一个新的可压缩性修正的k-ε模型

考虑结构可压缩性修正的影响,发展了一个同时考虑结构可压缩性修正和膨胀可压缩性修正的k-ε湍流模型,新模型包括Chang可实现性、Heinz湍流动能产生项以及Sarkar可压缩性三部分修正.新模型扩宽了以往发展的可压缩性修正模型的适用范围,适用于高超声速(M>5)复杂湍流流动中.通过对多个复杂超声速横侧射流工况的计算,验证了新模型的预测效果.与实验结果相比表明,几个工况下新模型的预测精度都显著高于标准k-ε模型.流体分离强度越大,新模型的修正效果越显著.与标准k-ε模型相比,新模型计算结果与实验更加接近.     

大气紊流环境下的纵向飞行参数响应分析

介绍了大气紊流环境下飞机纵向飞行参数的计算与匹配方法.由于紊流的随机性特征,飞行参数响应分析采用时域与频域相结合的方式,对某型飞机紊流环境下的纵向飞行参数进行了计算分析.采用所介绍的方法,可以获得比较合理的飞行情况来指导型号设计工作.     

桥梁断面静风荷载的格子Boltzmann方法数值计算

通过引入反映湍流涡粘性的湍流松弛,得到了模拟高雷诺数湍流的BGK方程.在速度相空间、物理空间和时间上对BGK方程进行离散得到了三维十九速离散速度模型;结合分区计算技术,设计了格子Boltzmann并行算法;根据亚格子Smagorinsky模型,提出了直接从粒子分布函数计算湍流松弛时间的方法.用开发的并行计算程序对分体双箱截面和闭口箱梁截面的静风荷载进行了数值识别,得到的静力三分力系数和流场压力分布与风洞试验结果及CFD宏观方法计算结果吻合,并从表面压力分布入手分析了两种桥梁截面的绕流特点.     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

基于自然激励技术的颤振边界预测

为了预测紊流激励条件下机翼的颤振边界,基于自然激励技术提取紊流响应的自由衰减信号,采用矩阵束方法识别模态参数,最后通过Z-W(Zimmerman-Weissenburger)方法计算稳定性判据,拟合判据变化曲线并外推颤振边界.对平板机翼模型进行了数值仿真分析,对单独机翼模型风洞颤振试验数据进行了计算.结果表明:采用自然激励技术与矩阵束方法能够较准确地识别紊流激励响应的模态参数,频率识别误差小于6％,阻尼比识别误差小于30％,结合Z-W方法能够在较低风速较早地预测颤振边界,有助于提高试验的安全性.     

吸气式高超声速飞行器机体推进一体化技术研究进展

吸气式高超声速一体化飞行器最显著的特点是子系统之间的耦合较其他类型飞行器更加强烈,这使得其设计具有挑战性。所有的子系统之间部件相互干涉,包括:气动、推进、控制、结构、装载和热防护等,特别是机体与超燃冲压发动机之间的耦合最为突出。飞行器的前体和后体下壁面既是主要的气动型面,又是超燃冲压发动机进气道外压缩型面和尾喷管的膨胀型面,在产生推力的同时也产生升力和俯仰力矩。机体与发动机的强耦合作用对飞行器的推力、升力、阻力、俯仰力矩、气动加热、机身冷却、稳定性和控制特性有直接的影响。本文介绍了国内外机体推进一体化技术的研究进展,重点介绍了中国空气动力研究与发展中心(CARDC)的相关研究工作,包括:密切曲锥曲面乘波进气道和基于双激波轴对称基准流场内转式进气道设计方法、独创的大尺度脉冲式燃烧加热风洞一体化飞行器带动力试验技术和高超声速内外流耦合数值模拟技术等。对高速飞行中激波边界层相互干扰、流动分离机理、可压缩湍流转捩及其控制、超燃冲压发动机燃烧流动机理等相关基础问题也进行了研究,强调了对高效高精度计算方法的迫切需求。     

含激波的湍流流动高精度大涡数值模拟方法

针对采用亚格子模型进行含激波的湍流流动模拟时会面临激波附近的精度损失问题,考虑从通过亚格子模型以及数值模拟方法两方面的改进来实现湍流流动大涡模拟的精度提高.大涡模拟采用了Yee及Sj(o)green (2009)提出的高阶低耗散方法.该方法采用自适应的流场探测器以控制计算中所需区域的数值耗散,并考虑对动力学模型采用在激波位置使用Sagaut和Germano(2005)提出的单边亚格子过滤器和(或)直接禁用亚格子项等方法加以改进.对于标准的马赫数1.5和3条件下的激波-湍流干扰问题,上述新方法相较于全区域采用亚格子模型的方法均表现出了相似的精度提升.同时实现的数值精度改进方案采用了Harten的亚单元分辨过程来定位和锐化激波,并在精确激波位置附近的网格点处采用了单边测试滤波.     

串列叶栅尾迹特性的实验研究

对串列叶栅的尾迹特性进行了实验研究。研究表明：串列叶栅尾迹平均速度和单列一样满足相似律。串列叶栅的尾迹平均速度和紊流脉动速度沿额线方向比单列叶栅变化平缓。对于串列叶栅在叶栅下游相对距离Ｘ／ＣＴ＞０．７时，尾迹纵向紊流度基本不随距离Ｘ／ＣＴ变化。