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

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

飞行仿真气动力数据机器学习建模方法

基于机器学习思想,提出了一种大空域、宽速域的气动力建模方法。该方法利用飞行仿真弹道数据辨识的气动力数据,采用人工神经网络技术,实现了对高度、速度、姿态和舵偏角等多维度强非线性特性的全弹道气动力数据的高精度逼近。首先,分析了神经网络层数、隐含层神经元个数等对建模误差的影响,通过对典型弹道气动数据的神经网络建模计算,确定了较合适的神经网络层数和较优的隐层神经元个数。进而,利用飞行仿真的弹道数据辨识出沿弹道的气动力,采用神经网络建立了包含多个弹道融合的气动力模型,输出量分别为三轴气动力系数和力矩系数。最后通过气动模型输出量与原样本数据的对比,以及4条未参与训练弹道气动数据的预测,验证了该气动力建模方法具有较高的精度。建模结果表明:采用神经网络方法建立的飞行器气动力模型,对拟合多源耦合输入全弹道非线性气动力是可行的和有效的,在样本覆盖的高度、速度、姿态和控制舵偏角范围内,气动力拟合能力较强,并具有一定的外推性。该项研究可以为基于飞行试验数据的气动建模提供新的方法,并且能为飞行器气动力数据挖掘、飞行仿真和总体性能分析提供参考。     

宽弦风扇叶片颤振预测的工程研究

为研究宽弦复合弯掠风扇叶片的颤振问题,以及适用于工程设计的颤振预测方法,分别采用经验法和数值模拟方法,对大涵道比风扇性能试验件叶片进行颤振预测,并将预测结果与试验结果进行对比,考察颤振预测方法的准确性和实用性。对比过程中,分析了颤振发生的机理,提出了抑制颤振的手段。结果表明:两种方法都准确地预测了风扇转子叶片发生的颤振,与试验现象吻合,计算量在可接受范围内,具备工程实用价值。     

考虑空腔的高超声速多流动区域同步数值模拟

先进的高超声速飞行器具有薄壁空腔结构,在飞行过程中受热会产生空腔内气体流动现象,从而影响流场和结构的温度分布。采用数值方法准确模拟高超声速流场、结构温度场和空腔内流动对热结构分析是很有必要的。以研究空腔流动对结构温度分布影响为目的,发展了一种适用于多流动区域流场/结构温度场耦合问题的同步计算方法,并以高超声速带空腔结构物体为例,数值研究了其外部气动热/结构热传导引起的空腔热对流问题。以已发展的高超声速外流场/结构温度场同步计算方法为基础,为了进一步考虑空腔内低速流场,采用了预处理矩阵方法。在流场与结构温度场的交界面两侧分别引入虚拟单元,从而高效地实现相邻场之间物理信息交换。首先通过标准算例验证了方法在求解单独气动热/结构热传导问题以及空腔自然对流问题中的准确性。进而对封闭和带有开孔的两种高超声速运动圆环分别进行多流动区域同步数值模拟。计算结果表明,由于结构温度不均匀引起的空腔内热对流反之也会对结构温度场分布产生轻微的影响。在空腔内气体流动的影响下,封闭圆环的前缘温度在35s内最多下降0.8%左右。对于带开孔空腔的圆环,其孔壁周边温度在0.5s内能够超过外流前缘驻点温度。     

高超声速全机外形气动加热与结构传热快速计算方法

发展了一种无黏流场解与工程计算方法相结合的高超声速全机外形气动加热与结构传热快速计算方法。该计算方法结合了三维块结构网格无黏流场数值计算技术可处理复杂外形流动的优点与工程计算方法效率高的特点,将气动热的计算简化为绕飞行器的无黏外流(边界层以外)数值解和边界层内热流求解两个部分,同时耦合了防热结构传热计算模型、高温化学非平衡热效应估算方法以及弹道状态动态插值方法,可用于快速计算与分析三维复杂外形高超声速飞行器在弹道飞行状态下全机热环境参数、防热结构内温度场等随飞行时间的变化特性。以RAM-CⅡ、类Ⅹ-37B等典型高超声速飞行器为研究对象,在设定的飞行条件及热防护方案下,进行了气动加热与结构传热问题的求解,给出了全机表面热流密度与防热结构材料温度的时变特性。结果对比表明,所发展的方法具有快速、高效的特点,且计算精度可满足工程设计初期选型需求,可为高超声速飞行器的热防护系统初期设计及热环境特性快速计算分析提供技术支持。     

空间电动力绳系电流峰值点变化规律研究

