翼型积冰的数值模拟

本文采用计算流体力学的方法对NACA0012翼型前沿的Rime ice积冰进行了预测.在流场计算中,用有限体积(FV)法对二维定常不可压粘流的时均N-S方程进行离散求解;采用四阶龙戈-库塔(RK)法求解水滴运动方程;假定水滴在与翼型相碰撞的点处完全凝结,并且冰沿着与当地翼型表面法向一致的方向增长,以此来预测积冰的形状.文中同时对迎角为4°时结冰翼型的气动特性进行了分析.     

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

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

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

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

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

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

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

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

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

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

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

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

计入离心力影响的直升机旋翼翼型结冰数值模拟

建立了一套计入离心力影响的直升机旋翼翼型结冰的数值模拟方法.首先生成围绕翼型的贴体正交网格,然后用Navier-Stokes（N-S）方程求解黏性绕流流场.在此基础上,利用拉格朗日法建立水滴运动方程.其中,为提高计算效率,提出了结合位移矢量的水滴所处单元寻觅方法.最后,结合桨叶工作特点,发展了一种计入离心力影响的三维结冰模型.通过与桨叶结冰实验的对比,验证了本文结冰预测方法的可靠性.对比常规结冰模型,桨叶结冰量减少22.3%;若考虑桨叶的挥舞运动影响,桨叶结冰量进一步减少,表明了离心力及桨叶运动在结冰数值模拟中的重要性.通过不同剖面间的结冰量和冰形对比,分析并获得了桨叶结冰特征.结果表明离心力的影响程度随径向位移的增加而增加,下翼面结冰量随挥舞角的增加而减少.      

低速两段翼型的数值模拟

在多段翼型风洞实验中,很难观察翼型不同迎角下主翼、襟翼上的绕流及缝道流动。使用了O-H型混合的结构化网格,采用k-ε二方程湍流模型求解可压的N—S方程来模拟二元风洞中两段翼型的流场。计算结果与实验进行了比较,结果表明：计算与实验结果吻合良好,说明本方法可以较好地模拟两段翼型的绕流,并能很好地显示流场的变化情况。另外,从结果中可以看出：随着迎角的变化,主翼尾流和缝道间的流动对襟翼附面层分离有很大的影响作用。     

粘性流绕俯仰振动翼型流动的数值模拟

本文直接从N-S方程出发,利用LU-ADI格式和Baldwin-Lomax代数湍流模型研究了粘性流绕俯仰振动翼型的流动。数值实验表明,本文的数值结果同实验吻合较好。通过数值模拟手段研究了振动对流场中激波和分离这两大主要特征的影响,结果发现;(1)激波滞后于攻角的变化,如当NACA0012翼型在负攻角状态下上翼面存在激波而下翼面无激波。(2)由于翼型的振动,分离被减轻升力相对提高。(3)随翼型攻角和振幅的变化,翼型振动的升力回线走向可不同。     

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

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

大粒径过冷水滴结冰模拟破碎模型研究

目前的民机结冰适航规范中,大粒径过冷水滴(supercooled large droplet, SLD)环境尚未得到规定,是一个重大的飞行安全隐患。由于SLD的粒径较大(大于50μm),在研究其结冰时需要考虑水滴动力学过程。针对SLD成冰模拟中破碎模型的不足,根据对破碎机理的认识提出了多因素破碎(multifactor effected breakup, MEB)判定模型,并结合积冰算法对SLD环境中的结冰进行了模拟,验证了模型的有效性,证明其能提高对SLD结冰的模拟能力。从翼面收集率及冰型生长两方面对SLD的结冰过程和水滴动力学过程产生的效应进行了分析。     

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

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

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

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

机翼积冰数值模拟

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

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

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

简化的积冰热力学模型及其应用

依据飞行器表面液态水凝结过程的热力学理论分析,提出了一种适用于不同环境温度的翼型结冰数值模拟方法.采用中心型有限体积方法求解N-S方程来获得翼型的绕流流场,结合给出的结冰数值模拟方法对两种典型温度下的翼型积冰情况进行了模拟.模拟所得的积冰冰形与文献中的实验数据吻合良好,表明该方法的可行性,并进一步给出了翼型表面的冻结系数分布.