结冰风洞中过冷大水滴云雾演化特性数值研究

为明晰结冰风洞中过冷大水滴（SLD）云雾演化特性,发展了基于欧拉法的SLD液滴运动、传热和传质耦合计算方法,针对3 m×2 m结冰风洞主试验段水平收缩构型,分析了SLD云雾沉降收缩特性、动量平衡特性和热平衡特性,探索了液滴变形破碎的影响,评估了构型出口处SLD液滴动量平衡和热平衡状态。研究结果表明：直径超过250 μm的SLD液滴在构型内会发生显著形变,液滴尺寸越大则变形程度越强,尤其在160 m/s工况下,当液滴直径超过750 μm后,SLD液滴甚至会发生破碎;液滴变形破碎效应会增大液滴加速度和液滴温度下降率,促使SLD液滴趋近动量平衡和热平衡状态;SLD云雾（最大液滴直径小于1 000 μm）在构型出口处会出现显著的粒径浓度分层、动量分层和热分层现象,其中直径小于100 μm的小尺寸液滴速度快、温度低且不断凝结,趋于平衡态,但直径超过500 μm的大尺寸SLD液滴速度慢、温度高且不断蒸发,则显著偏离平衡态;增大试验段气流速度尽管会减弱SLD云雾粒径浓度分层程度,但会增强动量分层和热分层程度,尤其在160 m/s工况下,SLD云雾会均匀分布在构型出口中心区域内（-0.75 m < Y < 0.75 m且-0.5 m < Z < 0.5 m）,与其平衡态间的最大速度差和温度差将分别超过18 m/s和20℃。     

飞机结冰中水滴撞击特性欧拉法的网格影响分析（英文）

过冷水滴撞击特性计算研究是飞机结冰预测与防除冰系统性能分析的基础,欧拉法是过冷水滴运动及撞击特性计算的常用方法之一,分析欧拉法在计算复杂表面水滴运动撞击时的准确性具有重要意义。以NACA0012翼型、冰风洞风道、S形进气道及三段翼为研究对象,利用欧拉法采用不同网格计算获得了水滴运动及局部水收集系数,结果表明：当水滴运动没有受到上游部件的导流或者遮挡时,欧拉法计算结果具有网格无关性;当水滴运动受到上游部件影响发生轨迹偏转时,欧拉法计算结果具有网格相关性,不同网格计算得到的液态水含量及局部水收集系数不一致,需要进一步考虑水滴运动受到上游效应影响时网格对欧拉法计算结果的影响。     

基于蒙特卡洛方法的冰晶撞击特性计算

航空发动机摄入冰晶易发生内部结冰、造成发动机推力损失,严重时引发航空事故。冰晶结冰数值模拟是研究冰晶结冰的重要手段,而撞击收集系数计算是冰晶结冰数值模拟的关键环节。基于NNWICE平台开展冰晶结冰前期的阻力系数、运动轨迹和撞击收集系数的数值模拟与计算研究,实现了基于蒙特卡洛方法的冰晶撞击收集系数的准确计算;在此基础上发展了拉格朗日框架下冰晶粒子运动、传热传质耦合数值计算方法;通过完成二维NACA0012翼型的相关算例,对比分析了不同阻力模型和冰晶相变对撞击收集系数的影响。结果表明阻力模型对冰晶撞击收集系数影响较小,冰晶相变过程会对冰晶运动轨迹及撞击收集系数产生显著影响,为进一步发展冰晶结冰数值模拟计算方法奠定了基础。     

面向过冷大水滴适航验证的结冰探测技术研究

过冷大水滴（Supercooled large droplet,SLD）结冰超出了常规防除冰系统的防护范围,是一种更为严重的结冰情形,极大地影响着飞行安全。在美国联邦航空管理局14 CFR 25.140条款中明确提到为保证飞机在SLD结冰条件下的安全运行,首先且必须要对SLD结冰气象环境进行探测。由于水滴破碎会改变水滴运动轨迹和表面撞击水分布,使水滴撞击极限变小,而水滴飞溅对结冰极限位置影响不大,因此本文通过研究大、小粒径水滴在飞机上收集范围不同的特点,提出一种可以满足SLD结冰适航符合性的结冰探测技术。该结冰探测技术采用常规结冰探测器,根据水滴收集范围将其布置于不同敏感结冰部位,可以实现全剖面飞行的结冰环境探测。     

折叠机翼变体飞机非对称变形控制效率分析