为研究电流峰值点在几种影响因素下的变化规律,建立一种能够捕捉绳系电子发射与收集自洽平衡过程的新算法——电路空间耦合算法。为验证该方法的计算精度,以1.35kW霍尔推力器为等离子体源,在真空舱内开展绳系的电荷收集试验,在电路参数方面,计算精度约为12.8%,在场参数方面,计算精度约为3.6%。在此基础上,针对不同偏置电压、绳系长度以及绳系直径,对绳系的绳上电流分布、电势分布以及空间电势分布等参数进行数值计算。结果显示:电动力绳系的电流峰值点会随着偏置电压升高、绳系长度增加及绳系直径增大而发生比例上的向阳极端漂移,揭示了电子轨道运动限制的机制在各类壁面电荷输运机制中占优,导致电流峰值点漂移的产生。     

基于正交优化设计的诱导轮空化性能研究

为了研究不同结构因素对诱导轮空化性能的影响及主次顺序,基于正交优化设计,采用数值模拟的方法,通过诱导轮可视化试验进行验证,对某型诱导轮及其改型进行了分析。仿真计算得到的不同空化数下空穴形态与实验结果吻合较好,断裂空化数的计算结果与实验结果误差为0.3%。仿真及实验结果表明,诱导轮发生扬程断裂不仅与空穴区占流道体积的大小有关,还与空穴区在流道中分布的位置有关,空穴区位置越靠后,越容易发生扬程断裂;对诱导轮空化性能及水力效率影响最大的结构因素为入口安放角,其次为入口修圆包角;入口等螺距段长度和变螺距段长度较小时可以提高诱导轮的空化性能。     

涡轮叶片尾缘凹坑/凸起结构气膜冷却特性研究

为探究吸力面凹坑和凸起结构对涡轮叶片尾缘气膜冷却特性的影响,在吹风比M=1.1时(雷诺数Re=2.5×10~5),采用数值模拟方法,通过在叶片尾缘吸力面上加入凹坑或凸起,对涡轮叶片尾缘的冷却性能和流动机理进行了详细分析。结果表明:与原始结构相比,叶片尾缘凹坑和凸起结构提高了劈缝出口下游远距离端X/H6 (H为劈缝宽度,为4.8mm)区域气膜冷却效率,对下游的X/H6区域气膜冷却效率影响较小;三种叶片尾缘结构,沿着流向方向会产生由二维展向涡到发卡涡,再到流向涡的变化过程,凹坑和凸起结构通过抑制流体的扰动,改变流体流动情况,提高了劈缝出口下游远距离端气膜冷却效率。     

背压对撞击式喷嘴雾化特性影响研究

为研究背压对撞击式喷嘴雾化特性的影响,将压力的变化等效为气体密度的变化,基于一种树形自适应加密算法,通过直接数值求解不可压Navier-Stokes方程组实现了不同背压条件下射流撞击雾化的数值模拟。首先将数值模拟结果与试验数据进行对比,验证数值模拟的有效性,在此基础上开展了高背压条件下雾化过程的数值模拟。结果表明,随着背压的提高,气动力相应增强,液膜的破碎更加剧烈,一次雾化区域的液滴数密度增大,雾场由稀疏向稠密发展;液膜在向下游运动过程中波动速度的幅值逐渐增大,并且随着背压的提高,液膜的波动由线性向非线性转变;背压增大导致破碎长度减小,经过参数修正得到了液膜破碎长度的经验公式,并与试验数据进行了对比;背压对液滴尺寸分布规律没有显著影响,但随着背压的提高,同一时刻大液滴所占的比例提高,整个雾场的Sauter平均直径有增加的趋势,当背压从0.1MPa增大到1MPa时,雾场的Sauter平均直径由155.5μm增大到166.9μm;背压增大,液滴粒径分布的均匀度指数减小,液滴尺寸分布更加不均匀。     

涡桨飞机发动机进气道排除异物特性数值方法研究

根据适航规范针对常见的两种外来异物：冰雹与冰块进行本体特性研究,确定两种外来异物的几何、质量以及初始姿态特性。耦合计算流体力学(CFD)以及六自由度方法(6DOF),开展涡桨发动机进气道以及旁通道中的异物运动排除特性数值模拟,分析其气动特性以及运动轨迹。针对运动中可能存在的碰撞现象,进一步联合LS-DYNA软件进行仿真分析,建立合理有效的涡桨发动机进气道外来异物排除数值模拟方法。以某型涡桨发动机为例,两种异物的计算结果表明,一旦异物与壁面发生碰撞,碎裂成若干很小的碎块,能量损失,对发动机威胁较小,无论是进入主发动机或者旁通道,均可认为排除。但是,由于下壁面结冰区的冰块容易直接进入主发动机,有可能造成严重影响,需要重点关注。