叶轮机械三维不可压Euler方程的直接解法

本文用近似因式分解交替方向隐式格式直接求解叶轮机械中原参数三维不可压Euler方程组。在连续方程中引入“拟压缩性”之后,将原来类型不确定的控制方程变成双曲型方程,给定初边值之后用时间推进法求得稳定解,对一个单级压气机转子内部流场进行了计算,结果与实验值比较符合。     

混合解析/数值方法及其在湍流数值模拟上的应用

针对带小时间尺度的源项的方程描述的流动问题,提出了混合解析/数值方法。混合解析/数值方法的基本思想是:分裂原始方程组为对流-扩散部分的偏微分方程和源项的常微分方程。偏微分方程采用合适的数值方法求解,而常微分方程采用解析方式积分。模型方程的理论误差研究表明,混合方法提高了源项处理的精度,降低了混合方法的整体数值误差。分析同时表明,基于时间分裂的算法在求解含源项双曲系统的定常类型问题,会存在数值振荡。为此发展了非分裂方式的混合解析/数值方法,在湍流模型数值计算中提高了数值稳定性,而且加快了计算的收敛速度。     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

环下润滑结构对径向收油环收油效率影响的数值计算研究

为研究环下供油润滑系统中的流动过程和收油环收油效率的影响因素及其影响规律,基于VOF(Volume of fluid)两相流模型和滑移网格模型的数值计算方法,建立了环下润滑系统内部流场的计算模型。通过数值计算,获得了收油环的收油效率并验证了数值计算结果的有效性,分析和探讨了不同结构参数对收油效率的影响。计算结果表明:利用本文数值计算方法得到的收油效率和流动现象与已有实验和数值计算结果均一致;供油喷嘴角度和收油环外径均主要影响滑油射流冲击收油叶片的飞溅量和收油叶片切割滑油量,收油效率取决于二者的相对大小,θ=0°,θ=9°和r0=2mm,r0=6mm时收油效率取极大值;收油叶片数量和供油喷嘴与收油叶片间距离分别影响收油叶片切割滑油射流的时间间隔和滑油射流在到达收油叶片前的流动,不影响滑油射流冲击收油叶片的飞溅量和收油叶片切割滑油射流的量,对收油效率的影响较小。     

高超声速飞行器宽速域翼型/机翼设计与分析

对于现代高超声速飞行器的设计而言,除了需要保证高超声速的性能外,还必须兼顾满足工程需求的亚跨超声速特性。首先,采用雷诺平均Navier-Stokes（RANS）方程流动求解器,结合基于Kriging模型的代理优化算法,开展了高超声速飞行器宽速域翼型的优化设计研究,设计出了一种下表面具有双"S"形特征的新翼型。综合性能评估结果表明,该翼型相比于常规的高超声速翼型,在跨声速和高超声速下具有更加优良的气动特性;其跨声速状态下的升阻比达到78.9,高超声速状态下的升阻比达到5.94,能够实现宽速域内良好的综合气动性能。其次,开展了仿德国"桑格尔号"（SANGER）空天飞机运载机机翼的气动特性研究,对配置宽速域翼型与常规高超声速翼型的机翼进行了气动力特性综合对比分析。结果表明,配置新翼型的机翼在宽速域范围内整体气动性能更优,说明所设计的宽速域翼型在三维机翼上也具有一定的实用价值。     

基于非结构网格流场超大规模并行计算

大规模并行的计算流体力学已成为现代航空工业研发的核心手段之一。基于非结构混合网格和有限体积法,发展了适用于工业级复杂外形气动计算的并行流动数值模拟方法。文中首先介绍了紧致数值离散格式、基于Metis的分布式多核系统网格分区技术、并行边界虚拟单元技术和MPI并行实现等相关算法。采用网格量相对较小的旋成体构型绕流模型对比分析多核并行计算结果与单核计算结果以验证并行计算的正确性,比较了不同并行规模下并行效率和残差收敛情况。然后通过对上亿网格单元的运输机复杂构型绕流模拟,开展并行效率的测试,结果表明,本文方法并行加速性能高,直到多达18816核并行效率都保持在80%以上。     

大型水陆两栖飞机侧风起降的滑流影响分析

