旋转反压与离心式压气机叶片相互作用数值模拟研究

由于旋转爆轰燃烧室具有自增压、热效率高的特点，将航空涡轮发动机的燃烧室替换为旋转爆轰燃烧室可进一步提高发动机性能。本文基于非稳态雷诺时均Navier-Stokes方法，采用k-ε湍流模型，针对旋转反压与离心式压气机叶片的相互作用开展三维数值模拟，研究旋转反压形成的前传压力波与扩压器及叶轮叶片的相互作用，分析旋转反压传播方向对前传压力波运动过程的影响。结果表明：旋转反压会在压气机内形成前传压力波，其与扩压器和叶轮的叶片相互作用，形成复杂的波系结构。当旋转反压传播方向与叶轮旋转方向相反时，前传压力波的压力峰值在扩压器内先下降后上升，在进入叶轮通道后迅速下降，无法到达叶轮入口。当旋转反压传播方向与叶轮旋转方向一致时，压力波波形发生明显变化，其与扩压器叶片相互碰撞导致强度明显下降，向上游传播距离更短。     

探针耙对跨声速环形涡轮叶栅流场的影响

采用数值模拟方法研究叶栅出口旋转总压探针耙对跨声速环形涡轮叶栅流场的影响,其中对比分析了探针耙在不同周向位置与无探针耙时的流场及叶栅性能参数。结果表明:常规总压探针耙会造成其相邻及上游叶栅通道流动减速、下游叶栅通道流动加速;探针耙尾缘燕尾型激波及其上游的正激波,与相临及下游叶片的尾缘燕尾型激波系互相干涉,形成了复杂的激波系和膨胀波系,造成相邻及上游若干叶栅通道堵塞,导致各叶片的载荷分布均不相同;叶栅各通道的叶片载荷、出口总压、端壁静压、出口气流角等参数偏离无探针耙情形,测得的叶栅出口流场与无探针耙时相差较大。     

尾缘冷却跨声速涡轮气动特性的数值模拟

为了明确跨声速涡轮中尾缘冷却措施的引入对其气动性能的影响,针对跨声速涡轮叶栅中压力面半劈缝尾缘冷却方式和尾缘全劈缝冷却方式条件下,叶栅性能和尾缘激波系结构的变化进行了数值研究。结果显示,两种尾缘冷却措施都降低了叶栅能量损失系数,压力面半劈缝尾缘冷却方式效果更好,最佳情况损失系数下降达到48%。冷却流量对作用效果有一定影响,存在最佳值。冷却措施的引入显著改变了尾缘激波系结构,尤其对PSEJ冷却方式,将尾缘激波系压力面分支由原本一道强激波分成三道弱激波。     

基于OpenFOAM的三维H2/Air连续旋转爆轰流场数值模拟

为进一步研究旋转爆轰流场特征,基于开源计算流体动力学软件OpenFOAM,采用9组分19步的基元化学反应模型,对H2/Air连续旋转爆轰流场进行了三维数值模拟,得到了旋转爆轰波稳定传播时燃烧室内部流场的详细结构,研究了燃烧室头部激波的传播特性,分析了旋转爆轰燃烧室的压力增益性能。结果表明:旋转爆轰波后的第一道反射激波在由燃烧室外壁面向内壁面传播过程中反射激波的高度增加并在靠近内壁面附近与滑移线交汇形成局部高温高压区域;旋转爆轰波在外壁面位置处相位约落后于内壁面0.003rad～0.15rad,其相位差随燃烧室曲率差的增加而增大;燃烧室头部反射激波数目受到曲率差和进气总压的影响,燃烧室曲率差增大,反射激波数目减少,进气总压增大,反射激波数目增多;燃烧室压力增益保持在0.3以上,在进气总压一定的条件下,压力增益随着燃烧室曲率差的增大有增加的趋势。研究结果揭示了三维旋转爆轰流场的精细结构和燃烧室头部激波的传播规律。     

旋转爆轰燃烧室高度对流场影响

