基于火焰面模型的超声速湍流燃烧数值模拟研究

在非结构混合网格有限体积框架下,将RANS方程和层流火焰面模型相结合,开展超声速湍流燃烧流动数值模拟研究,时间和空间离散分别采用LU-SGS格式和HLLC格式。用SST k-ω模型模拟湍流,FlameMaster 3.9生成火焰面数据库,并采用β-PDF和δ-PDF分布来考察湍流和化学反应的相互作用。采用该数值方法,对DLR氢燃料超燃燃烧室和北航乙烯燃料气动斜坡超燃燃烧室开展了数值模拟研究,并将数值计算结果与实验数据进行对比,二者符合较好。研究结果表明:在主流中心支板顺流喷射超燃燃烧室内,燃烧主要发生在支板后的燃料/空气混合层内,并且火焰核心落在恰当混合物分数Zst附近区域。      

不同PDF模型在超燃冲压发动机数值模拟中的应用

湍流与燃烧相互作用对复杂超声速燃烧流场产生影响,基于设定型PDF(ProbabilityDensityFunction)模型对其不同分布模型进行研究。考虑温度和组分联合概率密度函数,采用β-PDF和δ-PDF方法对HF2直连式发动机模型进行模拟,并与该发动机风洞试验数据进行对比。结果表明:无论是定量还是定性分析,两种PDF方法的预测结果与实验结果吻合好,说明湍流与化学反应相互作用不能忽略;采用β分布的算例壁面压强分布结果与实验结果吻合更好,燃烧效率更高,更能够捕捉湍流对于化学反应的影响。     

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

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

超声速H_2/Air湍流扩散燃烧RANS数值模拟

为了研究修正的火焰面反应进度变量燃烧模型在超声速湍流扩散燃烧问题中的适用性,对德国宇航中心(DLR)超声速燃烧室开展RANS数值模拟。基于Open Foam软件平台中密度求解器分别对三维冷态场和燃烧场进行模拟分析。将网格自适应加密技术用于流场的计算;燃烧场计算中,通过分析不同压力下层流火焰面数据库,引入了反应进度变量源项的压力修正系数,压力修正系数α等于2.2。计算结果表明,冷态场中压力分布、波系分布、速度分布以及燃烧场中波系分布、速度分布、温度分布结果均与实验值符合较好。压力修正方法能够较好地解决超声速湍流扩散燃烧问题。湍流Schmidt数敏感性分析表明,湍流Schmidt数Sct对湍流火焰结构有较大影响,文中Sct等于0.7时能得到与实验值较为一致的分布。     

乙烯超声速燃烧部分预混火焰面模型数值研究

湍流涡结构的作用导致以扩散形式喷射燃料的超声速燃烧过程通常以部分预混的方式进行,为了使火焰面模型对这一燃烧过程描述的更准确,在模型中对扩散燃烧和预混燃烧状态均需要予以合理考虑。基于低Ma条件多区域火焰面模型(MRF)思想,发展得到了适用于超声速条件的部分预混火焰面模型。以碳氢燃料超声速双燃烧室结构作为验证算例,采用k-ωSST湍流模型、部分预混火焰面模型和乙烯(C2H4)的28组分72步化学反应机理对超声速湍流燃烧流场进行了数值研究。计算结果与实验数据吻合良好,验证了模型的准确性。加权系数计算结果表明,在中心剪切层和下壁面附近区域流场主要由预混燃烧控制,而其他大部分区域则由扩散燃烧控制。     

湍流扩散火焰局部熄火现象的大涡模拟研究

为了研究湍流扩散火焰的燃烧特性和局部熄火现象的机理,使用火焰面/反应进度模型对Sandia湍流值班射流火焰E进行大涡模拟。火焰面/反应进度模型适用于非预混燃烧,基于守恒标量和反应标量构建火焰面数据库,将详细的化学反应过程映射到这两个特征标量的平面上。数值模拟结果与实验数据符合得较好,火焰面/反应进度模型能较好地模拟湍流扩散火焰和预测火焰的不稳定性。研究发现,火焰在轴向距离与中心射流直径比值为7.5至15的截面处局部熄火概率较大,这与该位置标量耗散率偏大有关;火焰在轴向距离与中心射流直径比值为30以后的流场重新点燃,原因是混合物分数和当量比变小,火焰更多处于贫燃状态。     

