螺旋波等离子体推力器放电过程的流体模拟

针对建立三维数值模型对螺旋波等离子体推力器放电进行数值模拟计算量过大的问题,使用COMSOL多物理场软件建立了二维轴对称结构,采用漂移-扩散流体模拟方法,分别改变工质气体种类、气体压强及射频功率,模拟了螺旋波等离子体放电过程,分析了不同参数条件下放电室中电子数密度、电子温度、碰撞功率损耗分布情况,结果显示,在一定条件下气体压强越大,电子数密度越高,电子温度越低,碰撞功率损耗越大,在1600～2400 W范围内,随着射频功率增大,电子数密度和碰撞功率损耗增加,电子温度变化较小,电子数密度在10~(18 ) m~(-3)左右,为螺旋波等离子体推力器的设计与实验研究提供了参考。     

离子推力器放电腔内等离子体流动规律的全粒子模型

为了详细描述等离子体在放电腔内的产生、演化过程,建立了放电腔内部等离子体流动的二维轴对称全粒子数值模型。在考虑壁面二次电子发射影响的前提下,模拟得到了稳态后放电腔内电磁场、电子与离子数密度,电流密度等一系列参数分布,且与相关实验数据进行了比较,吻合较好。模型得到了实验中难以观测到的放电腔内等离子体的产生、演化过程以及稳态下的分布规律,合理解释了放电腔工作的基本原理。     

LIPS-200环型会切磁场离子推力器热模型计算分析

为了建立国内自行研制的20cm口径LIPS-200环型会切磁场离子推力器放电室的热模型,研究了放电室内等离子体的产生过程,得到了二次电子的温度、离子密度以及电子密度分布规律,在此基础上得到放电室各个关键部件的电流沉积和能量沉积热模型。以热模型计算结果为依据,进行了推力器稳态工作下的有限元热分析以及热平衡验证试验。结果显示:推力器处于稳定放电时,放电室内电子温度分布范围为2~4e V;电离基本发生在放电室轴线附近,此处电离产生率和电子温度最高,并且整个放电室内离子密度约为1017/m3;放电室的内表面能量沉积主要来自二次电子及Xe离子。     

非麦氏电子分布对霍尔推力器等离子体与壁面相互作用的影响

为深入分析霍尔推力器放电通道的非麦氏电子分布对等离子体与壁面相互作用的影响,采用一维非稳态鞘层动力学模型,统计了等离子体与壁面相互作用的重要物理量。结果表明,非麦氏电子分布函数和麦氏电子分布函数下等离子体与壁面相互作用存在很大差异,电子服从非麦氏电子分布时入射电子在壁面上的能量沉积,以及二次电子对主流区电子的冷却作用都明显弱于电子服从麦氏分布的情形。     

等离子体气动激励抑制翼型失速分离的仿真研究

通过求解表面放电的二维流体体力模型,建立了翼型等离子体流动控制的数学模型,得到等离子体气动激励诱导的体力和热量分布,与Navier-Stokes方程耦合求解.进行了低雷诺数条件下,等离子体气动激励抑制NACA0009翼型失速分离的数值仿真研究,研究了等离子体激励的强度、激励电极数目和激励位置对流动分离抑制和翼型升阻特性的影响.在雷诺数为58000、攻角为24°的情况下,施加等离子体激励后,升力系数由0.7449增大到1.2404;阻力系数由0.4012减小到0.3503.      

电子回旋共振等离子体推力器放电机理数值模拟研究

电子回旋共振等离子体推力器(ECRPT)是一种高比冲、高效率且结构简单的新型电磁式推力器。为了研究推力器的放电原理和工作机制,采用漂移-扩散流体模拟方法,仿真模拟了微波等离子体放电过程。仿真结果表明,电子数密度达到10~(16)~10~(17)m~(-3)数量级,氙气的电子数密度比氩气高50%;电子数密度、碰撞功率损耗均随着计算域内压强的增大而增大,电子温度随压强的增大而减小;电子数密度、碰撞功率损耗随着入射微波功率的增大而增大。在未来ECRPT的实际应用中,可以通过使用氙气,适当增大推力器腔内压强以及入射微波功率,使其具有最佳的推力、比冲和工作效率。     

