基于变分分割模型的结冰冰形测量方法

结冰冰形是结冰风洞试验中关注的要素之一,针对其冰形测量问题,提出了一种基于变分图像分割技术的非接触测量方法。基本思路是对翼型结冰俯视图像进行图像分割,将获得的曲线标定后与翼型数模结合,得到最终结冰冰形。为了排除结冰图像中不相干物体和因素的干扰,更精确获取结冰冰形,提出了选择分割的思想,构造了区域特征函数,建立了新的能量函数。采用所提基于变分分割模型的测量方法,对结冰风洞试验冰形进行了测量,并定量分析了方法的误差,测试了方法的抗噪声能力。结果表明该方法对结冰冰形的测量是准确可行的,且具有较高的精度,另外,其对噪声具有鲁棒性。所提方法可推广至其他与物体形状相关的测量中。     

风力机叶片翼型的结冰数值模拟研究

研究了风力机叶片翼型的结冰数值模拟方法,及翼型结冰后对其气动性能的影响。求解雷诺平均N-S方程,引入k-ωSST湍流模型封闭方程,获得风力机叶片流场;采用拉格朗日法计算风力机叶片翼型周围的水滴运动轨迹,同时考虑了多尺寸水滴分布的影响,获得翼型表面的局部水收集系数分布;根据质量守恒和能量守恒原理,计算翼型表面的各项热流,获得翼型表面的结冰速率和结冰冰形;考虑到翼型结冰是时间的动态函数,采用多时间步长法完成结冰数值模拟。计算了风力机叶片翼型在不同环境条件和气象条件下的结冰冰形,同时模拟了风力机翼型结冰后周围流场的变化,并与干净翼型的气动特性进行了对比。结果表明,环境温度较高时形成的明冰对翼型气动性能的影响较大,结冰导致翼型升力下降,升阻比减小,最大减小幅度达到61%,同时结冰后的翼型会提前进入失速区,导致桨叶气动性能恶化。     

二元翼型结冰数值模拟研究

对二元翼型前缘结冰数值模拟进行了研究。采用中心有限体积法求解N—S方程；四阶龙格一库塔法求解水滴运动轨迹方程；逐次二分法寻找水滴撞击极限，然后计算局部水滴收集率；分析翼型表面的传质、传热项，建立相应质量、能量方程并迭代求解，得到翼型表面结冰质量；根据角度二分法生成结冰冰形。预测的NACA0012结冰冰形与实验及Lewice模拟结果吻合很好，结冰过程中的特征参数与冰型及相应冰形特征吻合。     

民用飞机翼面冰形的工程修正

翼面结冰是威胁飞行安全的重要因素之一。以数值模拟结冰为基础,进行了翼型结冰的工程修正方面的研究。运用结冰分析软件Sadrice,对不同状态下的不同翼型进行表面冰形数值模拟。结合结冰风洞试验结果,对影响冰形的参数进行了变参研究,得到了最有效的修正方法。对现有对流换热系数经验公式进行线性修正,使修正对流换热系数经验公式后的预测冰形和冰风洞试验结果之间能够较好地吻合,最终总结出了一套工程修正数值模拟冰形的方法。     

二维翼型结冰数值计算

翼型结冰计算是随着结冰时间步长变化的准定常过程,因此在每个时间步长后都要更新各个模块的计算工作,如此反复直到达总结冰计算时间.在各模块计算中,分别采用了泊松方程法生成网格、求解雷诺时均方程计算流场、拉格朗日追踪法求解水滴轨迹、基于Messinger模型计算冰形、采用变距滤点法重构边界.引入动态链接库技术将各模块集成为结冰计算平台,并将三种冰形(楔形冰、槽状冰及混合冰)计算结果与国外软件及试验数据进行了对比.      

基于多模态融合的任意对称翼型结冰预测方法

为解决目前绝大多数神经网络冰形预测方法只能针对特定翼型且不具备面向多翼型特征的普适性的问题，采用基于多模态融合的深度神经网络方法，以翼型截面图像与结冰工况参数作为输入，以二维冰形曲线傅里叶级数拟合参数作为输出，建立深度神经网络预测模型，实现了对任意对称翼型结冰特征的预测能力。结果表明：提出的模型可以准确地预测任意对称翼型几何特征条件下的结冰外形，冰形面积与最大冰厚等冰形主要参数预测误差均保持在10%以下。     

临界冰形确定方法及其对气动特性影响研究

