二维翼型结冰数值计算

翼型结冰计算是随着结冰时间步长变化的准定常过程,因此在每个时间步长后都要更新各个模块的计算工作,如此反复直到达总结冰计算时间.在各模块计算中,分别采用了泊松方程法生成网格、求解雷诺时均方程计算流场、拉格朗日追踪法求解水滴轨迹、基于Messinger模型计算冰形、采用变距滤点法重构边界.引入动态链接库技术将各模块集成为结冰计算平台,并将三种冰形(楔形冰、槽状冰及混合冰)计算结果与国外软件及试验数据进行了对比.      

基于准三维冰形算法的结冰预测技术及应用研究

机翼结冰是威胁民用飞机飞行安全的主要因素之一,准确预测机翼的结冰过程和特点,可为防除冰设计提供依据和支撑。提出了利用二维模拟三维的准三维冰形算法,该算法在保证结冰预测精度的同时,极大地降低了型号结冰计算成本;给出了顺气流和垂直于前缘两种典型民用飞机后掠翼机翼的二维剖面截取方法,并将采用该方法的计算结果与三维后掠机翼的结冰预测结果进行了对比分析。结果表明,相对于顺气流方式,垂直于前缘方式的准三维冰形算法与三维后掠机翼的结果更为接近。     

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

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

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

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

基于鲁棒伺服LQR的结冰飞机纵向控制律重构方法

由于防/除冰系统总会出现故障或者除冰不彻底,因此仅依靠防/除冰系统实现结冰条件下的安全飞行并非完全可靠,研究结冰后飞机控制律重构对飞机操纵安全和飞行安全极其重要。针对飞机的纵向运动建立了结冰影响模型和纵向动力学模型,采用鲁棒伺服线性二次型调节器(LQR)最优控制设计了飞机结冰后空中飞行纵向控制律,模拟了飞机在俯仰姿态保持模式下遭遇不同严重程度结冰后的动态响应特性,并与常规PID控制进行对比。结果表明,所设计的控制律能够有效改善结冰飞机的飞行性能和飞行品质,准确跟踪给定的俯仰角指令,且抗干扰能力、动态性能以及鲁棒性均优于常规PID控制。为飞机结冰后的重构控制问题和自动飞行控制,提供了新的思路。     

待机状态下机翼结冰的快速计算方法

飞机在飞行过程中遭遇结冰气象,尤其是过冷大水滴结冰条件所导致的机翼结冰,将严重影响飞机的气动性能与操纵品质,从而导致飞行故障或飞行事故,是飞机飞行安全的重要影响因素,是飞机研制中必须重点解决的难题。虽然针对结冰已建立有多种计算仿真程序,但快速获得结冰模拟冰形,包括过冷大水滴条件下的结冰模型是工程师们一直追求的解决方案。利用前向多层神经网络,尝试建立了一种针对待机状态下的机翼结冰模型。此方法基于坐标转换原理,将获取的采用直角坐标表达的标准翼型的机翼结冰数据转换为对应的极坐标表达的数据,以飞行参数、气象参数和极坐标角度作为输入,极坐标模值作为输出训练构建出此种机翼结冰模型。使用所建立的模型,仿真预测结果表明该计算方法具备快速性和精确性,其计算精度可满足实际要求。     

计算流体力学中的网格生成方法及其应用

讨论了多种网格生成方法及其流场的数值模拟。给出了一种生成现代双立尾战斗机外形的分区对接网格和Ｅｕｌｅｒ方程分区解。给出了一种带副翼偏转的三维机翼分区网格及绕流的Ｅｕｌｅｒ方程解。以二维多段翼型为对象分别讨论了分区重叠结构网格方法、非结构网格方法和自适应的笛卡尔网格法。     

一种交互式分区网格的生成方法和软件

介绍了一个具有交互性能的网格生成系统--BIGG系统。它是采用Windows环境下的C++语言与Fortran语言混合编程技术开发的,具有界面友好、显示灵活、交互性能强等特点。网格生成实例表明,采用BIGG软件可以缩短网格生成的周期,提高网格生成的质量。     

斜网格修正的低速计算方法及应用

采用无坐标变换的一种有限体积法对RANS方程进行离散求解,离散方法考虑了对网格非正交情况的二阶修正,通过构造一个动量插值实现了同位网格下的压力修正算法,并发展了相应的分区算法.利用以上方法对轿车绕流和三峡坝区主流摆动问题进行了数值计算,计算结果与实验吻合较好.     

基于网格缩小的非结构网格梯度重构及制造解精度测试与验证

网格收敛性研究是验证与确认中的重要手段,本文针对一般非结构网格难以自相似加密的问题,实现了基于网格缩小的精度测试方法。在典型的各向同性和各向异性拉伸网格上,利用网格缩小精度测试方法分别考察了梯度重构精度以及制造解流动模拟精度,将网格缩小精度测试结果与传统的网格加密精度测试结果以及理论分析结果进行对比,验证了网格缩小精度测试方法与网格加密精度测试方法及理论分析的一致性。最后将网格缩小精度测试方法应用到各向异性拉伸网格粘性制造解精度测试中,得到了预期中的结果,并初步考察了梯度重构方法、网格类型对数值模拟精度的影响,显示出网格缩小精度测试方法在一般非结构网格精度测试上的优势,具有一定应用前景。     

