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

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

飞机机翼结冰计算方法与风洞试验研究概况

介绍了飞机结冰的数值计算方法,包括绕流流场计算、过冷水滴轨迹方程计算和结冰量计算,并给出冰层增长的两种模型;针对结冰风洞试验,指出液滴粒径和液态水含量的测量是试验的难点和重点,并总结上述两个结冰参数的测量方法;给出两种常规结冰缩比方法,并简单介绍混合缩比模型的设计方法。     

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

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

混合相/冰晶条件下的结冰研究综述

针对混合相/冰晶结冰问题，阐述了与过冷水结冰的差异及其危害。从混合相/冰晶结冰地面模拟试验设备、试验现象以及数值模拟3个方面详细论述了混合相/冰晶结冰的研究进展。其中，地面模拟试验设备以Cox冰风洞、NRC高空试验设备、NASA推进系统试验室、德国布伦瑞克冰风洞为例，总结了各自的冰晶生成方法和试验能力；分别介绍了以上4个地面模拟设备中开展的混合相/冰晶试验，分析了其试验现象的异同和产生的原因；针对混合相/冰晶结冰数值模拟中面临的4个关键问题：冰晶运动融化相变、冰晶黏附、侵蚀、结冰过程模拟，总结了现有模型和模拟方法，提出了不足之处。最后，对混合相/冰晶结冰有待进一步重点关注的研究方向进行了展望。     

霜状冰结冰数值模拟研究

提出了一种三维霜冰结冰的数值模拟方法。基于欧拉-拉格朗日法计算空气-水滴两相流流场,获得了表面水滴撞击特性以及撞击量,求解质量守恒和能量守恒方程获得网格单元内的结冰量,三维霜冰的结冰计算程序采用用户自定义（UDF）编写,FLUENT提供的动网格功能实现了结冰后的网格重构。计算了4°攻角,液态水含量为1 g/m3,水滴直径20μm,结冰时间6 min,温度分别为-19.4°和-28.3°时NACA0012翼型的结冰,并与相关文献结果进行了对比,计算结果与文献实验结果吻合较好,验证了方法的合理性和可靠性。分析了不同的粒径、液态水含量、来流速度等对撞击特性以及结冰的影响,随着粒径、来流速度的增加,水滴收集系数、撞击区和结冰厚度增大,随着液态水含量的增加,水滴收集系数和撞击区是不变的,但是结冰厚度增加。     

飞机机头附近液态水含量分布规律的数值研究

在飞行器飞行过程中,安装在机头上的各类传感器结冰会对飞机的飞行安全产生不利影响。而液态水含量是影响飞机结冰的重要参数,为了更好地进行传感器结构及布局设计,对机头附近液态水含量分布规律进行数值研究。采用计算流体力学方法对飞机飞行时的空气-过冷水滴的稀疏两相流场进行数值模拟,通过求解N-S方程获得空气流场,采用欧拉法求解水滴运动场得到水滴运动轨迹。分析在不同来流马赫数、来流迎角及液滴直径等条件下,机头关键截面上广义水滴遮蔽高度的变化及机头关键位置上液态水含量的变化规律,研究不同状态参数对液态水含量分布的影响。结果表明:水滴遮蔽高度及浓度增加区范围随马赫数增大而增加,迎角的变化对机头上下表面液态水含量分布产生相反的影响。     

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

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

一种基于水滴收集量代理模型的结冰试飞空域确定方法

针对自然结冰试飞空域确定时使用的结冰指数只能给出结冰概率和结冰等级的问题，提出一种新的方法。通过对美国联邦航空条例（FAR）25部附录C连续最大结冰条件采样，对采样点进行空气流场和水滴撞击特性求解，获得不同工况的水滴收集量；基于POD (Proper Orthogonal Decomposition)和Kriging构建水滴收集量代理模型；使用WRF(Weather Research and Forecasting)对目标区域进行气象模拟，获得温度以及液态水含量分布；使用代理模型对目标区域内水滴收集量进行预测，以中度结冰强度对目标区域进行划分；最后，针对2种飞行速度对试飞空域的影响进行研究。结果表明：代理模型能够很好地预测温度、液态水含量、水滴中值体积直径、高度以及速度对水滴收集量的影响；WRF获得的目标区域的温度、液态水含量与观测值符合良好；基于代理模型可快速获得目标区域水滴收集量分布及随时间的变化，还可获得适合自然结冰试飞的目标区域及结冰速度；飞行速度的增加使得水滴收集量增加，进而引起试飞空域的变化。本文对结冰试飞空域确定具有一定参考意义。     

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

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

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

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