化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。     

局部催化特性差异对气动热环境影响的计算分析

高温气体非平衡效应及其壁面催化效应对高超声速飞行器气动热环境造成显著影响,是当前高超声速飞行器气动热环境预测和热防护设计的关键问题之一。考虑高温空气离解与电离等化学反应、气体分子热力学激发、流动中的非平衡效应和壁面催化效应,通过数值求解三维热化学非平衡Navier-Stokes方程和壁面处质量、能量平衡关系,完善了高温气体热化学非平衡流场有限催化气动热环境数值计算方法和计算程序,采用典型算例进行了考核验证。在此基础上,开展了不同条件下高超声速飞行器热化学非平衡流场气动热环境数值模拟,分析局部催化特性差异对气动热环境的影响。研究表明：所建立的高超声速飞行器热化学非平衡流场有限催化气动热环境数值计算方法及程序,其数值模拟结果与飞行试验、文献符合;局部催化特性差异会导致热流跳变,其热流跳变量与催化特性差异量、材料分布方式等有关;催化特性差异较大时,局部区域热流可能远远高于飞行器全表面完全催化的热流结果,此时将飞行器在全表面完全催化（FCW）和完全非催化（NCW）条件下的数值模拟结果作为实际飞行过程中表面热流的上、下限这一简化处理方式,是不可取的。     

基于嵌套网格的CHN-T1标模数值模拟

为响应第一届航空CFD可信度研讨会(AeCW-1),对CHN-T1标模翼身-尾翼组合体构型采用两种RANS求解器进行了基于结构化嵌套网格的网格收敛性研究和抖振特性计算,并与对接网格结果进行了对比分析。参考AeCW-1网格生成指南自主研发了一族以约3.3倍规模增长的稀-中等-密-特密嵌套网格系列。相比于分区点对接网格,融合了贴体网格和笛卡尔网格优势的嵌套网格从根本上降低了网格生成难度,单块网格拓扑形式更为合理且正交性较好,同时能够有效地平衡近场和远场的网格量。对于跨音速运输机构型,同等规模下嵌套网格的物面网格密度较对接网格更大,能够对激波和分离等复杂流动进行更理想的预测。对比计算过程中,自研求解器OFS3D表现出了较高的嵌套网格计算效率和可靠性,而NASA CFL3D求解器的嵌套网格计算能力有待进一步挖掘和验证。计算还发现,不同湍流模型得到的表面压力分布、分离区大小以及最大升力系数存在差别,而风洞模型支撑和机翼静气弹变形对飞行器力矩特性预测影响较大。     

航空橡胶密封件力学与密封性能检测技术

介绍了航空橡胶密封件的力学、泄漏、损伤、接触变形等特性的无损检测技术及应用情况。其中,橡胶密封件泄漏压差测量方法,用于获得不同压差下的气体泄漏率;橡胶密封件接触应力超声检测技术,实现了橡胶密封接触应力的准确测量;橡胶密封件损伤与泄漏的红外热成像测试技术,可以得到泄漏点的定位与泄漏量的精确测量;橡胶密封件接触变形数字图像相关技术,获得了织物接触界面变形场与应力场的演化规律等。同时,对航空橡胶密封件的非接触测量技术进行了展望。     

基于人机闭环稳定性的舵机速率限制边界

舵机速率限制是造成电传操纵飞机人机耦合的主要原因。利用描述函数法对舵机速率限制非线性进行建模,分析舵机速率限制非线性、人机耦合发生频率及人机闭环稳定性间的关系,并基于此提出舵机速率边界的确定方法;以典型放宽静稳定性飞机为例,基于最优McRuer驾驶员模型,确定人机闭环稳定性所需的最小舵机偏转速率;基于开环起始点(OLOP)准则对所确定的速率限制边界进行验证。结果表明:本文提出的舵机速率边界的确定方法最小成本地避免了人机耦合;所确定的舵机速率限制边界与OLOP准则边界对应的舵机速率基本吻合,即所建立的舵机速率限制边界确定方法合理。     

