二维翼型大尺寸过冷水滴撞击特性及冰形分析

基于现有大尺寸过冷水滴(supercooled large droplet, SLD)动力学特性,分析水滴变形对阻力的影响.并根据几种典型的反弹/飞溅模型,分析了SLD的阻力变化、反弹、飞溅等对水滴撞击特性的影响,采用软件FENSAP-ICE的飞溅模型和LEWICE 2.0的反弹模型研究了反弹及飞溅现象对冰形的影响.计算结果显示:水滴撞击前的破碎现象对水滴尺寸分布有较大的影响,进行撞击特性以及冰形计算的时候需进行考虑;SLD破碎、飞溅、反弹降低了局部水收集系数、减小了水滴撞击范围;飞溅现象主要发生在机翼前缘附近区域,反弹主要在撞击边缘区域;水滴直径增加,飞溅现象逐渐减弱,但边缘位置的反弹现象一直很明显.SLD变形带来的阻力影响对冰形及结冰区域影响很小;与未考虑飞溅及反弹现象得到的冰形比较,考虑飞溅及反弹得到冰形前缘区域形状变化不大,但是整体结冰区域减小.      

变形破碎特性对SLD结冰过程影响

基于结冰过程数值模拟,针对过冷大水滴（SLD）条件下的变形破碎特性及其对成冰过程影响进行了研究。利用动态阻力计算模型,分析了水滴变形对水滴运动轨迹影响,采用基于泰勒类比理论的计算模型,研究了水滴破碎过程、子液滴粒径分布特性等。完成了NACA 0012翼型典型SLD结冰算例的数值分析,与参考文献计算结果和实验数据进行对比。结果表明:对于SLD结冰,变形和破碎主要改变了水滴运动轨迹和表面撞击水分布,使水滴撞击极限变小,在加入变形破碎特性的计算模型后,可以较好预测结冰极限位置和冰形轮廓,说明该方法对SLD变形破碎效应及其影响的模拟是可行和正确的。      

过冷大水滴变形与破碎的影响因素

过冷大水滴（SLD）结冰过程中,影响水滴变形与破碎及其导致的效应的因素众多,多种影响因素耦合作用于水滴变形与破碎的规律还未得到清晰认识。采用数值计算的手段,模拟了SLD运动过程中的变形与破碎现象,分别研究了水滴直径、流场速度及绕流物体尺寸变化对水滴破碎与变形及其效应的影响。通过研究,提出了包含SLD的典型流场区域划分方法,得到了不同因素对水滴破碎区域面积、物面的水滴收集特性、破碎子水滴分布特性的影响规律,发现水滴直径、流场速度、翼型前缘半径与变形和破碎影响撞击特性的程度呈正相关关系。在此基础上,提出了一个综合影响参数X,线性回归分析的结果初步说明,X是表征水滴变形与破碎及其效应的合适参数。相关工作可以为SLD结冰的计算、试验以及防/除冰设计提供参考。     

旋转部件复杂表面水滴撞击计算

为模拟旋转体结冰问题的过冷水滴运动及撞击过程,基于欧拉方法及单旋转坐标系模型,建立了三维旋转水滴运动模型,并提出了相应的数值求解方法。采用单旋转坐标系对空气及水滴两相流动过程进行处理,通过引入惯性力,将惯性系下的周期转动边界转换为定常流动边界;利用欧拉方法,使用单向耦合形式描述空气—水滴流场;在单旋转坐标系下,向控制方程内引入科里奥利加速度及牵连加速度,进行非惯性系下欧拉方程的修正,从而描述水滴运动过程;采用有限容积求解器对空气及水滴运动的控制方程组进行求解,通过引入源项定义单旋转坐标系下的惯性力,得到空气流场及水滴场的速度、体积分数分布,进而得到表面水滴撞击特性。采用上述方法对旋转帽罩与叶片模型进行算例分析,结果表明所建立的旋转水滴计算方法有效,对比静止状态表面的结果,旋转对帽罩的水滴撞击特性影响甚小,而对桨叶存在显著影响;由于帽罩具有中心对称的特性,因而旋转带来的切向速度变化对其水滴撞击特性影响不明显;桨叶表面水滴收集系数随旋转角速度增大而增大,同时收集系数在表面的分布会向迎风方向偏移,较大的角速度对应了更为显著的收集系数增幅与偏移现象。     

