离子推力器欠聚焦冲击电流的数值模拟

离子光学系统的离子束引出过程是离子推力器重要的物理过程,该过程直接关系到推力器的推力、比冲、效率等参数。为研究离子在离子推力器光学系统中的运动特性,使用了基于IFE-PIC（Immersed Finite Element Particle-In-Cell）的离子推力器光学系统离子束引出过程的三维数值计算模型,计算了栅极间电场分布、电荷密度,栅极冲击电流及欠聚焦极限。计算结果表明,当屏栅极电压不同时,发生欠聚焦的等离子束电流也不同。在欠聚焦工况下,一部分离子与栅极碰撞,产生冲击电流。冲击电流随电离室等离子体数密度增加而增大。     

离子发动机栅极系统电子回流数值模拟

以离子发动机栅极系统为研究对象,基于网格质点法(PIC),建立了二维数值计算模型.利用模型能够较好地模拟离子在栅极系统中的运动,得到栅极系统周围的电势分布以及离子和电子的密度分布等.进一步的计算结果表明,在几何参数一定的条件下,加速栅极电压以及加速栅极孔径对加速栅极下游的电子密度分布有重要的影响.选取合适的加速栅极电压和孔半径,可以阻止电子回流的产生,使离子发动机获得较佳的性能.      

考虑电荷交换的三维离子发动机栅极粒子模拟

为了对离子发动机光学系统进行数值分析，开发了考虑电荷交换的三维离子发动机栅极粒子模拟程序。模拟程序采用单元内粒子（PIC）方法，光学系统为双栅极结构，推进剂采用氙，总加速电压为1600V，栅极间距为1mm，粒子之间的电荷交换碰撞模型采用蒙特卡洛方法。采用链式存储结构存储粒子信息，可以实现数组式存储方法所不能实现的功能。使用粒子模拟程序对离子发动机光学系统进行了粒子模拟，模拟得到了栅极间的电场分布和氙离子空间分布。模拟结果表明，在所取工作参数下氙离子全部通过了栅极孔，未撞击到栅极上，验证了参数的合理性。     

离子发动机加速栅极腐蚀深度的DFF测量与数值模拟

使用聚焦深度表面测量(DFF)方法对加速栅极下游表面腐蚀深度进行了测量,并将测量结果与数值模拟结果进行了比较,所使用的数值方法为PIC-Monte Carlo方法.利用数值模拟程序对离子发动机栅极腐蚀进行了数值模拟.以氙为推进剂,栅极材料为钼.用蒙特卡罗方法模拟了氙离子与中性氙原子之间的电荷交换碰撞.模拟得到了加速栅极下游表面离子溅射腐蚀的深度分布,腐蚀模式与"Pits and grooves"模式相吻合.      

离子发动机栅极系统中束流离子的三维模拟

应用网格质点法,对离子发动机栅极系统进行了三维数值模拟.模拟得到了离子发动机栅极系统中的电场分布、离子运动轨迹、离子在计算域内的空间分布和相空间分布等.计算结果表明,屏栅的离子透明度为69.3%,大于其物理透明度67.1%,加速栅的离子透明度为100%,远大于其物理透明度31.1%,说明所选模拟工况参数合理.模拟结果可以为离子发动机栅极的初步设计提供理论参考.      

离子发动机加速栅极孔扩大腐蚀的粒子模拟

为研究离子发动机加速栅极孔扩大腐蚀速的特征及影响规律,采用单元内粒子——蒙特卡罗碰撞(PIC-MCC)方法,结合溅射腐蚀模型,模拟了栅极间束流离子及交换电荷(CEX)离子的分布情况,得到离子发动机加速栅极孔扩大腐蚀情况及平均腐蚀率,为0.0176μm/h。进一步研究了工作参数和几何参数对其影响,并建立了适用于工程应用的加速栅极孔扩大腐蚀预测模型。结果表明,在影响孔扩大腐蚀率各种因素中,束流电流密度和中性气体密度的影响最大。因此,在离子发动机直径一定的情况下,提高放电室出口等离子体密度的均匀性及推进剂利用率,是提高栅极寿命的根本途径。      

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

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

