《1 前言》
1 前言
考虑空间发展战略的需求, 我国加快了空间开发利用的进度, 相继开始了微小卫星群、深空探测等系列空间计划的论证和实施。微小卫星具有体积小、重量轻、转动惯量小的特点, 用于卫星轨道与姿态控制所需的推力和力矩就小, 要求精度高。在空间开发领域我国已经启动了以“探月工程”为标志的深空探测系列计划。由于深空探测的飞行器飞行距离长, 为了提高有效载荷率, 就需要高比冲的推进器。因此, 从我国空间发展战略的需求来看, 研制高比冲 (1 000 s~3 000 s) 、高效率 (50 %~70 %) 、长寿命 (>10 000 h) 、小推力 (几十毫牛顿) 的稳态等离子体发动机 (SPT) 将具有重要的战略意义。
如图1所示, 在稳态等离子体推进器通道内存在由内外线圈产生的大致沿径向的磁场和阴阳电极之间产生的大致沿轴向的电场。工质氙气注入通道后被由磁场束缚的电子碰撞电离, 产生的离子在电场加速下以大于15 km/s的高速喷出产生反作用力, 为小卫星等飞行器提供动力
一样, 为了提高电离率, 增大等离子体密度, SPT需要设计磁场来约束等离子体。所以, 磁场位形的设计对于控制SPT通道的等离子体运动参数, 进而提高SPT的工作性能和技术水平有着至关重要的作用。纵观SPT的发展历史就会发现每一次SPT通道磁场设计的改进都会带来性能上的大幅度提高, 磁场位形特点也成了每一代SPT的最主要标志特征之一。
《2 SPT磁场设计发展简要回顾》
2 SPT磁场设计发展简要回顾
《2.1 径向磁场的建立标志着SPT概念设计的形成》
2.1 径向磁场的建立标志着SPT概念设计的形成
前苏联从20世纪60年代起开始着手研制SPT, 从最初的概念设计起磁场设计就一直在SPT的设计中起着至关重要的作用。如图1所示, 离子的加速力是由安培力F=IJ×B提供的, 要增大该加速力, 就需要大的加速电流IJ, 但同时放电电流还要尽量小, 以减少焦耳耗能。为了解决这一难题, 第一代SPT发明人A.I.Morozov提出了磁场方向沿环形通道径向, 以与轴向电场配合来产生电子沿周向的闭环霍尔漂移的设计思路。在该径向磁场位形下, 电子周向霍尔漂移运动形成的电流可以很大, 以满足离子加速的安培力需求, 同时, 由于电子的轴向运动被径向磁场约束, 沿轴向的放电电流就可以很小, 减小了焦耳耗能。因此遵循这一设计原则, 在轴对称的结构下, SPT在磁路结构上采用了4个缠在圆柱铁芯上的外线圈和1个内线圈来产生这种径向磁场位形
《2.2 正梯度磁场强度的建立标志着第一代SPT的真正诞生》
2.2 正梯度磁场强度的建立标志着第一代SPT的真正诞生
但是, 在上述正交电磁场位形下, 由等离子物理中的Bohm传导公式, 可预计到SPT通道中由于等离子体振荡而产生很高的轴向放电电流而无法正常工作, 因此, 在SPT研制初期遭到了很多物理学家的反对。后来, 这一问题的解决仍然得益于磁场设计的改进。A.I.Morozov教授经过实验和理论推导证实了设计磁场强度沿通道轴向以正梯度变化就会抑制这种高放电电流
其中:n为沿磁力线的等离子体密度,
H为该磁力线处的磁场强度。
又由等离子物理学中磁流体动力学可得
第一代SPT从前苏联设计成功起已经在100多颗卫星上使用
《2.3 弯曲磁场的提出是第二代SPT诞生的重要标志之一》
2.3 弯曲磁场的提出是第二代SPT诞生的重要标志之一
到了20世纪90年代, 在法国SEP公司资助下, A.I.Morozov和A.I.Bugrova教授又于1996年发布了一种新设计的SPT—ATON, 被称为第二代SPT
ATON发动机弯曲磁场设计带来性能大幅提升的主要原因在于凸向阳极的弯曲磁场会产生出近似凸向阳极形状的弯曲等电势面, 进而产生出指向通道中心区域的径向电场, 并使打向通道壁面的离子向通道中心迁移, 从而减少了壁面的离子溅射腐蚀, 提高了效率, 减小了羽流发散角。但是俄罗斯目前在弯曲磁场的设计中只是提出了凸向阳极的定性要求, 然后通过大量正交实验来优化曲率, 直到等离子射流满足要求为止。可见, 其磁场设计是建立在大量实验基础之上的, 还远没有形成相对比较完整的设计理论体系来指导通道的二维磁场设计。这不仅在工程应用中带来不便, 而且也很难形成磁场控制等离子体平行射流的具有普适性的科学结论, 从而为其他等离子源的设计提供一般化理论指导。所以, 在此基础上, 我们意图提出一套较完整的SPT通道二维磁场设计的理论化体系, 以指导SPT的磁场设计。
表1 SPT技术指标的发展简表
Table 1 The development process of some important parameters about SPT
《表1》
磁场位形的特征 |
效率/% | 比冲/s | 羽流发散角/ (°) |
径向磁场位形 |
45 | ~1 600 | ~45 |
正梯度磁场位形 |
50 | ~1 700 | ~45 |
弯曲磁场位形 |
68 | ~2 000 | <10 |
《3 SPT弯曲磁场设计理论体系的提出和构成》
3 SPT弯曲磁场设计理论体系的提出和构成
《3.1 SPT内离子径向二次流概念的提出及产生原因分析》
3.1 SPT内离子径向二次流概念的提出及产生原因分析
如图1所示, 在理想情况下, 希望离子以轴线方向平行喷射, 这时推进器具有最佳的比冲和效率, 羽流发散角几乎为零。但是在实际推进器中, 发现在沿径向的直磁场位形下, 离子在沿轴向形成射流的同时还存在从通道中心区域向两个壁面沿径向的串流 (见图2) 。在叶轮机械中, 由气体的三元流动理论可知, 气体的流动可以分解为沿叶片通道方向的主流以及相对于主流而言分别沿叶片径向和叶轮周向的二次流运动。主流运动是产生功能转换的主方向运动, 而气流的二次流运动则会带来能量的损失
将式 (2) 沿通道轴向分解, 可得式 (3) 和式 (4) , 这两个方程分别描述了离子的轴向主流运动和径向二次流运动。可见, 矢量方程的可分解性使得我们可以将离子径向二次流运动从主流运动中分离出来单独考虑其控制手段和方法。
从离子径向二次流运动方程 (4) 可以看出, 离子径向二次流运动的产生原因可以从离子的初速度和通道内的电场分布这两方面来分析。
首先, 从离子的径向初速度来源分析。离子的初速度继承于被电离的原子, 如果原子具有打向壁面的宏观定向速度, 那么离子也必然会具有指向壁面的径向分速度。第一代SPT的工质氙气流是通过均布在SPT阳极径向中心处的若干小孔进入放电通道的 (如图1所示) , 所以中性原子气体在通道的径向方向会存在密度差异, 从而会产生径向的定向扩散运动。可见, 在第一代SPT的结构下, 中性气体密度沿径向的不均匀性会使得离子具有径向初速度。其次, 从SPT通道中的电场分布分析, 认为外加电极在SPT通道内只产生轴向电场, 对径向离子的运动不产生影响。径向电场的产生将是等离子体自身自恰产生的一种结果。如果忽略电子的惯性力, 由电子沿径向磁力线的运动方程式 (5) 可以看出电子沿径向的压力梯度会自恰产生出径向电场Er。在认为电子温度沿磁力线恒定的前提下,
电子的径向压力梯度是由电子的径向密度梯度引起的。从电子运动的控制方程而言, 沿着磁力线的电子运动是自由运动, 因此密度差异不是控制方程作用的结果, 这种结果可能来源于2种物理效应。一是电子产生源的影响, 在电子产生的电离区里, 电子密度沿径向就可能是不均匀的, 这样在没有外加力扰动的情况下, 这种不均匀性会一直维持。另外一种情况就是壁面鞘层的影响, 为了满足SPT通道壁面鞘层形成条件, 在预鞘层区域需要有电势降以加速离子向壁面运动, 这个径向电势降的形成需要由电子密度沿径向分布不均匀来产生。例如在无碰撞的近似条件下, 形成加速离子的预鞘层电势需要等离子体的密度降约为60.5 %
《3.2 弯曲磁场控制离子径向二次流的基本原理》
3.2 弯曲磁场控制离子径向二次流的基本原理
离子径向二次流运动会使得高能量离子打向壁面, 不仅引起强烈的壁面溅射腐蚀现象, 降低SPT寿命, 而且高能离子在壁面处复合成中性原子, 造成能量损失。另外喷射出去的离子由于径向二次流运动会使得SPT羽流发散角增大, 通常可以达到30~40° (半角) , 这对推进器的效率和比冲以及星载器件均具有显著影响。在叶轮机械中同样有此问题存在。叶片通道中的气体二次流运动是低能流体, 在端壁处容易产生分离, 增加损失。为了控制该端壁损失, 王仲奇院士首先提出了弯曲叶片的成型技术, 成功控制了端壁低能流体向中心区域迁移, 降低了端壁损失, 提高了叶轮机械的效率
设计弯曲磁场以控制离子径向二次流运动的基本原理是基于磁力线就是等电势面的假设。如果忽略电子的质量惯性力, 电子的运动方程可以写作:
又由
再令热化电势
从式 (9) 可以看出描述离子径向二次流运动的方程 (4) 中的径向电场Er直接受到了轴向磁场分量Bz的控制。因此, 通过设计Bz就可以改变径向电场Er的大小和方向, 从而可以控制离子的径向二次流运动, 减小离子与壁面的溅射腐蚀, 满足近似平行喷射的设计要求。
《3.3 弯曲磁场工程化设计的几个问题》
3.3 弯曲磁场工程化设计的几个问题
如前所述, 目前SPT通道弯曲磁场位形的设计是在磁力线与等电势面重合的假设基础之上通过大量的正交实验来实现的, 这种方法显然大大制约了SPT磁场的工程化设计, 也无益于继续研究SPT性能提高的途径。为了建立起较为完整的弯曲磁场工程化设计体系方法, 作者从几个方面给出建立SPT弯曲磁场设计体系的构想:
首先, 从数学上需要给出描述弯曲磁场位形的特征参量以便于磁场优化。从几何构型的角度来讲, 当提出磁力线沿SPT通道的弯曲要求时, 除了磁场梯度这样一个表征磁场大小的特征参量以外, 还需要从数学上给定描述磁场弯曲程度的特征参量。为便于实际应用, 我们提出表征弯曲磁场的一个基本特征参数弯曲角α。其定义为α=180°-β (如图4所示) 。其中β角是弯曲磁力线在通道内外壁面处切线所形成的夹角。可见, 弯曲角α基本表征了弯曲磁场的弯曲程度, α越大表明磁力线向阳极的弯曲程度越大。可以由弯曲磁场控制离子径向二次流的基本物理过程得到α角与离子径向二次流运动的关系。然后在设计弯曲磁场时, 可以由所要控制的离子径向二次流运动的行为选定弯曲角α, 再由磁场位形自身的特性构造弯曲磁场的具体弯曲型线类型。这样就可以将用数学上描述出的弯曲磁场位形的特征参数α与离子径向二次流运动直接联系起来, 方便于根据工程设计的具体要求通过优化磁场位形的弯曲角α来设计出符合要求的最优弯曲磁场。
其次, 需要设计弯曲角α沿轴向的分布。弯曲磁场的设计主要是为了实现通道内离子的平行射流。在磁力线与等电势面重合的假设前提下, 磁力线弯曲程度沿轴向的设计需要根据离子在通道内的理想运动轨迹来实现。理想的离子运动轨迹虽然最终表现为在出口处平行射流, 但是在通道内的轨迹却是随着离子产生区域距离出口的不同而有所不同的。若离子产生于电离区之前, 那么其在后续运动中受到电场力的作用时间较长, 而在电离区以后产生的离子在后续运动中受电场力作用时间则短, 显然要使得出口处的不同离子的运动方向都能近似平行于射流两则的磁力线而所需的弯曲角α应该是不同的, 否则可能会有“过聚焦”的情况出现 (如图5 的a处所示) , 即离子可能会从一个壁面沿径向迁移而打向另外一个壁面, 这样仍然会带来壁面溅射损失。因此, 需要根据离子动力学的要求设计磁场弯曲角α沿通道轴向Z的最优分布。近阳极区磁场较弱, 基本可以忽略;从电离加速区开始, 离子大量产生且在该区急剧加速, 能量较高, 需要设置较大的弯曲角;在接近出口处校正离子近平行射流, 弯曲角基本为零。即沿通道轴向弯曲角α的定性变化基本要满足图6的要求。
最后是热化电势的干扰分析以及抑制措施设计。上述的分析是基于磁力线与等电势面重合的假设, 实际上它们二者之间并不相同, 差别在于等离子体密度沿磁力线分布不均匀所带来的干扰项, 如式 (7) 所示。分析干扰项具体表达形式为
《3.4 弯曲磁场设计的效果评估》
3.4 弯曲磁场设计的效果评估
由于磁场是SPT的关键设计因素, 因此设计了弯曲磁场位形后, 需要从几方面来评估弯曲磁场设计的效果和影响。
首先, 分析发动机的性能影响。设计了优化的弯曲磁场位形以后, 可以通过控制径向电场分布来控制离子的径向二次流运动, 使之近似平行喷射, 这样就可以大大减小SPT的羽流发散角, 如第一代SPT羽流发散角近45°, 而采用了凸向阳极的弯曲磁场以后, 第二代ATON型稳态等离子体发动机的羽流发散角<10°。随着弯曲磁场弯曲角α沿通道的分布设计体系的提出, 优化设计后的磁场位形将会使羽流发散角进一步减少, 从而将会大大降低对星载器件的干扰。羽流发散角的减少还会使得出口离子径向速度减少, 在相同的电压降加速离子的前提下, 离子出口的轴向速度就会增加, 因此表现为SPT的比冲将会进一步提高。另外, 弯曲磁场对离子径向二次流的控制使得在通道内打向壁面的离子数目将会减少, 从而大大降低离子壁面溅射腐蚀效应, 减少端壁损失, 提高SPT的效率和寿命。如第一代SPT的效率为50 %, 而ATON发动机的效率已达近68 %, 效果显著。
其次, 分析发动机的工作过程影响。磁场对电子运动有着强烈的约束作用, 在外加电磁场作用下电子整体将作霍尔漂移运动。当鞘层和壁面与霍尔漂移的电子发生相互作用时, 在破坏和重新建立霍尔漂移运动的过程中电子就会有一定的轴向位移产生, 表现为电子轴向电导机制。因此, 与径向磁场的情况相比, 在弯曲磁场下, 电子霍尔漂移运动的方向将会有所改变, 与壁面和鞘层发生相互作用时, 电子的轴向导电情况也会有所不同, 具体影响还有待于进一步研究。
《4 SPT磁场设计问题的进一步讨论和展望》
4 SPT磁场设计问题的进一步讨论和展望
由于磁场是约束等离子体运动, 提高SPT性能的最关键设计因素, 因此, 磁场设计将在未来不同用途的新型SPT设计中起到至关重要的作用, 这里我们试图提出一些在未来发展中的磁场设计问题。
1) SPT宽通道效应时的磁场设计问题
由于未来深空探测, 空间轨道转移等任务对于推进器都提出了大功率、高效率、高比冲的要求, 所以要满足这些空间任务的需求, 就需要研制大功率的SPT。功率的提高必然伴随着通道质量流率的加大, SPT的通道宽度也必然随之加大, 此时沿径向的环形通道宽度也将成为不可忽略的设计因素。通道的环形几何形状决定了沿着径向磁场强度是不均匀的, 内壁面的磁场强度大, 外壁面磁场强度小。在该径向磁场梯度的作用下, 电子和离子都会出现磁镜效应, 因此, 引起离子径向二次流运动的原因也变得复杂起来。对于要控制离子径向二次流运动的磁场设计而言, 就不仅要考虑弯曲磁场的弯曲角在通道轴向上的分布要求, 而且需要考虑弯曲角沿通道宽度方向上的设计要求。
2) SPT非定常工作情况下的磁场设计问题
稳态等离子体发动机的工作电源虽然是直流, 但是在实验中发现它在很宽的频带范围内都存在着电流的自发振荡, 大幅值振荡的发生频率为几kHz到几十kHz之间
3) 弯曲磁场设计应用于磁流体发电机 (MHD)
SPT是用于将电能转换成等离子流机械能的装置, 弯曲磁场的设计可以使得这种能量转换的效率提高;而如图7所示的MHD是将等离子流机械能转换成电能的装置, 磁场设计在其中仍然起着至关重要的作用。目前, MHD装置在近壁面处的流体速度变化剧烈, 流体切割磁力线产生的反电动势明显小于中心流场, 边界层内的电流密度较大, 发热量大, 温升高, 这种现象被称为哈特曼效应。哈特曼效应的存在严重地影响了通道壁面处的摩擦系数, 增大了加速通道中的摩擦损失, 成为磁流体加速器设计的技术壁垒