基于Myers模型的三维结冰数值仿真

三维结冰表面上的水膜流动和结冰增长是结冰计算模型应考虑的核心内容,其中广泛应用的是Myers模型。Myers模型考虑了空气剪切力和空气压力对结冰表面水膜流动的影响,以及冰层、水膜和空气之间的导热与对流传热对结冰速率的影响。本文在使用Myers模型进行结冰预测时,发现Myers模型对霜冰转化为明冰的判断标准存在缺陷,会在结冰极限处产生不合理的冰角。因此对Myers模型的结冰类型判断标准进行了修改,对机翼表面的结冰过程进行了更加准确的模拟,并应用了有效的离散算法计算水膜流动和结冰增长过程。对比了二维NACA0012翼型的单步法、多步法计算结果和实验结果。明冰结冰温度较低时,本文计算结果与实验结果吻合很好,明冰结冰温度较高时,本文对上冰角的计算与实验结果有一定差距。本文提供了三维GLC-305后掠翼的结冰计算结果和实验结果的对比,冰角厚度的计算结果略小于实验结果,但整体趋势一致。     

机翼三维结冰数值模拟

采用欧拉-拉格朗日法计算空气、水滴两相流动,根据质量守恒和能量守恒计算结冰量和溢流水量,溢流水的流动方向及不同方向的流量分配由当地空气的速度矢量、网格边界的法向量计算得到。采用非稳态方法、非结构网格计算了NACA0012平直翼、截面为GLC-305的后掠翼的结冰,并与文献中的实验结果进行了对比,结果显示最大结冰厚度与实验吻合较好,预测的冰形与实验冰形发展趋势一致。计算结果表明,提出的三维结冰计算方法是有效、可行的,为三维结冰计算提供了一条可选的途径。     

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

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

基于准三维冰形算法的结冰预测技术及应用研究

机翼结冰是威胁民用飞机飞行安全的主要因素之一,准确预测机翼的结冰过程和特点,可为防除冰设计提供依据和支撑。提出了利用二维模拟三维的准三维冰形算法,该算法在保证结冰预测精度的同时,极大地降低了型号结冰计算成本;给出了顺气流和垂直于前缘两种典型民用飞机后掠翼机翼的二维剖面截取方法,并将采用该方法的计算结果与三维后掠机翼的结冰预测结果进行了对比分析。结果表明,相对于顺气流方式,垂直于前缘方式的准三维冰形算法与三维后掠机翼的结果更为接近。     

翼身融合布局与传统布局飞机结冰特性对比分析

翼身融合（BWB）布局是提升未来民机综合性能的重要布局方式，其结冰飞行安全问题逐渐受到人们重视。针对翼身融合布局飞机与常规布局飞机的结冰问题，本文通过数值模拟的方法来开展结冰特性研究。本文提出基于Navier-Stokes方程对空气流场进行求解，并应用欧拉法计算水滴撞击特性，之后采用Shallow-Water结冰热力学模型的结冰计算方法。首先，通过将翼身融合布局飞机和传统布局飞机的计算结果与风洞试验数据进行对比，验证空气流场计算的正确性，并将两者进行对比分析；其次，数值预测飞机表面冰形的特征，将两种布局飞机结冰特性进行对比，结果表明，结冰对两种布局飞机气动外形的破坏程度从后掠翼翼根至翼尖逐渐变大，但传统布局飞机结冰只发生在机翼前缘和机头处，而翼身融合布局飞机前部几乎都发生了结冰，可为相关的结冰特性研究及防除冰设计提供技术参考。     

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

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

机翼积冰数值模拟

采用计算流体力学（CFD）方法对三段翼和MS-317后掠翼进行了积冰数值模拟研究,基于欧拉两相流理论,应用分区算法,对空气流场和水滴流场进行了数值求解,得到水滴收集率.求解3-D积冰模型得到积冰量.将MS-317的冰形计算结果与实验数据以及软件LEWICE3D的计算结果进行了对比.结果显示:三维积冰模型具有较好的精度,虽然预测的冰形与实验数据有一些差异,最大厚度方面,误差在13.1%以内,但是积冰的类型与体积和实验数据基本一致.      

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

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

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

