碳基复合材料超高温热化学烧蚀过程的数值模拟

对超高温热化学烧蚀的计算方法进行了研究,对碳基复合材料在不同环境下热化学烧蚀过程进行了数值模拟,对烧蚀产物进行了计算.计算表明,碳基复合材料在烧蚀过程中依次出现氧化、升华过程.10atm下,氧化、升华过程交替温度为330K.压力升高,交替温度也依次升高.氧化产物为CO,而升华产物依次为C3、C1.在3300K以下,烧蚀主反应为2C O2→2CO.在3300K以上,烧蚀主反应为:2C O2→2CO,2C N2→2CN,i(C)→Ci(g)i=1,3.并用哈尔滨工业大学研制的"星形"等离子发生器及国产的等离子体火炬进行了碳基复合材料烧蚀实验.实验光谱检测结果与超高温热化学烧蚀计算的烧蚀产物一致,验证了超高温热化学烧蚀计算的正确性.     

存在烧蚀的高超声速压缩拐角流动数值研究

耦合求解热化学非平衡流控制方程和烧蚀壁面边界条件,进行存在石墨烧蚀的压缩拐角流场数值模拟。流场化学反应采用16组元(N2,O 2,NO,N,O,NO +,N+2,O +,N+,CO,CO 2,C,C2,C3,CN,e-)29个反应的非平衡模型,热力非平衡的双温度模型下,不同反应采用不同控制温度。石墨材料表面反应包括碳的氧化反应、碳催化的O 原子复合反应和碳的升华反应。对15°、18°、24°压缩拐角模型,在自由流 Ma =10~30,总焓值6~55 MJ/kg 范围,分别进行无烧蚀的壁面催化与非催化条件和石墨烧蚀条件下的流场计算,分析各类条件下的流场结构、流动分离特性以及流场热化学参数分布特点,研究壁面条件对流动特性的影响。结果表明：流动分离可能性和分离区范围随着压缩拐角斜面倾角增大而增大,随来流马赫数增大而减小;相对于低壁温条件,无烧蚀的辐射平衡壁温和壁面烧蚀条件下流动分离区增大,斜面上压力、摩阻和热流峰值点也有所后移。     

碳/碳化锆复合材料烧蚀机理和计算方法研究

研究了碳/碳化锆复合材料氧化烧蚀机理,发现它们与传统的硅基和碳基材料烧蚀有很大差别。基体ZrC氧化后在表面形成一种膨松状的多孔固态抗氧化膜,氧化膜能有效阻止材料进一步氧化,使烧蚀量大大降低。研究了烧蚀过程抗氧化膜的形成、演化和流失行为,研究了氧化膜中氧气的扩散机制(包括分子扩散和Knudsen扩散),研究了材料原始层表面可能存在的化学反应,建立了分析碳/碳化锆复合材料烧蚀响应的物理数学模型。     

气相环境下EPDM绝热材料双区体烧蚀模型

针对气相环境下EPDM绝热材料的烧蚀行为建立了双区体烧蚀模型,着重考虑了沉积反应、气流剥蚀和膨胀现象。多孔介质区和固体区分开求解,采用交界面温度耦合的数值处理方法。对气相环境下的烧蚀实验开展了数值计算,分析了影响炭化层孔隙结构、质量烧蚀率的主要因素,计算结果与实验结果吻合较好。计算还获得了炭化层内烧蚀气体流速的量级是毫米到厘米,相对压强的量级是103～104Pa以及相对压强的分布情况。     

基于热解过程的EPDM变热物性数值传热

针对三元乙丙橡胶(EPDM)类绝热材料热载荷作用下的热解炭化过程,基于热解动力学和多孔介质传质传热理论,建立了芳纶/EPDM绝热材料热物性参数随温度和时间变化的变热物性模型,并通过与实验的对比,验证了模型的准确性与可靠性。随后对热载荷作用下的烧蚀热响应开展了数值计算,结果表明:热载荷作用初期,材料表面升温迅速,随着能量不断传递,温度推进速率明显降低,炭化层厚度增长减慢,部分材料仍为原始状态；温度对热解反应速率的影响呈指数级,距表面越近反应速率越快,反之则慢。所提出的变热物性模型对绝热材料的烧蚀研究具有一定的参考价值。      

