CN104406761B - 霍尔推力器低频振荡时间尺度内羽流发散角测量方法 - Google Patents
霍尔推力器低频振荡时间尺度内羽流发散角测量方法 Download PDFInfo
- Publication number
- CN104406761B CN104406761B CN201410578264.5A CN201410578264A CN104406761B CN 104406761 B CN104406761 B CN 104406761B CN 201410578264 A CN201410578264 A CN 201410578264A CN 104406761 B CN104406761 B CN 104406761B
- 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.)
- Active
Links
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 (3)
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＝0.12，并判断下一时刻的时间值是否大于预设周期的时间值，如果判断结果为是，则执行步骤五；如果判断结果为否，则返回执行步骤三；
步骤五、将获得的所有时刻的羽流发散角进行拟合，从而完成霍尔推力器低频振荡时间尺度内羽流发散角的测量。
2.根据权利要求1所述的霍尔推力器低频振荡时间尺度内羽流发散角测量方法，其特征在于相邻两个测量点的距离为0.2cm。
3.根据权利要求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 CN104406761A (zh) | 2015-03-11 |
| CN104406761B true 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) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105574883A (zh) * | 2015-12-30 | 2016-05-11 | 哈尔滨工业大学 | 一种图像法评估霍尔推力器羽流发散角的方法 |
| CN106772543B (zh) * | 2016-12-07 | 2018-11-02 | 哈尔滨工业大学 | 霍尔推力器羽流发散角测量修正方法 |
| CN106802125B (zh) * | 2017-01-24 | 2020-02-11 | 上海空间推进研究所 | 航天器用电推力器推力矢量偏心的确定方法 |
| CN109342861B (zh) * | 2018-12-10 | 2020-12-22 | 哈尔滨工业大学 | 一种霍尔推力器周向辐条特性测量方法及系统 |
| CN115682919B (zh) * | 2022-10-28 | 2023-08-11 | 哈尔滨工业大学 | 一种基于曲面积分的霍尔推力器推力矢量偏心计算方法 |
| CN115808315B (zh) * | 2022-11-30 | 2024-01-23 | 兰州空间技术物理研究所 | 一种基于霍尔推力器低频振荡性能评估的装置及方法 |
| CN117451238B (zh) * | 2023-12-19 | 2024-04-09 | 哈尔滨工业大学 | 基于神经网络的推进器推力波动在轨光学检测方法及装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6448721B2 (en) * | 2000-04-14 | 2002-09-10 | General Plasma Technologies Llc | Cylindrical geometry hall thruster |
| US6777862B2 (en) * | 2000-04-14 | 2004-08-17 | General Plasma Technologies Llc | Segmented electrode hall thruster with reduced plume |
| CN101334352A (zh) * | 2008-07-30 | 2008-12-31 | 哈尔滨工业大学 | 霍尔推进器寿命的估算方法 |
| CN101969737A (zh) * | 2010-08-27 | 2011-02-09 | 哈尔滨工业大学 | 一种实现霍尔推力器等离子体束聚焦的外回路控制装置及方法 |
-
2014
- 2014-10-25 CN CN201410578264.5A patent/CN104406761B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6448721B2 (en) * | 2000-04-14 | 2002-09-10 | General Plasma Technologies Llc | Cylindrical geometry hall thruster |
| US6777862B2 (en) * | 2000-04-14 | 2004-08-17 | General Plasma Technologies Llc | Segmented electrode hall thruster with reduced plume |
| CN101334352A (zh) * | 2008-07-30 | 2008-12-31 | 哈尔滨工业大学 | 霍尔推进器寿命的估算方法 |
| CN101969737A (zh) * | 2010-08-27 | 2011-02-09 | 哈尔滨工业大学 | 一种实现霍尔推力器等离子体束聚焦的外回路控制装置及方法 |
Non-Patent Citations (1)
| Title |
|---|
| 霍尔推力器中励磁与放电耦合振荡研究;魏立秋;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20110615(第6（2011）期);第19-28页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104406761A (zh) | 2015-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104406761B (zh) | 霍尔推力器低频振荡时间尺度内羽流发散角测量方法 | |
| CN102540127B (zh) | 低轨道航天器空间电位探测器校准平台 | |
| MY183579A (en) | Wafer producing method and processing feed direction detecting method | |
| Chattopadhyay et al. | Analysis of stator current of induction motor used in transport system at single phasing by measuring phase angle, symmetrical components, Skewness, Kurtosis and harmonic distortion in Park plane | |
| Wang et al. | Convex combination based target localization with noisy angle of arrival measurements | |
| Jin et al. | Fault location method for distribution lines with distributed generators based on a novel hybrid BPSOGA | |
| CN102811419A (zh) | 一种基于迭代的最小二乘定位方法 | |
| CN104796116A (zh) | 一种基于参数拟合的冲击波形峰值测量方法 | |
| Kong et al. | A new calibration model for pointing a radio telescope that considers nonlinear errors in the azimuth axis | |
| CN103954789B (zh) | 离子速度分布函数瞬时测量装置及方法 | |
| Cai et al. | Lightning electric‐field waveforms associated with transmission‐line faults | |
| CN103414184A (zh) | 一种频率变化情况下序分量的计算方法 | |
| Dhara et al. | Methods for localization of partial discharge sources within air insulated electrical substation | |
| Tong et al. | High current field test of impulse transient characteristics of substation grounding grid | |
| Gasch et al. | Uncertainty evaluation for the impact of measurement accuracy on power quality parameters | |
| CN103245876A (zh) | 利用电压降实部特性实现线路相间故障单端测距方法 | |
| CN103234457B (zh) | 基于数字成像的多光束位移测量方法 | |
| CN108196168A (zh) | 一种局放脉冲检测方法、系统、终端及可读存储介质 | |
| CN106018943B (zh) | 一种接地极电位分布测量方法 | |
| CN107834678A (zh) | 一元线性回归测频算法在直流充电机中的应用 | |
| KR101446083B1 (ko) | ExB 탐침의 신호 처리 방법 | |
| CN103293440A (zh) | 利用序分量实现线路单相接地故障单端测距方法 | |
| Zhao et al. | Simulation of downburst in a multiple fan wind tunnel and research on its load on high-rise structure by wind tunnel experiment | |
| Franci et al. | Bar-mode instability suppression in magnetized relativistic stars | |
| RU150629U1 (ru) | Устройство для управления разгрузкой турбин энергоблоков электростанций |
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 |