临界冰形是指在适航规章结冰包线内,每个可用飞行构型下,对飞机操纵性和稳定性影响最严重的冰形,临界冰形分析是飞机适航取证中的重要工作。对临界冰形确定方法,及临界冰形对气动特性的影响进行了研究。发展了基于 CFD 方法计算临界冰形的一般方法,包括临界冰形分析状态、敏感性分析截面确定、结冰参数敏感性分析、临界结冰条件确定、临界冰形确定等。流场计算采用中航工业空气动力研究院气动力计算平台(UNSMB),基于 Jameson 中心格式的有限体积法求解 N-S 方程;水滴撞击特性计算采用 Eulerian 方法求解水滴轨迹运动方程;结冰计算采用经典的 Messinger 热力学模型。选取 CRM 飞机为研究对象,以机翼外翼50%展长处为敏感性分析截面,在典型飞行条件下,分析了结冰对环境温度、水滴直径、飞行速度、飞行迎角等参数的敏感性。利用“几何外形敏感性分析方法”,即通过对比冰形的上下冰角角度和冰角厚度等冰形几何参数来确定最严重冰形,得到了CRM 飞机的临界结冰条件和临界冰形,其中敏感性分析截面在水滴直径为30μm 时上冰角厚度和下冰角厚度最大,冰角最大厚度约41 mm。计算了结冰后的气动性能衰减规律,临界冰形对飞机气动性能影响严重,导致升力降低6.7%~23.8%,阻力增加17%~70.9%。发展了45min 待机临界冰形确定方法,基于几何外形敏感性分析方法进行环境温度、水滴直径、飞行条件等各类参数的结冰敏感性分析,得到飞机的临界结冰条件和临界冰形,对于民用飞机设计和适航取证具有一定的工程应用价值。     

冰形表面粗糙度对翼型的失速特性影响分析

结冰翼型气动影响分析是飞机结冰问题重要的研究内容之一,可帮助进行飞机结冰安全分析和容冰设计。为填补现阶段关于冰形表面粗糙度对翼型气动影响的研究空缺,以冰形的粗糙度参数(冰形表面粗糙高度和冰形表面粗糙分散度)为依据,提出了一种冰形粗糙度叠加方法;并通过有限体积法求解空间离散的RANS方程,结合S-A湍流模型,分别针对流线冰型和角状冰型,量化分析了不同微观特征参数的冰形对翼型失速特性的影响规律,为研究飞机结冰问题提供了思路和基础。结果表明,对于流线结冰翼型,冰形表面粗糙度会明显影响翼型的失速性能;而对于角状结冰翼型,两个凸起的羊角状外形才是导致气动力损失的主要因素,粗糙度的影响往往可不考虑在内。     

翼型表面结冰准定常数值模拟

对翼型表面的结冰过程进行了准定常数值模拟。在每个时间步长内完成网格随着壁面边界的移动而更新、周围流场和水滴撞击特性重新计算、冰形计算及壁面边界的重构工作,如此循环直至所需的结冰计算时间。采用拉格朗日轨迹追踪法获得水滴撞击特性,应用考虑壁面粗糙度影响的边界层积分法计算壁面的对流换热系数,在此基础上求解Messinger结冰热力学模型,冰层始终沿壁面外法线方向生长。对NACA0012翼型在不同环境下生成的三种典型冰形进行了预测,并与实验结果进行了比较,表明本文所述方法的有效性。     

基于卷积神经网络的结冰翼型气动参数预测

针对传统CFD计算获取结冰翼型气动参数过程繁琐的问题,提出一种基于数据转换和深度学习的图形化预测方法,开展结冰翼型气动特性参数的预测研究。该方法利用卷积神经网络可有效提取图形特征的特点,通过数据转换提取翼型结冰图形,建立结冰翼型气动参数预测模型,通过训练网络建立冰形与气动参数之间的非线性映射关系。最后,以NACA0012翼型为例开展仿真试验,证明了方法的可靠性和高效性。研究结果表明,利用该方法预测气动参数平均相对误差在5%以内,预测时间在毫秒量级,可为后续飞机结冰飞行安全分析和防/除冰设计提供参考。     

过冷大水滴条件下结冰相似准则

虑到中国现有结冰风洞在过冷大水滴(SLD)结冰条件模拟能力方面的局限,发展SLD结冰风洞及其相应的缩比方法迫在眉睫。研究了大水滴运动机理,并依据绕流流场相似、水滴运动和撞击特性相似、结冰过程热力学特性相似原理,提炼出SLD结冰条件下的6种相似缩比准则,利用结冰数值模拟软件FENSAP-ICE对所提相似缩比准则进行筛选与比较,从而基于结冰相似百分比量化判断准则找出给定条件下的最优准则法。结果表明:考虑SLD运动特性的撞击参数与ONERA法相耦合的方法以及撞击参数与Ruff法相耦合的方法的计算结果与实际成冰构型吻合良好；此外,撞击参数与ONERA法相耦合的方法在计算直径为171 μm的SLD时优于其他。      

运输机翼型结冰的计算和实验

