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

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

变形破碎特性对SLD结冰过程影响

基于结冰过程数值模拟,针对过冷大水滴（SLD）条件下的变形破碎特性及其对成冰过程影响进行了研究。利用动态阻力计算模型,分析了水滴变形对水滴运动轨迹影响,采用基于泰勒类比理论的计算模型,研究了水滴破碎过程、子液滴粒径分布特性等。完成了NACA 0012翼型典型SLD结冰算例的数值分析,与参考文献计算结果和实验数据进行对比。结果表明:对于SLD结冰,变形和破碎主要改变了水滴运动轨迹和表面撞击水分布,使水滴撞击极限变小,在加入变形破碎特性的计算模型后,可以较好预测结冰极限位置和冰形轮廓,说明该方法对SLD变形破碎效应及其影响的模拟是可行和正确的。      

二维翼型大尺寸过冷水滴撞击特性及冰形分析

基于现有大尺寸过冷水滴(supercooled large droplet, SLD)动力学特性,分析水滴变形对阻力的影响.并根据几种典型的反弹/飞溅模型,分析了SLD的阻力变化、反弹、飞溅等对水滴撞击特性的影响,采用软件FENSAP-ICE的飞溅模型和LEWICE 2.0的反弹模型研究了反弹及飞溅现象对冰形的影响.计算结果显示:水滴撞击前的破碎现象对水滴尺寸分布有较大的影响,进行撞击特性以及冰形计算的时候需进行考虑;SLD破碎、飞溅、反弹降低了局部水收集系数、减小了水滴撞击范围;飞溅现象主要发生在机翼前缘附近区域,反弹主要在撞击边缘区域;水滴直径增加,飞溅现象逐渐减弱,但边缘位置的反弹现象一直很明显.SLD变形带来的阻力影响对冰形及结冰区域影响很小;与未考虑飞溅及反弹现象得到的冰形比较,考虑飞溅及反弹得到冰形前缘区域形状变化不大,但是整体结冰区域减小.      

结冰风洞水滴直径标定方法研究

过冷水滴粒径大小是重要的结冰云雾参数,获知结冰风洞中的水滴直径,是得到定量结冰风洞实验结果的基础。对于结冰风洞内水滴直径单一或者分布比较集中的情况,提出了一种采用数值计算和结冰风洞实验相结合的手段标定水滴直径的方法。该方法首先采用拉格朗日法数值计算水滴运动轨迹,得到撞击极限随水滴直径变化的关系曲线,在此基础上,进行结冰风洞实验,测量实验得到的水滴撞击极限,通过在撞击极限与水滴直径关系曲线上进行插值,进而得到实验水滴直径大小。采用该方法对0.3m×0.2m结冰风洞内的水滴直径进行了标定,分别计算和测量了25m/s和35m/s两种速度条件下的水滴撞击极限,得到的水滴直径值相差不超过1μm,初步说明该方法的合理性。同时,对于结冰风洞内水滴粒径多尺寸分布的情况,还提出了相应的标定其容积平均直径MVD的方法,该方法在计算水滴收集率的基础上,通过测量驻点处的结冰厚度,实现对MVD的测量。采用本文提出的两种方法进行结冰风洞水滴粒径标定,只需要一般的长度测量工具即可进行,操作方便,成本低廉,克服了常规的水滴直径测量或标定需要专门设备的不足。     

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

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

MVD对翼型冰形及气动性能影响数值模拟

结冰问题严重威胁现代飞行器的飞行安全。当前,数值模拟结冰过程是一种高效的结冰问题研究方法。本文首先对OAXXX翼型生成六面体计算网格,然后对该翼型的流场、水滴撞击、结冰过程进行了数值模拟研究,重点研究了不同的MVD(平均水滴直径)对于水滴收集效率、冰形的影响,并将不同MVD情况下结冰后的翼型气动特性和结冰前进行了对比。     

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

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

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

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

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

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

三维翼面结冰过程及其影响数值研究

为得到三维机翼及翼身尾构型在结冰条件大气中的结冰过程和影响,采用结构化网格对M6机翼和翼身尾构结冰过程进行数值模拟,对飞机结冰过程的数值模拟采用欧拉两相流法求解水滴撞击特性,基于Messinger热力学模型模拟冰形增长过程,对其水滴撞击特性、结冰形状及结冰后对气动特性影响进行研究.分析结果可得:机翼展向,随着弦长变小,结冰量逐渐增多,尤其对大展弦比机翼更为明显.同时机身前端迎风面也是比较容易积冰的区域.     

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.     

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

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

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

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

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

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

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

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

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

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

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

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

三维复杂表面水滴撞击特性计算

为预测多段翼和发动机进气道等三维复杂形状表面在结冰气象条件下的水滴撞击特性,提出了一种基于欧拉两相流模型求解水滴收集系数的方法.空气相和水滴相认为是单相耦合的,空气流场由Euler或Navier-Stokes(N-S)方程独立求得.给出了一种水滴相撞击壁面边界条件的处理方法.为避免因局部水滴容积分数异常而导致计算发散,提出了一种自适应的数值扩散模型以增加高阶计算的稳定性.验证计算了圆柱、MS(1)-0317翼型、球面等典型二三维物面的收集系数,并计算了水滴直径呈某种分布时的情形.与公开发表的试验数据对比表明,计算方法准确有效,能很好的应用于三维复杂表面的水滴撞击特性预测.