电推进羽流粒子束溅射的蒙特卡罗模拟

为了研究电推进羽流对推力器及航天器造成的溅射腐蚀作用,采用蒙特卡罗方法对羽流等离子体对材料的溅射进行模拟,研究了溅射的级联碰撞机理,并开展了Xe离子正入射、斜入射Cu材料的溅射率、溅射产物的能量分布以及空间角分布的研究.结果显示,蒙特卡罗方法对溅射率和能量分布的计算结果与试验结果吻合较好,误差<30％,但是斜入射溅射率...     

霍尔推力器放电室壁面溅射产额研究

霍尔推力器放电室壁面经受低能量(不大于300e V)离子的溅射是影响其寿命的关键因素之一。为了获得放电室壁面材料(BNSi O2)的溅射产额随离子入射角度和能量的变化规律,采用真实霍尔推力器提供275e V的氙离子在真空舱内轰击靶材,利用称重法获得实验参数下的溅射产额。为了克服单纯依靠实验测量耗时耗钱且更低能量的离子溅射实验测量误差会陡然增大的缺点,采用前面实验结果修正了基于蒙特卡罗(MC)方法的SRIM软件溅射产额计算参数,并采用文献实验结果对不同能量下的模拟结果进行验证。在此基础上,用SRIM软件较为详细地考察了入射离子能量低于300e V时入射角度和能量对霍尔推力器放电室壁面材料溅射产额的影响规律。结果表明,溅射产额随离子入射角度先增大后减小,而随入射离子能量则呈现增大的趋势,但当能量小于100e V时,溅射产额逐渐趋于一个非常小的数值。     

三栅极离子推力器电子反流失效影响参数的敏感性研究

电子反流失效模式是离子推力器关键失效模式之一,决定推力器工作寿命。为明确各参数对电子反流失效模式的影响程度,确定加速应力,为地面加速寿命实验验证方案和长寿命优化设计提供数据支持,采用Hybrid-PIC-MCC (Particle in Cell-Monte Carlo Collision)方法,构建了三栅极系统数值仿真模型。采用模型研究了地面真空舱本底压力、屏栅电压、加速栅电压、屏栅与加速栅间距、屏栅上游等离子体密度和放电室工质利用率等参数的影响敏感度对比。研究结果显示,真空舱本底压力可以作为加速寿命试验的首选加速应力,在推力器结构和工作本征参数中工质利用率为最敏感应力,其次是屏栅电压、屏栅上游等离子体密度、加速栅电压、屏栅和加速栅间距。     

基于加速栅溅射腐蚀失效的离子推力器寿命预测

离子推力器加速栅溅射腐蚀失效是制约离子推力器寿命的关键失效模式之一.针对离子推力器长寿命、多功率条件下运行的特点,基于坑和凹槽的溅射腐蚀数据,建立模型对其进行寿命预测.通过研究离子推力器加速栅中心凹槽腐蚀深度在不同功率段下随工作时间的变化规律发现:运行功率顺序对加速栅凹槽腐蚀率影响较小,进而采用累积损伤理论建立离子推力器多功率段下运行的寿命预测模型.最后, 对美国的NASA's Evolutionary Xenon Thruster(NEXT)进行了寿命预测,预测结果寿命为46041h,与试验结果符合较好.      

三栅极离子推力器电子反流失效影响参数的敏感性研究（空间推进专刊）

针对电子反流失效模式主导的三栅极离子推力器加速寿命试验加速应力选择及长寿命优化,需要开展影响参数的敏感性对比研究,采用Hybrid-PIC-MCC（Particle in Cell- Monte Carlo Collision）方法,构建了三栅极系统数值仿真模型。采用模型研究了地面真空舱本底压力、屏栅电压、加速栅电压、屏栅与加速栅间距、屏栅上游等离子体密度和放电室工质利用率等参数的影响敏感度对比。研究结果显示真空舱本底压力可以作为加速寿命试验的首选加速应力,在推力器结构和工作本征参数中工质利用率为最敏感应力,其次是屏栅电压、屏栅上游等离子体密度、加速栅电压、屏栅和加速栅间距。研究结果为三栅极离子推力器地面加速寿命试验验证方案设计和长寿命优化设计提供了数据支持。     