射频自偏压效应离子推力器束流特性实验研究

为了获得不同推力器工况对射频自偏压效应离子推力器的自偏压幅值和束流特性的影响，本文通过地面实验研究了栅极射频功率、线圈放电功率、工质种类对自偏压幅值和羽流区等离子体参数的影响，同时对直流栅极工况和射频栅极工况下的束流特性进行了对比。研究结果表明：栅极射频功率的增大会提高自偏压幅值并提升束流强度，但在较高栅极射频功率下，栅极下游区域将发生自持放电并形成等离子体；放电腔内放电模式转换会通过改变等离子体阻抗的方式大幅影响栅极直流自偏压幅值和栅极电压的射频分量，进而影响羽流区等离子体参数；与直流栅极工况相比，射频栅极能同时引出并加速离子和电子，并在栅极下游实现自中和，且在Ar,Kr,Xe三种工质下均具有自中和能力。     

不同环境温度对30cm离子推力器三栅极组件的离子刻蚀速率影响分析

为了对不同环境温度造成的30cm离子推力器三栅极组件离子刻蚀速率的影响进行分析,采用有限元仿真与试验验证相结合的方法,计算并试验验证了不同环境温度下的三栅极组件热平衡温度以及栅极间的相对位移变化,采用流体方法模拟了不同环境温度(20℃,-70℃,-120℃和-170℃)对三栅极组件的刻蚀影响,并结合短期寿命试验结果进行验证。结果显示:随着环境温度的降低,屏栅达到温度平衡的时间无变化,而加速栅温度平衡所需的时间则明显延长,20℃下的屏栅和加速栅热仿真结果与室温下推力器热平衡试验结果比对误差分别为7%和5%;其次环境温度的降低,会导致屏栅与加速栅的中心间距和边缘间距均缩小,而加速栅和减速栅的边缘间距却逐渐拉大,仿真结果与栅极热间距摄像测量结果符合性较好;根据三栅极组件的栅孔径扩大率随环境温度变化的计算结果来看,加速栅中心和减速栅边缘是离子刻蚀的主要位置,轰击至加速栅中心区域的离子数速率约是边缘的3倍,而轰击至减速栅边缘区域的离子数速率是中心的2.5倍,且环境温度的降低对加速栅中心区域离子刻蚀的影响更为强烈;经2100h的寿命试验验证,仿真结果与试验结果基本符合,误差经分析认为主要来自于流体方法的参数设置过程以及栅孔壁面均匀刻蚀的计算假设。     

LIPS-200离子推力器放电室出口离子密度分布研究

为了明确国内200 mm口径离子推力器放电室出口(即栅极上游附近)离子密度径向分布,采用实验与数值仿真相结合的方法对LIPS-200推力器放电室出口离子密度进行研究。应用法拉第筒分别测试推力器栅极下游50mm和100mm位置处束流特性,结合经验模型计算出栅极出口(z=0mm)束流离子径向分布。在此基础上,通过栅极数值模拟仿真,分析出栅极系统透过率随栅孔电流变化关系,进而反推计算出放电室出口离子密度径向分布。结果显示:放电室出口离子密度平均值约为9.0×10~(17)m~(-3),最大值约为1.54×10~(18)m~(-3),最小值约为4.6×10~(17)m~(-3);离子密度径向分布具有较好的中心轴对称性,离子密度从中心处沿着径向先缓慢减小,在径向位置约为50mm时出现快速下降;对比放电室出口与栅极出口离子密度径向分布发现,中心位置两者相差最大,边缘处相差最小。     

Numerical Simulation of Characteristics of CEX Ions in Ion Thruster Optical System

