一种霍尔推力器离子能量分布测量方法及系统

摘要

本发明公开一种霍尔推力器离子能量分布测量方法及系统。该方法通过测量多栅探针不同制止电压下的离子电流随时间变化曲线，获取同一时刻不同制止电压对应的电流值，得到该时刻电流与电压函数关系，然后求得该时刻的离子平均速度，求得该时刻离子能量；以设定时间间隔获取振荡周期内其它时刻的离子能量，将不同时刻的离子能量按照时间顺序拟合，得到低频振荡周期内离子能量分布随时间的变化曲线。本发明实现了测量羽流区离子能量分布的动态特性，获得羽流区离子能量分布在低频振荡周期内的变化曲线，解决了现有的测量方法无法获得几十微秒量级羽流区离子能量分布随时间动态变化的难题，为研究霍尔推力器羽流区离子能量分布变化提供有效技术途径。

说明书

技术领域

本发明涉及离子体推进领域，特别是涉及一种霍尔推力器离子能量分布测量方法及系统。

背景技术

霍尔推力器自1970年成功研制以来，因其具有高效率、高比冲及高可靠性等优点，被广泛应用于卫星位置保持和轨道转移等空间推进任务。霍尔推力器羽流一般由电子、高能离子电荷交换产生的低能离子以及中性原子组成。其中高能量的离子与航天器表面碰撞使其产生溅射和侵蚀，溅射的物质会对航天器上的太阳能帆板等敏感表面产生污染，进而使其性能下降或失效。同时高能离子直接与敏感表面碰撞，使敏感表面被破坏而性能降级或失效，从而导致航天器无法完成预定任务。

霍尔推力器的羽流离子能量分布和通道内的电离过程是一一对应的，因此可以通过离子能量分布分析推力器中性原子的电离机制以及磁场和电势的分布情况，进而可以分析通道内离子束流的加速特性，也被用来评价推力器对飞行器表面相互作用的一个重要参数。因此对霍尔推力器羽流离子能量分布的评估测量是霍尔推力器及航天器总体设计研究的重要内容之一。

霍尔推力器的低频振荡特性会使羽流区离子能量分布发生变化，目前的测量方法只能测量推力器在稳态条件下的羽流离子能量分布，没有实际应用价值，而在低频振荡周期内(通常为几十个微秒)，很难获得羽流中离子能量分布的动态特性。

发明内容

本发明的目的是提供一种霍尔推力器离子能量分布测量方法及系统，通过霍尔推力器低频振荡的调制作用，获得霍尔推力器的离子能量随低频振荡周期的动态变化特性曲线，为研究放电振荡与羽流区离子能量分布的关系了有效的技术途径。

为实现上述目的，本发明提供了如下方案：

一种霍尔推力器离子能量分布测量方法，所述方法包括：

获取多栅探针接收的离子电流随振荡周期变化的曲线L₁，L₂，L₃,......,L_n；所述多栅探针位于霍尔推力器通道上，所述多栅探针中心轴线上与所述通道中心轴线重合；所述离子电流随振荡周期变化的曲线L₁为所述多栅探针制止电压为U₁时得到的曲线；所述离子电流随振荡周期变化的曲线L₂为所述多栅探针制止电压为U₂时得到的曲线；所述离子电流随振荡周期变化的曲线L₃为所述多栅探针制止电压为U₃时得到的曲线；所述离子电流随振荡周期变化的曲线L_n为所述多栅探针制止电压为U_n时得到的曲线；

根据离子电流随振荡周期变化的曲线L₁，L₂，L₃，......L_n，分别获取振荡周期内时刻T₁的电压值U₁，U₂，U₃,......,U_n对应的离子电流值I₁，I₂，I₃,......,I_n；

根据时刻T₁的电压值U₁，U₂，U₃,......,U_n与对应的离子电流值I₁，I₂，I₃,......,I_n，获取时刻T₁的离子平均速度＜v(t₁)＞；

根据所述时刻T₁的离子平均速度＜v(t₁)＞，获得时刻T₁霍尔推力器的离子能量；

以设定时间间隔分别获取所述振荡周期内时刻T₂，T₃,......,T_n的所述霍尔推力器的离子能量；