基于半经验模型的霍尔推力器壁面腐蚀形貌预测方法研究

为了研究霍尔推力器通道壁面溅射腐蚀对于推力器寿命的影响，针对霍尔推力器通道壁面溅射腐蚀演化过程，建立了预测推力器腐蚀形貌的半经验模型。该模型根据试验测量形貌数据，反推离子源模型，结合腐蚀速率公式，对霍尔推力器的壁面形貌进行演化预测。分别以SPT-100，T-220及KM-45短时间的壁面形貌作为输入条件，对不同功率量级推力器的壁面腐蚀过程、预测误差和腐蚀速率进行了研究。数值模拟结果表明，对于10kW级T-220，1kW级SPT-100以及百瓦级KM-45推力器，壁面形貌的平均误差分别为2.65%，2.88%和3.64%。推力器壁面形貌的预测结果与实际测量值基本一致，该模型可用于霍尔推力器壁面形貌预测以及寿命预估。     

来流对吸气式激光推力器冲量耦合系数的影响

针对大气吸气模式激光推进的原理,采用二阶精度的Roe格式,在气体模型为完全气体时,分别对无来流以及来流马赫数为2.5和5的条件下,环聚焦激光推力器内外流场的演化过程进行了数值模拟。结果表明,无来流、来流马赫数为2.5和5时的冲量耦合系数分别为2.42×10-4N.s/J,1.00×10-4N.s/J和1.13×10-4N.s/J,无来流时激光推力器的冲量耦合系数高于同工况下有来流的情况,但来流马赫数为5时由于激波不会传播到推力器前体而产生负推力,因此冲量耦合系数较马赫数为2.5时高。研究结果可为激光推力器的构形设计以及实验模拟高空环境下推力器的飞行情况提供一定的理论指导。     

离子推力器羽流场模拟以及Mo+CEX沉积分析

离子推力器工作产生的羽流会对航天器产生影响,严重时甚至会造成航天器无法正常工作,为了精确评估离子推力器羽流特性及其对航天器的作用,采用基于粒子轨道理论(PIC,Particle-In-Cell)的模型对复杂的航天器的离子推力器羽流进行了数值模拟,并结合最近几年发展起来的浸入式有限元(IFE,Immersed Finite Ele-ment),采用结构网格准确计算复杂边界电场。通过模拟,获得了Mo+CEX离子在卫星表面的最大可能沉积分布,定量分析了卫星表面Mo+CEX离子的最大可能沉积率,表明在垂直于推力器主束流方向的卫星组件的表面上容易产生较大的Mo+CEX离子污染沉积率,而平行于推力器主束流方向上Mo+CEX离子污染沉积率相对较小。     

离子推力器加速栅寿命概率性分析

交换电荷离子对加速栅极的溅射腐蚀是离子推力器的关键失效模式之一,基于交换电荷离子对加速栅溅射腐蚀的物理机理,对离子推力器加速栅工作寿命进行了概率性建模。利用该模型对20cm Xe离子推力器加速栅寿命和其达到预期寿命的可靠度进行了评估。结果显示加速栅的寿命近似服从高斯分布,当推力器工作环境压力近似6.7×10-3Pa时,加速栅工作寿命达到3kh的可靠度为0.9352。     

基于Matlab/Simulink的离子推力器输入参数设计模型研究