临界冰形是指在适航规章结冰包线内,每个可用飞行构型下,对飞机操纵性和稳定性影响最严重的冰形,临界冰形分析是飞机适航取证中的重要工作。对临界冰形确定方法,及临界冰形对气动特性的影响进行了研究。发展了基于 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湍流模型,分别针对流线冰型和角状冰型,量化分析了不同微观特征参数的冰形对翼型失速特性的影响规律,为研究飞机结冰问题提供了思路和基础。结果表明,对于流线结冰翼型,冰形表面粗糙度会明显影响翼型的失速性能;而对于角状结冰翼型,两个凸起的羊角状外形才是导致气动力损失的主要因素,粗糙度的影响往往可不考虑在内。     

飞机机翼表面霜冰的三维数值模拟

基于欧拉两相流理论对三维情况下飞机机翼表面的霜冰进行了数值模拟.根据水滴拟流体模型建立三维水滴控制方程;提出一套水滴控制方程的数值求解方法;由三维水滴流场的求解结果计算机翼表面的水滴收集特性,提出一种三维积冰外形的生成方法,完成了对飞机机翼表面霜冰的三维数值模拟.对ONERA M6机翼在不同迎角下霜冰的积冰情况进行了数值预测,并分析了结冰条件对积冰的影响.      

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

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

模拟飞机迎风面三维积冰的数学模型

针对飞机迎风面的积冰建立了三维积冰模型,给出了模型的求解方法和求解步骤.该模型不但可以模拟明冰积冰、霜冰积冰和不结冰三种情况下的冰层生长,还可以模拟未凝结水膜在冰层表面的流动,而且在求解过程中可以自动判断过冷水滴撞击到迎风面以后的积冰形态.利用该模型数值模拟了翼型与平板相正交而形成的简单三维结构上的积冰过程,并将所得结果中翼型上三维冰形的二维截面与美国航空航天局(NASA)的冰形进行了对比,获得了较好的一致性,证实了所建立的积冰模型是合理的.      

发动机进气道短舱前缘结冰三维模拟研究

为计算发动机进气道短舱前缘的结冰冰形,在Messinger结冰热力学模型的基础上发展了一套纯三维的表面溢流水流动结冰模型。考虑空气摩擦力为表面溢流水流动的主要驱动力,用空气对进气道表面的剪切力来确定溢流水的流动方向及流量分配。为求解溢流水结冰模型,发展了一套重复查找表面控制体状态的方法,能够快速完成整个三维表面的计算。用该方法对某三维发动机进气道进行计算,得到了三维结冰冰形,并将计算结果与FENSAP-ICE计算结果进行了对比,结果显示两者的冰形轮廓基本一致,仅在冰角处存在差异,表明本文三维发动机结冰计算模型与计算方法是有效的,其计算精度与FENSAP-ICE结果相当。     

A spongy icing model for aircraft icing

Researches have indicated that impinging droplets can be entrapped as liquid in the ice matrix and the temperature of accreting ice surface is below the freezing point. When liquid entrapment by ice matrix happens, this kind of ice is called spongy ice. A new spongy icing model for the ice accretion problem on airfoil or aircraft has been developed to account for entrapped liquid within accreted ice and to improve the determination of the surface temperature when enter- ing clouds with supercooled droplets. Different with conventional icing model, this model identifies icing conditions in four regimes： rime, spongy without water film, spongy with water film and glaze. By using the Eulerian method based on two-phase flow theory, the impinging droplet flow was investigated numerically. The accuracy of the Eulerian method for computing the water collection efficiency was assessed, and icing shapes and surface temperature distributions predicted with this spongy icing model agree with experimental results well.