准直流放电等离子体对超燃燃烧室中碳氢燃料燃烧的影响研究

为研究低燃空比条件下,准直流放电等离子对超燃燃烧室中乙烯燃烧流场的影响,在凹腔上游以及底部前壁面处布置电极产生等离子体,通过数值模拟方法,分析了不同等离子体激励强度下,燃烧室凹腔后缘附近压力分布、燃烧室总压损失、乙烯燃烧效率和燃烧室中水的分布情况。研究结果表明:准直流放电等离子体激励强度越高,对凹腔后缘附近压力场稳定能力越强。等离子体的存在,使得燃烧室出口总压损失微弱增加,损失最大值增加1.9%。燃烧室中乙烯燃烧效率平均提高1.77倍,随着激励强度的提高,燃烧效率呈现先增高后降低的趋势。等离子体改善了燃烧室中水的分布,凹腔内部产物分布范围更广、燃烧更加充分。     

放电发射光谱在金属微弧氧化表面处理中的应用研究进展

液相火花放电现象是金属微弧氧化表面处理时典型特征,发射光谱(OES)技术是表征微弧放电光谱特征和探索微弧氧化机理的有效手段。本文综述了铝、镁、钛等金属微弧氧化过程中发射光谱的研究现状。介绍微弧放电区等离子体的电子温度、电子密度等特征参数计算原理。重点关注了不同金属基体材料、电参数及电解液组成条件下,等离子体放电行为对微弧氧化膜结构的影响规律,并比较不同放电模型的异同。基于OES谱线评估得到的各种金属微弧等离子体温度为3000~10000 K,为放电通道内快速熔化-凝固过程促进陶瓷膜生长机制提供证据。     

离子推力器放电室阳极壁面电流密度分布特性研究

为了准确掌握离子推力器放电室阳极壁面电流密度分布特性，并深入理解阳极壁面处等离子体运动特性，设计了近阳极壁面等离子体诊断的具体实施方案，并基于LIPS-200离子推力器开展了近阳极壁面处等离子体诊断试验研究，得到了主要磁极附近壁面等离子体参数，并得到阳极壁面吸收电流密度分布特性。试验结果表明：LIPS-200离子推力器阳极壁面处主要磁极附近的等离子体密度范围为，测试点的电子温度范围为，壁面电流密度范围为；柱段壁面电子温度相对锥段较低，但电流密度较大，尤其在中间极靴位置电流密度最大，约为阴极极靴处电流密度的3倍，约为屏栅极靴处电流密度的2倍，阳极电流主要在放电室中间极靴处发生损失。     

贫燃条件下当量比对等离子体辅助起爆的影响

为研究贫燃条件下当量比对交流介质阻挡放电等离子体辅助起爆的影响,采用松耦合方法仿真求解了当量比分别为1.0,0.8及0.6三种条件下氢氧预混气体起爆过程,对比研究了放电产物时空分布特性、不同阶段爆震管内流场参数分布、起爆时间和距离以及推力壁动态历程。结果表明放电产物在同一周期内的时空分布形态不随当量比改变而变化,但随着当量比下降,各主要活性粒子的密度增大,增幅则降低。随着管内流速上升,等离子体对流场结构的改变更加明显,加速了形成稳定爆震波的各个子过程;不同工况下等离子体都缩短了起爆时间和距离,但当量比越低,缩短程度越小。对推力壁受力的动态分析也显示等离子体对流场演化的加速效果随当量比减小而减弱,不过其受力大小与等离子体无关。在本文贫燃条件下出现低当量比不利于等离子体发挥作用的原因包含两点,一是燃料占比降低带来的不利影响大于等离子体助燃的效果,二是放电区域尺寸有限。     

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