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

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

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

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

30cm离子推力器三栅极组件设计参数对寿命的影响研究

为了研究30cm离子推力器三栅极组件设计参数对预估寿命的影响,在完成失效模式分析的基础上,通过PIC-MCC方法对离子推力器三栅极组件的离子溅射速率进行了计算,建立起栅孔二维寿命预估模型,并针对栅极设计参数对预估寿命的影响进行研究。结果显示:导致三栅极组件的主要失效模式为5kW高功率模式下的离子直接轰击所造成的栅极早期结构失效,且减速栅的过快离子溅射腐蚀成为影响三栅极组件寿命的关键,而不同工作模式不会产生新的失效方式,仅影响栅极的离子溅射速率以及寿命;在现有三栅极设计参数条件下,当推力器工作时,栅极引出的离子束流处于明显欠聚焦状态,且加速栅寿命预估值约为9062h,而减速栅约为2642h;通过PIC-MCC方法得到的栅极三个关键设计参数对寿命的影响模拟结果显示,降低加速栅电压对提升减速栅寿命的作用较小;缩小加速栅与减速栅冷态间距后,离子溅射速率会随着冷态间距的缩小逐渐降低,冷态间距由1mm缩小至0.6mm后,减速栅在5kW工况下的工作寿命可提升至10726h,且经试验验证该间距可满足推力器力学环境试验要求;缩小屏栅孔径对改变离子束流引出形状具有显著作用,单孔束流发散角度随着屏栅孔径的缩小出现了明显降低,且束流离子几乎不会再直接轰击至减速栅上游区域,当屏栅孔径由1.9mm缩小至1.6mm后,减速栅工作寿命可提升至9259h;分析结果对后续开展栅极组件的寿命优化设计提供了参考。     

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

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

射频离子推力器多元工质束流调节试验研究

为了获得射频离子推力器离子束流随放电参数的变化规律,采用试验研究的方法,就推力器引出束流与射频功率强度、工质种类、工质流量之间的调节规律开展了研究,搭建了射频离子推力器束流调节试验系统。研究结果表明:屏栅电压1200V,加速电压-250V,射频功率200W～700W,工质流量0.2mg/s～4.76mg/s,Xe,Ar,O_2,N_2四种工质下能够可靠放电并稳定引出,实现束流从54mA～467mA的调节,电离效率XeArO_2N_2,离子束流随射频功率和工质流量线性增加,在1.01mg/s的氙工质下,推力、比冲随射频功率从100W～400W线性增加实现推力7.35mN～27.5mN,比冲1191s～3696s大范围连续可调,工质利用率为21.1%～78.8%,并在射频功率为276W时工质利用率和功耗之间存在明显拐点,在应用中要根据任务选择最佳工作区间,合理控制工作参数可以提高推力器工作性能和效率。     

离子推力器栅极极限寿命优化

离子推力器的极限寿命最终取决于栅极的极限寿命。针对LIPS-200离子推力器延长寿命到20000h以上的工程应用需求,在分析离子推力器极限寿命所对应关键失效模式及磨损机理的基础上,基于加速电压能够有效调节关键失效模式发展进程的工作机制,提出了具有普适性的离子推力器栅极极限寿命优化的恒定加速电压方法和步进调节加速电压方法。结合LIPS-200离子推力器寿命试验的过程及最终结果数据,在完全继承推力器现有技术状态和成熟度的前提下,采用恒定加速电压方法可以将推力器的极限寿命从现有的14649h提高到17300h,采用步进调节加速电压方法可以将推力器极限寿命提高到20400h,从而实现LIPS-200延长寿命目标。     

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

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

两种会切场离子推力器对比研究

环型会切场离子推力器和柱型会切场离子推力器是当前广泛应用和研究的两种会切场离子推力器。基于30cm环型会切场离子推力器LIPS-300H和30cm柱型会切场离子推力器LIPS-300Z,对比研究了两类会切场离子推力器各自优劣及其机理。首先分析了两种会切场原理,总结给出了两种会切场差异,然后实验对比研究了两种会切场离子推力器束流均匀性、放电效率和寿命。实验结果显示：LIPS-300H相比LIPS-300Z在3kW和5kW工况下束流密度峰值分别降低25%和19%,放电电压分别降低7.8V和6.2V,放电损耗分别增加20W/A和32W/A,屏栅预测寿命分别增加6.7倍和3.2倍。试验结果表明：虽然LIPS-300Z比LIPS-300H具有放电损耗低的优点,但其较差的束流均匀性,较高的阳极电压和双荷离子比,使其在寿命和可靠性方面劣于LIPS-300H。     

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

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

700W功率HET-80霍尔推力器束流特性研究

