二维机翼混合相结冰数值模拟

近年来冰晶引发的飞机结冰问题逐渐引起人们的重视。针对冰晶和过冷水滴同时存在的混合相结冰问题,通过数值模拟的方法实现结冰冰形预测。空气流场和对流换热的计算中采用了转捩剪切应力输运（SST）湍流模型,基于欧拉法计算冰晶和水滴收集系数。在Messinger结冰热力学模型基础上进行扩展,分析了二维结冰部件表面在混合相气象条件下的传热传质过程,建立了适用于混合相的结冰热力学模型,同时考虑冰晶的黏附效应,添加黏附效率经验公式。利用FLUENT的用户自定义函数（UDF）编程求解混合相热力学模型,计算了霜冰和明冰条件下NACA0012翼型表面结冰情况,与国外风洞试验结果进行比较,验证了计算模型和方法的有效性。结果表明,冰晶黏附效应对混合相结冰量及冰形有很大影响,明冰条件结冰形状偏向楔形。     

飞机防冰与除冰的若干技术

阐述了飞机表面“水滴结冰”、“干结冰”、“升华结冰”等结冰现象、结冰形式及其影响因素；总结了现有的飞机液体、电热、气热防冰与气动带、电脉冲除冰技术；还对冰风洞、干空气飞行、模拟结冰飞行、结冰模拟以及自然结冰等飞机防冰试验进行了综述。     

民用飞机机翼结冰试验与数值预测

描述了某民用飞机在意大利IWT（icing wind tunnel）冰风洞的机翼结冰试验结果.根据冰风洞试验状态点,开展了基于拉格朗日方法和Messinger模型的冰形数值预测，并对试验冰形和数值预测冰形进行了对比分析.霜冰的数值预测结果与试验结果一致性良好,对于明冰,数值结果和试验结果存在一定的差异,原因在于:①明冰结冰过程中较强的绕流变化增大了冰积聚非定常效应对数值预测的影响；②采用的结冰热力学模型对明冰的模拟存在一定局限性；③明冰对于冰风洞设备工况和试验模型表面粗糙度的影响更为敏感.      

航空发动机冰晶结冰研究进展

飞机在高空区域巡航时可能遭遇冰晶气象，导致发动机吸入冰晶颗粒发生结冰故障，欧美等国已将冰晶结冰纳入适航 取证要求，中国也考虑在后续发动机适航取证中增加冰晶结冰考核项目。为了准确掌握国内外冰晶结冰研究现状，通过调研文 献，对比了冰晶结冰与过冷水滴结冰在结冰位置、结冰条件和结冰机理方面的区别，阐述了试验研究和数值仿真方面的重要进展， 其中，试验研究分别介绍了冰晶在静态部件和动态部件内的撞击与结冰过程，数值仿真按照空气流场计算、冰晶运动轨迹和撞击 特性计算及结冰计算的顺序，总结了欧拉法与拉格朗日法的区别、冰晶运动相变与黏附过程及结冰热物理模型，并介绍了冰晶结 冰软件在数值模拟中的广泛应用，提出了目前在冰晶碰撞、结冰热力学模型、结冰探测及风洞试验等方面需要进一步研究的内容， 以期为中国冰晶结冰研究提供参考。     

移动式冰风洞试验方法研究和应用

为了探索掌握移动式冰风洞校测、结冰和防/除冰试验的一般方法,使其可用于飞机进气系统地面防冰试验,开展了对移动式冰风洞模拟云雾参数的校测。通过金属圆柱管结冰对比确定水滴过冷距离,采用格栅测量云雾均匀性,利用机载雾滴组合探测器测量水滴直径与液态水含量,并对发动机短舱唇口模型、NACA23012翼型模型进行了结冰和防/除冰研究试验。研究表明虽然户外模拟结冰条件受环境因素影响较大,但移动式冰风洞喷雾性能良好,试验能够反映在各因素影响下的一般结冰规律,可以满足飞机进气系统防/除冰试验要求。     

基于欧拉方法的2 维翼型冰晶结冰数值计算

针对现有航空发动机冰晶结冰试验方法存在数据采集困难,且难以揭示冰晶结冰机理问题,对某楔形翼型进行冰晶结 冰数值计算。采用S-A 湍流模型获得流场结果,应用欧拉方法获得冰晶和液滴轨迹结果,基于Messinger 模型获得冰形,并以 NASA-NRC 的第139 号试验结果验证了数值方法的可行性。结果表明：在其他工况相同时,压力越低（飞行高度越高）,冰晶/ 液滴蒸 发越强烈,结冰情况越严重;分析不同液态水含量与总水含量比例（CLW/CTW）对结冰的影响可知,液态水过少或过多都不利于冰层形成。     

航空发动机内部冰晶结冰研究综述

为深入理解航空发动机内部冰晶结冰现象,把握冰晶结冰过程的特征规律,建立科学有效的冰晶结冰防护手段,对近年来冰晶结冰文献资料进行了调研和总结。回顾了冰晶结冰现象的发现与证实的过程,重点阐述冰晶结冰与过冷水结冰的区别、冰晶结冰模拟试验及计算方法,并对冰晶结冰问题有待进一步研究的方向做了展望。     