基于绳系并联机器人支撑系统的SDM动导数试验可行性研究

详细给出了在低速风洞中,采用绳系并联机器人(WDPR)支撑模型,用强迫振荡法进行标准动态模型(SDM)动导数试验可行性的研究。试验中将杆式六分量应变天平内置入模型中以测量模型的气动力和气动力矩,建立了适用于绳系并联机器人支撑系统的模型运动控制子系统和数据采集子系统。采用绳拉力作为参考信号,对气动力矩信号与位姿信号进行数据的同步处理,解决了绳系并联机器人支撑系统应用于动导数试验时所测力矩信号与位姿信号之间的相位差确定问题,给出了WDPR支撑下模型动导数的计算方法。整个试验样机置于某开口式低速直流风洞中进行了俯仰、带偏航角的俯仰以及升沉的动导数试验,通过测量和计算得到各动导数。试验结果与参考文献相比较具有合理的一致性。研究结果表明,采用绳系并联机器人支撑模型进行动导数试验是可行的,至少对于SDM是这样的结果;使用一套绳系并联机器人支撑系统,可以完成多套硬式支撑系统才能完成的动导数试验,从而提高试验效率,降低试验成本。     

防冰疏水微结构表面的设计

微结构表面具有出色的防水自净性能以及减阻性能,但目前超疏水微结构表面设计制备思路欠缺系统性,很多制备方法也比较复杂,在工程实践方面缺乏指导性。从工程制备的角度,概述了微结构表面设计的经典理论模型和国内外对微结构几何设计的理论研究,探究常用工程材料的理论本征接触角,通过力学平衡分析方法计算仿真确定疏水微观结构几何尺寸、几何截面外形,提出截面线长度及尺寸系数两个疏水微结构设计的参数,归纳了设计疏水微结构的流程,为制备疏水微结构或微结构模具提供理论基础。     

等温锻造对TA15合金组织性能的影响

研究发现,在(α+β)双相区内等温锻造工艺参数对合金组织有较大影响;TA15钛合金室温拉伸强度在950~990MPa范围内变化,室温冲击性能在45~59J/cm2范围内变化.随锻造温度的升高,室温拉伸强度有所降低,室温冲击性能有所升高;而随锻造速度的减小,室温抗拉强度和冲击性能均逐渐降低.     

转台电磁锁辅助脱锁电路

本文从实际需求出发,提出了一种转台电磁锁辅助脱锁电路,该电路利用反冲原理,可以使转台电磁锁快速可靠的脱开.通过对电路的原理分析、参数选择、关键参数的仿真及实例验证,最终设计出一种电路简单,快速可靠的电磁锁辅助脱锁电路,为某些转台电磁锁不能完全脱开问题提供了一个良好的解决方案.     

安装角对压气机叶栅气动噪声特性的影响

 基于大涡模拟(LES)和边界元方法对轴流压气机叶栅湍流流场以及流场诱导的噪声进行计算,在不同叶栅安装角下研究来流攻角和来流雷诺数对叶栅气动噪声产生、辐射的影响。研究表明:来流雷诺数不变时,同一安装角下,随着来流攻角从-5°~20°变化,叶栅监测曲线上的声压级先减小后增大,在0°来流攻角下声压级达到最小。安装角为45°时,外场总声压级随来流攻角的分布与30°安装角变化趋势相近。但安装角为60°时,总声压级的变化则明显变缓。在0°来流攻角下,总声压级比安装角为30°和45°时增加了近6 dB,但在其他正来流攻角下,变化并不明显。叶栅的最小声压值出现在弦线方向附近,安装角改变时,最小声压级出现的位置也不同。安装角不变,随着来流雷诺数的增大,叶栅表面的分离减小,损失降低。但叶栅表面的压力脉动随着来流雷诺数的增大而增大,使外场辐射噪声增加。