大型水陆两栖飞机AG600的动力装置为安装在机翼上的4台同向旋转涡轮螺旋桨发动机,针对1：15缩比模型带动力风洞试验显示的螺旋桨滑流对侧风起降状态的偏航力矩不稳定影响,对全机带动力风洞试验模型进行了大规模并行非定常数值计算,再现了风洞试验现象,通过流动机理分析明确其产生原因主要是左侧滑时右外翼分离和垂尾背鳍涡破裂,这些原因和数值模拟的准确性也为后期的风洞试验所证实。考虑到模型风洞试验中尺度限制造成的低雷诺数和高螺旋桨转速,为保证飞行安全,继续采用该非定常方法对全尺寸飞机真实侧风起降状态进行了详细数值分析和偏航稳定性评估。研究结果显示,在飞行雷诺数和螺旋桨转速下,相同侧风范围内风洞试验显示的流动不稳定因素基本消失,偏航稳定性允许的侧风范围明显增加。本研究实现了四发螺旋桨飞机起降状态横向气动特性的滑流影响非定常数值分析,建立了基于计算流体力学的风洞与飞行雷诺数效应的相互关系,进行了偏航稳定性的虚拟试飞评估,研究成果也为AG600飞机的首飞和飞行试验所验证。     

基于代理模型的大型民机机翼气动优化设计

先进的气动优化设计思想与方法,对于提升大型民机气动与综合性能具有至关重要的意义。探讨了大型民机超临界机翼气动优化设计的基本准则和要点,并结合代理优化算法,提出了一种面向工程应用的多轮次高效全局气动优化设计方法。首先,通过一系列解析函数/翼型/机翼优化测试算例进行了验证。其次,将所提出的方法与人工修型相结合,开展了针对宽体客机超临界机翼的两轮气动优化设计,使其气动性能得到显著改善。最后,采用不同的雷诺平均Navier-Stokes（RANS）方程求解器对安装优化机翼的全机巡航构型进行了典型状态的气动性能综合评估。研究结果表明,所提出的代理优化算法具有很高的优化效率、较强的约束处理和全局优化能力;将所发展的基于代理模型的多轮次气动优化设计方法与人工修型相结合,能够获得满足设计要求的气动外形,验证了该方法在大型民机超临界机翼气动设计中的有效性和工程实用性。     

翼身融合布局飞机机体-发动机气动干扰效应

翼身融合（BWB）布局飞机是一种相对较新的飞行器概念,具有商业运输飞机的潜在用途。研究表明,翼身融合布局客机可获得比常规布局客机更好的性能。但是,由于各方面限制,BWB飞机不宜使用传统的机翼安装或机身安装的发动机布局,发动机背部安装成为首选布局。然而,背部安装发动机容易产生激波、分离、进气畸变等空气动力干扰问题,发动机与机身一体化的气动设计已成为BWB飞机发展的关键技术。针对BWB布局飞机的机体和发动机之间的气动干扰进行了数值研究,结果表明,背部安装发动机对BWB布局飞机的全机气动特性和发动机本身的推力性能都会产生较为明显的影响。     

低声爆静音锥设计方法研究

 声爆精确预测及低声爆设计方法已成为新一代军民用超声速飞机研制过程中必须解决的关键难题之一。基于计算流体力学(CFD)、波形参数法以及MARK-VII方法构建了高精度声爆预测方法,利用该方法对低声爆静音锥的设计展开研究。研究结果表明,静音锥的设计存在临界长度,静音锥长度小于临界长度时,静音锥产生的激波在传播过程中与机头弓形激波合并,静音锥无法起到降低声爆的作用;静音锥长度大于临界长度时,声爆水平也会略有上升。静音锥临界长度随飞行高度和飞行马赫数的变化而变化,可以根据实际飞行状态采用可伸缩设计,达到最佳的低声爆效果。多级静音锥利用多道弱激波取代机头强弓形激波,其声爆水平较单级静音锥也更低,同样,多级静音锥每一级的长度也要达到临界长度。不同静音锥头部形状产生的脱体激波形状不同,脱体距离也不同,导致阻力系数以及静音锥壁面温度有所不同,但静音锥头部形状对远场声爆信号的影响并不明显。采用静音锥的低声爆方案与原始方案比较,声爆水平得到大幅降低,阻力系数略有上升。