翼面结冰过程中的冰晶运动相变与黏附特性

冰晶的运动、相变、黏附现象广泛存在于航空发动机内部结冰过程中，开展相应的计算方法研究，系统分析冰晶运动相变与黏附特性是研究发动机内部结冰、保证飞行安全的迫切需求。针对冰晶运动相变及黏附特性计算，建立了拉格朗日框架下冰晶运动-传热传质耦合的数值计算方法，分析了基于Monte Carlo方法的冰晶撞击与黏附收集系数计算方法，基于NNWICE平台开发了相应的计算程序。以NACA0012翼型为对象，计算研究了扁椭球状、六角平板状和长椭球状3种非球状冰晶颗粒与球状冰晶的运动、传热传质过程，系统分析了来流温度等参数对冰晶黏附特性的影响，得到了冰晶形状对冰晶运动轨迹和融化过程的影响规律及来流总温与液态水含量/总水含量（LWC/TWC）对冰晶黏附特性的影响规律。相关工作为进一步分析混合相结冰数值模拟计算奠定了工作基础。     

大型结冰风洞中冰晶热/力平衡特性数值研究

为明晰大型结冰风洞中冰晶热/力平衡特性,发展了基于欧拉法的冰晶运动和传热传质耦合计算方法,模拟了典型大型结冰风洞构型内冰晶运动和传热过程,从沉降收缩、动量平衡和热平衡三个方面,考察了颗粒形状和体积密度的影响.结果表明:颗粒形状和体积密度对小尺寸冰晶的热/力平衡特性无显著影响.降低颗粒球形度和体积密度会抑制大尺寸冰晶沉降...     

航空发动机冰晶结冰适航审定要点及符合性方法研究

为深入了解航空发动机冰晶结冰现象,开展冰晶结冰适航审定要点及符合性方法研究,回顾总结了冰晶结冰事件及相关研究内容,阐述了冰晶存在环境及对航空发动机的危害性。结合适航规章要求、冰晶结冰事件及已发布的适航指令,给出了冰晶结冰适航审定的关注要点。依据当前冰晶结冰计算模拟及地面试验研究现状,提出了一种表明航空发动机冰晶结冰符合性的方法。梳理了冰晶结冰符合性研究所面临的主要问题,并对后续研究方向进行了展望。     

飞机机翼结冰过程的数值模拟

提出了一种适用于各种结冰气象条件下飞机机翼结冰过程的数值模拟方法。结合欧拉坐标系下空气-过冷水滴两相流动控制方程的计算,对不同气象条件下飞机机翼的结冰过程进行了数值模拟,得到了4°攻角NACA-0012翼型表面霜状冰、混合冰及瘤状冰3种冰形。所得冰形与文献中的试验数据和计算结果吻合良好,表明本文方法有效。另外还计算了不同的冰形对机翼升、阻力的影响,结果表明瘤状冰对机翼气动特性的影响最大。      

基于旋转多圆柱的冰风洞水滴参数分析方法

阐述了基于旋转多圆柱测量仪的冰风洞水滴参数分析方法.对二维圆柱的水滴撞击特性进行了计算分析;对二维旋转圆柱的霜状冰结冰过程进行了数值模拟;根据水滴参数对旋转圆柱结冰量的影响, 结合数学逼近的方法, 编制了由旋转多圆柱测量结果分析计算冰风洞水滴参数的软件, 最后用国外试验数据对软件分析结果进行了验证, 结果表明:计算方法正确, 可以实现冰风洞水滴参数快速测量.      

3m×2m结冰风洞试验技术新进展（2020-2022年）

3 m×2 m结冰风洞是我国“十一五”国家重大科技基础设施，也是国际上尺寸最大的非季节性结冰风洞。自2013年建成以来，已经完成了70余项试验，有力支撑了我国飞机的自主研制和适航取证。本文首先介绍了3 m×2 m结冰风洞的组成和特点，其次重点阐述了2020年至2022年间风洞试验能力和试验技术的若干新进展，通过发展双闭环自适应温度控制技术、多路热气供气防除冰试验技术、冰形在线测量技术、发动机进气精确模拟技术、旋翼结冰与气动载荷同步测试技术等，使风洞的温度场模拟能力、热气防冰试验能力、冰形测量能力、进气模拟能力和直升机旋翼结冰试验能力得到增强，综合试验效率显著提升。最后，针对大型结冰风洞过冷大水滴试验面临的挑战，对下一步试验技术的发展进行了展望。     

爆轰驱动膨胀管性能研究

超高速流动一般指速度超过5 km/s的流动,由于流动具有高焓高速的特点,模拟超高速流动的地面试验设备面临极大挑战。膨胀管(风洞)是少数几种具备超高速流动模拟能力的地面试验设备之一。中国科学院力学研究所高温气体动力学国家重点实验室(LHD)通过将正向爆轰驱动技术和膨胀管结合在一起,建成了可实现最高速度10 km/s超高速试验气流的爆轰驱动膨胀管(JF-16),并开展了典型模型试验。在此基础上对JF-16进行了改造升级工作,为其设计喷管增加了膨胀风洞运行模式,对其性能进行了相关试验测试研究。同时,对膨胀管相关数值方法进行了介绍,并开展数值模拟对试验状态进行辅助诊断和分析。     

