自吸气脉冲爆震发动机超声速进气道气动性能的数值研究

采用有限体积法计算了某轴对称超声速进气道,在不同幅度和频率的出口扰动压强下的进气道内结尾正激波的运动情况,得出了进气道内结尾正激波运动特性和扰动压强的频率和幅度的关系。在计算中,本文采用了多块结构化网格,控制体积的界面无粘通量采用三阶迎风格式插值获得,同时采用了minmod通量限制器,以确保在激波处的解的物理特性;扩散通量采用二阶中心差分格式插值获得。定常计算采用当地时间步法,非定常计算采用双时间步法。离散的代数方程采用交替方向迭代法求解。     

超燃冲压发动机隔离段流动特性研究

隔离段是双模态超燃冲压发动机实现双模态和模态转换的一个重要部件,同时它把进气道和燃烧室隔离开,以防止燃烧室工作对进气道干扰,引起进气道不启动。由于隔离段内流场的复杂性及其广泛的工程应用前景,隔离段内的流动特性引起了人们广泛的关注。本文利用计算流体力学软件Fluent,采用k-ε湍流模型,近壁面采用非平衡壁面函数方法处理,控制方程的离散采用二阶迎风格式,对隔离段的流场进行了数值模拟,并分析了出口反压、来流马赫数等参数对隔离段内流场的影响。     

超燃冲压发动机隔离段流场数值模拟

隔离段是双模态超燃冲压发动机实现双模态和模态转换的一个重要部件,同时,它把进气道和燃烧室隔离开,以防止燃烧室工作对进气道干扰,引起进气道不启动。由于隔离段内气体流动的复杂性及其广泛的工程应用前景,隔离段内的流动特性引起了人们的广泛关注。本文使用隐式TVD格式,采用雷诺应力湍流模型,近壁面采用非平衡壁面函数方法处理,控制方程的离散采用二阶迎风格式,对隔离段内的流场进行了数值模拟,并分析了出口反压、来流马赫数对隔离段内流场的影响。     

超声速进气道湍流流场的数值模拟

采用有限体积法数值模拟了超声速进气道的湍流流场,改进了离散纳维尔－斯托克斯方程中粘性项的方法,使之具有较好的守恒性,研究使用了变系数的隐式残差光顺等加速收敛技术,并讨论了超音来流、收敛型进气道复杂流场初场的给定方法,所得到的数值解正确地展示了流场的复杂激波结构。     

二维超／高超声速进气道流场数值模拟

对超／高超声速三级压缩进气道流场进行了数值模拟,来流马赫数为4、6,进气道内流动为层流状态,根据二维Navier-Stokes方程,采用二阶精度Roe格式进行离散。按照流场特点,合理地设计网格分布及调整不同黏性范围的熵修正,防止了壁面附近过大的数值耗散,使计算结果更加合理。在进气道模型的各级压缩折转角处,获得了清晰的激波结构,在进气道内部的各种波系的相交、反射和激波诱导的边界层分离等现象都得到合理的描述。计算得到的压力分布,在各级压缩斜板上同简单波理论结果十分接近。用本文方法计算了另一个二级压缩进气道,沿上、下壁面的压力分布与试验比较符合得较好。     

基于SIMULINK的超声速进气道仿真方法

将进气道一维流动的偏微分方程离散为一系列可积分的线性常微分方程,基于此创建了超声速进气道仿真模型,提出了一种基于SIMULINK的进气道仿真方法,模拟了来流参数和燃烧室反压扰动时流场内各个参数的变化规律。对一个一维变几何超声速进气道进行仿真,获得了结尾正激波的稳态特性及动态特性。该模型定义并输出了激波位置,有效地将一维超声速进气道的分布参数问题转换为零维问题,解决了系统分析与控制系统设计中无法实现集总参数的困难。     

飞机／进气道系统内外流场计算

为了对飞机/进气道系统内外流场进行数值模拟,本文提出并研制了计算定常(或非定常)外压式飞机/进气道系统内外流场的计算方法和程序。它适用于任意几何外形的飞机/进气道在亚音速、跨音速或超音速来流条件下的内外流场和气动参数计算。求解的主管方程是在一般曲线坐标下的 Euler 方程、简化/雷诺平均的 Navier-Stokes方程。用隐式时间差分和近似因式分解方法离散和求解.针对飞机/进气道的复杂流场,本文采用了贴体网格生成技术和分区域方法.对于湍流流动,把 Baldwin-Lomax 湍流模型推广到飞机/进气道内外流场计算并使用了壁面律计算公式。应用上述方法对绕翼型的亚、跨音速、绕可压缩拐角的超音速湍流以及飞机进气道系统内外流场进行了计算,计算结果与实验数据和其它计算结果符合较好。     

弯曲进气道三维边界层预测

本文研究了弯曲进气道内三维可压缩湍流边界层的预测方法。采用非正交曲线坐标系中三维边界层积分方程,并引入掺混系数微分方程,用滞后掺混法求解,数值方法为预估-校正法,用三次样条函数解决了离散数据的插值、求导的精度问题。所编程序可在微机上运行。实例计算清楚地揭示了弯曲进气道内强横向流动对边界层发展的影响,计算所得的进气道总压恢复系数与实验值的误差为1～2％。     

锯齿形转捩片触发高超声速进气道边界层转捩的大涡模拟

为促进锯齿形转捩片在高超声速进气道中的应用，以地面风洞条件下的二元进气道为研究对象，采用高精度大涡模拟方法对锯齿形转捩片在三级压缩楔面上触发的边界层转捩现象开展了研究。数值方法基于隐式亚格子模型，空间离散采用高精度通量限制型紧致格式，时间推进采用显式Runge-Kutta方法。数值模拟清晰捕捉到了边界层转捩的空间发展演化过程，并获得了统计平均流场以及流场脉动特征。数值模拟结果表明转捩片能够有效触发进气道压缩面边界层转捩；通过与等熵压缩面及单楔面数值模拟结果的对比分析，获得了转捩片触发边界层转捩的内在机理，为后续研究工作奠定了基础。     

带单锥和双锥混压式进气道的冲压增程弹丸气动特性仿真分析

针对带单锥和双锥混压式进气道的两种冲压增程弹丸,对比分析了单、双锥混压式进气道的特点。采用仿真手段,对弹丸外部气动阻力和进气道工作特性进行了数值计算,分析了不同飞行马赫数条件下的弹丸外部气动阻力的变化、进气道临界状态下的总压恢复、捕获流量特性,并分析讨论了进气道的起动特性。研究结果表明:在飞行马赫数2.65到1.75之间,两种冲压增程弹丸的外部气动阻力几乎相同;在飞行马赫数1.75下,两种冲压增程弹丸进气道均能够起动;在低于设计马赫数时,进气道临界状态下,单锥进气道的捕获流量和总压恢复特性均优于双锥进气道,冲压增程弹丸采用单锥进气道将更为合适。     

ν-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。另外,基于该乘波前体和马赫数分布可控的内收缩进气道给出了一种双乘波的前体与进气道一体化设计方案,实现了内外流分别独立乘波,充分发挥了乘波前体和内收缩进气道的各自优势。