热容燃烧室基于单点温度的热流测量方法研究 *

为了研究推力室内壁面与燃气间的传热特性,获得多喷嘴燃烧室内壁热载分布,从而为推力室设计及热结构分析提供参考,开展了热容燃烧室壁面热流测量方法研究。以位于距内壁一定距离位置的单点测量温度作为输入条件,文章提出了三种热流计算方法,并进行火箭发动机试验研究,得到了对应位置燃气对内壁的热流密度。可以看出,运用单点方法得到的热流密度与传统两点法结果符合较好,两者计算结果误差在10%以内。结果显示:不同点火时长工况下在同一位置相同时间点上重复性很好,在燃烧室圆柱段热流密度随着点火时长的增加而减小,而在喷管段变化规律则相反,热流密度随时间增加而增加。在提出的计算平均热流的三种方法中,热积累方法的计算结果最高,热平衡方法次之,瞬态法最低。在火箭发动机热容试验中,本文提出的方法可以用于测量燃烧室内壁面热流密度。     

燃烧室中辐射热流分布的蒙特卡罗计算

本文用蒙特卡罗法计算了燃烧室壁面辐射热流密度分布。计算中考虑了燃气参数如:压力、温度、燃气组分及燃气的吸收和发射特性沿燃烧室轴向、径向和周向的变化。详细地说明了计算过程并进行了算例计算。计算所得热流密度值与分布规律是合理的。      

液体推进剂火箭发动机推力室再生冷却通道三维流动与传热数值计算

应用湍流模型对液体推进剂火箭发动机再生冷却推力室通道的流动与传热进行了三维数值模拟,冷却工质为氢气,其密度、导热系数、动力粘度随着温度和压力而变化,冷却剂比热容及金属固体物性随着温度而变化。计算采用标准k-ε双方程湍流模型及气-固耦合算法。结果表明:推力室燃气侧壁面的温度和热流密度的最高点均发生在喉部附近,喉部横截面固体区域最大温度梯度靠近燃气,喉部附近氢气在垂直主流方向的截面上产生了二次流。气固耦合面最大热流密度及最大对流换热系数同样位于推力室喉部附近。      

矢量喷管内燃气辐射与壁面温度的耦合计算

选取内、外调节片和隔热屏建立几何模型,基于封闭腔净辐射模型和壁面热平衡模型建立了燃气辐射与喷管壁面温度的耦合算法.波段为1~5μm的气体辐射采用窄波带模型计算,其他波段不考虑气体辐射,建立辐射净热流密度-有效辐射亮度-壁面温度的关联式求解燃气与壁面的辐射换热,采用牛顿-拉斐尔森迭代法求解壁面热平衡方程计算其温度.对某轴对称矢量喷管(偏转20°),计算了喷管壁面的红外光谱辐射和辐射净热流密度,以及各部分结构的温度.作为验证,还计算了文献中某液体火箭发动机轴对称矢量喷管壁面的辐射净热流密度,与文献的结果进行对比一致性较好.研究表明:轴对称矢量喷管偏转段沿周向的辐射热流密度和温度差异很大,沿偏转方向部位壁面的温度和辐射热流密度都较低,偏转方向壁面的温度比相反方向大约低10%,辐射热流密度大约低50%.      

气气喷注器推力室传热特性研究

设计了单喷嘴气气喷注器容热式推力室,进行了0.92～6.1 MPa范围内7个燃烧室压力工况,共17次热试车;采用壁面测温方法获得推力室沿轴向方向的内壁热流分布,得到各室压工况下燃烧室内壁热流分布曲线,不同工况的热流曲线显示出相似的分布;并从中分离出对流传热热流,得到对流传热热流与室压的关系。为拓展范围,采用多组分湍流N-S方程描述推力室内燃烧流动,采用6组分9步反应模型来描述氢氧反应,反应速率由Arrhenius公式计算,进行了5～20 MPa更高室压范围内的燃烧内流场的数值模拟,并耦合计算了各工况燃气与室壁之间的传热,获得了与试验相同规律的结果。     

轴对称喷管与圆转方喷管冷却换热特性的比较