为了研究旋转爆轰燃烧室几何结构对旋转爆轰波的影响规律,针对旋转爆轰燃烧室高度对爆轰流场结构和爆轰波传播特性问题,采用2维数值模拟方法分析了爆轰波在4种高度的燃烧室内的参数变化和波系特征。结果表明:在4种高度的燃烧室内爆轰波均能自持传播,且爆轰流场结构相似;当燃烧室高度为4 mm时,爆轰波后的反射激波压力较高;在一定范围内,爆轰波传播速度随燃烧室高度的增加而增大,在爆轰波高度为10、16和20 mm时,爆轰波波锋面的压力分别为3.93、4.61、4.69 MPa。     

大膨胀比跨声速涡轮流动结构及损失的数值研究

为了揭示跨声速大膨胀比涡轮损失的主要特点和两种不同尾缘冷却方式对损失的影响,以典型大膨胀比跨声速涡轮和跨声速叶栅为研究对象开展了数值研究。研究发现大膨胀比跨声速涡轮的主要损失是叶型损失,占到总损失的65%左右,尾缘激波损失是叶型损失的主要来源。尾缘全劈缝冷气入射通过提高尾缘基压区基压来减少尾缘膨胀波对气流的加速程度,从而降低最高马赫数和激波损失,尾缘压力面劈缝冷气入射通过改变叶片尾缘压力面激波波系结构,使原来的一道激波变成两道或者两道以上的弱激波,从而减少激波损失。两种尾缘冷气方式都有利于降低大膨胀比跨声速涡轮激波损失,但压力面劈缝冷气入射方式效果更为明显。     

跨声速涡轮叶栅激波损失控制方法

为了降低高负荷跨声速高压涡轮激波损失,发展了针对性的涡轮叶栅激波控制方法。针对吸力侧激波,提出可控膨胀设计概念,结合基于曲率的叶型设计方法,通过调整吸力面曲率分布以控制气流膨胀力度,减小了尾缘激波前马赫数,有效减弱了吸力侧激波强度和叶栅出口压力不均匀程度。针对压力侧激波,发展了消波设计方法,在吸力面的激波作用区域设计一鼓包型线,利用鼓包迎风面压缩波的预增压作用和外凸面膨胀波的消波作用,有效抑制了激波/边界层相互干扰,显著削弱了反射激波强度。可控膨胀设计和消波设计对叶栅尾缘两道激波的控制作用相互独立,可单独采用,当两种方法相结合时,吸力侧激波强度降低了29.66%,叶栅出口压力不均匀程度减小了29.28%,总压损失系数减小了12.11%。      

高负荷跨声速涡轮激波损失机理及控制技术研究

为了指导高性能、高负荷跨声速涡轮的设计,对其叶栅内的流场结构、尾缘波系结构、减小激波损失的机理及其控制技术进行了分析研究。结果表明:跨声速涡轮尾缘流场结构复杂,存在分离膨胀波、分离激波、基底区、再附激波、尾迹、吸力面反射波甚至激波边界层相互干扰等流动现象。通过采用收缩-扩张通道,喉道后采用直线型吸力面,减小吸力面尾缘弯折角、尾缘厚度和尾缘附近局部修型等措施,从而减弱激波强度,减小激波损失。     

大攻角范围超声压气机叶栅激波波系研究

对超声压气机叶栅的多攻角工况进行试验,利用纹影仪、油流试验及叶片表面等熵马赫数分布结果进行对比分析,观察到大攻角范围下叶栅激波波系结构发生了明显变化。为揭示激波结构变化原因,利用NUAA计算程序对叶栅进行仿真。研究发现,大攻角状态下叶栅通道中斜激波产生的原因,为前通道激波诱发附面层分离再附后,气流为沿叶片表面继续流动,从而形成斜激波;由于斜激波的增压降速,导致尾缘激波非常微弱甚至消失。     

大攻角范围超声压气机叶栅激波波系研究北大核心

