面向非结构混合网格高精度阻力预测的梯度求解方法

针对非结构混合网格的特点,通过改进传统的Green-Gauss梯度求解方法,提出了一种可提高非结构混合网格黏性计算精度的节点型Green-Gauss梯度求解方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析。改进后的梯度求解方法残差收敛更好,下降量级更多,阻力系数和试验吻合更好,激波区域压力分布和分离区域流场细节的模拟更精确,说明改进后的梯度求解方法有效提高了程序的鲁棒性和阻力预测精度,验证了方法的有效性。采用改进后的方法对第5届AIAA阻力预测研讨会的通用研究模型（CRM）进行了详细的模拟分析,结果表明：改进的梯度求解方法更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

基于有限时间理论的临近空间拦截器末制导律PWPF调节器研究

针对拦截临近空间高速机动目标,基于有限时间理论,设计了有限时间收敛末制导律和脉冲宽度脉冲频率（PWPF）调节器。首先,介绍了加入PWPF调节器的末制导数学模型及工作原理;其次,设计了有限时间收敛制导律,应用有限时间理论证明了收敛条件和性质,并推广到推力受限条件下;最后,考虑视线（LOS）角速率快速收敛和发动机最小工作时间要求,设计了PWPF调节器参数。仿真结果表明,视线角速率在制导结束前的有限时间收敛至零,对目标机动具有鲁棒性,具有较高的制导精度,避免了发动机频繁开启,节省了燃料。     

AZ31镁合金板热成形中的屈服和损伤:本构实现与数值分析

为了准确预测各向异性镁合金板的成形质量,将改进的GTN损伤模型与各向异性拉压不对称的CPB06屈服本构耦合,并考虑了屈服面形状随塑性应变累积的变化,得到了考虑本构参数随塑性应变演化的各向异性屈服CPB06-GTN损伤模型。基于该模型,在ABAQUS/Explict中编译得到了相应的VUMAT子程序,采用单个单元进行了单轴拉伸和压缩模拟,并通过与实验一致性对比验证了子程序的正确性。使用子程序不仅能够模拟镁合金的各向异性屈服及其不规则的硬化,同时也能够模拟镁合金的损伤破坏。此外,采用编写的子程序进行了热拉深成形的数值模拟,模拟预测与实验结果对比表明,采用各向异性损伤模型的计算结果能够准确预测镁合金的变形及其损伤破坏,模拟结果与实验数据吻合;采用合适的压边力和成形的温度条件及非等温成形方法能够提高镁合金的成形性。     

煤油液滴高压蒸发特性

采用真实流体模型描述高压下流体热物理性质的非理想性,并采用状态方程（EoS）法计算多组分高压气-液相平衡及环境气体溶解性,在此基础上建立包含亚临界和超临界两种不同机制的瞬态液滴高压蒸发模型。针对中国新一代高压补燃液氧/煤油发动机,详细研究了煤油液滴在超临界环境下的高压蒸发特性及各因素影响机理。结果表明：高压环境会显著加快液滴温升速率,但弱超临界环境下仍然为相平衡控制的亚临界蒸发状态;只有强超临界环境下才较容易发生扩散控制的超临界蒸发状态。在高压、高温环境下,忽略气相溶解性将导致液滴蒸发速率明显偏小。针对弱超临界环境,温度升高会使液滴蒸发速率单调增加;压力升高则在低温下降低蒸发速率,而在高温下加快蒸发速率。针对强超临界环境,温度升高只提升初始亚临界蒸发阶段的蒸发速率,而超临界蒸发阶段的蒸发速率与环境温度无关;压力升高则同样会提升初始亚临界蒸发阶段的蒸发速率,但会降低超临界蒸发阶段的蒸发速率,此时的总蒸发寿命随压力升高小幅下降。     

结冰飞机着陆阶段飞行安全包线确定及操纵应对策略

结冰会导致飞行安全包线收缩、严重威胁飞行安全,研究结冰后安全包线变化对于操纵应对策略设计及提高飞行安全具有重要意义。以美国国家航空航天局（NASA）的项目飞机GTM（Generic Transport Model）为研究对象,通过对飞机的气动参数进行多项式拟合,建立了结冰飞机纵向通道的动力学模型。为了得到能随结冰程度变化的安全包线,将可达性分析理论引入到对结冰飞机着陆过程的安全性分析。提出将正向可达集与反向可达集的交集作为飞行安全包线,其中可达集的确定是基于水平集方法求解哈密尔顿-雅克比方程的最优解。最后针对不同程度的结冰条件进行了操纵时域验证,并提出了相应的操纵控制策略。研究结果表明,轻度结冰对安全包线影响较小,整个着陆过程的飞行状态始终能在最优控制指导下保持在安全包线以内;但对于重度结冰,飞行安全包线收缩严重,常规操纵已经很难使飞机达到着陆要求,需要进行飞行状态改出处理。研究结果为指导飞行操纵及实时包线保护打下基础。     

动基座近舰面流场数值模拟

