同心旋流分层火焰的外激脉动特性统计学分析

为了更深入地理解民用航空发动机低排放燃烧室中的同心旋流分层火焰脉动机制，以丙烷为燃料，在100Hz频率的外激条件下对该类旋流火焰的CH*化学发光图像进行高频采样，重点通过本征正交分解法（POD）的统计学分析，提取其时间和空间域内的脉动特征。分析结果表明，旋流火焰中4个最主要的脉动模态分别是轴向振荡、径向振荡、火焰脱落和非对称螺旋运动，并且模态1轴向振荡和模态2径向振荡之间具有显著相关性。POD分析方法可作为中心分级贫油预混燃烧不稳定性机理探索的有效手段。      

同心旋流分层预混火焰的动力学模态分析

为了探究同心旋流分层火焰的不稳定燃烧问题,需要对火焰的动态脉动特性进行分析。采用高速摄像捕捉了自激和外激条件下同心旋流分层预混火焰的CH*化学发光动态图像,利用动力学模态分解(DMD)方法研究其主导脉动模态,提取出了相关模态的空间形态和脉动幅值,证明燃烧组织方式的改变会对火焰脉动的形态和规律产生影响。研究发现:模态分解可以反映分层旋流火焰的动态响应特征,并能解析脉动空间结构,验证火焰形态和外激频率对火焰动态响应的影响,其规律和火焰传递函数结果符合,且在分层比0.5～2内提高分层比使释热中心转移到预燃级和主燃级的剪切层内,300Hz～400Hz释热脉动将加强。验证了DMD方法适用于处理分析复杂分层旋流火焰的脉动特征,可为燃烧振荡控制提供参考。     

贫燃预混旋流火焰的模态转换燃烧不稳定特性分析

为进一步研究燃烧室燃烧不稳定特性,通过实验对贫燃预混旋流火焰的燃烧不稳定性模态转换现象进行了研究。在贫燃的情况下进行了甲烷-空气预混旋流燃烧实验,用动态压力传感器测量火焰场动态压力,用光电倍增管捕捉CH~*化学发光,用4kHz频率的高速摄像记录火焰图像。发现当量比渐增到0.72附近时,燃烧有明显的模态转换,压力脉动主频由264Hz突变为187Hz,脉动幅值大幅增加。采用本征正交分解法得到各模式空间分布,通过模态系数进行快速傅里叶变换(FFT)分析,得到流场时间分布特征。发现模态转换前平均火焰为喇叭形,转换后为"M"型,模态的转换与涡脱落模式有重要关系。     

当量比对分层旋流火焰燃烧不稳定性的影响

针对分层旋流燃烧室,实验测量了不同分级比和油气比下的动态压力。结果显示:在总燃烧当量比0.47~0.68,主燃级当量比0.48~0.65范围内,燃烧室内的压力脉动幅值随当量比增大而增大。基于层流火焰模型和局部火焰吹熄假设,公式推导了吹熄前后火焰厚度、火焰长度、层流火焰传播速度、气流速度间的关系,得出燃烧不稳定下的压力脉动幅值ps′与层流火焰传播速度SL0成正比,结合贫油状态下SL0随当量比增大而增大这一事实,得出ps′随当量比增大而增大,从而解释了实验结果。     

基于PaSR模型的低旋流燃烧大涡模拟

为了深入理解低旋流流场特征和燃烧稳定性,基于OpenFOAM平台,采用动态k方程模型和有限速率PaSR模型对甲烷/空气预混低旋流燃烧进行了大涡模拟,研究了气流入口速度、当量比和压力等流场参数对流场结构和燃烧非稳态特性的影响,分析了流场大尺度结构与火焰相互作用。结果表明,流场结构和火焰抬升高度受入口速度影响较小,流场和火焰形态能够保持自相似性;随着当量比和压力提高,流场扩张性增强并在燃烧区下游产生回流区,火焰稳定不依赖回流区,根部火焰锋面形状由U形转变为W形,火焰抬升高度降低。火焰锋面稳定在剪切层,剪切层产生的周期性有序涡结构引起当地流场速度脉动和火焰表面褶皱,反映了流场非稳态特性;通过剪切层监测点瞬时轴向速度分析,涡结构特征频率随速度增大而提高,由250Hz提高至300Hz,随当量比和压力提高而降低,由250Hz降低至125Hz。     

