高超声速喷管三维热化学非平衡流场的数值模拟

利用数值求解三维热化学非平衡Navier-Stokes方程的方法和多块网格技术,数值模拟了典型高焓激波风洞锥型喷管的三维非平衡流场,分析了热化学非平衡效应对喷管流动的影响,给出了平动温度、组分N2和O2振动温度、马赫数和组分质量分数在轴对称线上的分布,也给出了喷管出口附近的速度和温度沿径向分布等结果,获得了喷管流场的详细信息,并与已经存在的采用轴对称Navier-Stokes方程计算的结果,进行了比较,二者吻合得很好.计算结果表明:目前的代码是可以模拟多块结构网格下的热化学非平衡流动的.     

高焓激波风洞自由流参数测量的数值重建

高焓风洞中参数测量本身受到高温非平衡效应影响,单纯依靠实验测量不足以确定详细的风洞自由流参数.本文以JF-10高焓激波风洞为背景,通过喷管-试验段-测量仪模型非平衡流场的连贯计算,开展自由流参数测量的数值重建.在比较分析测量仪模型绕流场的计算结果和实验测量值的基础上,结合喷管流场的数值模拟结果,共同确定高焓风洞的详细自由流参数.     

高超声速喷管非平衡黏性流动的数值研究

分别采用平衡气体模型、化学非平衡气体模型和热化学非平衡气体模型,通过求解轴对称Navier-Stokes方程组,数值模拟了法国Marseille高焓激波风洞锥型喷管(H₀=10.3MJ/kg),分析了热化学非平衡效应对喷管流动的影响,给出了平动温度、振动温度、马赫数和组分质量数在轴对称线上的分布与喷管出口附近的速度和温度沿径向分布等结果。计算结果表明:化学反应速度加快,会导致喷管出口温度增加,振动能的冻结会导致喷管出口温度降低。     

高焓激波风洞喷管流场非平衡特性研究

高焓激波风洞是开展高超声速流动研究的重要地面模拟设备,但其产生的高焓气流在喷管中的膨胀过程是一种典型的热化学非平衡流动,试验段特征参数通过直接实验测量难以完全确定。本文通过求解耦合双温度模型的轴对称Navier-Stokes方程,研究了高焓激波风洞中典型状态下气流的热化学非平衡流动特性,分析了焓值对非平衡特性的影响规律。结果表明,喷管出口自由流均匀区域达到出口截面直径的75%以上,能够为实验提供足够的空间;喷管出口自由流处于热化学非平衡状态,在喷管喉道后约1/5喷管长度处气流即已处于冻结流状态,组分浓度和振动温度随气流流动基本不变;焓值在8.4MJ/kg~19.5MJ/kg之间变化时,非平衡程度随着焓值的增加而增强,但是低焓值时非平衡程度的增强更加剧烈。     

高焓激波风洞有效试验时间的诊断

高焓激波风洞能够产生模拟高马赫数飞行条件的气流总温,是研究高温真实气体效应以及再入物理问题的有效试验装备,但是激波风洞的试验时间较短,且随着气流焓值的提高大幅降低,仅为几毫秒,因此试验测试数据曲线中有效时间段的分辨十分重要,它直接影响到试验结果的可靠性及精度。鉴于此,采用压力测量、静电探针测量、非接触光学测量和热流测量的方式,针对中国科学院力学研究所JF-10高焓激波风洞16 MJ/kg总焓、7700 K总温的流场状态,对比研究了风洞喷管的起动时间以及有效测试时间。试验结果表明:静电探针测量方法最为有效地分辨了喷管起动时间段、有效试验时间段以及驱动气体的到达; JF-10高焓风洞在16 MJ/kg的状态下,喷管起动时间约为1.3 ms,风洞有效试验时间约为2 ms。     

不同热化学非平衡模型对高超声速喷管流场影响的数值分析

采用轴对称NS方程,数值研究不同热化学非平衡模型对高超声速喷管流场的影响,包括:(1)不同组元数的高温空气模型(5组元、7组元1、1组元)的比较;(2)热力非平衡(双温度)的化学动力过程与热力平衡(单温度)的化学动力过程的比较。计算结果表明,高焓风洞实验条件下喷管流场处于热力和化学都是非平衡的状态。在计算条件下,数值模拟以采用7组元或11组元的热化学非平衡模型为宜。     

