基于湍动能谱的自适应大涡模拟方法

为了在计算量有限的前提下实现对不同尺度上湍流燃烧结构的直接求解,提出了新的基于湍动能谱的自适应大涡模拟的概念。在这一概念下,利用湍动能谱和自适应网格技术建立了跨尺度算法,从而实现DNS和LES两种方法的耦合。在湍流非预混氢气射流火焰的模拟中,该方法既能捕捉大尺度湍流结构,又能分辨小尺度燃烧结构,计算出的物理量的瞬态和时均分布和DNS结果吻合得较好。计算偏差主要体现在DNS和LES区域的过渡处,这与亚格子湍动能在过渡处的变化不够平滑有关,表明跨尺度算法需要优化,以完善对不同尺度间湍动能传输的模拟。      

不同湍流预混燃烧模型在本生灯火焰中的比较

湍流预混燃烧的大涡模拟数值计算近来受到很大的关注,为了进一步发展湍流预混燃烧模型,利用正在发展的湍流预混燃烧模型等值面集合的G函数方法、层流预混火焰生成流型的化学热力学建表FGM方法、G函数/FGM耦合方法三种燃烧模型,对甲烷/空气湍流本生灯火焰进行了大涡模拟计算,并对上述燃烧模型进行了比较。大涡模拟结果表明,G函数方法能够较为准确捕捉火焰的传播特性,FGM方法能够捕捉到火焰和湍流的相互作用,但无法准确预测火焰的传播特性。而耦合的方法既可以捕捉火焰的传播特性,又能够捕捉到火焰和湍流的相互作用。     

输运概率密度函数中的小尺度标量混合建模

输运概率密度函数方法能精确地求解湍流燃烧中的有限化学反应速率,但小尺度上分子扩散引起的组分变化则需要通过标量小尺度混合模型来模拟。论文综述了当前湍流预混燃烧中混合模型研究的最新进展。首先介绍了现有标量混合模型在不同预混燃烧模式下的表现:对于火焰面燃烧模式,保持组分空间的临近性对于标量混合模型非常重要,标量混合频率受火焰结构影响显著;对于破碎反应区模式,组分空间保持临近性的重要性有所减弱,标量混合频率由湍流主导。接着介绍了最新基于直接数值模拟数据的标量小尺度混合机制和建模研究,包括湍流预混燃烧中反应标量小尺度耗散机制、反应标量混合时间尺度在不同燃烧模式下的演化规律、依据线性混合思路构建的标量混合时间尺度模型。针对不同湍流预混燃烧模式的模型验证表明:相比于已有湍流混合时间尺度模型,新模型显著提升了对湍流预混火焰中的标量耗散率和燃烧特性的预测精度。最后介绍了大涡模拟/输运概率密度函数方法在近极限湍流预混火焰,如甲烷预混值班火焰和湍流对冲预混火焰中的应用,验证了其对贫燃、高湍流度条件下的湍流-化学反应相互作用及局部熄火/再燃现象的预测能力。     

基于LES的非稳态火焰面/反应进度变量方法模拟部分预混抬举火焰

为了描述部分预混燃烧的火焰机制,捕捉燃烧现象中的非稳态效应,采用大涡模拟与非稳态火焰面/反应进度变量方法相结合来模拟湍流部分预混燃烧抬举火焰,并将计算结果与实验测量值进行了对比研究。结果表明基于大涡模拟的非稳态火焰面/反应进度变量方法捕捉到了湍流部分预混燃烧中的火焰抬举现象,计算得到的火焰抬举高度大致为35,和实验测量值完全一致。同时采用各热力学参数在混合分数空间的分布来研究部分预混抬燃烧中具有非稳态燃烧特性的局部熄火现象。各不同截面上的燃烧热力学参数与实验测量值吻合较好,说明采用该方法能比较真实地反映湍流部分预混燃烧的火焰机制。     

甲烷/空气预混火焰回火特性数值模拟

为了获得有效的预测回火极限的数值模拟方法,对甲烷/空气预混火焰进行了数值模拟。采用不同的反应机理和燃烧模型,针对一个二维轴对称的本生灯,分别对甲烷/空气预混层流火焰和湍流火焰进行了数值模拟,预测了甲烷/空气预混火焰的回火特性,并与Johnson的实验结果进行了对比。对层流火焰的模拟结果表明,采用Smooke-46反应机理能够较为准确地预测层流回火极限。而单步化学反应所预测的层流回火极限要比实验值低得多,这说明单步化学反应下的火焰更难以发生回火。对湍流火焰的模拟结果表明,涡耗散概念模型比有限速率模型更准确地预测了湍流回火极限。有限速率模型所预测的回火极限略高于实验值。     

