氢燃料双模态冲压发动机火焰结构及其稳定机制的LES研究

为深入理解双模态冲压发动机内流场结构和燃烧特性，本文采用精细的有限速率化学反应模型，对凹腔内氢气直喷式超声速燃烧室火焰进行了大涡模拟研究。发动机隔离段入口马赫数为2.0，滞止温度和压力分别为950 K和0.82 MPa。定性和定量的验证分析表明，计算结果良好符合试验所反映的物理规律，再现了两种典型的工作模态及其稳焰模式。当量比为0.1时，发动机处于超燃模态，为凹腔剪切层稳焰模式；当量比为0.3时，发动机处于亚燃模态，为凹腔辅助射流尾迹稳焰模式，分离涡的大尺度脉动及凹腔回流区的缺失致使火焰剧烈振荡。同时采用改进的火焰因子和过滤函数详细分析了局部火焰特征和流动模式，观察到了不同规律的局部熄火现象，并且剧烈的流动振荡对于局部火焰结构的稳定性有着不利影响。     

正癸烷在模型超燃燃烧室中的燃烧特性研究

为探究正癸烷在双凹腔结构的模型超燃燃烧室中的燃烧特性,试验采用蓄热式加热器提供高焓纯净空气,燃烧室进口的Ma数为2.03,来流总温在800K~1100K,常温液态的正癸烷经凹腔上游的直射式燃料喷嘴进入燃烧室。通过对试验过程中燃烧室壁面压力和流场中的CH*基分布的分析,发现了正癸烷的两种稳焰模式:双凹腔稳焰模式和单凹腔稳焰模式。双凹腔稳焰模式是通过射流形成的尾迹区和凹腔中的回流区共同作用实现稳焰;单凹腔稳焰模式则是通过主流中的激波与边界层干涉形成分离区与凹腔中的回流区实现稳焰。随着试验来流总温的降低,正癸烷的稳焰模式从双凹腔稳焰转变为单凹腔稳焰,直到稳焰失败,模式转变温度和临界稳焰温度分别为876K和842K。还利用一维分析方法对两种稳焰模式的燃烧效率、Ma数分布以及总压恢复系数进行了比较,结果发现双凹腔稳焰模式的燃烧效率和总压恢复系数均大于单凹腔稳焰模式。     

基于500 Hz OH-PLIF技术的超声速燃烧室火焰结构

针对超声速燃烧室火焰驻留时间短和高湍流度的特点,利用500 Hz高频OH-PLIF(planar laser induced fluorescence)技术,研究了氢燃料超声速燃烧室火焰结构,结合壁面压力,获得了点火、稳焰以及熄火过程的规律。结果表明：隔离段入口来流马赫数为2的情况下,氢气当量比0.1和0.3,在火花塞点火下,2 ms均能点燃并获得PLIF火焰图像,其中当量比为0.1凹槽内存在爆燃现象,火焰达到稳定燃烧状态较慢(约65 ms)。稳焰过程中,当量比为0.3燃烧相对不充分,燃烧区域更靠近下壁面,燃烧位置和火核位置变化较大。当量比为0.1和0.3,熄火时间均在6 ms内,其中当量比为0.1熄火时火焰是从凹腔中部传播到凹腔前缘位置熄灭并伴随有一个短暂火焰增强的过程,并且在传播至凹腔前缘时已接近熄灭,当量比为0.3在熄灭前则是慢慢变弱,最终在凹腔前缘至喷氢位置间熄灭。     

可变几何单凹腔驻涡燃烧室的贫油熄火性能试验

为了研究可变几何单凹腔驻涡燃烧室的贫油熄火性能,设计了一个可变几何单凹腔驻涡燃烧室,通过改变凹腔前进气槽的开度,实现对凹腔当量比的控制,并在此基础上进行了贫油熄火性能试验。试验在常压状态下进行,采用航空煤油作为燃料。试验中的主要研究参数如下:进口温度在287~487K变化;进口马赫数在0.2~0.3变化。研究结果表明:可变几何单凹腔驻涡燃烧室的贫油熄火油气比随着进口温度的增加而减小,随着进口马赫数的增加而增加;保持进口马赫数为0.25不变,当进口温度为287K时,贫油熄火油气比随着凹腔前进气开度的减小先增加后减小,而当进口温度逐渐升高到487K时,这一趋势逐渐转变为单调减小;保持进口温度为487K不变,贫油熄火油气比随着凹腔前进气开度的减小而减小。     