为了解和比较轴对称喷管与圆转方喷管不同的再生冷却换热特性,分别对轴对称喷管(推力室)与圆转方喷管(推力室)建立计算模型,通过数值模拟的方法重点研究和比较了轴对称喷管与圆转方喷管的流场、壁面热流密度和温度分布、冷却剂温升和冷却通道压降等换热特性.计算结果表明:圆转方喷管由于型面不连续,在转方位置后壁面出现了温度和热流密度的峰值,从而导致沿周向壁面温度和热流密度的分布也不均匀.      

液体火箭发动机推力室复合冷却流动与传热研究

为了预测液体火箭发动机推力室的复合冷却性能,建立了推力室再生冷却通道和超临界氢的三维仿真模型以及推力室内燃气和超临界氢膜的轴对称二维仿真模型。通过边界耦合发展了液体火箭发动机推力室复合冷却流动与传热的数值仿真方法。对航天飞机主发动机推力室内部燃气、超临界冷却膜、室壁和再生冷却剂进行了流动与传热耦合计算仿真研究。研究表明,仿真方法可较好地预测推力室燃气及再生冷却剂的流动和传热,计算得到航天飞机主发动机的燃气侧壁面最高热流密度为129MW/m2,最高壁温为885K,冷却剂温升为192K,压降为8.8MPa,结果与已有数据吻合较好。模型和仿真方法可用于液体火箭发动机推力室冷却系统传热计算和冷却结构的优化设计。     

不同材料发汗冷却结构下氢冷却性能

为了实现航空航天等领域高温大热流燃烧装置的有效冷却,研究了不同材料和工艺制成的发汗冷却结构在高温高热流密度下,氢的发汗冷却性能。模拟高压推力室的结构特点和高热流设计发汗冷却试验件,用电弧加热主流空气模拟高温燃气、以氢气为发汗冷却剂对多孔陶瓷、烧结多孔不锈钢和多孔层板材料进行了33次172 s热试验研究。试验的材料设计孔隙率为10%~40%,燃烧室压力为2.7~8.4 MPa,主流燃气温度约为3 600 K,主流空气流量为220~1 490 g/s,冷却氢气流量为9.6~57 g/s,注入率为0.005~0.029。试验结果表明:当冷却剂氢注入率为1%时,主流与多孔陶瓷材料壁面和粉末冶金多孔结构壁面之间的换热分别减少了30%和70%以上;当注入率为3%时,主流与光刻多孔层间结构壁面之间的换热也能降低60%。证明氢发汗冷却可以有效减小壁面与燃气之间的对流热流。最后还总结得出了常温氢气对高压大热流环境进行发汗冷却的性能关联式。      

再生冷却推力室耦合传热数值模拟

对大推力液体火箭发动机再生冷却推力室内部燃气、室壁和再生冷却剂进行了耦合传热数值计算.采用二维轴对称N-S方程描述推力室内部燃气的湍流流动与传热,对冷却剂流动采用简化的一维模型,通过室壁的偶合传热采用一维传热模型.N-S方程的求解采用贴体坐标系下的有限容积法,速度和压强的耦合采用可压缩的SIMPLE算法.湍流的模拟采用可压缩的标准κ-ε模型,辐射传热采用热流模型计算.研究表明,本文方法可较好地模拟燃气的二维流动,同时能快速计算壁面热流密度、壁温和冷却套温升,计算结果对大推力发动机推力室的设计具有一定的指导意义.     

高热流密度燃烧室壁面阵列空气射流-燃油组合冷却结构研究

针对高热流密度燃烧室壁面热防护需求，提出了一种空气阵列射流冲击和燃油冷却肋板的集成冷却方式，在射流平均雷诺数Re_j为1×10⁴~3×10⁴，燃油进口流速v_f为2.33～5.23m/s内，采用数值模拟方法对其传热特性进行了研究，并基于壁面加热侧当量对流换热系数的概念，分析了基准肋板以及燃油冷却肋板的传热增强作用。与无肋板靶面的阵列射流冲击相比，带肋板阵列射流冲击的面积平均当量对流换热系数是前者的1.6倍，压力损失系数相对提高了约25%；采用燃油冷却肋板，加热壁面综合传热能力进一步增强，在Re_j=1×10⁴时，采用燃油冷却肋板的面积平均当量对流换热系数是基准肋板的1.5倍以上，即使在Re_j=3×10⁴时，燃油冷却肋板的传热增强比也可以达到1.2；燃油冷却肋板的出口温度相对进口温度的提升在20~50K内，其提升幅度随着射流雷诺数或燃油进口流速的增大而减小。     

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