采用数值计算和结冰风洞实验相结合的手段,对某运输机机翼剖面缩比翼型的结冰特性进行了研究.结冰实验在气动中心0.3m×0.2m结冰风洞中进行,实验模型为弦长0.18m的层流翼型.翼型的绕流流场通过求解低速黏流的时均Navier-Stokes(N-S)方程得到,采用拉格朗日法计算水滴撞击特性,在此基础上,求解结冰热力学模型...     

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

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

引射式结冰风洞内圆柱结冰试验

设计并建造了用于结冰研究的小型引射式结冰风洞试验系统,对其结冰气象条件参数进行了标定,给出了试验段水雾均匀区范围,验证了试验冰形的可重复性.在霜状冰和瘤状冰条件下,对直径为25.4mm和50.8mm的圆柱在所建的结冰风洞内进行了结冰试验.试验结果表明:在所研究的结冰条件和结冰时间范围内,霜状冰与瘤状冰的结冰厚度均随时间推进而线性增长;在满足结冰相似准则的情况下,冰形相似性结果良好,无量纲厚度的最大偏差约为10%.      

航空发动机支板热滑油防冰性能试验

在冰风洞内开展了结冰条件下涡轴发动机进气支板的热滑油防冰系统的防冰性能试验研究。试验设计加工了滑油电加热系统,采用可编程逻辑控制器（PLC）监控滑油的温度和流量。在冰风洞中采用全尺寸模型开展滑油防冰性能试验,所开展的涡轴发动机支板热滑油防冰试验参数包括：来流温度为-10,-5℃,来流速度为40 m/s,液态水含量为0.5,1.0 g/m³,过冷水滴平均体积直径为20 μm。试验开展了不同结冰气候条件下、不同滑油通道位置滑油防冰进气支板防冰效果的研究,记录了支板表面温度的变化和结冰情况。试验同时得到了支板防冰能力不足时支板表面的结冰冰型和结冰环境下发动机支板热滑油防冰的特点。     

涡扇发动机短舱结冰试验相似方法

通过分析短舱表面结冰缩比试验相似要求,在保证绕流流场、水滴撞击特性和结冰冰型相似的基础上,分别采用基于气流雷诺数,气流韦伯数和水滴韦伯数3种速度选取方法建立了试验相似准则。利用该准则获得了明冰和霜冰条件下涡扇发动机短舱1/2缩比模型的结冰试验参数,开展了短舱进气道表面结冰冰型预测,对参考模型和缩比模型表面过冷水滴撞击特性、结冰特性、溢流水分布特性进行了分析。结果表明:采用基于气流韦伯数和水滴韦伯数速度选取方法的相似准则,在明冰和霜冰条件下均能保证缩比模型与参考模型表面水滴撞击特性相似,且相似性不受结冰温度影响；溢流水变化趋势与参考模型相似,流动极限相近；缩比模型表面冰型均与参考模型相似,该方法能够为短舱结冰试验参数选取提供依据。      

飞机结冰的不确定性量化研究进展

为深入认识影响飞机结冰的不确定性及其研究方法,从自然条件结冰、冰风洞试验、数值模拟等方面,介绍了飞机结冰不确定性来源,以蒙特卡洛法、多项式混沌法和随机配置法等为例,系统分析了各种不确定性量化方法在计算能力、求解精度等方面的优劣。考虑飞机结冰不确定性量化处于起步时期,重点综述了数值模拟中结冰条件不确定性对冰形和气动特性的影响。从结合单步法和多步法确定最佳结冰时间步长来提升结冰计算精度和效率、对其他关键部件进行结冰不确定性量化从而为更精细化的防/除冰系统的设计提供支撑,以及通过建立高精度的代理模型来代替原本复杂的数值模拟系统以应对考虑多个不确定性因素共同作用所带来的计算挑战等多个方面,全面展望不确定性量化方法及其在飞机结冰应用中的发展方向。     

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

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

结冰云雾参数对冰与固壁间剪切强度影响的初步研究

基于结冰风洞实验段,建立了冰与固壁间黏附界面剪切强度测量实验装置,获得了合理的剪切强度范围。在0.3 m×0.2 m结冰风洞实验段内开展了冰与固壁间黏附界面剪切强度的测量实验。掌握了结冰环境温度、平均水滴直径等结冰云雾参数对冰与固壁间黏附剪切强度的影响规律。研究表明:结冰环境温度为-15~-10 ℃间出现了剪切强度极大值。来流速度越大,水滴的惯性力越大,水滴间的间隙越容易被填充,冰与固壁间剪切强度也就越大。平均水滴直径(MVD)为35 μm附近出现了剪切强度极小值,水滴直径越小,水滴之间的间隙就越小,结冰就越致密。水滴直径越大,水滴越容易发生平铺,水滴之间的间隙就容易被填充,反而造成剪切强度增大。