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

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

三维笛卡尔网格的结冰数值模拟

<正>飞机结冰现象会对飞机安全飞行带来严重威胁,因此对飞机结冰现象研究的开展具有必要性。目前研究结冰的方法主要有两种:试验法和数值模拟法。数值模拟法,则是利用计算机建立机翼几何模型、划分计算网格,计算水滴撞击特性,进行机翼表面的传热、传质计算得到冰形。数值模拟不需要在现实中去模拟真实的飞机飞行情况,仅仅需要在数学上构建相关模型,这在成本上相对低廉,同时也无需承担试验风险,目前在结冰研究中成为了不可或缺的重要方法。国外三维结冰数值模拟研究开     

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

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

旋转机翼飞机旋翼/机身干扰流场数值计算分析

基于动量源方法,建立了旋转机翼飞机的旋翼/机身组合流场的数值模拟方法.首先通过Caradonna-Tung旋翼算例计算,验证了本方法的有效性;然后采用该方法对某旋转机翼飞机的悬停及低速前飞状态下的旋翼/机身干扰流场特性进行了详细的计算与分析.结果表明,该方法能够用于计算旋翼与机身的干扰流场,为旋转机翼飞机的设计提供了参考.     

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

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

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

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

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

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

机翼结冰分析与防除冰系统设计验证

飞机在结冰条件下飞行时可能发生结冰,飞机一旦结冰,会对安全飞行带来较大的隐患,如何降低飞机结冰带来的危害已成为飞机设计研究的重点内容.通过FENSAP-ICE对机翼进行数值模拟,并通过改进Messinger结冰热力学模型模拟更加真实的飞行情况;分析不同飞行环境下,飞机结冰前后机翼气动特性的变化,同时针对机翼设计一套防除冰系统并验证其可行性.结果表明:飞行速度越大,机翼表面的局部水收集系数越大;环境温度会影响机翼结冰的类型和结冰厚度,机翼发生结冰时,其升力系数减小、阻力系数增大,机翼的气动特性受到严重的影响;设计的电热防冰系统可以有效地预防机翼表面结冰,也可以进行周期性除冰.     

机翼结冰过程的数值模拟

为预测飞机机翼结冰的形状,建立了一套采用时间步进法进行结冰过程数值模拟的方法.开发了一个软件框架用于使结冰模拟和网格更新全自动化.基于Messinger的质量和能量守恒方程建立了结冰热力学模型.采用了基于Spalart-Allmaras湍流模型在粗糙壁面的扩展模型来计算表面传热系数.提出了一种新的结冰生长方式.能预测霜...     

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

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

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

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

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

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

基于水膜流动与耦合传热的翼型防/除冰计算

针对飞行器防/除冰过程中翼面上空气-水膜-冰层-机翼之间的耦合传质传热现象,建立了一种基于水膜流动与耦合传热模型的翼型防/除冰数值模拟方法。基于Myers水膜流动模型建立了防/除冰热载荷作用下翼面溢流水流动、积冰及内部温度分布的数值计算理论。对于翼型及冰层内的传热现象,利用焓理论及有限体积法建立了复杂多层结构传热的数值模拟方法,对于冰层相变过程,提出了一种基于焓理论的相变修正方法以考虑相变潜热对温度变化的影响。最终实现了翼型防/除冰过程的耦合计算,结果表明:通过结合不同界面处的传热边界条件和考虑了相变潜热效应的焓理论对水膜流动与翼型/冰层传热模型进行耦合求解,能够对翼型/冰层内温度分布进行准确计算,可实现对翼型防/除冰过程中溢流水流动及积冰特性的有效预测与分析。      

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

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

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

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

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

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

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

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

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

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

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

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