带凹腔火焰稳定器的固体火箭超燃冲压发动机燃烧室试验研究

以含硼贫氧固体推进剂为燃料,对带凹腔火焰稳定器的固体火箭超燃冲压发动机燃烧室构型首次开展了地面直连试验研究。试验模拟了23km,5.5Ma的飞行工况,通过测量压强、推力和流量等参数,得出了燃烧室性能。试验结果表明：一次富燃燃气在燃气发生器喉部沉积导致燃气流量持续提高,试验过程中当量比由0.44逐渐增加至0.54;本文所研究的凹腔稳焰结构提高了富燃燃气中气相可燃组分的燃烧效率,但对于凝相颗粒燃烧的促进效果不明显;试验工况下发动机总燃烧效率约48%,高空比冲约为423s,高于文献中所报道的中心支板喷射稳焰的固体火箭超燃冲压发动机试验比冲性能。     

中心分级低污染燃烧室贫油熄火模型

为了准确掌握中心分级燃烧室火焰稳定边界的影响因素,建立中心分级燃烧室贫油熄火边界预测模型,对中心分级燃烧室的熄火过程进行了试验与理论研究。研究获取了燃烧室结构、雾化和工况参数对燃烧室贫油熄火边界的影响规律,建立并验证了中心分级燃烧室熄火半经验预测模型。结果表明：相比反向涡流器,同向涡流器下游具有更大的回流区、更低的回流速度和更长的停留时间,从而减弱了主燃级与值班级之间的湍流交换,导致同向旋流火焰的贫油熄火性能明显优于反向旋流火焰的熄火性能。中心分级燃烧室熄火边界预测模型对单头部和全环燃烧室熄火性能预测的最大误差为20%,满足燃烧室工程设计需要。     

导流片结构参数对先进旋涡燃烧室性能影响研究

为了研究带导流片的先进旋涡燃烧室燃烧及流动性能,对导流片不同结构参数下燃烧室的流场分布、温度分布、总压损失、燃烧效率及污染物排放进行了数值模拟。结果表明,带导流片的先进旋涡燃烧室性能表现出明显的优越性,能够在不依赖凹腔射流的情况下形成理想的双旋涡结构,有利于增强凹腔稳焰及燃气掺混,大幅度提高燃烧效率,改善出口温度分布,降低NO排放。导流片结构参数对燃烧室性能影响呈现出一定的规律。当导流片伸入凹腔的长度与前钝体高度之比(a/B)为0.2,导流片到前钝体上端面的距离与燃烧室进气通道高度之比(b/H)为0.4,导流片到前钝体后端面的距离与凹腔长度之比(c/L)为0.1~0.2时,燃烧室在贫燃条件下可以实现低压降、低污染的稳定燃烧,而且出口温度分布均匀。     

低马赫数下多凹腔燃烧室非稳态燃烧过程

作为稳定火焰的有效手段之一,凹腔构型在冲压发动机燃烧室研究中占有重要地位。在对以煤油为燃料的多凹腔燃烧室冷/热态流动特性分析的基础上,重点研究低进口马赫数条件下燃烧室点火起动初期非稳态过程。结果表明：上游凹腔内大涡结构有助于提高燃料的驻留时间,未燃混气被高速主流带入下游凹腔内继续反应,进一步提高燃烧效率;燃油喷射速度决定被卷吸进回流区的燃油质量分数的大小,进而影响燃烧效率高低;燃烧室点火起动初期出现了主流熄火、火焰逆流传播以及主流再着火等复杂现象,火焰逆流传播现象是在上游凹腔内燃料自燃与下游燃烧释热压缩来流两种机制共同作用下完成的。     

