间隙高度对超声速膨胀器流场及性能的影响

采用三维雷诺平均N-S方程和标准k-ε湍流模型,对不同间隙高度超声速膨胀器的流场和性能进行了数值研究,结果表明:间隙高度显著影响三维流道内的局部流动特性和超声速膨胀器的整体性能,随间隙高度增加,气流最高相对马赫数降低,高速区范围逐步缩小;泄漏涡增强,尺度变大,横向和径向运动明显,泄漏损失增加,但激波及激波附面层相互作用的损失降低;超声速膨胀器的膨胀比先增大后减小,等熵绝热效率持续降低.下端壁、吸力面附近低能流体之间以及与壁面的摩擦损失和间隙泄漏损失是有间隙超声速膨胀器三维流道内损失的主要来源,超声速膨胀器的间隙高度宜在0.9%h₀~1.5%h₀之间选取.      

向心涡轮进气结构的气动性能及损失机理

针对某向心涡轮,采用二维流动分析方法设计矩形截面蜗壳,同时采用商用计算流体动力学软件CFX对带蜗壳的向心涡轮流动损失进行数值计算.将计算结果与原型设计中带有集气室的向心涡轮计算结果进行对比.结果表明:蜗壳内流动损失要小于集气室内的流动损失,向心涡轮采用蜗壳后流道内流场有明显改善,效率有所提高.采用设计的矩形截面蜗壳,向心涡轮的功率提高1.7%.通过内部流场的分析,揭示了内部流场结构和损失机理,为向心涡轮的设计和优化提供了一定的参考.      

两种空气/煤油燃气发生器富油燃烧组织方案对比

根据空气/煤油富油燃烧的特点,提出了两种空气/煤油燃气发生器富油燃烧组织方案,设计了采用钝体稳定火焰和旋流空气、二次喷注空气稳定火焰的两种燃气发生器.为了对比两种方案的点火和燃烧特性,对两种燃气发生器进行了一系列热试,结果表明余氧系数是燃气发生器最重要的工况参数.随着余氧系数的增加,燃气发生器的状态逐渐从启动失败变为中途熄火,最终呈正常启动状态.采用钝体稳定火焰的燃气发生器稳定工作的余氧系数边界为0.518,采用旋流空气和二次喷注空气稳定火焰可将该边界延伸到0.237,极大地扩大了燃气发生器的工作范围.与钝体稳定火焰的燃气发生器相比,旋流空气和二次喷注空气稳定火焰的燃气发生器的富油燃烧的燃烧效率提高了20%.两种方案结构复杂性相当,旋流空气和二次喷注空气稳定火焰的燃气发生器不需要冷却火焰稳定器,可提高燃气发生器的工作时间.      

慢车燃油分级对高温升燃烧室燃烧效率的影响

在双级旋流多点喷射直接混合(TAMDIM)高温升燃烧室单头部试验件上进行了燃烧效率试验,对比了两种不同主副级喉道间距尺寸下的慢车主、副分级供油方案及副油路单独供油方案的燃烧效率,分析了不同主、副级当量比(副油路供油比例分别为40%,56%,65%,100%)及喉道间距对燃烧效率的影响.试验表明:主副级采用分级供油策略时,喉道间距对燃烧效率有显著影响,喉道间距设计为19.3mm,副油路供油比例控制在40%~56%,可以获得比单独副油路供油更好的燃烧效率.在扇形试验件上进行了试验验证,进一步验证了慢车分级供油策略在TAMDIM高温升燃烧室上的可行性.      

用于压气机流动计算的3种模型比较

为了更好地对压气机流动进行模拟,在课题组自行开发的结构化有限体积解算器上实现了用于压气机流场计算的混合平面法、谐波平衡法及相滞后法.以NASA Stage 35为例,对3种方法的计算结果进行了比较分析.结果表明:相滞后法的计算精度最高,混合平面法的计算精度最低;相滞后法与半环的双时间推进法结果相近,计算速度提高了20倍;相比混合平面法,谐波平衡法能准确地模拟动静叶间的非定常干涉及进出口参数变化;在谐波阶数达到5阶后,谐波平衡法计算得到的结果不随阶数变化,且与相滞后法的结果基本吻合;混合平面法的计算效率远高于另外两种方法,相滞后法与谐波平衡法在谐波阶数为5阶时的计算效率相当.      

凹腔/支板结构亚燃冲压燃烧室性能

为了避免基于凹腔火焰稳定器的亚燃冲压燃烧室壁面喷注时燃料与主流空气掺混非均匀性问题和提高燃烧室的性能,提出在亚燃冲压燃烧室中使用支板喷注代替壁面喷注的方案,数值模拟了凹腔/支板结构亚燃冲压燃烧室中燃料分布及流场结构,并分析了支板结构对燃料空气混合及燃烧室性能的影响。研究表明:支板虽然使燃烧室出口的总压恢复系数相对于壁面喷注方式下的降低了63%,但能使燃料均匀分布于整个流道内,增强了燃料与空气掺混,使燃烧室出口的混合效率和燃烧效率分别提高了21.4%和20.5%。燃烧效率的提高弥补了采用支板导致的燃烧室内气流的额外总压损失所带来的机械能损失,使得支板喷注时燃烧室出口的比冲提高了39.6%。因此,在亚燃冲压燃烧室中设置凹腔/支板结构,有利于提高燃烧室整体性能。      

离心压气机失速模式及自循环机匣处理的作用机制

通过对离心压气机不同转速下实壁机匣和处理机匣两种状态流场的数值模拟,研究了失速模式随转速的变化及不同失速模式下机匣处理作用机制的转变.结果表明:随转速的降低,失速模式由100%设计转速时的扩压器失速逐步向导风轮失速过渡,导风轮失速也经历了由80%设计转速时的端壁失速和叶片前缘失速向50%设计转速的通道大范围完全失速的发展过程.自循环机匣处理对离心压气机失速裕度和效率的影响都与失速模式紧密相关,其扩稳作用仅对80%设计转速时的叶片前缘失速有效.在80%左右设计转速,机匣处理使主流效率有朝大流量工况点变化的趋势,从而使效率降低.当转速低至50%设计转速时,主流通道形成大范围回流区,机匣处理的二次回流可以减少叶轮对通道回流区气流的做功量,从而提高效率.      

应用傅里叶时间谱方法求解二维跨音速流动问题解的精度研究

为了从时间离散精度的角度评估傅里叶时间谱方法的实际计算效率,对用该方法求解二维跨音速流动问题的解进行了精度研究。通过求解Euler方程模拟俯仰振动NACA0012翼型的周期性非定常流动。结果表明:傅里叶时间谱方法可以使用很大的时间间隔较为准确地模拟有激波存在的周期性跨音速流动;对于此类流动问题,相比二阶向后差分公式法,傅里叶时间谱方法在预测翼型表面压力系数积分量特别是升力系数的时间变化规律方面具有明显的计算效率优势;用于模拟翼型表面压力系数自身的时间变化规律时,傅里叶时间谱方法具有不低于二阶向后差分公式法的计算效率。     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。