根据所述振荡周期内时刻T₂，T₃,......,T_n的所述霍尔推力器的离子能量，获得振荡周期内所述霍尔推力器的离子能量分布曲线。

可选的，所述根据时刻T₁的电压值U₁，U₂，U₃,......,U_n与对应的电流值I₁，I₂，I₃,......,I_n，获取时刻T₁的离子平均速度＜v(t₁)＞，具体包括：

根据时刻T₁的电压值U₁，U₂，U₃,......,U_n与对应的离子电流值I₁，I₂，I₃,......,I_n得到时刻T₁电流值与电压值的函数关系；

根据时刻T₁所述电压值与电流值的函数关系获得时刻T₁的离子速度分布函数g(v(t₁))；

获得T₁时刻的离子平均速度＜v(t₁)＞。

可选的，所述根据所述时刻T₁的离子平均速度＜v(t₁)＞，获得时刻T₁霍尔推力器的离子能量，具体包括：

利用公式计算时刻T₁霍尔推力器的离子能量E₁，其中m为离子质量，＜v(t₁)＞为离子平均速度。

可选的，所述根据时刻T₁所述电压值与电流值的函数关系获得时刻T₁的离子速度分布函数g(v(t₁))，具体包括：

获取所述电流值与所述离子速度分布函数的关系表达式为：其中e为单位电荷电量，s为多栅探头接受面积，m为离子质量，I(t₁)表示T₁时刻电压为U(t₁)时的电流值；

设定时刻T₁电压值与电流值的函数关系I(t₁)＝h(U(t₁))，将所述电流值与所述离子速度分布函数的关系表达式变形，获得变形后表达式为：即

根据离子速度与电压值关系表达式获取时刻T₁的离子速度分布函数g(v(t₁))为：

一种霍尔推力器离子能量分布测量系统，所述系统包括：

曲线获取模块，用于获取多栅探针接收的离子电流随振荡周期变化的曲线L₁，L₂，L₃,......,L_n；所述多栅探针位于霍尔推力器通道上，所述多栅探针中心轴线上与所述通道中心轴线重合；所述离子电流随振荡周期变化的曲线L₁为所述多栅探针制止电压为U₁时得到的曲线；所述离子电流随振荡周期变化的曲线L₂为所述多栅探针制止电压为U₂时得到的曲线；所述离子电流随振荡周期变化的曲线L₃为所述多栅探针制止电压为U₃时得到的曲线；所述离子电流随振荡周期变化的曲线L_n为所述多栅探针制止电压为U_n时得到的曲线；

电流获取模块，用于根据离子电流随振荡周期变化的曲线L₁，L₂，L₃，......L_n，分别获取振荡周期内时刻T₁的电压值U₁，U₂，U₃,......,U_n对应的离子电流值I₁，I₂，I₃,......,I_n；

平均速度获取模块，用于根据时刻T₁的电压值U₁，U₂，U₃,......,U_n与对应的离子电流值I₁，I₂，I₃,......,I_n，获取时刻T₁的离子平均速度＜v(t₁)＞；

离子能量获取模块，用于根据所述时刻T₁的离子平均速度＜v(t₁)＞，获得时刻T₁的霍尔推力器的离子能量，还用于以设定时间间隔分别获取所述振荡周期内时刻T₂，T₃,......,T_n的所述霍尔推力器的离子能量；

离子能量分布曲线获取模块，用于根据所述振荡周期内时刻T₂，T₃,......,T_n的所述霍尔推力器的离子能量，获得振荡周期内所述霍尔推力器的离子能量分布曲线。

可选的，所述平均速度获取模块具体包括：

函数关系获取单元，用于根据时刻T₁的电压值U₁，U₂，U₃,......,U_n与对应的离子电流值I₁，I₂，I₃,......,I_n得到时刻T₁电流值与电压值的函数关系；

速度分布函数获取单元，用于根据时刻T₁所述电压值与电流值的函数关系获得时刻T₁的离子速度分布函数g(v(t₁))；

平均速度获取单元，用于利用公式获得T₁时刻的离子平均速度＜v(t₁)＞。

可选的，所述离子能量获取模块具体包括：

离子能量获取单元，用于利用公式计算时刻T₁霍尔推力器的离子能量E₁，其中m为离子质量，＜v(t₁)＞为离子平均速度。

可选的，所述速度分布函数获取单元具体包括：

