大水滴撞击壁面的动态特性数值模拟

采用 VOF 方法建立了大水滴撞击壁面的计算模型,模拟了大水滴以不同直径、不同速度撞击光滑壁面的动态撞击过程和撞击特性。计算结果表明：在大韦伯数情况下,水滴撞击光滑壁面会在铺展过程中发生边缘水滴飞溅;在水滴撞击壁面的收缩阶段,随着水滴直径的减小和水滴速度的提高,会发生液膜缓慢收缩、边缘液环和液膜分离、中心处部分液膜和边缘液膜分离、液膜完全破裂等不同结果;当水滴直径和撞击速度增大时,同一时刻水滴的铺展半径、最大铺展半径、最大铺展系数增大;水滴在壁面上达到最大铺展系数所用的无量纲时间随水滴直径增大而增大,同一直径水滴不同初始速度下达到最大铺展系数所用的无量纲时间变化较小。     

考虑低温增压真实气体效应的运输机气动特性数值模拟研究

为了给低温风洞试验数据修正提供参考,本文利用数值模拟手段研究了低温真实气体效应相比较于完全气体对飞行器气动特性的影响,以及该影响与雷诺数影响相比所占比例的大小等问题。文章应用 Aungier-Redlich-Kwongz 方程,发展了适用于模拟氮气低温高压真实气体效应的 RANS 求解软件。与 NIST 数据的对比表明,该状态方程在5倍大气压下,比热比等参数误差在0.3%以下。同时,标模测试结果表明本文软件计算精度与国外软件相当。应用本文方法研究低温增压风洞中氮气真实气体效应：以典型运输机构型 DLRF6为模型,分别计算了高速、低速状态下各种不同温度和压力工况下的流动。计算结果表明,在低温增压情况下,真实气体效应引起的气动力差异很小,升力、阻力、力矩最大相对误差均在0.3%以下,与雷诺数效应引起的偏差相比可以忽略不计。因此,以氮气为介质的低温风洞试验研究可以采用完全气体假设。     

存在烧蚀的高超声速压缩拐角流动数值研究

耦合求解热化学非平衡流控制方程和烧蚀壁面边界条件,进行存在石墨烧蚀的压缩拐角流场数值模拟。流场化学反应采用16组元(N2,O 2,NO,N,O,NO +,N+2,O +,N+,CO,CO 2,C,C2,C3,CN,e-)29个反应的非平衡模型,热力非平衡的双温度模型下,不同反应采用不同控制温度。石墨材料表面反应包括碳的氧化反应、碳催化的O 原子复合反应和碳的升华反应。对15°、18°、24°压缩拐角模型,在自由流 Ma =10~30,总焓值6~55 MJ/kg 范围,分别进行无烧蚀的壁面催化与非催化条件和石墨烧蚀条件下的流场计算,分析各类条件下的流场结构、流动分离特性以及流场热化学参数分布特点,研究壁面条件对流动特性的影响。结果表明：流动分离可能性和分离区范围随着压缩拐角斜面倾角增大而增大,随来流马赫数增大而减小;相对于低壁温条件,无烧蚀的辐射平衡壁温和壁面烧蚀条件下流动分离区增大,斜面上压力、摩阻和热流峰值点也有所后移。     

机翼电脉冲除冰效果的仿真分析

电脉冲除冰系统是保障飞机在结冰气象条件下安全飞行的一种机械除冰系统,它具有高效节能、稳定性好、通用性强等优点。国外研究机构从20世纪70年代起就针对电脉冲除冰系统开展了大量的研究工作,国内关于该方面的研究起步较晚。为研究不同因素对电脉冲除冰系统除冰效果的影响,本文分别建立了简单翼型前缘结构以及 NACA2414、NACA2428、NACA2407等多种 NACA 翼型前缘结构的机翼电脉冲除冰系统的有限元模型,并提出了基于除冰率的系统除冰效果评判准则。基于机翼电脉冲除冰系统有限元模型,采用数值仿真,分别研究了不同幅值脉冲力作用下,脉冲力作用位置、翼型几何形状以及脉冲力作用序列对系统除冰效果的影响。通过对不同系统模型的仿真结果进行分析,可以发现：为达到相同除冰效果,采用对称双脉冲作用所需的幅值明显小于单脉冲作用;相同工况下,相对厚度较小的翼型前缘除冰效果较好,且同一翼型的下翼面的除冰效果略好于上翼面的除冰效果;通过适当调整脉冲触发顺序,可以达到或接近通过单纯增加脉冲载荷所达到的除冰效果。本文通过仿真计算,提炼了机翼电脉冲除冰系统除冰效果随各影响因素的变化规律,指出了机翼电脉冲除冰系统的优化方向,为机翼电脉冲除冰系统的设计奠定了基础。     

低速风洞引射短舱动力模拟技术新进展