折叠机翼变体飞机变形量大,变形引起的气动参数变化显著,提出一种将非对称变形作为操纵输入的控制方案,研究非对称变形的控制效率和有效区间。首先建立能够完整描述变形过程的非线性动力学方程和气动力模型;然后基于非对称变形控制方法建立一种非对称变形操纵模型;最后通过与常规操纵面效率对比和仿真的动态响应总结出非对称变形操纵的最大变形操纵有效区间。结果表明:在较低飞行速度下非对称变形操纵效率高,非对称变形操纵能够在基准折叠角度90°附近提供最高的滚转操纵效率。     

结冰飞机着陆阶段飞行安全包线确定及操纵应对策略

结冰会导致飞行安全包线收缩、严重威胁飞行安全,研究结冰后安全包线变化对于操纵应对策略设计及提高飞行安全具有重要意义。以美国国家航空航天局（NASA）的项目飞机GTM（Generic Transport Model）为研究对象,通过对飞机的气动参数进行多项式拟合,建立了结冰飞机纵向通道的动力学模型。为了得到能随结冰程度变化的安全包线,将可达性分析理论引入到对结冰飞机着陆过程的安全性分析。提出将正向可达集与反向可达集的交集作为飞行安全包线,其中可达集的确定是基于水平集方法求解哈密尔顿-雅克比方程的最优解。最后针对不同程度的结冰条件进行了操纵时域验证,并提出了相应的操纵控制策略。研究结果表明,轻度结冰对安全包线影响较小,整个着陆过程的飞行状态始终能在最优控制指导下保持在安全包线以内;但对于重度结冰,飞行安全包线收缩严重,常规操纵已经很难使飞机达到着陆要求,需要进行飞行状态改出处理。研究结果为指导飞行操纵及实时包线保护打下基础。     

正常类飞机自然结冰试飞适航审定技术

民用飞机如果申请已知结冰条件下飞行（FIKI）,应按照除防冰相关适航规章条款要求进行表明符合性适航验证,其中自然结冰试飞是必须完成的重要环节。以Y12F飞机自然结冰试飞为例,以中国民用航空局（CAAC）和美国联邦航空管理局（FAA）发布的已知结冰条件下飞行相关指导文件为基础,结合Y12F飞机自然结冰试飞适航审定过程,分析总结了除防冰系统、测试设备及安装、试飞空域、结冰大气条件、试飞内容、试飞状态、试飞程序、试飞结果和分析以及结冰探测、结冰大气探测、冰积聚情况监测等方面的适航审定要求和关键技术,并对试飞过程中发现问题的分类、解决措施、设计改进方案以及完成的补充验证工作进行了阐述。经CAAC和FAA同步审查,Y12F飞机获得CAAC和FAA已知结冰条件下飞行的批准。构建的自然结冰试飞适航审定方法,成为FIKI适航验证重要实践指导性材料。     

带压缩性修正的离散涡方法在结冰数值模拟中的应用

离散涡方法（DVM）是一种无网格的涡运动算法,适用于解决易产生分离流的非定常问题,将其应用于结冰过程中的流场求解,在有效模拟分离流动的同时能避免冰形尖角对网格质量的影响。但离散涡方法基于不可压 N-S 方程,无法应用于预测可压缩流动下的结冰过程。本文在离散涡方法的基础上添加普朗特—格劳尔特压缩性修正,进行基于离散涡方法的可压缩流动下的数值模拟,并将其应用于翼型结冰预测;对流场分布、结冰冰形和结冰模型计算过程等仿真结果进行对比和分析。结果表明：引入压缩性修正后的离散涡方法能较好地模拟可压缩流动,与实验值相比,基于该方法得到的结冰数值模拟结果符合良好,对结冰数值模拟在工程上的应用提供了一定的参考。     

翼型测力天平的研制

为研究翼型结冰状态下的载荷变化规律,对翼型进行测力实验.研制一台三分量翼型测力天平.天平设计的突出难点为设计载荷极不匹配,阻力测量载荷很小.设计时,针对翼型载荷的特点,采用片梁组合方式测量阻力,四柱梁测量法向力及俯仰力矩,并采用有限元方法对结构进行了验证计算.主要介绍天平的设计、校准及应用情况.     

非对称喷管超声速流场的PIV实验与数值模拟

选取了跟随性能以及散射性能都上佳的二氧化钛粒子,利用了粒子图像测速仪(particle image velocimetry,PIV)对不同落压比下的非对称喷管(single expansion ramp nozzle,SERN)流场进行了测量,成功获取了相应流动结构,以及沿程上壁面压力分布和速度场等.在测得的图像中,可以清晰地观察到由于过膨胀产生的激波与膨胀波系,并进一步获得出口速度分布.运用计算流体动力学数值模拟方法,对实验状态下的非对称喷管进行了模拟.对比结果发现:PIV实验测得流场结构能很好地与其吻合,这互相验证了CFD计算和PIV实验的正确性及可靠性.