基于热力学的涡扇发动机神经网络建模方法

由于无法掌握不同发动机的真实部件特性，传统热力学模型对在翼涡扇发动机的建模存在较大的建模误差；同时，热力学模型在特性图边界线附近迭代时，容易迭代到特性图之外，造成迭代过程的不收敛。针对上述问题，本论文提出基于热力学过程的涡扇发动机神经网络建模方法，在神经网络模型的训练过程中充分考虑对部件共同工作热力学约束的优化，提高发动机建模的准确性。通过构建部件级网络结构、部件共同工作损失函数及融合训练过程，将基于部件特性图的传统热力学模型迭代过程转化为部件级神经网络的多目标优化与训练过程，提高了模型的收敛性及建模准确性。模型在26 970条发动机实际飞行数据上进行了训练及测试，结果表明，在相当宽松的准稳态数据下，论文提出的建模方法最大误差可以达到7%左右，比基于部件特性图的热力学模型低5%左右。     

软着陆验证试验中月面特性模拟方法

月面特性的模拟是月面软着陆验证试验的重要设计因素。文章结合探测器的设计状态和试验需求,论述了月面特性模拟的具体要求;基于月面地形统计分布规律,实现了月面原始形貌的模拟;通过高程剔除设计和模块化组合设计,实现了对多种典型月貌的快速模拟;最后通过月表反射特性的模拟,全面满足了试验要求,并取得了良好的试验效果。文中所述的模拟方法可为我国后续行星表面探测器及着陆技术的验证提供借鉴。     

基于隐式温度修正的二维热线风速仪校准方法

热线风速仪主要用于湍流流场测量。热线测量的主要误差来源是由环境温度变化导致的热线校准参数的变化。Brunn提出的传统温度效应修正方法在环境温度变化较大的情况下会带来明显的误差。本文经过理论推导,提出了一种新的基于隐式温度修正的二维热线风速仪校准方法,并开展了热线校准实验研究。结果表明:当环境温度变化在4℃范围内时,本文提出的方法与Brunn方法的测量精度相当;当环境温度变化超过4℃时,本文提出的方法的测量精度远高于Brunn方法。     

发动机喷流对飞行器底部流动影响数值模拟

针对发动机燃气喷流对底部流动的影响开展研究。建立冷喷与热喷计算方法，与经典的高压空气尾喷管喷流试验数据进行了对比，验证了本文建立的三维喷流方法的可靠性。对本文选用的飞行器外形采用冷喷与热喷方法开展了对比计算并与飞行试验值进行比较，分析了两种方法结果的差异。采用热喷方法对来流马赫数 2.5 ，不同飞行高度及喷管进口总压开展计算，研究飞行高度及喷管进口总压对发动机喷流及底部流场的影响。结果表明，保持飞行高度、来流马赫数不变，喷管进口总压增加，底部压力系数逐渐提高。燃气质量浓度最大值位于底部空腔的壁面处，且保持一个恒定值。保持喷管进口总压、来流马赫数不变，飞行高度增加，喷流高速区向后移动且中心区最大马赫数增加。在一定飞行高度下，底部压力系数由负转正，即飞行器底部会出现正推力，这对飞行器的射程会产生重要影响，需要提前评估。     

基于析取关系直接变换的冲突解脱方法

针对自由飞行条件下的飞行冲突解脱问题,提出了一种基于析取关系直接变换的联立求解方法(DDTSA)。首先在不引入额外辅助变量的情况下,将基于析取关系的安全边界条件变换成易于联立求解的约束条件,使问题的变量和方程个数大幅度减少,从而降低问题的求解难度。然后为保证求解的精度和稳定性,采用基于Radau配置点的拉格朗日插值多项式对微分代数方程组进行离散化处理,继而求解离散化得到的非线性规划问题。最后分别对2架、3架以及4架飞机的冲突解脱问题进行仿真,仿真结果表明与相关文献中的方法相比,该方法具有更高的计算效率并且能够得到更优的结果。     

基于WENO-分段线性格式的旋翼流场数值模拟

为了降低二阶JST(Jameson-schmidt-turkel)格式导致的数值耗散,发展了一套适合于格心格式的基于加权本质无振荡(Weighted essentially non-oscillatory,WENO)-分段线性格式的旋翼流场数值模拟方法。采用运动嵌套网格方法生成围绕旋翼的网格系统,主控方程选择Navier-Stokes(N-S)方程。为了更加有效地对桨尖涡等流动细节进行捕捉,在笛卡尔背景网格上采用七阶Roe-WENO格式计算对流通量;在桨叶贴体非结构网格上采用二阶精度的Roe-分段线性格式计算对流通量。时间离散采用了高效的双时间隐式LU-SGS(Lower upper symmetric Gauss-Seidel)方法进行时间推进。最后,应用上述方法对悬停状态的C-T(Caradonna-tung)旋翼和Helishape 7A旋翼进行了数值模拟,将数值计算结果与实验数据进行了对比,计算值与实验值吻合较好;并将桨尖涡模拟效果与二阶JST格式的模拟效果进行了对比。对比结果表明:本文方法能有效对旋翼流场进行计算,且在相同计算条件下,WENO-分段线性格式能够更有效地捕捉旋翼涡流场的流动特性,表明在计算旋翼涡流场时WENO-分段线性格式相比传统二阶JST格式具有更低的数值耗散。     

高温铝合金电磁超声检测回波特性及因素分析

针对高温铝合金在线检测条件下，温度对铝合金电磁超声检测回波特性的影响规律尚不明确、高温检测时缺陷定量/定位补偿困难这一难题，以螺旋线圈电磁超声换能器（EMAT）为例，建立了高温铝合金EMAT检测过程的场路耦合有限元模型；研究了温度对EMAT激励/接收换能效率、EMAT激励/接收电路的功率分配特性、超声传播过程中的扩散/介质衰减特性、回波幅值和超声声速等因素的影响规律；研制了耐高温EMAT探头，对20~500℃高温铝合金试样进行了检测实验，并测定了高温铝合金的超声介质衰减系数和超声声速。在仿真和实验相结合的基础上，分析了高温检测时超声回波幅值变化特性及其影响因素。结果表明：对于铝合金这类非铁磁性金属材料，导致高温时超声回波幅值下降的主要原因是超声介质衰减系数随着温度的升高而增大，其次为高温时EMAT激励/接收电路的功率分配特性的改变。在激励EMAT在试样表面形成的洛伦兹力不变的条件下，其所激励的超声波回波幅值具有随着温度的增加而增加的特点，可以有效减缓超声回波幅值下降的趋势。     

地铁车站屏蔽门风压特性试验研究

针对A+型标准地铁列车模型（缩尺比1:20），通过动模型弹射试验研究了地铁列车过隧道、过站台、跟随工况下的隧道内风压、屏蔽门风压分布与变化规律。研究发现:过隧道工况下，列车经过区间泄压井时，会产生与列车进入隧道时类似的压力波，但风压极值略小。过站台工况下，受前方传来的压力波影响，屏蔽门风压出现一极大值；列车通过站台时，屏蔽门风压出现另一极大值，并在车头经过后立刻达到极小值。这些极值风压决定了屏蔽门的强度设计标准。跟随工况下，当前后方列车尚有一定距离时，屏蔽门受压力波影响而出现风压极值，随后压力略为减小并持续一段时间，这一侧向压力是导致屏蔽门无法正常开闭的主要原因。     

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

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