翼型表面粗糙度对结冰的影响分析

飞机的结冰表面出现微小的凹凸不平,即形成所谓的粗糙度,对飞机气动性能产生一定的影响。在考虑表面粗糙度时,针对对流换热系数的计算建立了热力学模型,同时对只有单个粗糙微元的表面进行了流场计算和表面对流换热效果分析;然后利用边界层积分方法,对某一翼型光滑表面和不同程度的粗糙表面分别进行了流场计算、对流换热系数计算以及结冰冰形的模拟仿真。结果表明：结冰表面粗糙度很大程度增强了表面的对流换热效果,导致在翼型前缘位置结冰厚度更大、冰角突出更为明显并且距离驻点区域更近。     

飞机机翼积冰研究方法和进展分析

通过对目前国内外研究飞机机翼积冰三种方法的分析,指出了它们各自优缺点,重点介绍了现行国内外研究的进展和新趋势,提出了研究方法的改进方案和进一步研究的建议。     

应用混合网格分析复杂积冰翼型气动性能

通过发展一种点-点对接的结构/非结构混合网格生成方法,避免了复杂积冰外形难以划分结构化网格的问题,并对复杂积冰翼型的气动性能进行了分析计算.在靠近积冰边界的内层采用了非结构网格以拟合复杂的边界,在非结构网格外采用结构化网格以节省存储空间和计算时间.为了检验网格对选用湍流模型的影响,整个流动区域分别采用SST和SA湍流模型,求解了雷诺平均N-S方程.数值计算结果表明SST模型更适于模拟复杂分离流动,积冰对翼型的气动性能造成了严重的影响.     

机翼积冰数值模拟

采用计算流体力学（CFD）方法对三段翼和MS-317后掠翼进行了积冰数值模拟研究,基于欧拉两相流理论,应用分区算法,对空气流场和水滴流场进行了数值求解,得到水滴收集率.求解3-D积冰模型得到积冰量.将MS-317的冰形计算结果与实验数据以及软件LEWICE3D的计算结果进行了对比.结果显示:三维积冰模型具有较好的精度,虽然预测的冰形与实验数据有一些差异,最大厚度方面,误差在13.1%以内,但是积冰的类型与体积和实验数据基本一致.      

机翼曲面成型技术对气动特性影响的研究

增升减阻一直是空气动力设计工作者追求的目标。机翼设计是飞机气动设计的核心,而机翼控制剖面间过渡方式是机翼设计中十分有效的增升减阻措施。首先对飞机翼根两剖面没有添加约束的机翼外形进行了详尽的气动分析计算,针对分析结果出现的问题,重新确定了设计思路和方案,探讨了在翼根两剖面翼型和位置都不变的情况下两翼型之间的过渡对机翼气动性能的影响,研究了如何过渡才能改进飞机的气动特性,最后将三种配置的设计结果进行了对比,获得了翼根曲面过渡几何参数影响的初步规律。对比结果表明,以加入过渡约束的上凸与线化过渡的结果取得了十分理想的改进结果。     

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

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

三角翼大迎角绕流特性数值模拟的网格处理技术研究

研究网格对三角翼大迎角绕流特性数值模拟的影响.结果表明,网格的生成需要结合流动现象,对网格拓扑结构和网格点分布进行选择与搭配.C-H型网格适宜模拟尖前缘分离涡流态,法向网格在一定范围内应等距增长,沿流向逆压梯度较大的区域内适当增加网格点,尾迹区网格则应做上翘处理.Euler方程具有模拟三角翼旋涡及预测涡破裂特性的能力,但对二次涡等粘性引起的流动细节把握能力不足.利用层流假设的N-S方程,通过合适的网格,也可得到满意的计算结果,但对涡破裂后的强烈非定常湍流流动模拟能力不足.采用旋涡螺旋度可准确反映主涡与二次涡流动,描述旋涡的破裂现象.用轴向速度迅速减小并小于来流速度的点作为涡破裂判据似应更合理.     

基于结构动网格的无人机地面效应研究

为了更好地分析地面效应对无人机气动特性的影响,采用一种新的结构动网格生成技术,运用CFD手段对某无人机的地面效应进行了研究。在一套静态网格CFD计算基础上,运用结构动网格生成技术,获得非定常的网格移动,得到一系列的不同攻角、地面效应高度的CFD网格及初始流场。在此基础上进行定常流场求解,计算得到地面效应的定常CFD结果。最后对某双尾撑布局无人机的地面效应进行了CFD计算研究,对相应的全机气动特性的影响进行了分析。计算结果表明：地面效应对平尾的影响大于对机翼的影响并使平尾的失速迎角提前;该方法是可行的,具有一定的工程应用价值。     

NACA4412翼型分离涡计算的混合网格策略研究

以NACA4412翼型为例,针对低速大攻角含有分离涡的复杂流动,提出了一种采用推高局部附面层网格的方法来改善分离涡模拟的新网格生成策略。数值算例表明,这种推高各向异性单元区域总层高的混合网格能够较合理地捕捉到分离流动涡区的位置和大小,数值实验和分析表明,这种网格生成策略是合理有效的。