航母尾迹流场对舰载机的着舰有较大影响,所以需要对其流场特点进行研究,分析不同状态下舰载机气动特性的变化。采用嵌套网格技术对航母处于垂荡状态下无人机的着舰进行了模拟。首先,利用SFS2舰船进行数值计算,验证了舰船流场的数值模拟方法。然后,对比了单独无人机定常与非定常计算结果,表明所建立的嵌套网格适用于无人机流场的模拟。接着,对航母单相流和两相流的流场结果进行了分析,结果显示甲板下方的流动对甲板上方流场没有大的影响。因此,忽略了水的影响只对航母在空气流场中的特性进行研究,结果表明航母尾迹非定常特性明显,静止航母下滑轨迹上的速度均处于周期性波动状态,且波动幅值随着远离航母而逐渐衰减;而在垂荡情况下航母尾迹变得更加紊乱,水平方向速度波动的周期性减弱,但垂向速度的波动幅值进一步增大。对于静止航母,无人机在不同时刻着舰气动特性的变化也存在差异;当航母处于垂荡状态时,无人机的升力和俯仰力矩在短时间内会有更大的波动。     

带压缩性修正的离散涡方法在结冰数值模拟中的应用

离散涡方法（DVM）是一种无网格的涡运动算法,适用于解决易产生分离流的非定常问题,将其应用于结冰过程中的流场求解,在有效模拟分离流动的同时能避免冰形尖角对网格质量的影响。但离散涡方法基于不可压 N-S 方程,无法应用于预测可压缩流动下的结冰过程。本文在离散涡方法的基础上添加普朗特—格劳尔特压缩性修正,进行基于离散涡方法的可压缩流动下的数值模拟,并将其应用于翼型结冰预测;对流场分布、结冰冰形和结冰模型计算过程等仿真结果进行对比和分析。结果表明：引入压缩性修正后的离散涡方法能较好地模拟可压缩流动,与实验值相比,基于该方法得到的结冰数值模拟结果符合良好,对结冰数值模拟在工程上的应用提供了一定的参考。     

NNW-TopViz流场可视分析系统

流场可视化技术采用图形图像直观地表现CFD数值模拟的计算结果,使用户能够方便地对这些数据进行分析、比较和研究。然而,CFD数值模拟的流动复杂,其产生的流场数据规模巨大、数据类型复杂、特征提取困难,传统的串行可视化软件效率低、交互手段单一,难以满足数据分析的需求。国家数值风洞（NNW）工程研制了一套流场数据处理可视化软件系统（NNW-TopViz,简称TopViz）,具有对流场数据处理与特征提取、几何图形绘制等可视化与交互功能。根据可视分析效率需求,TopViz实现了线程并行,在多核计算环境下有效提高了可视化计算和交互效率;针对流场特征提取困难、常规方法效率低的问题,TopViz实现了基于卷积神经网络的流场旋涡特征提取方法,提升了特征提取准确率和效率;为提高软件交互效率并提供便捷的交互方式和体验,基于头戴式显示设备和体感控制器构建沉浸式虚拟显示与交互平台,TopViz实现了手势和眼球凝视2种交互方法,提供沉浸式环境下多视图、多角度流场探测方式。     

阻旋栅对密封静力与动力特性影响的数值分析与实验研究

密封动力特性对旋转机械转子系统稳定性影响较大,在密封入口端部设置阻旋栅是提高密封稳定性的有效方法。设计加工了无/有阻旋栅共5种密封结构,从数值分析与实验研究两个方面研究阻旋栅对密封静力与动力特性的影响规律。建立阻旋栅密封静力特性CFD理论模型,数值分析阻旋栅对密封泄漏量、切向速度以及周向压力分布的影响;应用不平衡同频激励法实验研究阻旋栅对密封动力特性的影响。研究结果表明:阻旋栅可降低密封的泄漏量,减小密封内流体的切向速度,进而降低密封内流体的周向压力差,且随着阻旋栅周向稠度与径向长度的增加,这种作用逐渐增大,这是阻旋栅抑制气流激振力的主要原因;预旋是密封产生交叉刚度的重要因素,密封的交叉刚度随进出口压比与转速的增加而增大;阻旋栅可有效降低密封的交叉刚度,增加密封的主阻尼,提高密封的稳定性。本文研究揭示了阻旋栅抑制密封气流激振力的机理,为设计阻旋栅密封提供理论依据。     

轴向间距对压气机通道堵塞及总性能的影响

为了研究转静子叶片排之间的轴向间距对压气机内部流动堵塞及气动性能的影响,选取某单级轴流压气机为研究对象,采用多通道非定常数值计算方法对其5种不同轴向间距下的内部流场进行了全三维数值模拟。结果表明:在每一种轴向间距下,当压气机节流至某一工况之后,压气机通道内的流动堵塞区主要集中在转子叶顶间隙区域和动叶吸力面尾缘附近以及静叶吸力面轮毂角区内;在同一流量下,随着轴向间距的减小,转子叶根吸力面尾缘处的流动堵塞区有所扩大,但转子叶顶间隙区域及静叶吸力面轮毂角区内的流动堵塞区体积却不断减小,压气机通道内回流区的总体积也随之减小,其结果是压气机的静压升能力和流动稳定性增强且效率增大。通过进一步研究发现:在同一流量下,当轴向间距减小时,转子叶顶间隙区域内的主流轴向动量增大且泄漏流的轴向动量减小,其结果是转子叶顶间隙区域内流动堵塞区的体积减小。