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

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

改进的翼型积冰数值模拟方法

基于微团概念,以雨滴运动的基本假设出发,对翼型积冰的拉格朗日方法进行了改进.新方法克服了原方法计算效率低、冰形同雨滴分布密度及网格密度相关的缺点,继承了原方法的鲁棒性,不存在欧拉方法的收敛性问题.为验证方法的正确性,计算了NACA0012翼型在0°和4°攻角下的霜冰积冰情况、表面水收集量和局部撞击系数.结果表明新方法减少了计算量,预测的冰形和实验结果一致,并分析了积冰对气动性能的影响.      

离心力对于旋转叶片表面积冰影响的数值模拟研究

为了加深对旋转叶片表面积冰现象的认识,发展了一种适用于旋转表面积冰预测的数学模型,并采用该模型研究了离心力对于旋转叶片表面积冰的影响。通过对积冰控制体内的由于撞击水滴形成的薄水膜流动的质量、动量和能量守恒进行分析,发展了旋转贴体非正交曲线坐标系下的积冰模型并给出了相应的计算方法。所发展模型的动量方程中考虑了离心力对水膜流动的影响,并且在能量方程中考虑了净流出水带走的能量。采用文献中旋翼表面的积冰试验结果对所发展的模型和计算方法进行了验证,所计算得到的冰厚与试验结果吻合较好,计算得到的冰厚精度相比于LWICE软件有所提高。计算结果还表明,离心力增大会导致水膜速度在叶片展向上的分量变大,致使流出控制体的水质量增加,最终导致冰厚在驻点附近略微变小。     

飞机机翼积冰研究方法和进展分析

通过对目前国内外研究飞机机翼积冰三种方法的分析,指出了它们各自优缺点,重点介绍了现行国内外研究的进展和新趋势,提出了研究方法的改进方案和进一步研究的建议。     

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

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

应用混合网格分析复杂积冰翼型气动性能

通过发展一种点-点对接的结构/非结构混合网格生成方法,避免了复杂积冰外形难以划分结构化网格的问题,并对复杂积冰翼型的气动性能进行了分析计算.在靠近积冰边界的内层采用了非结构网格以拟合复杂的边界,在非结构网格外采用结构化网格以节省存储空间和计算时间.为了检验网格对选用湍流模型的影响,整个流动区域分别采用SST和SA湍流模型,求解了雷诺平均N-S方程.数值计算结果表明SST模型更适于模拟复杂分离流动,积冰对翼型的气动性能造成了严重的影响.     

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

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

复杂积冰翼型气动性能分析

为了克服复杂外形难以划分结构化网格的问题, 采用结构/非结构混合多块网格对复杂积冰翼型的气动性能-进行了分析计算.在靠近积冰边界的内层采用了非结构网格, 在非结构网格外采用结构化网格以节省存储空间和计算时间.不同类型网格之间采用了点点对接的方式以简化通量重构的处理方式.整个流动区域采用SST湍流模型, 求解了雷诺平均N-S方程.数值计算结果表明结构/非结构混合多块网格在复杂积冰翼型气动性能计算上的可行性, 积冰对翼型的气动性能造成了严重的影响.      

锥角对旋转整流罩积冰影响的模拟实验

采用缩比模型对4种不同锥角的旋转整流罩进行了冰风洞积冰模拟实验.推导出旋转表面积冰相似参数,并根据重要参数的匹配确定缩比整流罩模型的积冰实验参数.在实验中采用高速摄像系统记录冰生长过程及最终冰形.结果表明:4种锥角的整流罩表面积冰均由初期生成、分布连续的明冰和后期增长迅速的白色霜冰构成;锥角小于等于74°的整流罩表面霜冰为针状或粒状,积冰厚度较小;锥角大于80°的整流罩表面霜冰为羽毛状,积冰厚度较大,并伴随冰脱落现象,锥角较大的整流罩冰脱落位置向下游移动.      

翼型结冰过程数值模拟验算与分析

应用FENSAP-ICE结冰计算软件,对NACA0012翼型进行了流动特性、水滴撞击特性以及冰型生成过程的计算;同时,对结冰前后的翼型进行气动力特性计算对比分析,其中包括升力特性对比、阻力特性对比、流场细节分析以及压力系数分布对比。计算结果表明：翼型前缘结冰后,导致翼型前缘气流提前分离,最大升力系数、失速攻角大幅减小,...     

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

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

三维机翼结冰模拟

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

飞机机翼结冰过程的数值模拟

提出了一种适用于各种结冰气象条件下飞机机翼结冰过程的数值模拟方法。结合欧拉坐标系下空气-过冷水滴两相流动控制方程的计算,对不同气象条件下飞机机翼的结冰过程进行了数值模拟,得到了4°攻角NACA-0012翼型表面霜状冰、混合冰及瘤状冰3种冰形。所得冰形与文献中的试验数据和计算结果吻合良好,表明本文方法有效。另外还计算了不同的冰形对机翼升、阻力的影响,结果表明瘤状冰对机翼气动特性的影响最大。      

瘤状冰结冰过程的数值模拟

提出了一种可考虑粗糙度影响的瘤状冰结冰过程数值模拟方法。结合欧拉坐标系下空气-过冷水滴两相流动控制方程的计算,对飞机机翼在不同温度下的瘤状冰结冰过程进行了数值模拟,得到了NACA0012翼型4°攻角下结冰表面的瘤状冰冰形。通过与实验数据对比,表明该的方法是有效的。      

尾部可调的发动机进口支板结冰数值模拟研究

基于模块化思想,在FLUENT软件基础之上开发了发动机进口支板的结冰计算程序,针对尾部可调支板的结冰问题进行了计算研究。通过分析发动机进口支板结冰的物理过程,将支板结冰计算分成四个模块,分别为两相流流场计算、水滴收集计算、热力学分析和表面结冰量计算、动网格计算。首先将各模块单独进行研究,通过FLUENT提供的UDF函数及宏函数逐一实现各模块功能,最后将各模块整合,实现了发动机进口支板的结冰计算。通过计算获得了调整角度为30°时,支板的液态水含量分布、不同时间下支板结冰的区域以及不同温度下支板结冰的特点,给出了可调支板应重点防冰区域。数值模拟方法为工程应用提供了一条可选的途径。     

基于CFD的机翼结冰过程分析

采用一种适用于分析机翼结冰过程的数值模拟方法,该方法在飞机结冰气象条件下采用计算流体力学（Computational Fluid Dynamics,CFD）对机翼流场进行计算,获得流场中二维机翼翼型的水滴撞击特性,并遵循质量和能量守恒方程,计算出机翼表面的结冰量,预测出机翼的积冰形状。所得冰形与文献中的试验数据吻合良好,表明数值模拟方法有效。另外还分析了结冰对机翼升力、阻力和压力系数的影响,其结果可以指导飞机机翼防／除冰系统的设计研究     

飞机机翼表面霜状冰结冰过程的数值模拟

提出了一种处理结冰后结冰表面固壁区域移动的移动边界技术,结合欧拉坐标系下空气-过冷水滴两相流动控制方程的计算,对霜状冰的结冰过程进行了数值模拟,得到了NACA0012机翼在0°和4°攻角下机翼表面结冰后形成的冰形,与文献中的实验数据对比,表明本文的方法是可行和有效的。