模型燃烧室中值班火焰稳定器的燃烧不稳定性研究

针对值班火焰的燃烧不稳定性问题，通过试验得到了不同当量比工况下的火焰结构、压力脉动信号和热释放脉动信号。对脉动信号进行傅立叶分析，获得了脉动信号的频率和振幅。结果表明，随着当量比增加，火焰形态发生变化，燃烧室内发生了164Hz三阶模态向109Hz二阶模态的转换。对火焰平均结构、火焰瞬态结构和火焰正交分解（POD）结果进行分析后发现，漩涡脱落频率和燃烧室某一阶声学模态耦合是维持燃烧不稳定性的机理，耦合的强弱和漩涡脱落的尺度大小决定了热声振荡的振幅大小，火焰形态的变化导致热释放中心位置的变化是引起模态转换的机理。     

不同气量分配下声激励对中心分级旋流火焰动力学的影响

采用高速摄像拍摄火焰化学自发光和本征正交分解（POD）法研究了多激励频率下中心分级旋流分层火焰的宏观结构和动力学特性。结果表明：随着气量分配比的增加，火焰宏观结构逐渐沿径向方向扩张，火焰分层现象变得更加明显。激励频率的变化对火焰宏观结构的影响较小。POD分析结果表明：外激频率为200 Hz时火焰脉动最明显。随着主燃级气量占比的增加，释热脉动在激励的作用下沿轴向方向减弱，在径向上出现极值交替且振荡变强。     

雾化气流脉动下乙醇喷雾旋流火焰特性

针对雾化气流脉动条件下乙醇喷雾旋流非稳态燃烧问题,开展了不同雾化流量下乙醇喷雾旋流火焰对雾化气流脉动响应特性的实验研究,揭示了不同脉动频率和脉宽时长下火焰长度和抬升高度的响应特性,获得了乙醇喷雾燃烧在雾化气流脉动条件下的火焰结构、温度分布等动态规律.结果表明:随着雾化流量增加,旋流火焰由附着演化为抬升;此外,雾化气流脉...     

OH与CH2O双组分平面激光诱导荧光对旋流燃烧室火焰结构与脉动特征的研究

介绍了利用平面激光诱导荧光(PLIF)技术对航空发动机的旋流燃烧室模型在贫燃状态下工作特性的研究。通过对OH与CH2O双组分进行同步PLIF测量,获得了不同工况下燃烧室反应区以及预热区的瞬态结构信息。应用本征正交分解(POD)方法对OH PLIF的图像进行处理,得到了旋流火焰的主要脉动模态,并通过扩展本征正交分解(EPOD)方法计算出了相应POD模态的CH2O荧光信号分布。实验结果表明:随着燃烧室热功率的增大,火焰的整体结构、脉动模式均出现了明显的变化。在火焰高度增加的同时,轴向不稳定性逐渐增强,涡核旋进(PVC)的脉动特征相对减弱。在较大的热功率下,在燃烧室的外回流区(ERZ)出现未燃烧的燃料。      

贫燃预混旋流火焰PIV和OH-PLIF瞬态同步测量

为了研究贫燃预混旋流火焰流场与火焰结构特性,设计了强受限预混旋流燃烧器,搭建了PIV(粒子图像测速)和OH-PLIF(平面激光诱导荧光)瞬态同步测量系统。针对典型贫燃工况开展了瞬态同步测量研究,获得了多个截面的瞬态和时均流场速度和OH分布信息。研究表明:低当量比不稳定燃烧工况下,伴随着角涡回火现象,燃烧室底部角落存在明显火焰反应区;由于旋流拉扯效应导致OH分布结构边缘产生褶皱,并出现多个OH“孤岛”;旋流火焰外沿位于喷嘴射流的内侧剪切层;增加当量比,火焰流场结构与反应区分布与低当量比类似,但燃烧更稳定,火焰高度有所增加;内部回流区的经典双涡结构分裂为多个旋涡,存在旋涡进动、破碎的现象。     

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