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

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

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

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

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

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

热气防冰参数对机翼表面溢流结冰的影响研究

机翼热气防冰数值模拟中,根据N-S方程求解流场,用Euler法获得水滴撞击特性,通过蒙皮导热将求解得到的内外流场进行耦合传热并达到稳定后,开始模拟结冰,来进行机翼热气防冰及形成溢流结冰的数值计算。计算结果表明,热空气防冰数值模拟是可行和合理的。采用数值模拟方法对机翼热空气防冰过程进行了模拟,得到了由于引气温度不足或机翼热空气保护面积不同而导致的不同溢冰高度和位置。分析了供气温度、防冰区域对翼面溢流冰形成的影响。结果表明:供气温度直接正向影响热交换后蒙皮表面温度。供气温度越低,溢流冰形高度越高,对气动特性影响也越大;热防护区域范围对溢流结冰结果也会产生影响,热防护区域越大,冻结位置距离驻点越远,而且冰形高度越低,对气动特性影响也越弱。     

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

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

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

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

不同来流条件对机翼结冰及气动特性影响研究

针对三维机翼结冰问题,基于Navier-Stoke控制方程,采用多时间步长的计算方法,通过求解气液两相流的Eluer法得到空气中水滴的轨迹,再根据Messinger模型计算结冰冰形来实现对机翼结冰的数值模拟。首先以三维NACA0012机翼为研究对象,进行结冰数值模拟方法验证,随后计算不同迎角和不同来流速度对三维ONERA M6机翼结冰冰形的影响,最后分析了机翼结冰对气动特性的影响规律。结果表明,结冰后机翼的阻力系数增大,升阻比减小;来流速度对机翼前缘的结冰类型影响较大,当机翼前缘的结冰类型由霜冰变为混合冰时,将使机翼的气动特性恶化。     

非结冰气象条件下机翼热气防冰系统数值模拟

为了在自然结冰飞行前预估飞机防冰系统的性能,减少飞行风险,开展了非结冰气象条件下飞机机翼热气防冰系统的数值模拟研究。数值模拟采用了计算流体力学方法;外部对流传热系数的计算采用附面层积分方法。机翼热气防冰系统表面温度的计算是将蒙皮外部散热热流、内部热气加热热流以及蒙皮导热三者进行热耦合,并进行了改进,改进的方法不要求防冰腔内、外两套网格在重合面网格处一致,是通过双向线性插值将一方面网格信息插值传递到另一方表面网格上,提高了计算效率。     

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

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

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

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

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

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

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

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

冰形表面粗糙度对翼型的失速特性影响分析

结冰翼型气动影响分析是飞机结冰问题重要的研究内容之一,可帮助进行飞机结冰安全分析和容冰设计。为填补现阶段关于冰形表面粗糙度对翼型气动影响的研究空缺,以冰形的粗糙度参数(冰形表面粗糙高度和冰形表面粗糙分散度)为依据,提出了一种冰形粗糙度叠加方法;并通过有限体积法求解空间离散的RANS方程,结合S-A湍流模型,分别针对流线冰型和角状冰型,量化分析了不同微观特征参数的冰形对翼型失速特性的影响规律,为研究飞机结冰问题提供了思路和基础。结果表明,对于流线结冰翼型,冰形表面粗糙度会明显影响翼型的失速性能;而对于角状结冰翼型,两个凸起的羊角状外形才是导致气动力损失的主要因素,粗糙度的影响往往可不考虑在内。     

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

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

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.     

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

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