有蒸发稀液雾冷态射流的大涡模拟研究

为研究气相湍流场中液滴弥散和蒸发过程的特点和影响因素,加深对气液耦合作用的认知,用欧拉-拉格朗日相结合的办法数值模拟了稀液雾冷态受限射流和自由射流。欧拉坐标系下的气相湍流场使用大涡模拟技术模拟,离散液相则使用拉格朗日颗粒轨道模型进行描述,考虑相间质量、动量、能量交换,采用相平衡蒸发模型模拟液相的蒸发效应。模拟得到的液相速度场、粒径分布和质量流率等均和实验数据较好吻合。有回流的受限射流中初始惯性对Stokes数大的液滴弥散影响较大,而Stokes数在1附近的小液滴弥散更受气相流场结构影响。液雾单位体积蒸发速率和液滴弥散分布密切相关。自由射流中液滴受初始惯性和受气相流场作用大体一致,大小液滴的弥散分布基本满足自相似性。     

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

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

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

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

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

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

环形燃烧室两相喷雾燃烧的大涡模拟

在三维任意曲线坐标系下采用欧拉-拉格朗日两相大涡模拟方法对环形燃烧室火焰筒气液两相紊流瞬态反应流进行数值模拟.采用椭圆偏微分方程生成三维贴体网格,计算中所采用的数学度模型有:k方程亚网格尺度模型估算亚网格紊流粘性;亚网格EBU燃烧模型估算化学反应速率;热通量辐射模型估算辐射换热.并在非交错网格体系下气相采用SIMPLE算法和混合差分格式求解,液相采用随机离散模型(Stachasttc Separated Flow,简称SSF),在拉格朗日坐标系下追踪各油珠群沿各自轨道运动、质量损失及能量变化.通过计算结果与实验数据相比较,表明在三维贴体坐标系下对燃烧室火焰筒两相紊流油雾燃烧流场进行大涡模拟,采用欧拉-拉格朗日两相大涡模拟方法能反映两相紊流化学反应流流动及实际燃烧过程.     

湍流两相燃烧不稳定性分析

采用数值模拟的方法对液体火箭发动机燃烧室内湍流两相流动与燃烧过程进行了分析,同时研究了燃烧室内高频不稳定燃烧现象。气相控制方程用欧拉坐标系下的Navier-Stokes方程组描述,液相控制方程在Lagrangian坐标系下描述,湍流模型采用高雷诺数的k-ε双方程模型,并采用高压蒸发燃烧模型。通过数值模拟分析了当燃烧室处于稳态燃烧时从喷注面到喉部之间两相混合物的变化规律;当发生高频不稳定燃烧时所研究的控制体内混合气体和喷雾的量发生振荡,并且具有相同的频率和半个周期的相位差。     

湍流贫燃预混射流火焰的大涡模拟

为了研究高Karlovitz(Ka)数下湍流预混火焰的驻定机制和火焰结构,采用化学热力学建表方法结合假定概率密度函数(PDF)模型,对值班预混射流火焰(PPJB)进行大涡模拟(LES)研究。选取自点火模型作为建表模型,使用假定PDF模型考虑湍流和化学反应之间的相互作用,其中假定双混合物分数的概率密度分布为Dirichlet分布,反应进度变量的概率密度分布为β分布。模拟结果表明,能较好预测PPJB火焰的温度和主要组分的分布。使用中间组分(CO和OH)的质量分数来定性描述PPJB火焰结构,结果表明高Ka数下火焰结构受到湍流涡旋的显著影响。     

大涡模拟模型燃烧室燃烧性能计算

对带双级扩压器的模型燃烧室气液两相瞬态喷雾燃烧过程,在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法进行数值研究,同时采用多维经验分析法预估燃烧性能.采用 k 方程亚网格尺度模型模拟亚网格湍流黏性;亚网格EBU(eddy-break-up)燃烧模型预估化学反应速率;多维经验分析法计算燃烧性能;并在非交错网格体系下气相采用SIMPLE(semi-implicit method for pressure-linked equations)算法对控制方程进行求解,液相采用随机离散模型,两相之间的耦合采用PSIC(particle-source-in-cell)算法.通过大涡模拟瞬态及时均计算结果表明:与粒子图像测速仪(PIV)测量的瞬态速度场、出口温度分布试验数据吻合,表明在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法,数值模拟模型燃烧室两相喷雾燃烧流场,所采用的亚网格模型可以用于燃烧室气液两相喷雾燃烧流场的大涡模拟;燃烧性能计算结果与试验测量结果基本一致,说明所采用多维经验分析法可以用来数值模拟航空发动机燃烧室燃烧性能的计算,特别是污染物的预估,为设计低污染高性能航空发动机燃烧室提供有用的设计依据.      