对超声压气机叶栅的多攻角工况进行试验,利用纹影仪、油流试验及叶片表面等熵马赫数分布结果进行对比分析,观察到大攻角范围下叶栅激波波系结构发生了明显变化。为揭示激波结构变化原因,利用NUAA计算程序对叶栅进行仿真。研究发现,大攻角状态下叶栅通道中斜激波产生的原因,为前通道激波诱发附面层分离再附后,气流为沿叶片表面继续流动,从而形成斜激波;由于斜激波的增压降速,导致尾缘激波非常微弱甚至消失。     

ν-gap度量及其在飞行控制律评估中的应用

 传统控制律评估方法主要用于单输入单输出（SISO）系统，且对模型参数摄动考虑不够全面,针对这些不足，研究了ν-gap度量方法。在介绍系统广义稳定裕度相关概念的基础上，给出了ν-gap度量的定义、特点和性质以及近似摄动模型的计算，提出ν-gap度量评估控制律的步骤。实例结果表明，该方法不仅克服了上述传统评估方法的缺陷，而且还有根据所求的各摄动影响情况忽略影响小的元素，以减少计算量及可以找到最坏情况下的参数摄动组合等优点。     

基于飞机油箱模型形状特征油量测量切片步长选择方法研究

 在分析飞机数字式油量测量过程中目前广泛使用的切片法油量测量原理的基础上,针对现有的定步长切片法无法得到准确、可靠的燃油质量特性数据库的缺陷,结合对飞机油箱模型形状特征的分析,提出了基于飞机油箱模型形状特征的油量测量切片步长选择方法。此方法包括切片步长整体和局部选择两个过程,整体选择以实现相邻两切片平面所夹油箱模型体积近似相等为目的来确定切片步长,以体现油箱模型截面整体变化规律;局部选择以设计切片平面与截面突变平面重合或尽可能接近的方式,突出油箱截面的局部变化特征。实验结果表明：该切片步长选择方法较定步长方法能够建立更为合理、可靠的燃油质量特性数据库,从而提高了油量测量精度。     

Integrated Entry Guidance for Reusable Launch Vehicle

A method for the implementation of integrated three-degree-of-freedom constrained entry guidance for reusable launch vehicle is presented. Given any feasible entry conditions, terminal area energy management interface conditions, and the reference trajectory generated onboard then, the method can generate a longitudinal guidance profile rapidly, featuring linear quadratic regular method and a proportional-integral-derivative tracking law with time-varying gains, which satisfies all the entry corridor constraints and meets the requirements with high precision. Afterwards, by utilizing special features of crossrange parameter, establishing bank-reversal corridor, and determining bank-reversals according to dynamically adjusted method, the algorithm enables the lateral entry guidance system to fly a wide range of missions and provides reliable and good performance in the presence of significant aerodynamic modeling uncertainty. Fast trajectory guidance profiles and simulations with a reusable launch vehicle model for various missions and aerodynamic uncertain-ties are presented to demonstrate the capacity and reliability of this method.     

二维翼段颤振的μ控制

 采用超声电机作为作动器来实现含控制面的翼段颤振鲁棒抑制。针对作者设计的二维翼段颤振主动抑制系统，通过理论与实验相结合的方法，建立了考虑沉浮方向阻尼和作动器模型参数不确定性的控制系统模型，设计了μ控制器，并对控制器做了降阶处理。数值仿真和风洞试验表明，μ控制器可有效地抑制颤振的发生，将颤振临界速度提高23.4%。相对于H_∞控制器，μ控制器的控制效果和鲁棒性更好。     

The effect of inlet conditions on the flow and heat transfer in multiple rotating cavity with axial throughflow

