导航
流体力学与飞行力学栏目所有文章列表
Select 1. 基于超声导波的积冰量辨识及定位方法
张晏鑫, 张鸿健, 熊建军, 赵照, 冉林, 易贤
doi: 10.7527/S1000-6893.2023.29293
2024, 45 (16): 129293-129293.
PDF [6632KB](124)
为实现飞机结冰区域化定量探测能力,提出了一种基于Lamb波的冰层定位和定量辨识方法。首先,建立了压电耦合的波动力学仿真模型,通过三维有限元数值计算分析了Lamb波在铝板、冰层中的传播特性及模态转换过程,选取S0/B1模态作为结冰监测模态,并建立了一种基于小波变换的信号分析方法以提取S0/B1模态波的时间延迟特征量。进一步,研究分析了不同冰层结构尺寸对Lamb波信号的影响规律。结果表明,随着冰层长度(波导传播方向上的积冰量)的增大,B1模态波的时延比线性增大;冰层宽度(波导传播展向垂直方向上的积冰量)对B1模态波时延存在一个影响区域,在影响区域内(冰层宽度不大于40 mm),时延比与冰层宽度成线性关系;冰层厚度与B1模态波时延比总体上具有正相关变化趋势,冰层厚度在1.5 mm以内时,其与时延比之间具有较好的线性关系。基于此,建立了Lamb波压电阵列定位成像的数值计算模型,并提出了一种自适应阈值概率重构算法,进一步采用时延比来表征探测区域内的积冰程度,最终验证了冰层的定位成像、结冰定量辨识的可行性,提升了超声导波结冰在线快速评估的能力,为超声导波结冰区域化定量探测奠定了基础。
Select 2. 基于合成双射流的翼型阵风载荷减缓
王浩, 罗振兵, 邓雄, 周岩, 张鉴源, 赵志杰
doi: 10.7527/S1000-6893.2023.29660
2024, 45 (16): 129660-129660.
PDF [5875KB](95)
为研究基于合成双射流的阵风减缓策略,针对NACA0012翼型在其上表面阵列式布置合成双射流激励器,通过反向射流方式主动诱导流动分离,并详细研究分离涡的演化及动态涡脱落过程,为开展基于零质量射流的载荷控制提供了新的理念。结果表明,随着动量系数Cμ 的增大,合成双射流对于阵风载荷控制能力逐渐提高且没有造成升力系数CL 较大的脉动,当Cμ =0.033时,后缘分离区发展到翼型中部位置,CL,max响应幅值降低了47.9%。由于后缘处分离区较大,关闭射流激励器后分离涡在上翼面的卸载过程会对CL 造成较大影响,在气流冲击分离涡的过程中,上下翼面压差迅速增大进而造成CL 迅速减小;针对这一过程研究Cμ 阶梯递减控制,设置Cμ 区间式变化,相较于持续施加恒定Cμ 的情况,这一控制方法在达到相同减缓幅值的情况下能耗更低,且翼型处于高载荷冲击的时间也更短。
Select 3. 进气畸变条件下非轴对称风扇设计及扩稳机理
余军杨, 傅文广, 孙鹏, 张韬, 王春雪, 赵伟
doi: 10.7527/S1000-6893.2024.29725
2024, 45 (16): 129725-129725.
PDF [10080KB](141)
动力系统的进气畸变问题导致风扇性能下降、引发气动不稳定,从而对飞行器的安全构成潜在威胁,在军、民航空领域已受到广泛关注。为确保进气畸变条件下推进系统的可靠运行和飞行安全,本研究基于数值方法对某小涵道比涡扇发动机开展抗畸变风扇设计及扩稳研究。针对其畸变影响特征,采用非轴对称静子(NAS)布局来消除其不利影响,研究表明:NAS能够有效改善进气畸变条件下的风扇性能。相较于原型,在同流量工况下,绝热效率和总压比可分别增加1.37%和1.26%,稳定裕度可增加8.31%,实现了该NAS风扇的性能提升及扩稳目的;NAS设计直接作用于畸变区,有效改善了畸变区内的静叶冲角,抑制了静叶通道涡系结构的发展和迁移,较大程度减小静叶的扩压因子,从而有效消除了静叶的角区分离,使得静叶冲角周向分布相对均匀,有效改善风扇静叶内部流场,整体提高风扇的气动性能。此外,NAS设计还具有改善风扇后部内涵流场的潜力,使内涵出口的不均匀流场朝着周向均匀化方向发展,并且在有效减弱动态总压畸变强度、改善流场均匀性的同时,还兼具改善总温畸变的能力。
Select 4. 催化床结构对HAN基单组元发动机性能的影响
孙得川, 张国强, 姚天亮, 关亮
doi: 10.7527/S1000-6893.2023.29818
2024, 45 (16): 129818-129818.
PDF [2896KB](40)
针对硝酸羟胺(HAN)基推进剂单组元发动机,研究催化床长度和空腔缺陷对发动机工作过程的影响。采用多孔介质假设、Volume of Fluid(VOF)模型以及推进剂分解模型对300 N发动机稳态工作状态和起动过程进行数值仿真。研究结果表明:增加前床长度有利于提高性能,将前床长度由30 mm增加到40 mm,真空比冲提高2.7%。但是,增大前床长度会使响应变慢,不仅使燃气充填阶段末的压强值和温度值分别减小6.7%和12.5%,而且使催化床升温引起的升压时间增加。空腔缺陷会显著降低发动机的性能,前床长度为40 mm时,5 mm的空腔就会使真空比冲降低15%。
Select 5. 黏性对液滴偏心撞击旋转壁面铺展特性的影响
杨文, 张鑫玺, 汪翔宇, 刘川扬, 张云波
doi: 10.7527/S1000-6893.2024.29847
2024, 45 (16): 129847-129847.
PDF [5241KB](54)
液滴撞击旋转壁面广泛存在于工程实践中,如在航空领域发生的发动机吸雨吸冰现象,这些现象会导致发动机性能下降甚至引发安全问题,因此研究液滴撞击旋转壁面具有重要意义。搭建并开展液滴撞击旋转圆盘实验并用高速相机采集图像,在保持液滴尺寸、撞击速度、表面张力和撞击位置不变的情况下,研究了转速和黏性对液滴偏心撞击后铺展规律的影响。通过后处理液滴铺展动态图片,定量获取液滴沿径向和切向铺展数据。实验结果表明,液滴径向-切向铺展动态模式在旋转邦德数和旋转雷诺数下划分出4个区域:撞击惯性力-黏附力、撞击惯性力-黏附力/黏滞力过渡、撞击惯性力/离心力过渡-黏附力、撞击惯性力-黏滞力。径向上撞击惯性力主导模式下,最大径向铺展因子随奥内佐格数和旋转邦德数增大而减小,且在回缩阶段径向铺展因子与无量纲时间t/(Oh·T)满足-1/3幂次标度律。切向上黏附力主导模式下,液滴切向铺展速率与无量纲时间t/(Oh·T)满足2/3幂次标度律;切向上黏滞力主导模式下,切向铺展存在最大铺展因子。
Select 6. 飞翼布局高低速一体化气动优化设计
赵轲, 邓俊, 黄江涛, 陈树生, 高正红
doi: 10.7527/S1000-6893.2023.29367
2024, 45 (15): 129367-129367.
PDF [11262KB](297)
飞翼布局由于在气动、隐身和结构等方面的综合优势,是未来最有潜力的飞行器气动布局形式。关于飞翼布局的气动设计研究,学术界开展了大量的设计分析,研究发现低速起降性能和操控是飞翼布局性能和安全的主要难点。对于飞翼布局无人机和作战飞机,由于隐身设计的影响,使得其低速起降性能、操控特性更加严苛,以往的研究主要集中在高速巡航的设计,关于低速设计研究很少,因此基于高性能离散伴随优化设计平台,开展了某飞翼布局无人机的高低速综合设计研究,分析对比了不同低速设计模型对飞机低速特性的影响,在此基础上,建立了高低综合设计模型,全面提升了飞机的高低速性能,经过对结果的分析,总结了高低速一体化设计的要点和规律,为飞翼布局的气动设计提供了有力的设计模型和参考。
Select 7. 基于空间嵌套径向基函数的高效并行网格变形方法
路宽, 宋文萍, 郭恒博, 叶坤, 王跃, 韩忠华
doi: 10.7527/S1000-6893.2023.29433
2024, 45 (15): 129433-129433.
PDF [6772KB](176)
高效的网格变形方法可以大幅提高流固耦合数值模拟、基于高精度计算流体动力学(CFD)分析的气动外形优化设计等问题的计算效率。常规的基于径向基函数(RBF)的网格变形方法在变形控制点增多时,会引起网格变形计算量激增,而采用减少变形控制点的方法减小计算量,会带来拟合精度的损失。针对已有RBF网格变形方法计算效率与计算精度不能两全的现状,提出了一种空间嵌套径向基函数模型(SN-RBF),发展了基于SN-RBF模型的高效网格变形方法,在保证网格变形精度的同时大幅度提高了网格变形效率。发展的SN-RBF模型采用多个物理空间相互重叠的子模型代替样本点较多的径向基函数模型,大幅度缩短了径向基函数法网格变形的建模时间。鉴于其便于并行的优良特性,还发展了基于该网格变形方法的并行建模方法与并行网格变形方法,使得网格变形效率得到进一步提升。测试算例表明,建模样本点越多,网格变形效率提升越显著。其中CRM翼身组合体算例,建模效率提升高达16 947倍,网格变形效率提升高达5 218倍。
Select 8. 基于Liutex的数据驱动湍流模型修正
龙家俊, 刘陈飘, 秦飞, 张加乐, 徐圣冠, 高宜胜
doi: 10.7527/S1000-6893.2023.29579
2024, 45 (15): 129579-129579.
PDF [6468KB](66)
现有广泛应用的湍流模型对于大分离流动问题的计算结果往往与实验结果存在着较大偏差,这主要是由于大分离流动中存在的强逆压梯度导致湍流模型的基本假设不再成立。为了提高大分离流动问题的计算精度,需要对湍流模型进行修正。采用流场反演和数据驱动相结合的方法,提出了基于Liutex的Spalart-Allmaras (S-A)一方程模型生成项的神经网络修正方法。首先采用基于离散伴随的流场反演获得S-A一方程模型生成项的修正系数;由机器特征选择方法表明在翼型附近大分离流动区域本文选取的Liutex变量相对于修正系数具有最高的相关性,适合作为神经网络输入;再将Liutex等变量作为输入构造神经网络近似S-A模型生成项修正系数,建立生成项修正系数的神经网络结构,通过S809和S814翼型分离流动结果验证了该修正方法能够有效改善大分离流动的计算精度。
Select 9. 多步法冰晶积冰数值模拟
卜雪琴, 黄平, 林贵平, 楼燕霞
doi: 10.7527/S1000-6893.2023.29308
2024, 45 (14): 129308-129308.
PDF [3594KB](65)
近年来,冰晶结冰问题逐渐成为飞机/发动机防除冰领域研究的热点。为了实现对冰晶积冰过程的准确预测,建立了冰晶结冰计算模型,包括考虑了不规则形状的冰晶运动轨迹模型、冰晶结冰质量动量守恒模型、黏附模型、侵蚀模型。采用时间多步法模拟积冰过程,在每个时间步长内更新空气流场、冰晶粒子运动轨迹、结冰热力学计算结果,并利用动网格技术更新冰形的几何边界。研究结果表明,利用时间多步法进行计算,能够实现稳定冰形的模拟。稳定冰形形成的原因在于:一方面,积冰减少了表面冰晶的黏附量;另一方面,积冰增大了对已形成积冰的侵蚀量,最终使黏附积冰过程与侵蚀过程达到动态平衡。本文模型和计算方法适用于低融化率情况,能够对三维结构冰晶结冰数值模拟提供指导,但对于高融化率情况还需要进一步考虑积冰向下游流走的情况。
Select 10. 针对高阶间断伽辽金数值格式的Gibbs现象智能去噪方法
刘嘉文, 王明振, 欧阳文轩, 虞建, 刘学军, 吕宏强
doi: 10.7527/S1000-6893.2023.29323
2024, 45 (14): 129323-129323.
PDF [9705KB](120)
在使用高阶间断伽辽金方法进行高速可压缩流场计算时,激波附近会出现影响数值精度甚至导致计算失败的非物理数值振荡,这类似于图像处理领域不断堆积的Gibbs噪声。如何抑制激波振荡或消除Gibbs现象并确保计算过程稳定,已经成为了高阶间断伽辽金方法研究领域的一个挑战。针对这一问题,利用机器学习技术,提出了一种由图注意力机制和图卷积网络构成的Gibbs现象智能去噪模型,该模型能够抑制间断伽辽金方法计算中激波附近的振荡,在确保间断伽辽金方法计算顺利进行的同时提升了捕捉激波的效果。该模型使用间断伽辽金方法计算中产生的Gibbs噪声数据构造训练数据集,在图卷积滤波器的指导下进行图神经网络训练。对跨声速和超声速来流条件的NACA0012翼型进行了数值模拟,结果表明在间断伽辽金方法计算过程中嵌入所构建的Gibbs现象智能去噪模型,能够消除Gibbs现象,有效抑制激波振荡。
Select 11. 多段翼构型结冰计算方法及结冰影响分析
任靖豪, 王强, 陈宁立, 刘宇, 易贤
doi: 10.7527/S1000-6893.2023.29328
2024, 45 (14): 129328-129328.
PDF [8359KB](66)
为了揭示不同来流条件下多段翼构型表面结冰特性,针对30P30N增升构型开展了结冰数值模拟研究。在计算方法方面,通过引入最小壁面距离,实现了基于拉格朗日方法的高效液滴收集率计算。采用局部网格拓扑重构技术改善了复杂冰形下网格重构稳健性差的问题,实现了更加鲁棒的多时间步结冰过程自动化计算能力。基于上述算法,开展了结冰特性及带冰后气动特性模拟,计算结果显示,多段翼结冰区域主要分布在前缘缝翼迎风面、主翼段及襟翼下翼面。在上游部件尾迹的影响下,多段翼表面水滴收集特性呈现显著的动态分布特性。前缘缝翼表面形成的冰角结构及其缝道附近发生的积冰堵塞是导致多段翼构型气动性能退化的主要因素,相比之下襟翼表面结冰对气动特性的影响程度较小。
Select 12. 液体火箭发动机高频纵向燃烧不稳定现象的随机性
汪广旭, 杨宝娥, 谭永华, 高玉闪, 李斌
doi: 10.7527/S1000-6893.2023.29485
2024, 45 (14): 129485-129485.
PDF [2518KB](112)
围绕自燃推进剂液体火箭发动机燃烧室高频纵向燃烧不稳定现象的随机性进行了深入探讨,分析了结构振动数据在表征燃烧室压力振荡方面的有效性,对比了理论声模态和实测振型,基于蒸发过程作为速率控制的假设,对该燃烧室的稳定性进行了分析,并提出了造成不同供应系统下稳定性差异的2种假设模型。研究表明:由结构振动传感器测得的是燃烧室结构对燃气压力振荡场的响应结果,后者存在一定的频率选择性和非线性;实测突频幅值分布与该频率下的理论声模态基本吻合,确定该燃烧室发生了高频一阶纵向燃烧不稳定现象;线性模型结果表明,该燃烧室具有一定的稳定性裕度,但不能解释系统2条件下燃烧不稳定现象随机出现这一现象;系统2条件下该燃烧室的稳定性恶化或不同工况间的稳定性波动幅值增加,造成部分工况点落入非线性不稳定性区间,成为潜在的不稳定工况。
Select 13. 强旋翼尾迹涡干扰对直升机空速测量的影响
刘永华, 徐国华, 史勇杰, 杜振宇, 张慧鑫, 胡志远
doi: 10.7527/S1000-6893.2024.29536
2024, 45 (14): 129536-129536.
PDF [9956KB](46)
准确的空速测量对于直升机飞行安全至关重要。本文建立了基于非定常动量源理论的旋翼/机身气动干扰流场数值模拟方法,研究了旋翼尾迹涡干扰对直升机空速测量的影响。以某参考直升机为例,分析了不同旋翼拉力下,空速管布置区在全速度范围内空速值的变化特性。由于旋翼桨叶运动、机身阻滞等影响,不同位置空速管的静压相位和峰值存在差异,其中相位差异受到桨叶激励时刻影响、峰值差异受桨盘载荷与机身阻滞影响。研究也发现并解释了在特定速度段存在的由于旋翼桨尖涡通过空速管区域引起的“空速迟滞”现象。空速管的空速变化可分为下洗流主导区、迟滞区、来流主导区3个区域,对于常规构型直升机迟滞区仅在中低速度段。进一步的参数研究表明:空速管的位置越靠近机身前部、旋翼拉力越小,迟滞区的速度区间越靠前;当靠近机头时空速管并无明显的“迟滞”。最后,本文从气动干扰的角度给出了空速管的布置建议。
Select 14. 马赫数6飞行条件圆截面超燃冲压发动机流动燃烧特征分析
于江飞, 汤涛, 闫博, 汪洪波, 杨揖心, 熊大鹏, 孙明波
doi: 10.7527/S1000-6893.2023.29575
2024, 45 (14): 129575-129575.
PDF [4824KB](170)
本文基于RANS方法对马赫数6飞行条件下的乙烯燃料圆截面超燃冲压发动机燃烧流场开展了数值研究,计算中采用了压力相关的火焰面/进度变量模型,得到的数值结果与试验吻合,并据此分析了0.18和0.06当量比工况的火焰结构和燃烧特征。计算表明,0.18当量比工况释热强度高,在流场中形成了热壅塞,并且预混燃烧和扩散燃烧共同主导了整体的释热。而0.06当量比工况流道中心仍是超声速状态,扩散火焰主导了整体的释热。2个工况中,射流背风区和轴对称凹腔都是重要的反应区,其中射流背风区释热强度高,而凹腔起到了增强混合和提高燃烧效率的重要作用。
Select 15. 双油箱燃油热管理系统性能分析
魏宇豪, 葛玉雪, 赵倩, 裴扬
doi: 10.7527/S1000-6893.2023.29629
2024, 45 (14): 129629-129629.
PDF [2621KB](45)
随着飞机热载荷的增加,燃油热管理逐渐成为现代飞机设计和运行中的重要考虑因素。建立了考虑热损失的双油箱拓扑结构燃油热管理系统的动力学模型,给出了考虑热损失情况下的管理策略,并利用热续航时间、热沉利用效率和废热利用效率这3个指标对该系统的性能进行评价。结合典型示例,对所推导的动力学模型进行了验证。与传统单油箱系统相比,双油箱拓扑结构燃油热管理系统可以提升系统的热续航时间。对于本文所研究的双油箱燃油热管理系统,其增幅可达36.8%。在此基础上,分析了热管理系统设计参数对运行效果的影响。结果显示:总燃油质量流率的降低可以增强系统的热沉利用效率和废热利用效率。在油箱总容量不变的前提下,适当地缩小再循环油箱容量,则可以提高系统的热续航时间。对于本文算例,当总燃油质量流率为1 kg/s时,再循环油箱容量减少100 kg,热续航时间提高约109 s。本文研究可为飞机燃油热管理系统设计提供参考。
Select 16. S弯喷管正后向红外辐射特征快速预测模型
是介, 周莉, 史经纬, 王占学
doi: 10.7527/S1000-6893.2023.29639
2024, 45 (14): 129639-129639.
PDF [3905KB](90)
为了在航空发动机排气系统红外隐身设计初期快速计算不同遮挡特性S弯喷管正后向的红外辐射特征,发展了基于图像法和一维流场模型的S弯喷管红外辐射特征快速预测模型。该模型采用图像法将S弯喷管的几何遮挡关系转化为多层二维图像的像素运算,根据可压缩管流和射流理论建立S弯喷管一维流场模型,基于统计窄谱带模型计算燃气吸收发射特征,结合燃气可见区域长度和容积的拟合函数考虑遮挡特性对燃气辐射的影响。通过与离散传递法对比验证了模型的计算效率和准确性。结果表明:快速预测模型可将S弯喷管正后向红外辐射特征计算时长由数小时缩短至秒级,其计算结果与离散传递法在趋势和数值上均吻合较好。不同遮挡特性S弯喷管和不同喷管工况时总红外辐射强度最大相对误差仅为6.5%,壁面和燃气辐射的最大相对误差分别为4.1%和5.2%。快速预测模型针对轴对称和二元喷管也具有良好的泛用性。
Select 17. 两级入轨航天器级间分离姿态精确控制
姜权峰, 陈树生, 杨华, 李祚泰, 高正红
doi: 10.7527/S1000-6893.2023.29270
2024, 45 (13): 129270-129270.
PDF [4164KB](241)
两级入轨航天器(TSTO)级间分离过程中对轨道级、助推级的姿态精确控制尤为重要。数值模拟可以较为准确地预测级间分离这种飞行器多体分离问题,运用非结构嵌套网格技术及双时间步法耦合流动控制方程、刚体运动方程求解TSTO在攻角-2°、马赫数6、海拔高度30 km的级间分离过程,发展了一种优化轨道级、助推级在级间分离过程的运动姿态的技术。该技术通过松弛迭代法提取轨道级、助推级的无舵偏气动力数据,基于无舵偏气动力数据设计气动舵面、燃气舵所需要提供的俯仰力矩,使得轨道级、助推级的运动姿态更加接近设计目标。数值模拟结果表明,TSTO在分离过程中两级之间会存在强烈的气动干扰,严重影响两级的姿态控制,发展的技术可为TSTO分离过程中轨道级、助推级的舵面运动、燃气舵工作状态设计提供参考,提高TSTO级间分离的效率、安全性。
Select 18. 十字型柔性控制面气动控制力风洞试验验证
韩译辉, 胡俊, 于勇, 于剑桥
doi: 10.7527/S1000-6893.2023.29280
2024, 45 (13): 129280-129280.
PDF [7130KB](85)
现代战争中,战场环境复杂多变,需要弹箭具备快速大角度机动转弯能力。然而,传统的气动舵无法满足弹箭类飞行器超90°的高机动转弯需求。本文创新性地提出了一种新型柔性控制机构——大尺度可展开可变构型的十字型柔性控制面,用于产生超量气动控制力,实现弹箭的迅捷机动转弯;同时对一个缩比1∶5的可变构型十字型柔性控制面所产生的气动控制力进行了风洞试验。试验结果显示,通过在不同挂载点位收缩控制绳的方式下,相比于没有构型改变的对称构型十字型柔性控制面,改变构型的十字型柔性控制面可以基本稳定在收缩控制绳的挂载点位一侧,不再出现旋转、摆动,这一结果说明改变构型的十字型柔性控制面可以产生较稳定的径向控制力;对变构型的十字型柔性控制面的气动控制力测量值进行定量分析,进一步确认了柔性面气动控制力的指向性;在来流Re≈2.3×105,控制绳收缩比为0.9的情况下,产生了气动力系数为0.29~0.34的可用径向控制力。研究结果确认了变构型十字型柔性面控制力的可用性,为进一步研究柔性控制力的气动特性、变构型控制策略设计提供了基础。
Select 19. 等离子体合成射流激励器阵列破除翼前缘三维冰特性
程盼, 景向嵘, 罗振兵, 高天翔, 周岩, 邓雄, 孙乾
doi: 10.7527/S1000-6893.2023.29137
2024, 45 (12): 129137-129137.
PDF [4099KB](123)
无人机结冰将严重影响其安全性,机翼前缘更容易结冰及积冰,对无人机气动性能影响更为严重,而机翼前缘结冰形状具有三维特性,亟须发展有效破除机翼前缘三维冰的低能耗、快响应除冰技术。分别根据NACA0012、NACA0018及NACA0024 3种翼型前缘曲率制作了3种不同弧度的冰型,基于等离子体合成射流激励器响应快、能耗低、射流强度高等优点,开展了单个及阵列布置激励器破除三维形态冰特性研究,分析了45°和90°两种不同方向出口对破冰特性的影响,阐明了等离子体合成射流激励器的破除冰机理及规律。结果表明:放电总能量为16.66 J时,在单个激励器作用下,对于翼前缘类似NACA0012的曲率半径较小的冰,射流冲击应力易在前缘驻点处集中,在翼前缘驻点线处能产生长达20 cm的贯穿裂纹,从而实现有效破除冰;对于曲率半径较大的冰,破冰半径约5 cm。保持放电总能量16.66 J不变,激励器采用沿驻点线的“一字型”阵列布置时对3种弧度的冰都能产生贯穿裂纹实现破除冰,表明激励器阵列能有效拓展破除冰面积,具备低能耗、大面积破除冰能力。研究成果可为无人机低能耗、大面积破除冰提供理论及实践参考。
Select 20. NNW⁃ICE软件的三维结冰模型及精度验证
陈宁立, 易贤, 王强, 任靖豪
doi: 10.7527/S1000-6893.2023.29188
2024, 45 (12): 129188-129188.
PDF [3143KB](108)
结冰现象严重危害飞行器的飞行安全,而数值模拟是研究结冰的重要手段之一。NNW-ICE软件是一款由中国空气动力研究与发展中心开发的中国完全自主可控的结冰数值模拟软件。该软件中的三维结冰相变模型相比于传统的SWIM、Myers模型增加考虑了表面接触角对水膜流动的影响,且能量方程中增加考虑了伴有蒸发现象的霜冰结冰。数值计算方法方面,软件采用有限体积法对方程进行离散,拓展了模型对于复杂几何构型的非结构网格的适用性,并且集成了一阶显式、四阶四步龙格库塔-显式、一阶隐式、二阶一步龙格-库塔隐式共4个求解器。本文还采用了冰风洞试验结果、商业软件FENSAP-ICE计算结果对软件精度进行了验证。结果表明:NNW-ICE软件的计算结果与冰风洞试验结果吻合均较好,且精度总体优于FENSAP-ICE软件。
学术质量
下载排行