液体箭发动机燃烧不稳定性分析模型

通过一个简化的物理模型。运用计算流体动力学分析方法，研究了液体火简发动机的不稳定燃烧现象，验证了该模型及分析方法的简捷和有效性，描述发动机燃烧室内气相流动过气液两控制方程是非稳准一维欧拉方程组，液相则用拉格朗日方法描述，方程组中的源项反映了气液两相耦合及液雾蒸发燃烧的影响。针对不稳定燃烧的机理，考察了压力脉动随着时间的推进对声学扰动的敏感性，计算格式采用两步积分的预测校正法进行迭代求解。     

超声速气流中液体喷雾流动数值模拟研究

超声速气流中液体喷雾流动特性对超声速燃烧基础研究及其工程应用具有重要意义。为了定量探索超声速气流中液体横向射流雾化特性，本文对超声速气流中液体喷雾流动进行了数值模拟研究。数值模拟方法基于Eulerian-Lagrangian两相流计算架构，考虑气液双向耦合，采用KH/RT液滴二次破碎模型计算液滴雾化过程，采用大涡模拟计算气相流动。结果表明，该数值模拟方法所获得的液雾场突起结构、穿透高度、液滴平均速度分布等液雾特性均能与试验结果较符合；初始液滴直径分布对破碎后液滴平均速度影响较小而对破碎后液滴平均直径及液相平均体积分数影响较大，初始液滴直径分布需在后续的建模与模拟中进行更多研究。     

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

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

二维翼段颤振的μ控制

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

Abnormal Shape Mould Winding

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

基于H∞性能指标的质量矩拦截弹鲁棒控制

 以所建立的质量矩拦截弹数学模型为基础，利用微分几何的反馈线性化理论，得到一个解耦线性系统，考虑到拦截弹的鲁棒性要求和3个滑块的协调控制问题，提出采用双回路的设计方法，内回路采用线性二次调节器（LQR），外回路采用考虑混合灵敏度问题的H_∞控制设计。仿真结果证明了该方法动态性能满足设计要求，同时对系统参数的不确定性具有较强的鲁棒性。     

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

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

跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。     

基于马赫数分布可控曲面外/内锥形基准流场的前体/进气道一体化设计

提出了一种高超声速飞行器乘波前体的外锥形基准流场设计方法,在锥面马赫数分布规律给定的条件下,通过有旋特征线法实现反设计,提高了基准流场设计的灵活性。该基准流场通过锥形"下凹"弯曲激波和波后等熵压缩波系压缩气流,可以在较短的长度内完成高效压缩。基于反正切马赫数分布外锥形基准流场设计的乘波前体具有较高的容积率,乘波特性良好且出口均匀,设计点时有黏升阻比为1.89。另外,基于该乘波前体和马赫数分布可控的内收缩进气道给出了一种双乘波的前体与进气道一体化设计方案,实现了内外流分别独立乘波,充分发挥了乘波前体和内收缩进气道的各自优势。     

Potentials of manufacture and repair of nickel base turbine components used in aero engines and power plants by laser metal deposition and laser drilling

IntroductionExpensive turbine parts like HPT(HighPressure Turine)blades or vanes are replaced bynew parts in case of damage.For example theburn through of the inner side of a blade or vane(Figure 1)is a frequently appearing damage,which cannot be repaired…     

Enhancing the efectiveness of film cooling

Advanced gas turbine stages are designed to operate at increasingly higher inlet temperatures to increase thermal efficiency and specific power output.To maintain durability and reasonable life,film cooling is needed in addition to internal cooling,especially for the first stage.Film cooling lowers material temperature by forced convection inside film-cooling holes and by forming a layer of coolant about component surfaces to insulate them from the hot gases.Unfortunately,each cooling jet forms a pair of counter-rotating vortices that entrains hot gas and causes the film-cooling jet to lift off from the surface that it is intended to protect.This paper gives an overview of efforts to enhance the effectiveness of film-cooling.This paper also describes two new design concepts.One design concept seeks to minimize the entrainment of hot gases underneath of film-cooling jets by using flow-aligned blockers.The other design concept shifts the interaction between the approaching hot gas and the cooling jet to occur further above the surface by using an upstream ramp.For both design concepts,computational fluid dynamics results are presented to examine their usefulness in enhancing film-cooling effectiveness.      

Guide for Authors

<正>About Journal Chinese Journal of Aeronautics(CJA)is a comprehensive academic journal dealing with the fields of aeronautics and astronautics.It reports researches concerning the two fields in China and abroad to promote the academic exchange.Founded in 1988 and sponsored by the Chinese Society of Aeronautics and Astronautics and Beihang University,CJA publishes papers bimonthly,with issues released in February,April,June,August,October and December.