Charge exchange(CEX) ions could inflict severe damages on the ion thruster optical system. This article is aimed at investigating the characteristics of the CEX ions and their influences upon the optical system by means of particle-in- cell(PIC) ion simulation and Monte Carlo collision(MCC) methods. The results from numerical simulation indicate that despite the fact that CEX ions appear in the entire beamlet region near the ion optical system, the ones that present themselves downstream of the accelerator grid have good reason for attracting more attention. As their trajectories are significantly affected by the local electric field, a great number of CEX ions are accelerated toward grids resulting in sputtering erosion. When the influences of the CEX ions are considered in the numerical simulation, there could hardly be observed augments in the screen grid current, but the accelerator grid current increases from zero to 1. 4% of the beamlet current. It can be understood from the numerical simulation that the CEX ions formed in the region far downstream of the accelerator grid should be blamed for the erosion on the downstream surface of the accelerator grid.     

Ion thruster accelerator grid erosion mechanism under extreme conditions of cylindrical erosion and chamfer erosion

In this paper, the abnormal experimental phenomenon on barrel erosion under extreme working conditions in the ultra-long life experiment(>10000 h) of ion thruster ion optics is studied by the Immersed-Finite-Element Particle-In-Cell Monte-Carlo-Collision(IFE-PIC-MCC) method and the grid erosion evaluation model. The transport process of beam ions and Charge Exchange(CEX) ions in the grid system, and the characteristics and mechanisms of the aperture barrel erosion under extreme erosion condit...     

离子发动机羽流特性的数值模拟

离子发动机羽流中产生的交换电荷(CEX)离子返流会影响航天器的正常工作.建立离子发动机羽流模型,采用单元内粒子方法(PIC)对羽流场进行数值模拟计算.结合DS-1探测器飞行实验的测量结果,分析了卫星电势、电子温度、卫星几何尺寸、推力器工作特性等对相关因素对CEX离子返流特性的影响.结果表明:从推力器出口附近到卫星背面,CEX离子密度为108~1012m-3.当卫星电势从-15V变化到27V,测量点位置处CEX离子密度从0.65×1012m-3变化到1.5×1012m-3.羽流中CEX离子密度和电势结构随电子温度变化不大,但电势大小随电子温度成比例地变化.同一位置处不同工况下CEX离子的密度可根据CEX离子生成率与工作点参数间的关系式准确地估计.卫星安装推力器的表面起着对CEX离子返流屏蔽和降低的作用.      

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

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

Z源NPC三电平变换器恒功率并网控制策略研究

Z源中点钳位式(NPC)三电平变换器把逆变和升压两种环节结合在一起,结构简单、工作效率高,但是针对这种结构输出的无功功率研究较少。提出了一种Z源NPC三电平变换器的恒功率并网控制策略。首先对这种变换器的结构及原理进行分析；再把恒功率控制引入到这种变换器的并网控制中,并采用空间矢量脉宽调制(SVPWM)实现了变换器的并网控制。MATLAB/Simulink软件仿真结果证明了这种变换器能够提高电压升压比。该控制策略可使变换器输出的有功功率、无功功率稳定于设定值,且并网电流谐波较低。     

Interface flux reconstruction method based on optimized weight essentially non-oscillatory scheme

Aimed at the computational aeroacoustics multi-scale problem of complex configurations discretized with multi-size mesh, the flux reconstruction method based on modified Weight Essentially Non-Oscillatory (WENO) scheme is proposed at the interfaces of multi-block grids. With the idea of Dispersion-Relation-Preserving (DRP) scheme, different weight coefficients are obtained by optimization, so that it is in WENO schemes with various characteristics of dispersion and dissipation. On the basis, hybrid flux vector splitting method is utilized to intelligently judge the amplitude of the gap between grid interfaces. After the simulation and analysis of 1D convection equation with different initial conditions, modified WENO scheme is proved to be able to independently distinguish the gap amplitude and generate corresponding dissipation according to the grid resolution. Using the idea of flux reconstruction at grid interfaces, modified WENO scheme with increasing dissipation is applied at grid points, while DRP scheme with low dispersion and dissipation is applied at the inner part of grids. Moreover, Gauss impulse spread and periodic point sound source flow among three cylinders with multi-scale grids are carried out. The results show that the flux reconstruction method at grid interfaces is capable of dealing with Computational AeroAcoustics (CAA) multi-scale problems.     

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.