迎风面三维积冰过程中水膜流动的计算方法

对影响迎风面三维积冰过程的冰层表面水膜流动进行分析,建立描述其流动的数学模型,引入在流动前锋下游预设薄水膜的方法处理其在表面干区的流动推进过程,提出超过临界厚度的水膜被吹除的方法处理角状冰下游局部水膜被吹离表面的现象.通过对翼型-平板结构中翼型上的冰形对比和水膜流动规律分析,验证了该水膜流动数学模型的合理性和计算方法的可行性.研究表明:与Messinger模型相比,积冰模型中引入水膜流动后,可以较好地模拟三维明冰积冰,同时明冰极限的模拟结果也更接近试验值;冰层表面的水膜厚度约为10-5m,流动速度约为10-2m/s.      

考虑水膜蒸发的三维明冰积冰数值研究

明冰积冰过程中冰层表面会覆盖一层薄水膜,水膜在冰层表面流动和传热,水膜蒸发会引起质量和能量传递,对积冰相变产生影响。为了研究水膜蒸发对积冰的影响,在已有三维积冰模型的基础上,考虑了明冰积冰过程中水膜的蒸发,建立了考虑水膜蒸发的三维积冰数学模型,并开发了相应的积冰计算程序。选取NASA提供的典型算例,模拟了明冰积冰条件下NACA0012翼型的积冰过程,并将计算结果与NASA试验和数值结果进行对比。结果表明水膜的平均蒸发量为1.5g/(m~2·s),蒸发吸收的热量占总传热的35%。与不考虑水膜蒸发的计算结果相比,本文得到的积冰量和冰形都与试验结果具有更好的一致性,说明三维明冰计算中考虑水膜蒸发是有必要和合理的。     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法.外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格-库塔方法.采用四阶龙格-库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率.使用非结构动态网格技术对积冰后的翼型网格进行快速调整.在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性.     

非结构网格在积冰数值模拟中的应用

基于非结构网格对翼型表面的积冰进行了数值模拟研究.通过求解气体欧拉方程得到翼型的绕流流场,使用穿透型壁面边界模拟水滴对翼型表面的撞击,采用迭代法对水滴控制方程进行数值求解,得到翼型表面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对NA...     

直升机旋翼桨叶翼型积冰的数值模拟

本文对直升机旋翼桨叶翼型积冰进行了数值模拟.建立了空气、水滴两相流控制方程,通过求解气体欧拉方程得到空气流场,求解水滴控制方程得到翼型表面的水滴收集率.根据结冰表面的对流换热特性,使用积冰数值模拟的热力学模型确定翼型表面的积冰量,利用冰层时间推进法模拟积冰过程,采用积冰法向生长模型预测积冰形状;对NACA0012翼型在不同环境温度下的积冰进行了预测,最后对预测结果进行了分析.     

基于非结构动态网格的翼型积冰过程数值模拟

讨论了积冰气象条件下翼型积冰过程的非结构动态网格数值模拟方法。外部无黏绕流流场采用基于非结构网格的分布式并行计算技术,空间离散采用VanLeer迎风有限体积格式,时间推进为5步显式龙格 库塔方法。采用四阶龙格 库塔方法求解过冷水滴运动轨迹方程,使用了一种合理、简单的物面水积冰量计算方法,提高了计算效率。使用非结构动态网格技术对积冰后的翼型网格进行快速调整。在对NACA0012翼型积冰形状进行校核验证后,对某型飞机在积冰气象条件下飞行时翼型积冰过程进行了数值模拟,得到了与参考文献结果相符的冰形特性。     

基于神经网络的翼型积冰预测

采用神经网络方法对翼型表面的积冰进行了预测研究.建立了一种型面曲线坐标系,通过坐标变换将积冰外形转化为几何曲线.利用傅立叶级数对此曲线进行拟合,实现积冰外形的参数化.以结冰条件作为输人,冰形参数作为输出构建神经网络,对NACA 0012翼型表面的积冰进行预测.通过将预测结果与实验结果以及数值模拟结果进行对比,说明所提出...     

考虑传质传热效应的翼型积冰计算