电流值与速度分布函数获取子单元，用于获取所述电流值与所述离子速度分布函数的关系表达式为：其中e为单位电荷电量，s为多栅探头接受面积，m为离子质量，I(t₁)表示T₁时刻电压为U(t₁)时的电流值；

函数关系表达式变形子单元，用于根据设定的时刻T₁电压值与电流值的函数关系I(t₁)＝h(U(t₁))，将所述电流值与所述离子速度分布函数的关系表达式变形，获得变形后的表达式为：即

速度函数获取子单元，用于根据离子速度与电压值关系表达式获得：

根据本发明提供的具体实施例，本发明公开了以下技术效果：

通过霍尔推力器低频振荡的调制作用，获得霍尔推力器的离子能量随低频振荡周期的动态变化特性曲线，实现了测量羽流区离子能量分布的动态特性，解决了现有的测量方法无法获得几十微秒量级羽流区离子能量分布随时间动态变化的难题，为研究霍尔推力器羽流区离子能量分布变化提供有效技术途径。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明霍尔推力器离子能量分布测量方法流程图；

图2为本发明霍尔推力器离子能量分布测量系统结构图；

图3为本发明霍尔推力器离子能量分布测量系统安装装置图；

图4为本发明霍尔推力器离子能量分布测量系统多栅探针测量电路图；

图5为本发明霍尔推力器离子能量分布测量方法及系统实施例1的多栅探针制止电压下的放电电流曲线和所接受到的离子电流曲线；

图6为本发明霍尔推力器离子能量分布测量方法及系统实施例1的离子能量分布曲线。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

图1为本发明霍尔推力器离子能量分布测量方法流程图。如图1所示，所述方法包括：

步骤101：获取离子电流的变化曲线。具体是指获取多栅探针接收的离子电流随振荡周期变化的曲线L₁，L₂，L₃,......,L_n；曲线L₁为多栅探针制止电压为U₁时得到的曲线，用函数I₁(t)＝f(U₁(t))表示；曲线L₂为多栅探针制止电压为U₂时得到的曲线，用函数I₂(t)＝f(U₂(t))表示；曲线L₃为多栅探针制止电压为U₃时得到的曲线，用函数I₃(t)＝f(U₃(t))表示；曲线L_n为多栅探针制止电压为U_n时得到的曲线，用函数I_n(t)＝f(U_n(t))表示。

步骤102：获取T₁时刻离子电流值，具体是根据步骤101得到的离子电流随振荡周期变化的曲线L₁，L₂，L₃，......L_n，分别获取振荡周期内同一时刻T₁的电压值U₁，U₂，U₃,......,U_n对应的离子电流值I₁，I₂，I₃,......,I_n。

步骤103：获取T₁时刻离子平均速度，具体是根据步骤102获得的同一时刻T₁的电压值U₁，U₂，U₃,......,U_n和对应的离子电流值I₁，I₂，I₃,......,I_n，得到时刻T₁的电流值与电压值函数关系表达式I(t₁)＝h(U(t₁))；将电流值与离子速度分布函数的关系表达式变形，获得变形后表达式为：即其中e为单位电荷电量，s为多栅探头接受面积，m为离子质量，I(t₁)表示T₁时刻电压为U(t₁)时的电流值；根据离子速度与电压值关系表达式获取时刻T₁的离子速度分布函数g(v(t₁))为：利用公式获得T₁时刻的离子平均速度＜v(t₁)＞，将代入后，得到

步骤104：获取T₁时刻的霍尔推力器的离子能量。具体为利用公式计算时刻T₁霍尔推力器的离子能量E₁，其中m为离子质量，＜v(t₁)＞为离子平均速度。

步骤105：获取振荡周期内其他时刻的霍尔推力器的离子能量。具体为按照设定时间间隔(例如，时间间隔为2微秒)分别获取振荡周期内时刻T₂，T₃,......,T_n的电压值U₁，U₂，U₃,......,U_n对应的电流值，得到T₂，T₃,......,T_n时刻的电压值与电流值函数关系表达式，参照步骤103、104依次得到T₂，T₃,......,T_n时刻的霍尔推力器的离子能量。

步骤106：获取离子能量分布曲线。具体为根据步骤105中得到的T₂，T₃,......,T_n时刻的霍尔推力器的离子能量，生成离子能量与时间的关系分布曲线。

