姚露;杨文将;王宝军;汤海滨
【摘 要】由于传统提供附加磁场的液冷铜线圈存在体积大、质量大、电耗大等显著缺点,因此电流密度高、体积小、损耗低的高温超导(HTS)磁体在附加场磁等离子发动机(AF-MPDT)上的应用受到了广泛的关注.然而由于超导材料的临界特性,HTS磁体容易受到外部磁场的影响而发生失超.文章通过对MPDT放电回路的简化,采用仿真和试验的方法来分析简化回路通电时激发的磁场,并研究其对外围HTS磁体的影响程度.仿真结果表明,MPDT放电室内部磁场以周向磁场为主,且由于阳极壁面电流对内部磁场的屏蔽作用,使得外部几乎没有磁场分布.试验结果和仿真也取得了良好的一致性,因此在设计超导磁喷管时几乎可以不考虑来自MPDT放电室磁场的影响.该研究为HTS磁体在AF-MPDT上的应用提供了参考. 【期刊名称】《中国空间科学技术》
【年(卷),期】2018(038)005
【总页数】8页(P30-37)
【关键词】语音云平台磁等离子发动机;磁喷管;高温超导磁体;电磁场;临界电流
【作 者】姚露;杨文将;王宝军;汤海滨
【作者单位】北京航空航天大学 宇航学院,北京 100191;北京航空航天大学 宇航学院,北京 100191;北京航空航天大学 宇航学院,北京 100191;北京航空航天大学 宇航学院,北京 100191
【正文语种】中 文
【中图分类】V439+.4
磁等离子发动机(Magnetoplasmadynamic Thruster, MPDT)在推力和比冲方面具有其他发动机不可比拟的优势，是未来深空探测的理想推进方案之一[1]。自场磁等离子发动机(Self-Field Magnetoplasmadynamic Thruster, SF-MPDT)和附加场磁等离子发动机(AF-MPDT)的最初研究可以追溯到1963年。文献[2]很好地总结了之前MPDT的研究工作，说明AF-MPDT由于附加磁场的电磁加速作用，在中低功率范围内，性能较SF-MPDT更加优越。
MPDT的工作原理之一是基于洛伦兹力的作用加速等离子体来产生推力。然而，由于场强较小且能级较低，自身场加速将导致严重的粒子发散和巨大的等离子体能量损失。磁喷管可以限制压缩发散的粒子流，促进粒子动能向轴向转化，其作用效果如同增加了阳极的长度[3]。许多学者研究了磁喷管场强和形状对MPDT工作性能的影响[4-5]。研究表明，推力的增加与JTB呈线性关系，其中JT是加载在推力器的总电流密度，B是磁场的磁感应强度。另外，磁喷管有利于降低放电室阴极的腐蚀，从而延长MPDT的使用寿命[6]。
磁喷管磁场一般由电磁线圈或者永磁体产生，其磁感应强度通常在0.1～0.5 T。在大多数AF-MPDT的试验中，该磁场由水冷铜线圈提供。然而铜线圈却存在线圈质量过大，对电源和水冷系统要求过高等缺点。随着高温超导(High-Temperature Superconductor, HTS)技术的发展，HTS磁体的应用受到越来越广泛的关注。HTS磁体高电流密度、小质量、小体积、强磁场等优点，有利于AF-MPDT的小型化和工程应用。但由于HTS材料受到温度，电流和磁场的限制，带来诸如稳定性、低温条件等问题，因此需要对磁体的设计和工作环境等开展进一步的研究。
HTS磁体的性能容易受到磁场环境的影响，包括自身磁场[7]和外部磁场。本文主要对磁体
的外部磁场展开研究。目前，MPDT的功率为2 kW～4 MW，放电电流最高可以达到18 kA[2]。如此大的电流，在MPDT工作时会激发出强度不容忽视的磁场，可能对HTS磁体的工作性能产生影响。因此需要对该外部磁场的大小及影响展开研究。
本文首先简要说明了AF-MPDT的物理机制，介绍了磁喷管尤其是超导磁喷管在空间推进上的应用。考虑到HTS磁体对外部磁场的敏感性，本文对MPDT放电室的磁场进行了研究。文章首先采用合理方法简化了MPDT的放电回路，随后对放电回路模拟器进行了仿真计算，分析了它内部磁场的组成，以及沿轴向和径向的磁场分布。最后本文将放电模拟器与HTS磁体结合，进行了临界电流的测量试验。通过本文的研究，为应用于MPDT的超导磁喷管的设计提供了重要参考。
1 AF-MPDT的推力机制
图1是一般AF-MPDT的结构示意。核心部件是阴极和喷管形状的阳极。推进剂(He，Xe等)在放电电流的作用下电离生成等离子体，然后在放电室内经过高温加热、电磁场作用等复杂过程喷出形成推力。
图1 AF-MPDT结构示意Fig.1 Structural representation of AF-MPDT
根据经典的Tikhonov模型[8]，AF-MPDT的推力由气动推力FGD、磁喷管推力FH和自身磁场推力FSF三个部分组成，它们的表达式分别为：
(1)
FH=2KHIB0RA
(2)
(3)
式中： KGD为气动力系数；a0为阴极出口的音速；是质量流率；KH为磁喷管推力系数；I为放电电流；B0为附加场磁感应强度；RA为阳极半径；μ0为真空磁导率；RC为阴极半径。
Herdrich在Tikhonov模型的基础上，通过数学推导和不同研究机构的推力器试验数据拟合得到修正和提高的磁喷管推力模型[9]：
(4)
式中：是修正后的磁喷管推力系数，并对阳极半径引入了指数X的修正。在AF-MPDT样机的初步设计中，X可取5/3。
2 超导磁喷管在空间推进中的应用
磁喷管是AF-MPDT产生推力的主要部件之一。与拉瓦尔喷管的工作机理类似，磁场可以将高温推进工质的径向动能转化轴向动能。这里不同的是，磁喷管与高温等离子体几乎没有接触，而且可以方便地改变磁场的构型，从而调节推力[10]。
封装外壳不仅仅在MPDT中，磁喷管在螺旋推力器(Helicon Thruster)[11]、变比冲的磁等离子火箭(Variable Specific Impulse Magnetoplasma Rocket, VASIMR)[12]上均有应用。
国际上对超导磁喷管早有研究。2002年美国约翰逊空间中心应用HTS磁体替换了安装在VASIMR推力器上液冷铜线圈[13]。磁体采用Bi-2333超导线材，内径0.23 m，整体高度接近0.04 m，中心磁场在105 A的运行电流下可以达到0.28 T，总质量不超过5 kg。2008年，一个由Tai-yang研究公司和NASA先进空间推进实验室组成的团队设计、建造和测试了一个用于空间推进装置的HTS磁体和流动冷却系统[14]。磁体的内径0.16 m，高0.055 m，
在126 A工作电流下，可以达到0.5 T的中心磁场。这项研究证明了在空间推进装置中应用流动冷却HTS磁体的方案可行性，同时说明该方案具有结构紧凑、质量小、效率高和可靠性高等优点。
3 高温超导带材的特性
只有在达到对温度、电流和磁场的特定要求，HTS磁体才能实现并维持超导状态。图2是由函数f1(J,H,T=0)，f2(J,T,H=0)和f3(T,H,J=0)围成的HTS磁体的临界表面。只有在该曲面下HTS磁体才处于超导状态。其中3个最重要的参数是临界温度Tc、临界电流密度Jc和临界磁场强度Hc。这3个参数相互关联。当温度为常数时，根据Kim模型[15]，临界电流密度和磁场强度的关系为：
Jc≃
(5)
式中：αc和H0是常数，αc表示用于平衡洛伦兹力密度的趋向力密度；H为磁场强度。
半透明纸由式(5)可知，在一定温度下，HTS磁体的临界电流密度与外界磁场强度成反比。因此在设计HTS磁体时，不能忽略外界磁场的影响。对于应用于AF-MPDT的超导磁喷管，其外部磁场可能主要来自与放电室电流所激发的磁场，因此需要开展研究。
图2 HTS磁体的三个临界参数Fig.2 Three critical parameters of HTS magnets
4 MPDT放电回路模拟器的仿真研究
4.1 放电回路的简化
MPDT内部的放电机制十分复杂，至今仍没有建立出合适的理论模型。在这里，只考虑MPDT放电回路所产生的磁场环境，忽略产热和等离子体间的相互作用等因素。放电回路如图3所示，其中Ia为轴向电流，Ir为径向电流[16]。MPDT的放电电弧沿着等离子体的流动方向在阴阳极之间不均匀的分布，但有如下关系式：
(6)
式中：k为放电弧柱的数目。
图3 MPDT放电示意Fig.3 Discharge schematic of MPDT
根据式(6)，将各弧柱的电流叠加起来，可以等效在一个电弧面上。事实上，放电电弧存在一定的弧度，但是沿轴向和径向两个方向分解为Irx和Irz。其中Irz与阴级和阳极通过的电流Ia激发的磁场形状相当，通过叠加耦合可以用电流Ia来等效替代。这样由∑Irx就可以得到垂直于轴向z的电弧面。得到的简化放电回路模型如图4所示，之后的研究就是围绕该模型进行。
图4 简化的放电回路模型Fig.4 Simplified discharge circuit model
4.2 放电模拟器的有限元仿真
本文采用COMSOL Multiphysics 有限元软件进行仿真计算。根据简化的放电回路建立的三维仿真模型和网格结构如图5所示，几何形状参考某型AF-MPDT来确认，如图6所示。模拟器通过的电流分别以外电流密度(A/m2)的形式添加在阳极、阴极和电弧面上。其中电弧面的电流密度是在柱坐标中沿半径r的方向施加的。销子材料
图5 三维仿真模型网格结构Fig.5 Mesh structure of 3D simulation model
图6 放电回路模型几何尺寸Fig.6 Geometry of the discharge circuit model
计算得到的磁场分布云图如图7所示，图中模拟器通过电流为100 A。
图7 100 A时放电模拟器磁感应强度分布Fig.7 Discharge simulator magnetic flux density distribution at 100 A
4.3 仿真结果分析
醋酸乙酯生产由图8可知，放电模拟器内部磁场以周向磁感应强度By为主，轴向磁感应强度Bz和径向磁感应强度Bx十分微弱。该周向磁感应强度与径向放电电流作用产生的轴向洛伦兹力，是自身磁场推力FSF的主要组成部分。
在不同放电电流下在电弧面处沿径向的磁场分布如图9所示。由图可知模拟器内部磁感应强度与电流值成正比，且沿径向不断衰减，特别是在阳极壁面衰减速度迅速增加，在外围减小到0。阳极壁面电流发挥了类似于封闭筒的作用，将磁场束缚在了模拟器内部。坐标x=80 mm处，即设计的HTS磁体的内径位置，磁感应强度趋近于0，因此对附加线圈的工作性能几乎不会产生影响。
多媒体电教
图8 1 000 A时放电模拟器磁感应强度组成Fig.8 Discharge simulator magnetic flux density composition at 1 000 A
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议