This paper discusses experimental results from two different build configurations of a heated multiple rotating cavity test rig.Measurements of heat transfer from the discs and tangential velocities are presented.The test rig is a 70% full scale version of a high pressure compressor stack of an axial gas turbine engine.Of particular interest are the internal cylindrical cavities formed by adjacent discs and the interaction of these with a central axial throughflow of cooling air.Tests were carried out for a range of non-dimensional parameters representative of high pressure compressor internal air system flows(Re up to 5×106 and Rez up to 2×105).Two different builds have been tested.The most significant difference between these two build configurations is the size of the annular gap between the(non-rotating) drive shaft and the bores of the discs.The heat transfer data were obtained from thermocouple measurements of surface temperature and a conduction solution method.The velocity measurements were made using a two component,LDA system.The heat transfer results from the discs show differences between the two builds.This is attributed to the wider annular gap allowing more of the throughflow to penetrate into the cavity.There are also significant differences between the radial distributions of tangential velocity in the two builds of the test rig.For the narrow annular gap,there is an increase of non-dimensional tangential velocity V/Ωr with radial location to solid body rotation V/Ωr=1.For the wider annular gap,the non-dimensional velocities show a decrease with radial location to solid body rotation.      

Abnormal Shape Mould Winding

为解决网格化芯模的缠绕问题,本文提出了复合材料面片缠绕机理;接着详细分析了面片缠绕过程中的芯模凹曲面上纤维滑线和架空现象,应用微分几何曲面理论和空间几何理论,提出判据及其解决方案;最后,针对飞机发动机进气道的缠绕成型,编制缠绕控制程序并进行相应的实验,验证了面片缠绕方法的可行性。     

基于H∞性能指标的质量矩拦截弹鲁棒控制

 以所建立的质量矩拦截弹数学模型为基础，利用微分几何的反馈线性化理论，得到一个解耦线性系统，考虑到拦截弹的鲁棒性要求和3个滑块的协调控制问题，提出采用双回路的设计方法，内回路采用线性二次调节器（LQR），外回路采用考虑混合灵敏度问题的H_∞控制设计。仿真结果证明了该方法动态性能满足设计要求，同时对系统参数的不确定性具有较强的鲁棒性。     

航空发动机管路测量数据分割方法

 为了提高航空发动机管路测量数据的反求建模效率,提出了一种区域增长分割算法。该方法主要是利用管路表面在曲面索引系数映射、主曲率映射、高斯映射上的特性,并基于均值漂移算法和遗传算法提出了区域分割算法中种子区域的选择策略。然后利用区域增长分割方法实现了对管路测量数据的分割。经过仿真和实测数据验证,所提管路分割算法具有较好的分割质量和效率。     

航天器返回地球的气动特性综述

航天器返回地球的飞行过程中,气动特性是实现将宇宙飞行速度减到落地前速度、保证再入飞行得到有效控制以及再入防热安全可靠的关键因素。针对简单旋成体气动外形、半弹道式再入控制、烧蚀防热类返回航天器,综述了返回地球过程中变化的空气流域特性、航天器周围的气体绕流环境、空气与航天器作用产生的动力学与热效应等。系统地给出了该类航天器的再入气动特性参数与飞行性能的共性规律,包括:气动阻力与再入减速、气动升力与再入轨迹控制、配平攻角与飞行稳定性、气动加热与防热,以及再入过程中不同气动特性航天器、气象条件变化等对再入飞行性能的影响规律。为航天器开展返回飞行过程的跨流域气动性能工程研制提供设计参考。     

基于马赫数分布可控曲面外/内锥形基准流场的前体/进气道一体化设计

提出了一种高超声速飞行器乘波前体的外锥形基准流场设计方法,在锥面马赫数分布规律给定的条件下,通过有旋特征线法实现反设计,提高了基准流场设计的灵活性。该基准流场通过锥形"下凹"弯曲激波和波后等熵压缩波系压缩气流,可以在较短的长度内完成高效压缩。基于反正切马赫数分布外锥形基准流场设计的乘波前体具有较高的容积率,乘波特性良好且出口均匀,设计点时有黏升阻比为1.89。另外,基于该乘波前体和马赫数分布可控的内收缩进气道给出了一种双乘波的前体与进气道一体化设计方案,实现了内外流分别独立乘波,充分发挥了乘波前体和内收缩进气道的各自优势。