值班甲烷/空气射流火焰的小火焰模型模拟

将小火焰模型和假定的PDF方法相结合,模拟甲烷湍流射流扩散火焰,稳态火焰面结构由层流小火焰数据库得到,采用详细化学反应机理描述甲烷的氧化和NO以及碳烟的生成,数值模拟结果和实验数据吻合较好,说明小火焰模型可以较准确地计算出扩散湍流火焰中的燃烧过程。      

采用不同建表方法的火焰面模型在燃烧室中的应用研究

为加深对航空发动机燃烧室中湍流燃烧过程的理解,采用不同建表方法的火焰面模型对航空发动机模型燃烧室内的湍流燃烧过程进行数值模拟,包括层流火焰面数据库的构建和反应进度变量的PDF类型两个方面。其中,层流火焰面数据库的构造方法包括基于扩散火焰的FPV和基于预混火焰的FGM模型,反应进度变量的PDF类型包括δ和β分布。LISA和KHRT模型分别用于模拟液膜和液滴的破碎过程,非平衡Langmuir-Knudsen模型用于模拟液滴的蒸发过程。LISA模型得到的液膜破碎距离约为4.6mm,液滴直径在文氏管出口下游迅速减小到10μm左右,并在头部出口下游附近完全蒸发。通过与相干反斯托克斯喇曼散射(CARS)和可调谐二极管激光吸收光谱(TDLAS)测量温度的对比,验证了FPV和FGM模型的精度,并表明在流动变化较大的位置FPV模型具有更高的精度,而其他位置FGM模型具有更高的精度,采用β分布作为反应进度变量PDF的模型,可以有效提高温度的预测进度,而且主燃区内的误差基本都在5%以内。此外采用β分布作为反应进度变量PDF的FGM模型,可以更好地描述未燃混合物被回流燃气点火的过程,而且反应进度变量的PDF类型比层流火焰面数据库构建方法的影响更为显著。     

建库温度在超声速火焰面模型中的影响

在火焰面模型中,数据库计算结果的准确程度对之后湍流燃烧流场的预测精度有重要影响。为了研究建库温度在超声速火焰面模型中的影响作用,在建立火焰面数据库时增加了建库氧化剂温度这一维度,并提出了一种由流场各点平均温度确定其对应的建库氧化剂温度的查表方法。针对DLR结构的超声速燃烧室,以煤油作为燃料,在五个温度级别条件下建立火焰面数据库,考察了建库温度对一维火焰面数据库以及二维燃烧流场计算结果的影响。研究结果表明,建库温度对一维火焰面的温度和组分分布有较为明显的影响,增加建库温度这一查表维度有助于改善超声速湍流燃烧流场的预测精度。     

采用不同建表方法的火焰面模型在燃烧室中的应用研究（英文）

为加深对航空发动机燃烧室中湍流燃烧过程的理解,采用不同建表方法的火焰面模型对航空发动机模型燃烧室内的湍流燃烧过程进行数值模拟,包括层流火焰面数据库的构建和反应进度变量的PDF类型两个方面。其中,层流火焰面数据库的构造方法包括基于扩散火焰的FPV和基于预混火焰的FGM模型,反应进度变量的PDF类型包括δ和β分布。LISA和KHRT模型分别用于模拟液膜和液滴的破碎过程,非平衡Langmuir-Knudsen模型用于模拟液滴的蒸发过程。LISA模型得到的液膜破碎距离约为4.6mm,液滴直径在文氏管出口下游迅速减小到10μm左右,并在头部出口下游附近完全蒸发。通过与相干反斯托克斯喇曼散射(CARS)和可调谐二极管激光吸收光谱(TDLAS)测量温度的对比,验证了FPV和FGM模型的精度,并表明在流动变化较大的位置FPV模型具有更高的精度,而其他位置FGM模型具有更高的精度,采用β分布作为反应进度变量PDF的模型,可以有效提高温度的预测进度,而且主燃区内的误差基本都在5%以内。此外采用β分布作为反应进度变量PDF的FGM模型,可以更好地描述未燃混合物被回流燃气点火的过程,而且反应进度变量的PDF类型比层流火焰面数据库构建方法的影响更为显著。     

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