混合相态冰晶积冰的数值研究

在发动机内流的高温作用下,所吸入冰晶会部分融化为液态水,冰水混合相态条件下发动机内部表面会形成积冰,冰晶积冰会导致发动机喘振、熄火,甚至会由于冰脱落而造成内部结构损伤。为了对混合相态条件下冰晶积冰问题进行深入研究,以NACA0012翼型为对象,通过数值计算分析研究了环境温度、马赫数等参数对积冰形态、收集系数以及积冰生长率的影响,并分析了融化率对积冰过程的作用机制。结果表明:混合相结冰条件下若达到最大结冰厚度,需满足有足够的冰晶和液态水含量条件;环境温度直接影响了湿球温度变化,而随环境温度升高,液膜的厚度和润湿范围也随之增大。此外降低环境温度或增大马赫数,翼型前缘驻点处结冰量和积冰速率均有明显增加。     

民用航空涡扇发动机进气系统结冰高空模拟试验技术要求研究

高空模拟试验是民用航空发动机适航取证的符合性验证手段之一。基于中国民航规章CCAR-33部33.68条款“进气系统结冰”的要求，对进气系统结冰高空模拟试验验证开展研究，以表明符合性。本文提出了进气系统的结冰试验点的选取方法、分析了试验的流程、明确了试验结果评估方法，并针对试验设备提出了高空舱、测试装置和试验设备配备要求，介绍了典型符合性验证案例。本研究可为国内民用航空发动机开展进气系统结冰适航符合性验证工作提供技术支持。     

旋转帽罩电加热防冰瞬态过程研究

基于欧拉两相流理论,对旋转帽罩水滴撞击特性进行了数值模拟,提出了旋转部件表面对流换热计算方法,并基于改进Messinger结冰模型开发了旋转帽罩电热防冰计算程序,对旋转帽罩瞬态防冰过程进行了数值模拟研究.结果表明:旋转运动对旋转帽罩表面对流换热起到增强作用,且在供给相同加热热流密度时转速越大,防冰表面温度越低;防冰系统启动阶段,旋转帽罩表面会发生结冰和冰脱落现象,考虑结冰过程后系统的响应时间缩短;当电加热功率相同时,周期电加热防冰方式更为节能;当电加热能耗相同时,周期电加热方式系统响应更快.     

大型结冰风洞云雾场适航应用符合性验证

结冰风洞云雾场符合性是大型结冰风洞适航应用的基础。为验证3 m×2 m结冰风洞云雾场符合性，发展了基于SAE ARP5905的结冰风洞云雾场符合性验证方法，针对主试验段构型，开展了结冰云雾场符合性验证试验，获得了试验段内液滴尺寸和液态水含量拟合关系，考察了喷嘴水压、液滴尺寸、试验段气流速度和喷嘴数量对试验段液态水含量的影响，形成了主试验段结冰云雾控制包线。结果表明：主试验段内液滴尺寸分布具有显著的单峰分布特征，体积中值直径（MVD）模拟范围近似在10~75 μm之间；试验段中心处液态水含量随着喷嘴水压和MVD的增大而增大，同时近似与试验段气流速度成反比，与喷嘴数量成正比；增大喷嘴水压和喷嘴数量会提高试验段内云雾液态水含量空间均匀性，但是增大气流速度却会减弱试验段内云雾空间均匀性；3 m×2 m结冰风洞主试验段结冰云雾控制包线可以覆盖大部分适航条例25部附录C结冰气象条件，但对低液态水含量结冰条件的模拟仍存在局限。     

结冰风洞小流量发动机进气部件防冰试验技术

为满足国产航空发动机小流量进气防冰试验的需求,基于结冰风洞的发动机进气模拟系统,分析了典型小流量发动机进气部件的进气特征,开展了小流量发动机进气模拟方法研究。提出采用水环真空泵解决了大流阻进气导致吸气管道压力大幅度降低问题,在此基础上,优化了发动机进气模拟系统,建立了在结冰风洞开展小流量发动机进气部件防冰试验的流程,并应用于某型号小流量发动机进气部件防冰试验。结果表明：发动机进气模拟系统满足小流量发动机进气部件防冰试验要求,进气流量控制精度达到±0.05 kg/s（±0.33%FS）,试验过程中的流量变化可用于辨别防冰效果。相关研究为发动机进气模拟系统设计、优化和小流量发动机进气部件的防冰特性考核提供了参考依据。     

瘤状冰结冰过程的数值模拟

提出了一种可考虑粗糙度影响的瘤状冰结冰过程数值模拟方法。结合欧拉坐标系下空气-过冷水滴两相流动控制方程的计算,对飞机机翼在不同温度下的瘤状冰结冰过程进行了数值模拟,得到了NACA0012翼型4°攻角下结冰表面的瘤状冰冰形。通过与实验数据对比,表明该的方法是有效的。