发动机整流支板大尺寸过冷水滴撞击特性

目前航空发动机防冰系统通用设计方案是以过冷小水滴的撞击特性为依据,但最新的研究结果显示大尺寸过冷水滴(SLD)环境下的结冰会对飞行器带来更加恶劣的影响.基于此,对大尺寸过冷水滴环境下发动机进口支板的水滴撞击特性进行了研究.首先对水滴运动控制方程进行了修正,以模拟水滴的变形及破碎的影响;然后对水滴在壁面处飞溅、反弹现象进...     

过冷大水滴条件下结冰相似准则

虑到中国现有结冰风洞在过冷大水滴(SLD)结冰条件模拟能力方面的局限,发展SLD结冰风洞及其相应的缩比方法迫在眉睫。研究了大水滴运动机理,并依据绕流流场相似、水滴运动和撞击特性相似、结冰过程热力学特性相似原理,提炼出SLD结冰条件下的6种相似缩比准则,利用结冰数值模拟软件FENSAP-ICE对所提相似缩比准则进行筛选与比较,从而基于结冰相似百分比量化判断准则找出给定条件下的最优准则法。结果表明:考虑SLD运动特性的撞击参数与ONERA法相耦合的方法以及撞击参数与Ruff法相耦合的方法的计算结果与实际成冰构型吻合良好；此外,撞击参数与ONERA法相耦合的方法在计算直径为171 μm的SLD时优于其他。      

水滴撞击特性的重力影响分析

以NACA0012翼型为研究对象,采用拉格朗日法分析了重力对水滴撞击特性的影响,获得了过冷大水滴（SLD）的撞击特性。计算结果表明：当水滴直径超过50μm时,重力对水滴轨迹以及极限轨迹的初始位置影响很大;SLD的水收集系数和收集区域都较直径为10～40μm的水滴（美国联邦航空条例FAR25部附录C规定的防冰系统设计用水滴直径范围）有较大幅度的提高,该结论可为飞机防冰系统的设计以及冰风洞中喷嘴安装位置的确定提供参考。     

水滴撞击特性的数值计算方法研究

 对水滴撞击特性的数值计算方法进行了研究。在对防冰部件流场进行计算的基础上,通过对水滴所在单元的判断和单元内流场速度的插值求解以及水滴运动轨迹计算边界的判断,采用差分法对水滴运动方程进行了数值求解,得到了水滴运动轨迹,从而确定了水滴的撞击极限、总收集系数和局部水收集系数等水滴撞击特性参数,为飞机防冰系统的设计奠定了基础;此外,以发动机进气道的水滴撞击特性的计算为算例,研究了飞行高度、飞行速度及水滴半径对水滴撞击特性的影响。     

水滴撞击特性的高效计算方法

针对拉格朗日方法计算水滴撞击特性效率低、通用性差等问题,发展了一种水滴撞击特性的高效计算方法。在求解绕流流场的基础上,结合逐级结构化管理的边界信息存储方式,采用目标扩散追踪方法对水滴所在网格单元进行快速计算,并插值得到该点处的流场信息,逐个求解水滴运动方程得到各水滴的运动轨迹,从而确定水滴撞击极限、收集系数等撞击特性参数。通过对NACA0012翼型、GA-W(1)两段翼型和某三段翼型的计算得到不同状态下的水滴撞击特性,计算结果表明,该方法与传统方法相比具有计算效率高、结果可靠、通用性好等优点。     

动力对全机水滴收集率的影响计算

