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

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

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

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

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

为明晰结冰风洞中过冷大水滴（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℃。     

直升机旋翼翼型需求分析及技术发展展望

与机翼翼型相比,旋翼翼型在前飞时来流速度和迎角处于大范围变化之中,需要针对多种气动环境进行综合权衡设计.针对常规直升机和新构型高速直升机,综合考虑多种飞行条件,分析了旋翼翼型设计的性能需求.总结了先进旋翼翼型系列的发展及应用情况,以OA和TsAGI翼型为例分析了旋翼翼型的气动性能特点和发展趋势.从旋翼翼型指标分解、动态气动特性计算与试验、气动优化设计等方面介绍了旋翼翼型技术的发展现状,并对未来旋翼翼型技术的发展进行了展望.     

不同舵翼操纵面形式的水动力及艇后流场特性数值模拟研究

为研究共翼型舵和非共翼型舵两种操纵面形式的艇后流场特性及螺旋桨推进特性,基于SST(Menter)湍流模型建立了SUBOFF标准潜艇模型尾流场数值预报模型。经试验结果验证,所建立数值模型进行潜艇尾流场及螺旋桨推进性能预报具有较高精度。将SUBOFF潜艇模型的水平舵改进为共翼型舵及非共翼型舵,对艇后流场及水动力性能进行了预报。数值结果表明：在舵角小于10°时,共翼型舵使潜艇的俯仰力矩和垂向力相对非共翼型舵提升20%以上;在舵角超过10°时,共翼型舵的水动力优势随着舵角增大而减小。尾部流场预报显示：共翼型舵在小舵角时可以有效的消除舵翼结合处的涡流,同时共翼型舵可以有效的降低桨盘面伴流的不均匀性,对尾流品质的改善效果优于非共翼型舵。螺旋桨计算结果显示：共翼型舵螺旋桨推进系数在大部分舵角下都小于非共翼型舵,在舵角为20°时,共翼型舵相对非共翼型舵推力系数下降3.5%,扭矩系数下降2.4%;同时共翼型舵的桨盘面流场均匀度要优于非共翼型舵,舵角为5°时,共翼型舵桨盘面处流场不均匀度相对于非共翼型舵要降低7.1%,舵角为25°时则降低25.1%。     

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

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

跨声速层流翼型的混合反设计/优化设计方法

跨声速层流翼型设计须兼顾优良的超临界特性和自然层流特性,因而对设计方法提出了更高的要求。针对现有反设计方法和直接优化设计方法的不足,发展了一种适用于跨声速层流翼型的混合反设计/优化设计方法。该方法引入了基于经验的局部流场特征作为反设计目标,翼型性能指标作为直接优化设计目标,然后加权形成了混合反设计/优化设计总目标,并同时考虑了气动和几何约束。优化算法采用基于自适应并行加点技术的代理优化,流动数值模拟采用耦合基于线性稳定性理论的e^N转捩自动判定的雷诺平均Navier-Stokes（RANS）方程求解器。针对现代中短程民用客机需求,以NPU-LSC-72613翼型为基准,开展了层流翼型减阻的混合反设计/优化设计。分别将局部目标压力分布、总阻力作为反设计和直接优化设计目标,得到了较好的优化结果,验证了方法的有效性。经过2轮优化结果显示混合反设计/优化设计总目标显著下降。所设计翼型吸力面局部压力分布与目标压力分布基本一致,总阻力下降15.5%;吸力面和压力面层流范围均大于55%倍弦长,激波强度显著减弱,说明所设计翼型同时具有优良的超临界和层流特性。将所设计翼型配置到机翼上,通过三维数值模拟进行校验,结果显示所设计跨声速层流机翼升阻比提高了6.64%,在一定升力系数范围内,气动性能均有显著提高,验证了所设计跨声速层流翼型在机翼设计中的适用性。     

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

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

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

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

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.