三类气动导纳数值识别方法的适应性研究

针对气动导纳函数的数值识别方法,借助于CFD,在简谐脉动来流、湍流和竖向阶跃来流下对平板断面和箱梁断面的导纳函数函数进行研究。首先,在无断面存在的空流域内详细研究了简谐脉动来流、湍流和竖向阶跃来流的传播特性及其数值计算方法。其次,对有断面存在的情况进行了数值计算。最后,识别得到了平板断面和箱梁断面在三种不同来流下的气动导纳函数。结果表明:对平板断面,三种方法识别得到的气动导纳函数与Sears函数吻合良好,验证了三种数值计算方法的可行性;对箱梁断面,简谐来流和湍流下识别的气动导纳差别不大。相比之下,完全基于线性叠加原理的阶跃来流方法产生了实质性的偏差,表明该法不宜用于钝体断面。计算效率方面,湍流的计算效率适中且对任意断面适用;简谐脉动来流的计算效率最低,适用于气动导纳与风场无关和弱相关的断面;竖向阶跃流方法具有计算时间短的优势,但它仅能用于气动导纳与风场完全无关的断面。     

低温液氮贮箱增压及排气流量控制方法

在对以液氮为工质的低温贮箱进行热分析的基础上推导了增压速率计算模型,并通过两种增压工况进行了验证。主要开展了直排模式、排气模式及并行模式下的气枕压力及蒸气排气流量的对比实验,通过后两者模式下的蒸气排放质量与直接排放相比,分析了压力控制带内单循环各个模式下热力排气系统（TVS）的效能,结果表明:采用低温蒸气节流的并行模式在蒸气排放质量和气枕降压速率上具有优于其他模式的明显优势,气相并行模式比直排模式质量减少97%。      

基于气动外形优化和主动流动控制的减阻技术

增升和减阻是民用飞机设计的永恒目标。本文对民用飞机减阻技术最新进展及未来发展进行综述,在此基础上,报道了计算所团队通过数值优化设计及主动流动控制实现减阻的研究进展。为此,发展了基于壁面模化大涡模拟(WMLES)准确模拟湍流边界层的先进数值方法,以及利用伴随方程方法实现气动外形快速数值优化的技术。在具体实践中,通过组合基于自由变形(FFD)的几何参数化模块、基于线弹性体法的网格变形模块,以及基于伴随方程方法和序列二次规划(SQP)的优化算法模块,搭建了面向工程、适于超大规模变量优化设计的整机气动优化设计平台,在NASACRM模型上实现了设计变量超过600的大变量、跨声速气动数值优化设计,在合理的工程约束条件下,有效削弱了机翼表面激波,总体减阻率达到2%以上;通过对湍流平板边界层外层结构进行射流非定常控制,在中等雷诺数Reτ=4700条件下,实现了5%～6%的当地减阻率。研究证实,数值优化设计和基于外层射流控制的主动流动控制技术将是大型民用飞机减阻中非常有前景和值得重视的两种方法。     

离子推力器多模式化研究

为了实现离子推力器多模式化,分析了离子推力器功率宽范围调节限制因素,提出了两种宽范围调节策略;针对我国小行星探测任务,完成了30cm多模式离子推力器研制、功率宽范围调节限制条件确定、以及两种调节策略下多模式工作点设计及对比研究。结果显示,通过降低放电室磁场强度可延伸离子推力器最小稳定工作功率,提高束流均匀性,实现离子推力器更宽功率范围多工作点设计;功率宽范围调节主要是屏栅电压和束电流的宽范围调节,二者通过栅极导流系数限制和交叉限制而约束;推力随功率增加呈线性增加关系,比冲随功率的增加总体上呈先快速增加后趋于稳定的趋势;30cm多模式离子推力器在0.25kW～5kW内稳定工作,推力10mN～186mN,比冲1522s～3586s。     

钝体稳定射流火焰的火焰面模型数值研究

基于稳态火焰面(SLFM)和交互式非稳态欧拉颗粒火焰面(EPFM)模型对Sydney大学CH4/H2钝体稳定湍流扩散火焰行了数值研究,采用修正的雷诺应力模型(RSM),同时对两种不同规模的甲烷详细化学反应动力学机理进行研究,比较了燃烧模型和反应机理对湍流火焰结构、活性自由基以及氮氧化物预测精度的影响,与实验数据对比结果表明:两种反应机理得到的温度场和主要组分分布基本相同;SLFM模型能对速度场和标量场的分布进行较为准确预测,采用EPFM模型修正后,部分区域OH预测结果更加靠近实验结果;采用EPFM模型对SLFM模型耦合GRI-Mech 211的计算结果修正后,NO量级降低近2倍,预测精度明显改善,与实验结果实现较好的符合,验证了化学反应动力学机理以及非稳态效应对氮氧化物预测精度的影响。     