针对飞机在飞行中遭遇过冷水滴撞击并结冰现象,建立了适合于发动机带动力情况下结冰过程水滴收集率计算的三维数值方法和计算程序.其基本思路为:采用多块技术与SIMPLE方法计算空气流场,以流场分布的计算结果为基础,求解水滴相的控制方程,进而获得物体表面的水滴收集率.空气相控制方程和水滴相控制方程均写成典型输运方程的形式,采用一致的有限体积法离散求解,方便了计算程序的编制.对某型运输机巡航构型有/无动力条件的水滴收集率进行了比较计算,获得了不同直径水滴在飞机表面的撞击特征以及水滴收集率在飞机机翼、平尾、垂尾和发动机进气道唇口上的分布规律.研究表明:(1)发动机是否带动力对机翼、平尾、垂尾的水滴收集率基本无影响;(2)飞机带动力主要影响发动机进气道唇口处的水滴收集率,带动力后唇口的收集率比无动力情况高,水滴撞击范围增大,在进行防除冰研究和设计时需引起重视.     

运动油滴/固体壁面斜碰撞的状态辨识及特征分析

 航空发动机轴承腔油气两相流动与换热研究是轴承腔润滑与换热设计的重要基础工作,而运动油滴与轴承腔壁面的碰撞状态辨识及特征分析是轴承腔油气两相流动与换热研究的组成部分。基于VOF方法建立了运动油滴与固体壁面斜碰撞数值仿真模型,实现了运动油滴与固体壁面斜碰撞状态的数值模拟,分析了碰撞油滴直径、入射角和碰撞速度对沉积油膜铺展长度、油膜堆积厚度和溅射油滴数目等碰撞状态特征参数的影响规律,提出了运动油滴与固体壁面斜碰撞条件下的状态判断准则,并得到了文献的物理试验工作的支持。研究结果表明:随着油滴入射角减小,沉积油膜铺展长度和油膜堆积厚度增大,溅射油滴数目减少;随着碰撞速度和油滴直径的增大,沉积油膜铺展长度和溅射油滴数目均增大,但油膜堆积厚度对应前者呈减小趋势、对应后者仍呈增大趋势。     

过冷大水滴动力学特性对结冰影响数值研究

对比常规水滴,过冷大水滴结冰对飞机气动性能和飞行安全具有更大的危害性。在常规水滴结冰过程数值模拟基础上,针对过冷大水滴条件下各种动力学特性及其对成冰过程的影响进行了数值研究。采用结构化网格和中心有限体积法求解 N-S 方程获得空气流场,用拉格朗日法求解水滴流场,基于改进的 Messinger 热力学模型完成结冰过程模拟;对于过冷大水滴,采用泰勒类比理论和水滴碰撞模型,对比分析了水滴动力学特性及其对结冰过程的影响;通过数值模拟获得了水滴 Langmuir D 分布的撞击特性及结冰计算结果,研究了水滴多尺度分布对成冰过程的影响规律。通过一些典型结冰算例的数值模拟,并与参考文献结果及实验数据的对比分析,说明了本文计算方法及计算模型模拟过冷大水滴结冰是可行和正确的。     

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

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

浅析过冷大水滴规章对防冰系统设计的影响

阐述了FAA颁布过冷大水滴规章的背景,浅析了即将新增的FAR25.1420过冷大水滴结冰条件条款,以及其对机翼防冰系统和结冰探测系统设计的影响;相比14CFR附录C结冰条件,过冷大水滴结冰条件的水滴撞击区较大,更容易在结冰防护区后部形成严重的后溜冰,应增加上表面弦向防护区的范围;基于已有的结冰探测系统,在机翼上表面容易形成后溜冰的部位,装备能够齐平安装的结冰探测器作为补充,是一种有效探测过冷大水滴结冰条件的方式;浅析了FAR25附录O过冷大水滴结冰条件的分类;给出了修订后的25.1093（b）的防冰译文。     

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

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

过冷水滴撞击三维机翼的数值模拟

基于欧拉两相流理论,提出一种数值模拟过冷水滴撞击三维机翼的方法.根据水滴拟流体模型建立三维水滴控制方程;提出一种可穿透型壁面边界模拟水滴对机翼表面的撞击;由水滴流场的求解结果得到机翼表面的水滴收集特性.对ONERA M6机翼表面的水滴收集特性进行了研究,计算结果表明所提出的数值模拟方法是有效、可行的,可以为三维积冰的数值模拟研究奠定良好的基础.