CO_2激光烧蚀水所形成气蚀空穴机理研究

分析了激光烧蚀液体时形成表面空穴的机理。利用高速相机获得CO2激光烧蚀蒸馏水的流场照片,对空穴的发展过程进行了分析,发现空穴的快速发展期在最初的100μs内,之后在几百μs内保持平衡。用简单模型计算了不同能流密度下,空穴平衡时的形状,并分析了模型的局限性。研究结果对于揭示激光烧蚀水产生推力的机理有重要意义,有利于促进水工质激光推进的发展。     

EPDM的烧蚀模型

作为固体火箭发动机燃烧室壁面绝热层的ＥＰＤＭ是一种炭化型热防护材料,受热时形成炭化层和热解层,在燃气和粒子冲蚀下不断减薄并带走热量,有效地保护燃烧室壁。采用气动热化学烧蚀机理,扩散和化学动力学双控制机制,并计入气流与粒子的侵蚀,建立了ＥＰＤＭ的烧蚀模型;同时将烧蚀与移动边界下的传热相耦合进行烧蚀率预示。根据模型预示的烧蚀率与试验结果吻合很好。     

含碱金属杂质的高超声速烧蚀流场特性数值模拟

针对高超声速飞行器的烧蚀现象,采用Newton迭代的空间推进方法,数值求解化学非平衡的抛物化Navier-Stokes(PNS)方程.以黏胶基碳布-氨酚醛树脂为烧蚀材料,建立了含碱金属杂质的18组分33反应的化学反应模型,计算了小钝锥体的壁面温度,并和现有文献结果进行对比验证.在此基础上,对含碱金属杂质的高超声速烧蚀流场进行了模拟,深入分析了含碱金属杂质的烧蚀反应对流场温度、压力和电子数密度的影响.结果表明:烧蚀反应使得流场温度较非烧蚀情况增加了10%~15%;考虑碱金属杂质Na后,流场电子数密度增大了1~2个数量级.      

基于化学反应器网络模型法的预混气组分对NOx反应路径影响研究

摘 要：为探究预混气组分对不同反应路径下NOx生成影响,针对高压射流反应器,基于CFD流场特征等,构建优化了单PSR、2PSR、3PSR以及PSR+PFR四种化学反应器网络模型。在贫预混燃烧条件下,研究了纯CH4及其混合气燃烧时不同反应路径下NOx的生成情况,得出了燃烧室内不同反应区域每条反应路径NOx的生成量。结果表明,纯CH4燃烧产生的NOx主要来自于热力型、快速型和N2O-中间体型三条反应路径。绝热火焰温度的提高主要促进了热力型和N2O-中间体型NOx的生成。随着CH4或CO中加入H2摩尔分数的增加,NOx排放总量降低,快速型NOx生成速率降低,特别是火焰刷区域的。另外,混合气中CO摩尔分数的增加会导致NOx生成量增多。     

C/SiC复合材料烧蚀机理和通用计算模型研究

研究了C/SiC复合材料氧化烧蚀机理,发现它们与传统的硅基和碳基材料烧蚀有很大差别。C/SiC的烧蚀取决于氧的分压、表面温度和材料晶态结构及成份,可能出现活性氧化和惰性氧化两种破坏机制。研究了氧化膜的形成和破坏条件,以及氧化膜中氧气的扩散机制,建立了适用于C、Si、SiC和C/SiC(C和SiC可有不同混和比)四种材料烧蚀计算的通用物理数学模型。     

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

二维翼段颤振的μ控制

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

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.      

Abnormal Shape Mould Winding

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

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

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

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

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

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

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

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

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