离子推力器的输入参数设计直接决定其工作性能以及寿命,为了有效降低在前期设计阶段,由于输入参数设计缺陷导致的输出性能和预估寿命不一致性,以及后期优化改进所带来的风险及高昂成本,根据离子推力器工作过程的宏观描述以及寿命理论预估方法,基于Matlab/Simulink建立起离子推力器输入参数设计模型,实现了根据推力、比冲、效率等设计指标完成输入参数计算以及寿命预估的目的。结果显示:以30cm离子推力器为例,根据参数设计模型得到5kW额定工况下的关键输入参数与目前实际输入参数的比对误差均在5%以内,而模型得到的关键指标参数与冷启动试验测试值的比对误差同样在5%以内,3600h后的加速栅孔径刻蚀计算结果与试验结果比对误差在10%以内,证明了离子推力器输入参数设计模型的合理性,输入参数设计模型可用于后续同类型离子推力器的早期输入参数确定以及设计指标匹配性评价。     

Experiment of Hall thruster plume effects on different ground experimental conditions

Reducing the influence of back-sputtering effect can effectively improve the accuracy of the measurement of the Hall thruster plume effect. Quartz crystal microbalance (QCM) was used to measure the deposition and sputtering distribution of HET-40 thruster plume on two different experimental conditions: case 1, using liquid nitrogen heat sinks and case 2, without using liquid nitrogen heat sinks. Meanwhile, X-ray photoelectron spectroscopy (XPS) was used to analyze the composition of the QCM surface after two experiments. The results of the two experiments showed that the sputtering rate under the condition of case 1 was slightly higher than case 2. Especially within the range of 90° to 110° relative to the thruster axis, case 1 experiment result showed sputtering effect, while case 2 experiment showed deposition effect. Through analysis of the experimental results, it can be found that using liquid nitrogen heat sink to reduce the temperature of the inner wall surface of vacuum chamber can effectively adsorb the particles sputtered by the plume and reduce the concentration of back-sputtering particles, leading to the above phenomenon.      

电子回旋共振推力器C/C复合材料栅极的力学性能

为了研究电子回旋共振推力器在搭载火箭升空过程中,过载和振动环境对栅极带来的影响,采用有限元分析软件建立了10 cm C/C复合材料栅极的有限元分析模型,计算分析了不同C/C复合材料栅极的力学性能、频率响应和振动模态。结果表明:20g过载下,栅极的最大变形量为6.12×10-4mm;在5~800 Hz激振下,栅极频率响应的最大位移为2.3×10-6mm;过载对栅极的影响要大于振动对它的影响,但是这两个值都在安全门限内;综合分析栅极频响和模态分析的计算结果,栅极仅在其第一阶模态处于外界激振发生共振。     

一种基于螺旋桨部件特性的螺旋桨建模方法

通过缩比法,利用螺旋桨通用部件特性获得期望研究的螺旋桨部件特性,提出了飞行速度不为零条件下的螺旋桨数学模型建模算法,同时,借鉴缩比后螺旋桨部件特性、螺旋桨定桨叶角工作性能曲线以及螺旋桨空气动力学原理,分析了静拉力状态下的螺旋桨功率系数、拉力系数、桨叶角、螺旋桨静态推力进距比阈值以及螺旋桨几何设计参数的相互作用关系,提出了静拉力状态下的螺旋桨数学模型建模算法。所述算法与Gas Turbine Simulation Program (GSP)软件仿真数据进行了数字仿真对比验证。结果表明:所提出的螺旋桨建模算法具有有效性,在前进状态下,螺旋桨拉力相对误差最大不超过6.6059×10-6，需求功率相对误差最大不超过5.5098×10-6，效率相对误差最大不超过6.6955×10-6。      

High-temperature oxidation behavior and analysis of impedance spectroscopy of 7YSZ thermal barrier coating prepared by plasma spray-physical vapor deposition