氢气超音速燃烧非定常过程数值模拟

用有限差分法求解二维Ｎ－Ｓ方程和组分守恒方程，使用代数涡粘性湍流模型，两步化学反应模型和ＭａｃＣｏｒｍａｃｋ时间分裂法，模拟超音速气流中横喷氢气自动着火及火焰传播过程，清晰地描述了喷氢、火焰传播到主流、稳定燃烧过程中，温度、压力、速度和油气当量比的变化。计算结果表明，高温附面层和喷嘴前回流区能促使氢气自动着火，进口气流的静温是火焰向主流传播的重要影响因素。计算还证明了喷嘴后氢与空气两股流呈现平行流动，其间存在着强烈的湍流交换和化学反应边界层，超音速燃烧效率则主要取决于两股气流的混合速率。     

贫燃料湍流预混燃烧火焰稳定性的数值研究

对突扩燃烧室内甲烷-空气预混燃烧进行了数值模拟,时均控制方程组的封闭采用RNG 湍流输运模型和premixed combustion模型.通过模拟研究了不同当量比及进口速度对甲烷-空气预混燃烧效果的影响.CFD数值模拟结果给出了燃烧室内湍流预混反应流的速度场及温度场分布.模拟结果表明当量比为0.6、进口速度为30 m/s时燃烧室内产生稳定预混火焰.该结果对于航空燃气轮机低NOx排放燃烧技术的研究提供了一定的参考条件.     

反应机理对湍流预混值班火焰稳燃影响的研究

针对反应机理对湍流预混值班火焰稳燃的影响,使用大涡模拟(LES)耦合输运概率密度(PDF)方法模拟了湍流射流火焰的燃烧流场。化学反应过程采用被广泛应用的SMOOKE (16组分,46步)和DRM22 (22组分,104步)机理,概率密度输运的确定采用欧拉随机场来实现。计算结果表明,LES/PDF模型可以准确地模拟该类火焰的速度分布,火焰褶皱等结构。反应机理的点火延迟时间的准确与否是模拟值班射流火焰的关键,预混射流出口处的小尺度掺混引起的点火和火焰传播过程是值班射流火焰稳燃的关键。     

甲烷-空气贫燃料预混燃烧的数值模拟

针对燃气轮机燃烧室火焰筒中的甲烷-空气贫燃料预混燃烧问题,给出了基于八步化学反应动力学机理的数学模型。以某型航空发动机燃烧室火焰筒为例,对甲烷-空气贫燃料的预混燃烧进行数值模拟。研究结果表明,基于八步化学反应动力学机理的数值模拟方法,可以比较准确地反映燃烧产物的形成过程,在分析航空发动机燃烧室火焰筒内的贫燃料预混燃烧问题时具有较强的实用性。     

爆震波与非预混燃料射流相互作用的研究

为研究旋转爆震发动机（Rotating Detonation Engine, RDE）中燃料/氧化剂喷射和掺混对爆震波的影响及非预混环境下的爆震波的快速起爆与稳定传播,本文采用线性模型爆震发动机（Linear Model Detonation Engine,LMDE）来简化实际燃料喷射与爆震波相互作用的物理问题。通过RNG K-?湍流模型结合7步7组分氢气/空气机理的三维非定常反应流模拟方法,探究真实喷射条件下爆震波与混气相互作用、爆震波衰减及自持的特性。结果表明：氢气/空气的掺混均匀度至少要达到0.6才能使爆震波在非预混环境下传播;氢气孔与空气缝的入口压比需要满足爆震波进入燃烧室时,非均匀混气区恰好集中在氢气孔附近,燃料完全释放能量维持爆震波传播。     

ν-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_∞控制设计。仿真结果证明了该方法动态性能满足设计要求，同时对系统参数的不确定性具有较强的鲁棒性。     

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

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

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

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

Interfacial Properties Modification of Carbon Fiber/ Polyarylacetylene Composites

This work was dedicated to performing surface oxidation and coating treatments on carbon fibers (CF) and investigating the changes of fiber surface properties after these treatments, including surface composition, relative volume of functional groups, and sur- face topography with X-ray photoelectron spectroscopy (XPS) and atom force microscopy (AFM) technology. The results show that, after oxidation treatments, interfacial properties between CF and non-polar polyarylacetylene (PAA) resin are remarkably modified by removing weak surface layers and increasing fiber surface roughness. Coating treatment by high char phenolic resin solution after oxida- tion makes interface of CF/PAA composites to be upgraded and the interfacial properties further bettered.