本文采用计算流体力学的方法对NACA0012翼型前缘涉及传质传热效应的glaze ice积冰进行了预测.在流场计算中,用有限体积法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;通过求解在翼型表面控制体积内积冰过程中所遵循的质量和能量守恒方程,并且假定冰沿着与当地翼型表面法向一致的方向增长,来预测积冰的形状.文中同时对结有glaze ice冰翼型的气动特性进行了计算与分析.     

翼型电热防/除冰系统的数值模拟

飞机积冰是飞机在积冰气象条件下飞行时,大气中的液态水在部件表面冻结并积聚成冰的一种物理过程.开展飞机积冰研究,可以使人们更深入地认识积冰产生的机理,为飞机防/除冰装置的设计提供理论基础.采用数值模拟的手段,对二维翼型的结冰过程和防/除冰过程进行了系统的研究.根据经典Messinger模型对机翼的积冰过程进行了模拟,并在此模型的基础上计算了所需要的防冰热功率,同时采用焓方法对飞机机翼的融冰过程进行了模拟.     

机翼积冰数值模拟

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

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

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

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

描述了某民用飞机在意大利IWT（icing wind tunnel）冰风洞的机翼结冰试验结果.根据冰风洞试验状态点,开展了基于拉格朗日方法和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.     

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

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

民用飞机临界冰型确定策略浅析

分析了适航规章第25部121修正案中关于民机临界冰型的修订内容,提出大型客机结冰适航验证所用临界冰型的确定策略。针对飞机不同飞行阶段分析飞机结冰特征,对飞机不同部件开展结冰参数临界性分析,确定临界结冰条件,结合结冰数值预测和冰风洞试验手段提出相应的临界冰型确定方法;同时提出适航验证所需失效冰型、延迟冰型等特殊冰型的确定方法。根据性能、操稳和动力专业对于结冰临界性的不同要求,对比评估各种冰型对飞机影响的严重程度,最终确定飞机用于适航验证的临界冰型。本文可为大型客机的结冰适航验证工作提供参考。     

三维机翼结冰模拟

针对机翼结冰问题,利用欧拉法求解空气水滴两相流流场,通过Fluent软件用户自定义函数功能编程求解三维机翼表面水滴局部收集系数以及结冰情况.在计算中.将水溢流方向沿机翼弦向和展向分解,以此对Messinger结冰热力学模型进行改进,使之更符合三维情况,并在此基础上进行三维传质传热分析;修正传统的二维表面对流换热系数求解...     

飞机三维结冰模型及其数值求解方法

 对飞机结冰外形进行纯三维数值模拟,是目前飞机结冰预测中的一大难点。为了建立三维结冰数值模拟方法,基于Messinger的二维结冰模型,提出了一种考虑液态水溢流效应的三维结冰计算模型,并针对该模型建立了表面单元内溢流水流动的分配方案、发展了相应的迭代求解方法。采用本文方法对MS-317后掠翼结冰进行了计算,并与实验和Lewice3D的计算结果进行了对比。研究结果显示：本文方法具有较好的收敛性,计算的霜冰和明冰外形均与Lewice3D计算的冰形一致;对于霜冰和结冰时间较短的明冰,本文计算的冰形与实验吻合较好;对于结冰时间较长的明冰,本文计算的冰形与实验对比还有一定差异,但冰生长的总体趋势和大致体积与实验一致。同时,对比了MS-317无后掠翼与有后掠翼结明冰的外形,发现机翼后掠导致的三维溢流效应对结冰外形有明显影响,因此,对于三维结冰分析,如果用二维截面的结果来代替三维结果,其合理性还需验证。     

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

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

结冰条件下的飞机飞行动力学仿真

基于飞机结冰参量模型对飞机在结冰条件下的飞行动力学特性进行了仿真研究.提出一种飞机结冰参量模型的参数计算方法,并针对一定结冰气象条件下的不同结冰时间对飞机的结冰严重程度进行了评估,基于小扰动理论对飞机结冰前和不同结冰严重程度下的纵向、横侧动稳定性以及操纵响应进行了仿真计算,根据仿真计算结果分析了不同结冰严重程度对飞机造...