喷注方案对凹腔稳焰燃烧室低频振荡的影响研究

为了研究在入口来流马赫数2.52,总温1486K的超声速来流条件下,稳焰凹腔上游不同位置乙烯横向喷注对模型发动机燃烧室内低频燃烧振荡特性的影响,通过1kg/s直连式超燃试验平台,利用高频压力传感器、高速摄影相机等设备,对凹腔上游近距离、远距离喷注等方案的发动机内部压力与火焰动态特性进行了研究。试验结果表明:在当前当量比条件下,当稳焰凹腔上游近距离喷注燃料时,燃烧室存在较大范围亚声速区域,并出现由热声不稳定性激励的低频压力振荡,频率分布范围较宽(50～400Hz)且振幅较弱。对于燃料喷注位置到稳焰凹腔距离较远的情况,燃烧室内出现以火焰逆传和火焰吹脱为特征的周期性火焰振荡现象。分析认为较远喷注距离有利于燃料-空气充分混合并形成预混区,导致火焰快速逆传。火焰逆传与DDT (爆燃转爆震)中的火焰加速传播过程有关。周期性火焰逆传与火焰吹脱过程相耦合形成了具有特定主频(约85Hz)且振幅较大的低频压力振荡。     

轴向双级反旋预混喷嘴燃烧调整边界实验研究

干式低污染燃烧室为保证安全、稳定、高效的运行状态，通常需要进行燃烧调整。本文针对一种干式低污染燃烧室用轴向双级反旋预混喷嘴，在常温常压条件下开展燃烧调整边界实验研究。测量并探究了影响燃烧调整边界的因素，如贫熄火边界、污染物排放边界以及动态压力脉动边界，着重分析了喷嘴内的剪切层稳焰机理以及火焰根部跳动与壁面相干作用导致的热声振荡机理。研究结果表明，由NOx、CO、火焰状态、动态压力脉动和熄火边界共同决定的稳定运行边界在低负荷条件下为当量比0.576～0.714，随着热负荷升高，运行边界有所收窄。对比分析了纯甲烷以及掺50%H2两种燃料的燃烧调整边界，燃料加氢后边界整体向低当量比移动，运行窗口宽度得到拓展。     

ν-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.     

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

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

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.      

二维翼段颤振的μ控制

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

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

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

Abnormal Shape Mould Winding

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

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.     

Thermal Stresses and Cracks During the Growth of Large-sized Sapphire with SAPMAC Method

本文利用有限元法对冷心放肩微量提拉法(SAPMAC)生长大尺寸蓝宝石单晶过程中的热应力进行了模拟计算,应用临界缺陷理论解释了裂纹的萌生与扩展机理。研究表明,晶体内的热应力主要与晶体生长速率,环境温度和结晶取向有关;较大热应力多出现在结晶界面,放肩、收尾以及直径急剧变化等位置;最大热应力总是出现在籽晶与新生晶体的界面附近。裂纹将在临界缺陷位置产生,并在应力作用下沿a或m面扩展。计算结果与实验结果基本吻合,通过晶体生长系统和生长工艺的改进,对晶体的开裂问题得到了有效的抑制。     

Optimization Design System for Composite Structures Based on Grid Technology

To solve the topology optimization of complicated multi-objective continuous/discrete design variables in aircmit structure design, a Parallel Pareto Genetic Algorithm （PPGA） is presented based on grid platform in this paper. In the algorithm, the commercial finite element analysis （FEA） software is integrated as the calculating tool for analyzing the objective functions and the filter of Pareto solution set based on weight information is introduced to deal with the relationships among all objectives. Grid technology is utilized in PPGA to realize the distributed computations and the user interface is developed to realize the job submission and job management locally/remotely. Taking the aero-elastic tailoring of a composite wing for optimization as an example, a set of Pareto solutions are obtained for the decision-maker. The numerical results show that the aileron reversal problem can be solved by adding the limited skin weight in this system. The algorithm can be used to solve complicated topology optimization for composite structures in engineering and the computation efficiency can be improved greatly by using the grid platform that aggregates numerous idle resources.