再入体实验模型热化学非平衡绕流流场的数值分析

利用数值求解一个和多个振动温度热化学非平衡Navier-Stokes方程的CFD计算程序,对爆轰风洞球锥试验模型的热化学非平衡绕流流场进行了数值模拟,分析了分子组分振动温度在全流场(头身部流场和底部流场)中的分布规律.研究结果表明:(1)计算的压力、电子数密度以及流场的光辐射数据与试验数据符合较好;(2)在再入体身部,大多数分子的振动温度峰值大大高于平动温度的峰值;在再入体近尾出现振动温度冻结并大大高于平动温度的现象;CO2的两个振动态的振动温度分布非常接近平动温度分布.     

防热试验用高温超声速燃气流场热环境分析

对用于防热系统试验考核的高温超声速燃气流场的流场结构和热环境特性进行了计算和分析.将CFD模拟方法与工程计算方法相结合,对两类锥形喷管形成的燃气流场进行了数值模拟计算和流场结构分析,并与试验过程中高速拍摄的图片进行了比对.利用模拟得到的流场气动参数,计算了沿流场轴向位置球头驻点模型的热流和压力,并与试验测量数据进行了对比,结果吻合较好.锥形喷管形成压缩波和膨胀波交替的流场结构,导致流场气动参数与模型热流和压力形成正弦变化的特点.     

超高速高焓流动研究进展

本文介绍了中国科学院高温气体动力学重点实验室在超高速高焓流动模拟技术和试验方法方面取得的研究进展.文章主要包括三部分研究内容:第一部分是关于发展先进的超高速试验模拟技术,包括爆轰驱动高焓激波风洞和爆轰驱动高焓膨胀管.高焓激波风洞产生的超高速气流速度的范围是3.5km/s～6.0km/s,高焓膨胀管能够模拟速度为6.5km/s～10km/s的超高速气流.第二部分介绍高焓激波风洞喷管流场诊断结果,用来检验喷管产生的超高速流场的流场品质及其与飞行条件的差异.第三部分是关于超高速流动的试验方法和数值技术研究,包括高焓流动中真实气体效应对飞行器俯仰力矩变化的影响;热化学反应流动中表面催化效应诱导的气动热变化规律;喷管流场的气流非平衡效应对试验结果可能产生的影响.     

高焓流场球头外形气动热试验研究

再入地球大气层时,飞行器的再入速度极高,面临严重的气动加热问题。为了研究高焓流动导致的热化学非平衡现象,在高焓膨胀风洞FD-14X中开展了球头外形的热流测量试验以及CFD仿真模拟。FD-14X为中国空气动力研究与发展中心新建成的高焓膨胀风洞,速度模拟能力达到第二宇宙速度,总温模拟能力超过10 000 K,能够产生总焓70 MJ/kg的试验气体。试验来流总焓16.9～63.5 MJ/kg,球头直径20～50 mm,流场采用自发光拍照,同时CFD仿真采用Park双温非平衡模型计算球头绕流流场。试验与仿真结果表明：来流总焓大于5 MJ/kg时,球头绕流场存在显著的热化学非平衡现象;304钢模型壁面在来流总焓小于20 MJ/kg时表现为非催化壁面特性,在来流总焓大于30 MJ/kg时表现为催化壁面特性;当球头表面镀氧化锆膜、来流总焓49.5 MJ/kg时,球头壁面表现为非催化壁面特性。     

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

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

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

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

二维翼段颤振的μ控制

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

μ理论在电液负载模拟器中的应用

 研究μ理论在电液负载模拟器中的应用，提出电液负载模拟器的鲁棒控制策略。电液负载模拟器是一个复杂的机-电-液复合系统，且是一个强耦合、时变受控对象。综合考虑参数变化、模型变动和干扰等不确定性的影响，利用μ综合控制理论设计电液负载模拟器的鲁棒力控制器，并通过μ分析使用鲁棒力控制器时系统的鲁棒稳定性和鲁棒性能。给出使用鲁棒力控制器和经典力控制器时的实验结果，实验结果表明所设计鲁棒力控制器的有效性和优越性。     

Abnormal Shape Mould Winding

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

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

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

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

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

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

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

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.      

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…