引射短舱可以模拟发动机短舱的喷流影响,并部分模拟进气影响,能用于研究发动机短舱与机翼及增升装置的气动干扰特性,且具有研制周期短、造价低等特点,是在风洞中开展飞机/发动机一体化设计研究的一种重要试验技术。本文介绍了气动中心低速所在引射短舱设计技术和试验技术方面的新进展。采用商业软件对引射短舱进行了三维流场数值模拟,获得了引射短舱性能和三维流场信息。对引射短舱内部流场进行了分析和研究,对引射喷嘴数量、位置进行了优化,增加了引射短舱的进气流量,改善了尾喷口流场均匀度,明显提高了引射短舱性能。发展了空气桥技术,采用有限元方法进行了优化设计,对空气桥和天平进行一体化设计,并进一步发展了空气桥影响修正技术,解决了供气管路对天平测力的影响问题。发展了高精度流量测量控制技术,采用了数字阀、流量控制单元、短舱内部测量耙等技术,提高了流量的控制测量精度及测量不确定度,流量控制精度达到了0.1%,流量测量不确定度达到了0.3%,引射短舱落压比控制精度优于0.01。研制了短舱移动支撑装置,能够实现引射短舱的独立支撑,并实现短舱前后和上下位置的变化,用于开展短舱位置优化研究。最后,介绍了引射短舱的地面性能测试及风洞试验应用,给出了性能测试与数值模拟的对比结果和典型的风洞试验结果,试验结果表明动力影响使得飞机0°迎角升力减小,升力线斜率增大,失速迎角推迟。     

采用新型基准流场的高超内收缩进气道试验研究

由于新型变中心体基准流场具有压缩效率高、反射激波弱的优点,采用该基准流场设计了矩形转圆形内收缩进气道,在设计点马赫数Ma=6.0进行了风洞试验研究。试验中得到了进气道压缩面的沿程压力分布、隔离段出口皮托压分布等参数。通过和数值模拟对比分析,结果表明：进气道外压段的压力分布明显具有先增大后减小的特征,内压段的压力分布具有两级爬升的特点,且压升较小,流场结构较好。由于内压段流场激波强度弱,进气道总压恢复系数较高,达0.518,并产生了52倍的增压比,其抗反压能力在144倍以上。试验研究表明,采用新型变中心体基准流场能改善矩形转圆形内收缩进气道的内压段流场及隔离段流场,并能有效提高进气道的总压恢复系数。     

周向平均方法在某风扇/增压级分析中的应用

通过对三维（3D）Navier-Stokes方程进行周向平均,得到了通流模型的控制方程,对其采用时间推进有限体积方法进行数值求解。为实现风扇/增压级在设计初期的快速性能评估,考察了周向平均方法在风扇/增压级分析中的准确性。分别利用NUMECA三维数值模拟软件和周向平均通流模型（CAM）对某高通流风扇/增压级进行了性能分析,从对比结果来看,周向平均通流模型在近设计点给出了与三维数值模拟十分接近的特性参数,最大误差不超过2.0%。在风扇转子中,由于周向平均通流模型能捕获通道激波,其物理本质与三维平均结果有所区别,因此径向参数分布与三维有所差异。而在亚声速流动下的增压级及外涵道各叶片排出口参数的径向分布与三维数值模拟结果都能很好地吻合。     

纳秒脉冲等离子体激励器用于圆柱高速流动控制的数值模拟

利用数值模拟,研究了纳秒脉冲介质阻挡放电（NS-DBD）等离子体激励器在圆柱高速流动控制中的应用。首先,研究了单电极NS-DBD等离子体激励器在静止空气中放电后的流场特性。研究表明在介质阻挡放电形成的等离子体区域,有局部能量快速注入,放电结束5 μs后在上极板后端点位置形成了一个局部温度高达900 K的热点,由此引发很强的压力扰动,形成以上极板后端点位置为中心,扩散速度约为声速的半圆形压缩波。在此基础上,通过数值模拟研究了NS-DBD等离子体激励器布置在直径为6 mm的圆柱上,来流马赫数为Ma_∞=4.6时,对圆柱脱体激波的控制作用。研究表明介质阻挡放电形成的半圆形压缩波对于脱体激波有很强的干扰作用,激波距离增加了15.7%,激波强度也有相应的减弱,导致阻力减少了13%。     

几何尺寸对高超声速进气道气动性能的影响

为了探索模型缩尺比对高超声速进气道气动性能的影响,对不同缩尺比的二元高超声速进气道开展了数值模拟研究,结果表明：随着缩尺比的增大,进气道流量系数、隔离段出口总压恢复系数和马赫数均逐渐增大,而静压比逐渐减小,且来流马赫数越高,上述参数变化幅度越大。由理论与数值模拟分析可知,上述现象主要是由于不同缩尺比下,进气道当地雷诺数不同,导致进气道附面层相对厚度变化,进而影响进气道气动性能。理论分析了进气道总压恢复系数与缩尺比的定量关系,就进气道而言,进气道进口处附面层相对厚度减小1%,隔离段出口总压恢复系数提高约0.7%。     

湍流燃烧模型在燃烧室数值模拟中的对比分析

以某单头部矩形燃烧室为研究对象,采用多种湍流模型和燃烧模型进行组合计算,模拟燃烧室内部的速度场和温度场,并对计算结果进行对比分析。结果表明:湍流模型主要影响火焰筒内部主燃孔横截面上游的速度分布,Standard k-ε和Realizable k-ε模型的速度场计算结果差异相对较小;湍流动能预测受湍流模型的影响较大,并具有一定规律性;不同组合模型对燃烧室内部和出口温度分布的局部细节模拟差异较大,燃烧模型影响最大;PDF模型计算的温度值较合理,另外三种燃烧模型在单步完全反应燃烧机理下的计算值偏高,计算获得的OTDF也相差较大。