图2为本发明霍尔推力器离子能量分布测量系统结构图。如图2所示，系统包括：曲线获取模块201，用于获取多栅探针接收的离子电流随振荡周期变化的曲线L₁，L₂，L₃,......,L_n；

电流获取模块202，用于根据离子电流随振荡周期变化的曲线L₁，L₂，L₃，......L_n，分别获取振荡周期内时刻T₁的电压值U₁，U₂，U₃,......,U_n对应的电流值I₁，I₂，I₃,......,I_n；

平均速度获取模块203，用于根据时刻T₁的电压值U₁，U₂，U₃,......,U_n与对应的电流值I₁，I₂，I₃,......,I_n，获取时刻T₁的离子平均速度＜v(t₁)＞；

离子能量获取模块205，用于根据所述时刻T₁的离子平均速度＜v(t₁)＞，获得时刻T₁的霍尔推力器的离子能量，还用于以设定时间间隔分别获取振荡周期内时刻T₂，T₃,......,T_n的霍尔推力器的离子能量；

离子能量分布曲线获取模块206，用于根据振荡周期内时刻T₂，T₃,......,T_n的霍尔推力器的离子能量，确定离子能量随振荡周期变化的关系分布曲线。

图3为本发明霍尔推力器离子能量分布测量系统安装装置图。如图3所示，包括多栅探针301，霍尔推力器302，羽流区303，通道中心轴线304，真空罐305。进行霍尔推力器离子能量测量时，将多栅探针301放置在通道上，多栅探针301中心轴线与通道中心轴线304重合，高度与霍尔推力器通道中心平齐，距离通道出口截面有一定距离，例如距离可以为75cm。

图4为本发明霍尔推力器离子能量分布测量系统多栅探针测量电路图。如图4所示，401表示多栅探针接收到的离子电流，402表示多栅探针的制止电压。

本发明霍尔推力器离子能量分布测量系统及方法实施例1：

将多栅探针放置在通道中心轴线上，距离通道出口截面75cm，高度与推力器通道中心平齐。

调整多栅探针制止电压为U₁时，测得多栅探针接收的离子电流电流随振荡周期变化的曲线L₁和放电电流曲线；调整多栅探针制止电压为U₂时，测得多栅探针接收的离子电流电流随振荡周期变化的曲线L₂和放电电流曲线，直至将所选的n个多栅探针制止电压下的离子电流随振荡周期变化的曲线和放电电流曲线测量完毕。由于低频振荡周期几乎不变，在每个多栅探针制止电压下获得的放电电流的波动曲线相似，因此设为n个多栅探针制止电压下的放电电流曲线相同，设为L_d，如图5所示。

在L_d曲线上选择一个振荡周期，在该振荡周期选择一个时刻，依次获得曲线L₁，L₂，......L_n上该时刻的电流值，根据该时刻不同曲线获得的电流值与电压值生成该时刻的电流与电压关系函数，该时刻的电流与电压关系函数计算出该时刻的羽流区离子平均速度，进而利用公式计算时刻T₁霍尔推力器的离子能量E₁，其中m为离子质量，＜v(t₁)＞为离子平均速度。

再在L_d曲线上每隔2微秒取一个时刻，并进行与上述一样的离子能量获得方法，直至取完这一个振荡周期内的所有时刻，便可得到i个时刻的霍尔推力器的离子能量，将得到的i个离子能量值根据时间拟合，便可以得到在振荡周期内离子能量的分布曲线。如图6所示，图6为本发明霍尔推力器离子能量分布测量方法及系统实施例1的离子能量分布曲线。

本发明霍尔推力器离子能量分布测量系统及方法实施例2：该实施例与实施例1不同的地方是：在放电电流曲线L_d上以0.15微秒为间隔选取一个时刻值，其他的测量过程与实施例1相同。

本发明霍尔推力器离子能量分布测量系统及方法实施例3：该实施例与实施例1不同的地方是：多栅探针距离通道出口截面50cm，其他的测量过程与实施例1相同。

本发明霍尔推力器离子能量分布测量系统及方法实施例4：该实施例与实施例1不同的地方是：多栅探针距离通道出口截面100cm，其他的测量过程与实施例1相同。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。

本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。