基于连续损伤力学的轮盘超转破裂行为预测与验证

为有效预测轮盘的超转破裂行为,建立了1种基于连续损伤力学的预测方法。给出了该方法下基于光滑拉伸试验数据的直接时效GH4169合金弹塑性本构模型参数与损伤模型参数确定方法。对缺口半径分别为0.5、2和5 mm的3种缺口圆棒试样的拉伸响应曲线进行了有限元计算,并与试验结果进行了对比。结果表明:该模型可以较准确地预测缺口试样的极限拉伸载荷和拉伸断裂变形。采用弧长法对模拟盘的超转破裂过程进行了非线性有限元计算,使用模拟盘超转破裂试验进行了验证。该模型不仅可以较准确地预测模拟盘的破裂转速,还能有效预测模拟盘的超转破裂模式。对于实际轮盘破裂转速预测具有一定的指导意义。     

旋转弹体马格努斯效应数值模拟方法研究

旋转的弹体通常会存在马格努斯效应,严重影响飞行稳定性和弹道轨迹。提出基于旋转壁面法的多对称面弹体旋转效应数值模拟方法,无需借助动网格和多参考系模型,仅需要将旋转速度附加壁面,即可使用定常方法模拟马格努斯效应;同时,使用SOCBT标准弹体模型对该方法进行验证分析。结果表明:计算值与试验值吻合较好;由于旋转引起的附面层堆积以及大攻角分离效应、激波的干扰,计算时附面层的网格节点要足够密,并且应该针对不同的攻角选择精度与效率兼顾的计算模型。     

不同来流条件下旋转对气膜冷却的影响

大涡模拟考察了旋转状态来流条件对单孔平板气膜冷却的影响,气膜孔沿流向倾斜30°,气膜出流的雷诺数为2 600,吹风比为0.5,计算了静止和旋转数为0.02时气膜冷却的流动和换热,对比两种主流进口条件下旋转对气膜冷却的影响。计算结果表明：（1）均匀来流条件下,旋转主要使发夹涡结构产生非对称分布,裹挟气膜向高半径方向偏转;（2）充分发展的来流条件下,旋转使来流边界层内产生湍流结构,淹没了射流进入主流时产生的发夹涡,引起更强烈的射流扩散,增大了气膜覆盖范围,降低了气膜冷却效率的峰值;（3）旋转通过改变来流边界层内的湍流结构对气膜冷却的影响更显著。     

来流边界层拟序结构对超声速混合层流场结构的影响研究

超声速混合层涉及可压缩湍流的根本问题,具有重要的应用背景。通过设计超声速混合层实验装置、应用新近提出的高分辨率NPLS测试技术,拍摄了来流边界层分别为层流和湍流流态下混合层的流向和展向流动图像。根据流动图像的特征,分析了混合层的流向与展向流场切面中拟序结构的成因;深入讨论了来流边界层中拟序涡结构与混合层涡结构的相互作用问题;比较了层流和湍流来流条件下混合层拟序结构的异同及其对混合效率的影响。结果表明：当来流边界层为湍流时,对应的混合层具有较高的混合效率。     

Ti-46Al-6(Cr,Nb,Si,B)合金高温变形特性研究

采用Thermecmaster-Z型热压缩试验机,在900~1250℃温度范围内、和10-3~1s-1应变速率条件下对铸态和挤压态Til-46Al-6(Cr,Nb,Si,B)at%合金(以下简称G4合金)进行了热压缩模拟试验,建立了两种状态下G4合金的加工图.并以加工图为基础,结合组织观察,研究了高温下该合金的变形特性.结果表明:G4合金的高温变形性能受温度和应变速率强烈影响,并呈现不同特征;流变应力随变形温度升高而减小,随应变速率增大而增大;挤压态G4合金具有比铸态G4合金更好的稳定流变能力和更宽的可热加工窗口;动态再结晶(DRX)是导致G4合金流变软化和稳定流变的主要原因;铸态G4合金的最佳变形温度为1150~1200℃,应变速率为10-2.5～10-3s-1,挤压态G4合金的最佳变形温度为1050~1150℃,应变速率为10-1.5~10-2.5s-1;G4合金的主要失效模式包括表面开裂、局部流动和楔形开裂.