霍尔推力器的束流特性能够体现推力器的工作状态,为了了解额定功率1350W的Hall Effect Thruster-80(HET-80)霍尔推力器在700W功率下的工作状态,研究设计了一套等离子体诊断系统,包括了法拉第探针、平面型Langmuir单探针、阻滞能量分析仪(RPA)和发射光谱,对HET-80霍尔推力器在700W功率下的束流特性进行了诊断实验。实验结果显示:700W功率工况下HET-80霍尔推力器束流发散半角为45°,整体呈现出纺锤形分布,束流中心区域的离子电流密度分布呈现出圆环型霍尔推力器典型的双峰结构;由于等离子体的复合碰撞导致离子量变少,实验中法拉第探针测量得到的积分电流随着轴向距离增加而减小;在束流区边缘发生了较多的CEX碰撞,使该区域电子温度骤升;中心区离子能量接近247e V,相较于推力器放电电压310V,评估霍尔推力器加速效率达79.6%,加速效率较好,整体束流结果表明HET-80霍尔推力器在700W功率工况下表现正常。实验中发现法拉第探针结合RPA测量出来的离子数密度较Langmuir探针测量出的离子数密度更为准确;由于霍尔推力器只有一个空心阴极,发射光谱得到的出口附近电子温度和离子原子相对浓度分布均为空间非均匀状态。     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

一种三维结构网格并行算法

基于三维结构网格,以MPICH2作为并行开发环境,研究了对等模式下流场的并行计算问题。在三维结构网格已经生成的前提下,针对计算区域的划分,网格输出格式和数据通信,提出了一种并行算法。算法的核心思想是将需要通讯的其他区域的单元按区号从1到n排列,并获得本区域通信单元的I-J-K编号与其他区域通讯单元的排列序号的映射关系。最后,在IBM BladeCenter HS22刀片机平台上进行了算例验证,结果表明算法具有很好的可靠性与通用性。     

网格技术

网格计算是一种新兴的解决大规模复杂问题的应用技术，网格是Internet的发展方向。文章阐述了网格的定义、体系结构、以及实现网格计算的关键技术，介绍了国内外网格技术的研究现状、应用前景和发展趋势。     

飞机外形网格拓扑和多块网格拼接技术

多块网格方法成功的关键在于网格拓扑的选择和块与决之间的拼接方法。针对真实飞机外形多块网格生成,本文分析研究了翼类和体类部件的各种网格拓扑;论述了翼身组合体的多块网格非结构连续拼接技术;指出连续拼接方法的缺点,以及相应的解决方法;同时,还讨论了分块策略和确保块与块之间网格线完全连续拼接的方法;通过实例进一步阐述了本文工作的另外一些独到之处。     

二维及三维贴体网格的一种生成方法

本文讨论了二维及三维情况下应用Ｌａｐｌａｃｅ方程生成贴体网格的方法，并应用ＡＦ１格式求解方程以大大加速迭代收敛的速度，对于二维及三维情况，分别讨论了相应用Ｄｉｒｉｃｈｌｅｔ边界问题解法。对于二维情况，还讨论了翼型的Ｏ－、Ｃ－和Ｈ－型三种典型的网格生成方法。着重讨论了有关割缝上的网格点的处理方法，以保证网格的光滑性。三维情况下，以简化航天飞机外形为例，讨论了网格的生成方法，尤其着重讨论了边界点的处理     

复杂空气动力外形网格生成技术

通过对用于网格生成的Ｌａｐｌａｃｅ方程进行抛物化差分离散,发展了一种快速抛物型方程贴体及与边界正交网格生成方法,另外,采用源项修正技术求解用于网格生成的Ｐｏｉｓｓｏｎ方程,发展了一种可实现对网格线与边界间距及正交性进行双重控制的椭圆型方程方法,运用这种两种基本方法以及这两种方法的混合,分别生成了真实歼击机全机外形,带十字叉翼导弹以及双发,下单翼运输机外形的计算网格。并已成功地秀于跨声速得杂流场的娄     

基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。     

先进复合材料格栅结构制造工艺技术的研究进展

介绍了复合材料格栅结构的特点及发展现状,并回顾了传统复合材料成型工艺方法和格栅结构的早期制造方法。另外,介绍2种新型的格栅结构加工方法,即混合加工法、膨胀模块加工法以及改进的格栅结构:互锁复合材料格栅结构、加套管的槽形连接格栅结构和TR IG结构;并分别说明了每种结构的优点。其中重点介绍了混合加工法并从理论上推导其可行性。最后总结了复合材料格栅结构制造工艺的发展趋势。     

基于三维非结构混合网格的离散伴随优化方法

伴随优化方法在优化过程中不受设计变量个数的限制,有巨大的工程应用前景。基于非结构混合网格技术,采用雷诺平均Navier-Stokes方程、离散伴随优化方法以及自由形面变形（FFD）技术,建立了飞行器气动外形优化设计系统。离散伴随方程采用公式推导的方法直接获取,并运用LU-SGS迭代求解。通过与差分结果进行比较的方式,对通量的雅可比矩阵和目标函数敏感导数开展了确认,并对文中黏性通量项、限制器以及湍流黏性系数的近似处理方法做了分析对比,验证了近似处理方法的合理性和可行性。采用建立的优化系统,完成了ONERA M6机翼在跨声速条件下的减阻优化,并比较了有无容积约束下优化结果的差别。模拟验证结果表明,建立的飞行器气动外形优化设计系统具有较高的可靠性和有效性,可以用于三维飞行器外形的减阻优化。