Quasi-columnar structure 7YSZ(yttria stabilized zirconia) thermal barrier coatings(TBCs) are prepared by plasma spray-physical vapor deposition(PS-PVD) onto pretreated and un-pretreated bond coating, respectively. An isothermal oxidation experiment of 7YSZ TBCs is carried out in the atmosphere of 950 °C in order to simulate the high-temperature oxidation process of engine blades. The isothermal oxidation process of 7YSZ thermal barrier coatings is investigated systematically by impedance spectroscopy. The electrochemical physical model and equivalent circuit of columnar 7YSZ coatings are established. Results show that the isothermal oxidation kinetic curve of columnar 7YSZ thermal barrier coatings appears to follow the parabolic law. A pretreatment of bond coating can reduce the growth rate of the thermally grown oxide(TGO) layer, restraining the initiation and propagation of microcracks between YSZ and TGO layers. The oxidation rate constants of 7YSZ coatings with pretreated and un-pretreated bond coating are 0.101×10~(-12) cm~2·s~(-1) and 0.115 × 10~(-13) cm~2 ·s~(-1), respectively. Impedance analysis shows that the content of oxygen vacancies decreases and the density increases after the TGO layer is oxidized for 150 h. In addition, shrinkage microcracks formed by sintering during the oxidation process is the main reason for an increase of the capacitance and a decrease of the resistance in the grain boundary of YSZ.     

速率偏频激光陀螺偏频速率模糊PID自适应控制的设计与实现

应用模糊 PID自适应技术设计并实现了激光陀螺偏频速率的控制系统,结果表明 :偏频台速率平稳性优于 3.1 2× 1 0 - 4;角速率正反换向速度快,在± 60°/ s之间的转换时间约为 1 4 ms,而且换向过程速率超调很小。该控制系统在实际系统中得到了良好的应用     

R141b在矩形微尺度通道中的两相流传热特性

设计搭建水力直径分别为1mm和0.5mm的矩形微尺度通道实验台,研究了以R141b型制冷剂作为工质的两相流沸腾传热特性。实验取热流密度为1～16kW/m2、质量流速为111.1～333.3kg/(m2·s)和质量干度为0～1,分析了三者对平均传热系数的影响,探究影响换热的主导因素。结果表明:热流密度较高时,平均传热系数随热流密度增加而减小,流动换热主要受到沸腾传热的影响;当质量流速较大且热流密度较低时,平均传热系数随热流密度增加而有所增长;热流密度较低时,平均传热系数随质量流速变化明显,热流密度升高到一定值后,质量流速对平均传热系数的影响很小;当质量流速处于111.1～333.3kg/(m2·s)时,平均传热系数随质量干度的增加而减小。      

空间电动力绳系电流峰值点变化规律研究

为研究电流峰值点在几种影响因素下的变化规律,建立一种能够捕捉绳系电子发射与收集自洽平衡过程的新算法——电路空间耦合算法。为验证该方法的计算精度,以1.35kW霍尔推力器为等离子体源,在真空舱内开展绳系的电荷收集试验,在电路参数方面,计算精度约为12.8%,在场参数方面,计算精度约为3.6%。在此基础上,针对不同偏置电压、绳系长度以及绳系直径,对绳系的绳上电流分布、电势分布以及空间电势分布等参数进行数值计算。结果显示:电动力绳系的电流峰值点会随着偏置电压升高、绳系长度增加及绳系直径增大而发生比例上的向阳极端漂移,揭示了电子轨道运动限制的机制在各类壁面电荷输运机制中占优,导致电流峰值点漂移的产生。     

三浮陀螺仪沿输出轴误差系数DO的产生机理分析

目前,三浮陀螺仪误差模型包含沿输出轴的一次项误差系数DO,该误差项尚无明显的物理意义,但 DO的标定结果在0.0002(°)/(h·g)～0.015(°)/(h·g).针对三浮陀螺仪的结构特点和测试方法分析了DO的产生机理,通过调整安装误差角的试验说明安装误差角带来的天向地速和一次项分量会影响DO,因此系统在选用陀螺误差模型时,需关注由测试方法带来的误差.同时,陀螺在设计、生产和调试过程中对安装误差角的影响应予以重视.最后用统计检验的方法说明DO对误差模型的线性影响不显著,为系统选取陀螺误差模型提供依据.     

全离散孔气膜冷却实验研究

模拟加力燃烧室工况, 采用两种不同开孔率的离散孔板, 研究了在不同温度、压差等情况下的冷却效率和换热系数, 得到了一些有价值的冷却换热规律及相应的实验公式。