CN104406761A - 霍尔推力器低频振荡时间尺度内羽流发散角测量方法 - Google Patents
霍尔推力器低频振荡时间尺度内羽流发散角测量方法 Download PDFInfo
- Publication number
- CN104406761A CN104406761A CN201410578264.5A CN201410578264A CN104406761A CN 104406761 A CN104406761 A CN 104406761A CN 201410578264 A CN201410578264 A CN 201410578264A CN 104406761 A CN104406761 A CN 104406761A
- Authority
- CN
- China
- Prior art keywords
- low
- frequency oscillation
- hall thruster
- time
- plume
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Plasma Technology (AREA)
Abstract
霍尔推力器低频振荡时间尺度内羽流发散角测量方法,涉及等离子体推进领域。它是为了获取霍尔推力器低频振荡时间尺度内羽流发散角的动态特性。它通过探针测量每个测量点的低频振荡电流波形图和离子电流波形,选取一个时刻的低频振荡电流值并从离子电流曲线上找到与之对应的点,将采集的点进行拟合后得到沿径向的离子分布曲线,并计算得到该时刻的羽流发散角,以此类推,可以得到羽流发散角随时间的变化曲线。本发明实现了测量羽流发散角的动态特性,获得羽流发散角在低频振荡时间尺度内的变化曲线,为研究霍尔推力器羽流发散角变化提供有效技术途径。本发明适用于霍尔推力器低频振荡时间尺度内羽流发散角测量。
Description
技术领域
本发明涉及等离子体推进领域。
背景技术
近年来,霍尔电推进系统因其效率高、比冲适中和长寿命等特点,在卫星位置保持和轨道提升中获得成功应用,已成为航天推进领域的重要研究方向。霍尔推力器羽流为高速稀薄等离子体且较发散,对航天器表面、太阳能电池和通信设备等存在污染,因此对霍尔推力器羽流发散角的评估测量是霍尔推力器研究的重要内容之一。霍尔推力器的羽流发散角是指其喷出等离子体在羽流空间的束流扩张角度,因此通常用它来反映等离子体束的集中程度和聚焦特性,也是用来评价推力器对飞行器表面相互作用的一个重要参数。
霍尔推力器的低频振荡特性会使等离子体的发散角增大,增大羽流污染,而目前的测量方法只能测量推力器在稳态条件下的羽流发散角,在低频振荡时间尺度内(通常为十几个μs),很难获得羽流发散角的动态特性。
发明内容
本发明是为了解决现有的测量方法无法获得几十μs量级羽流发散角随时间动态变化的问题,从而提供一种霍尔推力器低频振荡时间尺度内羽流发散角测量方法。
霍尔推力器低频振荡时间尺度内羽流发散角测量方法,它由以下步骤实现:
在霍尔推力器的出口处设置N个测量点;N为大于2的整数;
该N个测量点位于一点直线上,且该直线沿霍尔推力器径向分布;
N个测量点等间隔设置;设r的初始值为1;
步骤一、将探针安装在第r个测量点处,采用所述探针探测第r个测量点处霍尔推力器的低频振荡电流和离子电流;获得第r个测量点处霍尔推力器的低频振荡电流随时间变化的曲线,以及离子电流随时间变化的曲线;
步骤二、将r的值加1,并判断r的值是否大于N,如果判断结果为是,则执行步骤三;如果判断结果为否,则返回执行步骤一;完成步骤二后,获得霍尔推力器径向方向所有测量点的低频振荡电流和离子电流,共获得N条低频振荡电流随时间变化的曲线,以及N条离子电流随时间变化的曲线;
步骤三、设每个时间点处N条低频振荡电流的波形相同,在步骤二获得的霍尔推力器的低频振荡电流随时间变化的曲线上选取第t时间点对应的低频振荡电流值,并在N条离子电流随时间变化的曲线上读取在该时间点对应的离子电流值,将读取的N个离子电流值进行拟合,获得该时间点的离子电流沿径向的分布曲线,进而获得该时间的羽流发散角;t的初始值为1;
步骤四、将t加Mμs作为下一时刻,M为正数,并判断下一时刻的时间值是否大于预设周期的时间值,如果判断结果为是,则执行步骤五;如果判断结果为否,则返回执行步骤三;
步骤五、将获得的所有时刻的羽流发散角进行拟合,从而完成霍尔推力器低频振荡时间尺度内羽流发散角的测量。
与现有技术相比,本发明的有益效果是:本发明解决了现有的测量方法无法获得几十μs量级羽流发散角随时间动态变化的问题,本发明通过霍尔推力器低频振荡的调制作用,实现了低频振荡时间尺度内羽流发散角的测量,避免了短时间尺度内羽流发散角测量所需的高速移动设备,为研究放电振荡与羽流发散角的关系提供了有效的技术途径。
附图说明
图1是为每个测量位置点的低频振荡电流曲线和离子电流曲线仿真示意图;
图2是沿径向的离子分布曲线仿真示意图;
图3是羽流发散角测量原理示意图。
具体实施方式
具体实施方式一、霍尔推力器低频振荡时间尺度内羽流发散角测量方法,它由以下步骤实现:
在霍尔推力器的出口处设置N个测量点;N为大于2的整数;
该N个测量点位于一点直线上,且该直线沿霍尔推力器径向分布;
N个测量点等间隔设置;设r的初始值为1;
步骤一、将探针安装在第r个测量点处,采用所述探针探测第r个测量点处霍尔推力器的低频振荡电流和离子电流;获得第r个测量点处霍尔推力器的低频振荡电流随时间变化的曲线,以及离子电流随时间变化的曲线;
步骤二、将r的值加1,并判断r的值是否大于N,如果判断结果为是,则执行步骤三;如果判断结果为否,则返回执行步骤一;完成步骤二后,获得霍尔推力器径向方向所有测量点的低频振荡电流和离子电流,共获得N条低频振荡电流随时间变化的曲线,以及N条离子电流随时间变化的曲线;
步骤三、设每个时间点处N条低频振荡电流的波形相同,在步骤二获得的霍尔推力器的低频振荡电流随时间变化的曲线上选取第t时间点对应的低频振荡电流值,并在N条离子电流随时间变化的曲线上读取在该时间点对应的离子电流值,将读取的N个离子电流值进行拟合,获得该时间点的离子电流沿径向的分布曲线,进而获得该时间的羽流发散角;t的初始值为1;
步骤四、将t加Mμs作为下一时刻,M为正数,并判断下一时刻的时间值是否大于预设周期的时间值,如果判断结果为是,则执行步骤五;如果判断结果为否,则返回执行步骤三;
步骤五、将获得的所有时刻的羽流发散角进行拟合,从而完成霍尔推力器低频振荡时间尺度内羽流发散角的测量。
具体实施方式二、本具体实施方式与具体实施方式一所述的霍尔推力器低频振荡时间尺度内羽流发散角测量方法的在于,相邻两个测量点的距离为0.2cm。
具体实施方式三、本具体实施方式与具体实施方式一所述的霍尔推力器低频振荡时间尺度内羽流发散角测量方法的在于,M=0.12。
具体实施方式四、本具体实施方式与具体实施方式一所述的霍尔推力器低频振荡时间尺度内羽流发散角测量方法的在于,探针端面距离霍尔推力器出口平面为3cm,探针轴线距离霍尔推力器中轴线10cm。
实验过程中可参考如下的公式:
式中:SP表示探针的面积,Jck+d表示定向电流和附加的杂乱电流,Jck表示杂乱电流。
测得电流密度后jd,电流为:
式中:r表示推力器径向积分面积的半径;
根据电流密度的测量结果,通过差值计算每隔0.5cm位置的电流j′i(插值),然后由:
式中:
当:j′i→0,即可停止计算,找到Ii/I≥95%的位置,I为总的离子电流,记此位置的半径为R,如图3所示,因此羽流发散角为:
式中:L为探针端面距推力器出口平面的轴向距离,r′为推力器陶瓷通道外半径。
Claims (4)
1.霍尔推力器低频振荡时间尺度内羽流发散角测量方法,其特征是:它由以下步骤实现:
在霍尔推力器的出口处设置N个测量点;N为大于2的整数;
该N个测量点位于一点直线上,且该直线沿霍尔推力器径向分布;
N个测量点等间隔设置;设r的初始值为1;
步骤一、将探针安装在第r个测量点处,采用所述探针探测第r个测量点处霍尔推力器的低频振荡电流和离子电流;获得第r个测量点处霍尔推力器的低频振荡电流随时间变化的曲线,以及离子电流随时间变化的曲线;
步骤二、将r的值加1,并判断r的值是否大于N,如果判断结果为是,则执行步骤三;如果判断结果为否,则返回执行步骤一;完成步骤二后,获得霍尔推力器径向方向所有测量点的低频振荡电流和离子电流,共获得N条低频振荡电流随时间变化的曲线,以及N条离子电流随时间变化的曲线;
步骤三、设每个时间点处N条低频振荡电流的波形相同,在步骤二获得的霍尔推力器的低频振荡电流随时间变化的曲线上选取第t时间点对应的低频振荡电流值,并在N条离子电流随时间变化的曲线上读取在该时间点对应的离子电流值,将读取的N个离子电流值进行拟合,获得该时间点的离子电流沿径向的分布曲线,进而获得该时间的羽流发散角;t的初始值为1;
步骤四、将t加Mμs作为下一时刻,M为正数,并判断下一时刻的时间值是否大于预设周期的时间值,如果判断结果为是,则执行步骤五;如果判断结果为否,则返回执行步骤三;
步骤五、将获得的所有时刻的羽流发散角进行拟合,从而完成霍尔推力器低频振荡时间尺度内羽流发散角的测量。
2.根据权利要求1所述的霍尔推力器低频振荡时间尺度内羽流发散角测量方法,其特征在于相邻两个测量点的距离为0.2cm。
3.根据权利要求1所述的霍尔推力器低频振荡时间尺度内羽流发散角测量方法,其特征在于M=0.12。
4.根据权利要求1所述的霍尔推力器低频振荡时间尺度内羽流发散角测量方法,其特征在于探针端面距离霍尔推力器出口平面为3cm,探针轴线距离霍尔推力器中轴线10cm。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410578264.5A CN104406761B (zh) | 2014-10-25 | 2014-10-25 | 霍尔推力器低频振荡时间尺度内羽流发散角测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410578264.5A CN104406761B (zh) | 2014-10-25 | 2014-10-25 | 霍尔推力器低频振荡时间尺度内羽流发散角测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104406761A true CN104406761A (zh) | 2015-03-11 |
CN104406761B CN104406761B (zh) | 2017-05-03 |
Family
ID=52644408
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410578264.5A Active CN104406761B (zh) | 2014-10-25 | 2014-10-25 | 霍尔推力器低频振荡时间尺度内羽流发散角测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104406761B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105574883A (zh) * | 2015-12-30 | 2016-05-11 | 哈尔滨工业大学 | 一种图像法评估霍尔推力器羽流发散角的方法 |
CN106772543A (zh) * | 2016-12-07 | 2017-05-31 | 哈尔滨工业大学 | 霍尔推力器羽流发散角测量修正方法 |
CN106802125A (zh) * | 2017-01-24 | 2017-06-06 | 上海空间推进研究所 | 航天器用电推力器推力矢量偏心的确定方法 |
CN109342861A (zh) * | 2018-12-10 | 2019-02-15 | 哈尔滨工业大学 | 一种霍尔推力器周向辐条特性测量方法及系统 |
CN115682919A (zh) * | 2022-10-28 | 2023-02-03 | 哈尔滨工业大学 | 一种基于曲面积分的霍尔推力器推力矢量偏心计算方法 |
CN115808315A (zh) * | 2022-11-30 | 2023-03-17 | 兰州空间技术物理研究所 | 一种基于霍尔推力器低频振荡性能评估的装置及方法 |
CN117451238A (zh) * | 2023-12-19 | 2024-01-26 | 哈尔滨工业大学 | 基于神经网络的推进器推力波动在轨光学检测方法及装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6777862B2 (en) * | 2000-04-14 | 2004-08-17 | General Plasma Technologies Llc | Segmented electrode hall thruster with reduced plume |
US6448721B2 (en) * | 2000-04-14 | 2002-09-10 | General Plasma Technologies Llc | Cylindrical geometry hall thruster |
CN101334352B (zh) * | 2008-07-30 | 2011-06-22 | 哈尔滨工业大学 | 霍尔推进器寿命的估算方法 |
CN101969737B (zh) * | 2010-08-27 | 2012-06-20 | 哈尔滨工业大学 | 一种实现霍尔推力器等离子体束聚焦的外回路控制装置及方法 |
-
2014
- 2014-10-25 CN CN201410578264.5A patent/CN104406761B/zh active Active
Non-Patent Citations (1)
Title |
---|
魏立秋: "霍尔推力器中励磁与放电耦合振荡研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105574883A (zh) * | 2015-12-30 | 2016-05-11 | 哈尔滨工业大学 | 一种图像法评估霍尔推力器羽流发散角的方法 |
CN106772543A (zh) * | 2016-12-07 | 2017-05-31 | 哈尔滨工业大学 | 霍尔推力器羽流发散角测量修正方法 |
CN106802125A (zh) * | 2017-01-24 | 2017-06-06 | 上海空间推进研究所 | 航天器用电推力器推力矢量偏心的确定方法 |
CN106802125B (zh) * | 2017-01-24 | 2020-02-11 | 上海空间推进研究所 | 航天器用电推力器推力矢量偏心的确定方法 |
CN109342861A (zh) * | 2018-12-10 | 2019-02-15 | 哈尔滨工业大学 | 一种霍尔推力器周向辐条特性测量方法及系统 |
CN109342861B (zh) * | 2018-12-10 | 2020-12-22 | 哈尔滨工业大学 | 一种霍尔推力器周向辐条特性测量方法及系统 |
CN115682919A (zh) * | 2022-10-28 | 2023-02-03 | 哈尔滨工业大学 | 一种基于曲面积分的霍尔推力器推力矢量偏心计算方法 |
CN115682919B (zh) * | 2022-10-28 | 2023-08-11 | 哈尔滨工业大学 | 一种基于曲面积分的霍尔推力器推力矢量偏心计算方法 |
CN115808315A (zh) * | 2022-11-30 | 2023-03-17 | 兰州空间技术物理研究所 | 一种基于霍尔推力器低频振荡性能评估的装置及方法 |
CN115808315B (zh) * | 2022-11-30 | 2024-01-23 | 兰州空间技术物理研究所 | 一种基于霍尔推力器低频振荡性能评估的装置及方法 |
CN117451238A (zh) * | 2023-12-19 | 2024-01-26 | 哈尔滨工业大学 | 基于神经网络的推进器推力波动在轨光学检测方法及装置 |
CN117451238B (zh) * | 2023-12-19 | 2024-04-09 | 哈尔滨工业大学 | 基于神经网络的推进器推力波动在轨光学检测方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104406761B (zh) | 2017-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104406761A (zh) | 霍尔推力器低频振荡时间尺度内羽流发散角测量方法 | |
CN103870654A (zh) | 基于并行矩量法与物理光学混合的电磁散射仿真方法 | |
CN103941090B (zh) | 基于谱线能量插值的谐波测量方法 | |
CN104796116A (zh) | 一种基于参数拟合的冲击波形峰值测量方法 | |
CN104597346A (zh) | 有界波环境下电子设备电磁环境效应实验方法 | |
CN103712606B (zh) | 一种海冰监测系统及其工作方法 | |
CN103440497A (zh) | 一种gis绝缘缺陷局部放电图谱模式识别方法 | |
Zhao et al. | Criterion to evaluate power system online transient stability based on adjoint system energy function | |
CN102157938A (zh) | 电力系统电压稳定薄弱节点在线识别方法 | |
CN104166793A (zh) | 一种基于二次曲线拟合的冲击波形峰值测量方法 | |
Jin et al. | Fault location method for distribution lines with distributed generators based on a novel hybrid BPSOGA | |
CN105223481A (zh) | 基于差值能量函数的局部放电特高频信号起始时刻确定方法 | |
CN103543026A (zh) | 基于振动传递率函数和支持向量机的结构损伤识别方法 | |
CN104101869B (zh) | 一种极坐标下的地波雷达运动目标仿真建模方法 | |
Shen et al. | Two‐step locating method for aircraft wake vortices based on Gabor filter and velocity range distribution | |
CN103530441A (zh) | 利用三维激光扫描技术对汽轮机组进行三维建模的方法 | |
RU2532760C1 (ru) | Способ определения места повреждения разветвленной линии электропередачи | |
CN112881812B (zh) | 基于机器学习编码的全闪实时定位方法及装置 | |
CN103592584A (zh) | 一种基于电磁波检测的变电站放电源定位方法 | |
CN103675777B (zh) | 基于拟合法的机载雷达杂波模拟方法及装置 | |
RU2632583C2 (ru) | Способ определения расстояния до места повреждения на линии электропередачи | |
CN104659767A (zh) | 一种利用序列重叠差分的三角形环网线路单元保护方法 | |
Tang et al. | Simulation of the Doppler echoes from wind turbine based on scattered electric field calculation | |
CN103559332B (zh) | 一种一阶螺旋杆频散特性的提取方法 | |
CN103278128B (zh) | 一种快速的kdp晶体最佳匹配角精确测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Wei Liqiu Inventor after: Ding Yongjie Inventor after: Yu Daren Inventor after: Li Jing Inventor after: Han Liang Inventor after: Yang Ziyi Inventor before: Li Jing Inventor before: Wei Liqiu Inventor before: Han Liang Inventor before: Yang Ziyi Inventor before: Yu Daren |
|
COR | Change of bibliographic data | ||
GR01 | Patent grant | ||
GR01 | Patent grant |