搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

纵磁作用下真空电弧单阴极斑点等离子体射流三维混合模拟

王振兴 曹志远 李瑞 陈峰 孙丽琼 耿英三 王建华

王振兴, 曹志远, 李瑞, 陈峰, 孙丽琼, 耿英三, 王建华. 纵磁作用下真空电弧单阴极斑点等离子体射流三维混合模拟. 物理学报, 2021, 70(5): 055201. doi: 10.7498/aps.70.20201701
引用本文: 王振兴, 曹志远, 李瑞, 陈峰, 孙丽琼, 耿英三, 王建华. 纵磁作用下真空电弧单阴极斑点等离子体射流三维混合模拟. 物理学报, 2021, 70(5): 055201. doi: 10.7498/aps.70.20201701
Wang Zhen-Xing, Cao Zhi-Yuan, Li Rui, Chen Feng, Sun Li-Qiong, Geng Ying-San, Wang Jian-Hua. Three-dimensional hybrid simulation of single cathode spot vacuum arc plasma jet under axial magnetic field. Acta Phys. Sin., 2021, 70(5): 055201. doi: 10.7498/aps.70.20201701
Citation: Wang Zhen-Xing, Cao Zhi-Yuan, Li Rui, Chen Feng, Sun Li-Qiong, Geng Ying-San, Wang Jian-Hua. Three-dimensional hybrid simulation of single cathode spot vacuum arc plasma jet under axial magnetic field. Acta Phys. Sin., 2021, 70(5): 055201. doi: 10.7498/aps.70.20201701

纵磁作用下真空电弧单阴极斑点等离子体射流三维混合模拟

王振兴, 曹志远, 李瑞, 陈峰, 孙丽琼, 耿英三, 王建华

Three-dimensional hybrid simulation of single cathode spot vacuum arc plasma jet under axial magnetic field

Wang Zhen-Xing, Cao Zhi-Yuan, Li Rui, Chen Feng, Sun Li-Qiong, Geng Ying-San, Wang Jian-Hua
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 真空电弧的特性直接受到从阴极斑点喷射出的等离子体射流的影响, 对等离子体射流进行数值仿真有助于我们深入了解真空电弧的内部物理机制. 然而, 磁流体动力学和粒子云网格仿真方法受限于计算精度和计算效率的原因, 无法有效地应用于真空电弧等离子体射流仿真模拟. 本文开发了一套三维等离子体混合模拟算法, 并在此基础上建立了真空电弧单阴极斑点射流仿真模型, 模型中将离子作宏粒子考虑, 而电子作无质量流体处理, 仿真计算了自生电磁场与外施纵向磁场作用下等离子体的分布运动状态. 仿真结果表明, 单个阴极斑点情况下真空等离子体射流在离开阴极斑点后扩散至极板间, 其整体几何形状为圆锥形, 离子密度从阴极到阳极快速下降. 外施纵向磁场会压缩等离子体, 使得等离子体射流径向的扩散减少并且轴线上的离子密度升高. 随着外施纵向磁场的增大, 其对等离子体射流的压缩效应增强, 表现为等离子体射流的扩散角度逐渐减小. 此外, 外施纵向磁场对等离子体射流的影响也受到电弧电流大小的影响, 压缩效应随电弧电流的增加而逐渐减弱.
    Vacuum arc is a special metal vapor discharge phenomenon, because its discharge medium totally comes from the evaporation and ionization of electrode materials. In the case of low current, the vacuum arc is completely composed of plasma jets emitted from discrete cathode spots on the cathode surface and the current carried by each spot depends on the cathode material. When the arc current exceeds a certain value, a certain number of cathode spot plasma jets will appear. Vacuum arcs play a very important role in some industrial applications such as vacuum circuit breakers, vacuum coatings and electric thrusters. As an important plasma control method, the external axial magnetic field (AMF) has an important influence on the macroscopic morphology and microscopic parameter distribution of the vacuum arc. Various studies of vacuum arc under AMF have been carried out and some progress has been made. However, the existing literature about the simulation research of vacuum arc is mostly concentrated in the case of large current, and less attention is paid to the case of small current. The reason is that the traditional methods, magneto-hydrodynamics or particle-in-cell, are limited by either accuracy or efficiency, and cannot be effectively applied to the low current vacuum arc plasma jet simulations. In this paper, we develop a fully three-dimensional hybrid plasma simulation algorithm to study the single cathode spot vacuum arc plasma jet under AMF. In this model, ions are modelled as particles while electrons are treated as massless fluid, and the self-generated magnetic field is also considered. To simplify the condition, the cathode spot in our model only exists as a plasma jet source, thus the detailed mechanism of producing plasmas is neglected. And the movement of the cathode spot is not considered either. The results show that the single cathode spot plasma jet diffuses into the interelectrode in a cone shape after leaving the cathode spot, and the ion density drops rapidly from cathode to anode. Under the simulation conditions in this paper (I ≤ 150 A), the self-generated magnetic field will not have a significant influence on the plasma jet itself in the case of low current. The external AMF has a compressive effect on the diffusion of the vacuum arc plasma jet. Under the AMF, the radial movement of the ions is suppressed, and the decrease of the ion radial velocity leads to a smaller diffusion radius of the jet. This compression effect of the AMF on the plasma jet is related to both the intensity of the external AMF and the magnitude of the arc current. In the case of a constant arc current magnitude, the compression effect gradually increases as the value of the AMF intensity gradually increases; in the case of a constant value of the external AMF, the compression effect gradually decreases as the current gradually becomes larger.
      通信作者: 王振兴, zxwang@xjtu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 51937009, 51807147)和陕西省自然科学基金(批准号: 2019JM-158)资助的课题
      Corresponding author: Wang Zhen-Xing, zxwang@xjtu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 51937009, 51807147) and the Natural Science Foundation of Shaanxi Province, China (Grant No. 2019JM-158)

    真空电弧是一种特殊的金属蒸气放电现象[1], 其放电介质来自电极材料的蒸发和电离, 在阳极未活跃时该放电通道完全由阴极表面离散的阴极斑点发射出的等离子体射流组成. 每个阴极斑点所能承载的电流取决于电极材料, 当电弧电流超过一定值时, 则需要一定数量的阴极斑点射流共同承担. 作为等离子体物理领域一种常见的放电形式, 真空电弧在真空断路器[2]、真空镀膜[3]和航天电推进器[4,5]等应用中都扮演着非常重要的角色.

    作为一种重要的等离子体控制手段, 外施纵向磁场对于真空电弧的宏观形貌及微观参数分布有着重要的影响[6]. 国内外研究者针对外施纵磁下真空电弧等离子体进行了大量的实验与数值仿真研究. Rondeel[7]通过实验和理论分析发现, 在磁场的作用下, 等离子射流中的电子受磁场作用做减速运动并产生径向电场, 离子的运动受该径向电场的影响; Keidar等[8]建立了大电流真空电弧模型, 仿真解释了外施纵向磁场与电弧电压之间的“L”形关系曲线; 王立军等[9]基于电子和离子双流体模型建立了大电流情况下真空电弧的magneto-hydrodynamics (MHD)模型, 仿真发现外施马鞍形分布的纵向磁场能够更有效地抑制大电流情况下真空电弧的收缩; Jia等[10]使用MHD方法建立了外施纵向磁场下考虑阴极处电流密度分布的大电流真空电弧模型, 考虑了阴极斑点分布不均对于真空电弧特性的影响; Wang等[11]实验研究了外施纵向磁场对真空电弧等离子体射流的影响, 发现随着纵向磁场强度的增加, 等离子体射流开始收缩并且其几何形状从圆锥形逐渐变为圆柱形.

    相比于实验研究方法, 数值仿真计算因其高效和经济的特点已经成为了研究等离子体物理问题的重要工具. 现有文献对于外施纵向磁场下真空电弧等离子体的仿真研究大都集中于多个阴极斑点同时存在的情况下[8-10], 而较少地关注单个阴极斑点等离子体射流情况下的物理问题. 其原因在于现有的仿真计算方法在处理单个阴极斑点等离子体射流时, 无法有效地平衡计算精度和计算效率二者之间的关系. 对于MHD仿真方法, 当电弧电流较小时, 极板间等离子体密度较低且极不均匀, 碰撞频率相对较低, 局部区域可能会偏离热平衡条件从而导致流体假设不再适用. 另一种常常采用的等离子体仿真方法是基于第一性原理的粒子云网格(particle-in-cell, PIC) 仿真方法[12]. 此种方法将电子和离子作为宏粒子考虑, 可以更精确地描述等离子体特性. 但是由于要追踪大量宏粒子的运动信息, 使得仿真计算对于计算机的性能要求较高, 耗时较长, 因此只适合将其应用于时空尺度较小或者等离子体密度较低的物理问题之中.

    最近越来越多的研究者开始关注混合等离子体仿真方法(hybrid plasma simulation), 以解决某些特定物理问题仿真中存在的计算精度与仿真效率之间的矛盾[13-15]. 混合仿真方法可以灵活地将不同类型的粒子建模为宏粒子或者流体, 还可以针对物理问题在不同区域采用不同的建模方法, 使得其可以在保留特定粒子动力学效应的前提下平衡仿真模型的效率和准确性. 但是, 此种仿真方法尚处于探索阶段, 并无成熟的商业软件可用, 为相关研究带来了困难.

    本文的目的是通过混合等离子体仿真方法研究外施纵向磁场对真空电弧单阴极斑点等离子体射流特性的影响. 第2节首先介绍了真空电弧单阴极斑点等离子体射流物理模型, 在模型中将离子作宏粒子考虑而电子作无质量流体处理, 同时考虑了电弧电流产生的自生磁场, 然后介绍了本课题组开发的一套等离子体混合仿真代码的数学模型及其推进步; 第3节展示了不同外施纵向磁场强度下的仿真计算结果; 第4节从能量守恒的角度理论分析了外施纵向磁场对于单阴极斑点等离子体射流扩散半径的影响, 验证了仿真结果的正确性; 第5节为本文的结论部分.

    模型的建立基于如下假定:

    1)等离子体射流完全电离, 极板间等离子体仅包括离子和电子两种带电粒子, 不考虑中性粒子的作用;

    2)等离子体是准中性的, 因此给定离子的电荷密度, 可以确定电子密度;

    3)在求解区域中, 等离子体参数满足le d, τe te, 其中le为电子平均自由程, d为电极间距, τe为电子平均碰撞时间, te为电子从阴极运动到阳极所需时间, 因此在仿真模型中用流体理论来近似描述电子输运过程;

    4)仿真模型中忽略了电子的惯性分量, 电子被认为是无质量流体, 并且其温度保持恒定[16];

    5)在本模型中“阴极斑点”不是真正的阴极斑点, 而是略微高于真实阴极斑点的等离子体射流的截面, 模型中其仅作为等离子体射流源存在, 不考虑近阴极斑点区域等离子体产生的细致过程, 下文为了描述方便, 仍称之为阴极斑点, 并且针对不同大小的电弧电流和外施磁场, 阴极斑点的大小均保持不变;

    6)不考虑阴极斑点的运动.

    真空电弧单阴极斑点等离子体射流模型如图1所示, 等离子体射流由阴极表面的阴极斑点喷射出的电子和离子组成. 电子和离子从阴极斑点发射出后扩散至极板间, 在向阳极运动的过程中同时伴随着沿径向的扩散运动. 根据实验测量结果[17], 在真空电弧中阴极斑点喷射出的离子的动能为40—60 eV, 而离子的温度只有0.5—1.5 eV, 所以离子的运动更多体现的是粒子特性. 而电子则刚好相反, 其热速度要远大于其漂移速度, 所以电子的运动更多体现的是气体特性. 因此在我们的仿真模型中, 将离子当作宏粒子处理, 而电子则当作流体处理.

    图 1 物理模型示意图\r\nFig. 1. The schematic of physical model.
    图 1  物理模型示意图
    Fig. 1.  The schematic of physical model.

    电子被当作无质量流体, 其动量方程用下式描述:

    tnemeve=0=ene(E+ve×B)Pe+enejσ,
    (1)

    式中ne, me, vee分别代表电子数密度、电子质量、电子速度和电子电荷; Pe是电子压力; j是电流密度; σ是等离子体电导率; EB分别为电场强度和磁感应强度. (1)式等号右边的最后一项代表电子和离子之间碰撞阻力的影响[18].

    电子被当作理想气体:

    Pe=nekBTe,
    (2)

    式中, kB是玻尔兹曼常数, Te是电子温度. 电子和离子之间的碰撞频率νei、等离子体电导率σ分别为[19]:

    νei=(lnΛ/10)Zine3.5×1010T3/2e,
    (3)
    σ=nee2meνei,
    (4)

    式中, lnΛ为库仑常数, Zi为离子平均电荷.

    离子作为宏粒子考虑, 其位移由运动方程决定, 加速度由洛伦兹力决定:

    dxidt=vi,
    (5)
    midvidt=qi(E+vi×B)qijσ,
    (6)

    式中, mi, xi, qivi分别代表离子质量、离子位移、离子电荷和离子速度. (6)式等号右边的最后一项和电子动量方程(1)式中的对应项平衡.

    电子密度和离子密度满足电中性条件:

    neZini.
    (7)

    电子和离子是电流的载体, 电流密度与电子和离子运动速度之间的关系为

    j=qinivieneve.
    (8)

    电流密度由安培环路定律给出, 其中位移电流项被忽略(达尔文静磁近似):

    j=1μ0×B,
    (9)

    式中μ0为真空磁导率. 磁场的演变遵循法拉第电磁感应定律:

    Bt=×E.
    (10)

    除了上面方程中的安培环路定律(9)式和电磁感应定律(10)式, 麦克斯韦方程组还另外包括磁通连续性原理和高斯定律两个方程. 在本文的仿真模型中, 磁场被储存在仿真网格的网格面上, 并且直接在网格面上计算其变化. 因此当磁场按照电磁感应定律计算时, 在数值上磁场将会保持为舍入误差内的无散度场, 满足了磁通连续性原理[20]. 由于电中性方程(7)式, 高斯定律满足D=0.

    (1)式—(10)式就是等离子体混合算法所包含的泛定方程, 在每个时间步结束后需要更新并代入下个时间步的变量包括离子的速度信息和位置信息(vi,xi)及磁感应强度(B). 代码推进采用蛙跳格式, 离子速度存储在时间整点上, 离子位置和磁感应强度则存储在时间半点上. 在t = n时间步时变量信息为(vni,xn1/2i,Bn1/2), 程序执行步骤如下.

    1)使用t = n时刻的离子速度vni推动离子位置从xn1/2ixn+1/2i:

    xn+1/2ixn1/2iΔt=vni,
    (11)

    同时基于离子的位置信息及电中性条件(7)式, 获得t = n时刻的离子数密度nn1和电子数密度nnet = n + 1/2时刻的电子数密度nn+1/2e.

    2)使用t = n – 1/2时刻的磁感应强度Bn1/2计算电流密度jn1/2:

    jn1/2=1μ0×Bn1/2.
    (12)

    3)采用预测-修正算法[21]分两步推动磁感应强度从Bn1/2Bn+1/2, 首先是预测步. 计算t = n时刻的电子速度vne(*号是由于此处使用了t = n – 1/2时刻的电流密度jn1/2代替 t = n时刻的电流密度jn, 从而导致等号左右两端时间步不匹配):

    vne=qinnivnijn1/2enne.
    (13)

    计算t = n时刻的电场强度En:

    En=vne×Bn1/2Peenne+jn1/2σ.
    (14)

    使用t = n时刻的电场强度En推动磁感应强度从Bn1/2Bn:

    BnBn1/2Δt/2=×En.
    (15)

    4)修正步. 使用t = n时刻的磁感应强度Bn计算电流密度jn:

    jn=1μ0×Bn.
    (16)

    计算t = n时刻的电子速度vne:

    vne=qinnivnijnenne.
    (17)

    计算t = n时刻的电场强度en:

    En=vne×BnPeenne+jnσ.
    (18)

    使用t = n时刻的电场强度En推动磁感应强度从Bn1/2Bn+1/2:

    Bn+1/2Bn1/2Δt=×En.
    (19)

    5)更新离子速度. 首先使用t = n+1/2时刻的磁感应强度Bn+1/2计算电流密度jn+1/2:

    jn+1/2=1μ0×Bn+1/2.
    (20)

    然后计算t = n+1/2时刻的电场强度E(n+1/2) (*号是由于此处使用了t = n时刻的电子速度vne):

    E(n+1/2)=vne×Bn+1/2Peenn+1/2e+jn+1/2σ.
    (21)

    使用Boris算法更新离子速度从vnivn+1i:

    mivn+1ivniΔt=qi(E(n+1/2)+vn+1i+vni2×Bn+1/2)ejn+1/2σ.
    (22)

    t = n时间步结束后, 离子速度、离子位置和磁感应强度的信息从(vni, xn1/2i, Bn1/2)更新到了(vn+1i, xn+1/2i, Bn+1/2), 仿真计算进入下一个时间步. 程序的执行框图如图2所示.

    图 2 程序执行步骤\r\nFig. 2. Execution steps of the program.
    图 2  程序执行步骤
    Fig. 2.  Execution steps of the program.

    对于阴极侧, 阴极斑点作为等离子体源, 其有效半径设置为R = 1 mm, 阴极侧除阴极斑点之外的区域设置为吸收边界. 计算区域阴极侧的电子温度设置为Te = 2.5 eV, 离子温度为Ti = 0.5 eV, 离子初速度在轴向上分量vi0, 离子初速度为1.0 × 104 m/s, 离子平均电荷取Zi = 1.85. 假定阴极斑点处电流密度j0和离子密度n0分布均匀, 满足如下条件:

    j0=I0πR2,
    (23)
    n0=γj0mivi0,
    (24)

    式中: I0是阴极斑点总电流; γ为电极侵蚀率, 取35 μg/C (电极选为纯铜触头); vi0为离子初速度在轴向上的分量. 对于磁感应强度, 由阴极侧电流密度分布计算磁矢位A来给出其第一类边界条件:

    2A=μ0j,
    (25)
    B=×A.
    (26)

    对于阳极侧, 假定其在现有条件下还不活跃, 没有向极板间喷射金属蒸气, 被看作是电流和粒子收集器. 整个阳极表面看作是一个等电位面, 在阳极和等离子体之间有一个阳极鞘层, 阳极鞘层的电压降用下式表示[22]:

    φsh=jσ+Peene+ve×B,
    (27)
    φsh=kBTeelnjejth,
    (28)
    jth=14ene8kBTeπme,
    (29)

    式中: φsh表示阳极鞘层电压; jth表示随机电子电流密度, 和电子温度相关; je表示电子电流密度.

    在本文中, 模拟区域大小为20 mm × 20 mm × 10 mm, 阴极和阳极分别位于z = 0和z = 10 mm处, 阴极斑点被设置在阴极中心, 外施纵向磁场在整个仿真区域中均匀分布.

    三维空间的离子数密度分布如图3所示. 电弧电流为30 A, 无外施纵向磁场(Bz = 0 mT). 从图3可以看到, 等离子体射流在离开阴极斑点后呈锥形扩散状, 离子数密度从阴极到阳极逐渐减小, 且射流中心的离子数密度大于射流边缘处的离子数密度. 对于单阴极斑点的等离子体射流来说, 由于电荷电流较小, 所以其自生磁场对等离子体射流的影响较弱. 所以等离子体射流的形状受离子初始喷射角的影响较大.

    图 3 三维空间离子数密度分布\r\nFig. 3. The distribution of ion number density in 3D space.
    图 3  三维空间离子数密度分布
    Fig. 3.  The distribution of ion number density in 3D space.

    由于在本文中只考虑了单个阴极斑点等离子体射流, 因此等离子体参数在x-z平面和y-z平面具有相同的对称分布, 因此在图4中仅展示了电流和自生磁场在x-z平面的分布. 在图4(a)中轴向电流的方向从阳极指向阴极, 由图3可知在等离子体射流从阴极向阳极运动的过程中, 作为电流载流体的等离子体横截面逐渐增大, 所以在图4(a)中可以观察到轴向电流密度从阴极到阳极逐渐减小. 值得注意的是, 等离子体射流电流密度分布可能与等离子体密度分布不同. 原因是等离子体密度分布主要和离子与电子的位置有关, 而电流密度的分布则主要受电子的漂移速度影响. 在真空电弧中, 由于电中性条件, 电子和离子的密度相同, 但是由于离子电流只占总电流的10%[23], 因此从(8)式中可以计算得到电子的漂移速度约为离子的11倍.

    图 4 (a) 轴向电流密度分布; (b)自生磁感应强度分布 I = 30 A, Bz = 0 mT\r\nFig. 4. (a) The distribution of axial current density; (b) the distribution of self-generated azimuthal magnetic field I = 30 A, Bz = 0 mT
    图 4  (a) 轴向电流密度分布; (b)自生磁感应强度分布 I = 30 A, Bz = 0 mT
    Fig. 4.  (a) The distribution of axial current density; (b) the distribution of self-generated azimuthal magnetic field I = 30 A, Bz = 0 mT

    图4(b)展示了方向与平板电极平行的等离子体射流自生环向磁场, 其方向遵守右手螺旋定则. 从图中可以看到磁场分布和等离子体射流的锥形分布相对应, 在轴线上(x = y = 0)自生磁感应强度为0, 并且磁感应强度从中轴线上到仿真区域的径向边界先上升后下降. 这是由于等离子体电流通道并未占据整个仿真区域, 所以自生磁感应强度的最大值出现在等离子体射流边缘处.

    图5展示了电弧电流为30 A时不同外施纵向磁场强度(0, 25, 50, 75 mT)下离子数密度在x-z平面的分布. 从图中可以看出, 在同等电弧电流情况下, 随着纵向磁场强度的增加, 离子密度分布在径向的扩散变得越来越小, 射流形状逐渐由锥形被压缩成圆柱形. 这说明外施纵向磁场对离子的径向扩散有抑制作用, 随着纵向磁场磁感应强度的增加, 纵向磁场对等离子体射流的压缩效应越来越强.

    图 5 在I= 30 A时不同外施纵磁条件下离子数密度分布 (a) Bz = 0 mT; (b) Bz = 25 mT; (c) Bz = 50 mT; (d) Bz = 75 mT\r\nFig. 5. Ion number density distributions under different external AMFs at I= 30 A: (a) Bz = 0 mT; (b) Bz = 25 mT; (c) Bz = 50 mT; (d) Bz = 75 mT.
    图 5  I= 30 A时不同外施纵磁条件下离子数密度分布 (a) Bz = 0 mT; (b) Bz = 25 mT; (c) Bz = 50 mT; (d) Bz = 75 mT
    Fig. 5.  Ion number density distributions under different external AMFs at I= 30 A: (a) Bz = 0 mT; (b) Bz = 25 mT; (c) Bz = 50 mT; (d) Bz = 75 mT.

    图6展示了电弧电流为30 A时, 不同外施纵向磁场强度下轴线上(x = y = 0)从阴极到阳极的离子数密度变化. 从图6可以看到, 不同外施纵向磁场情况下, 离子数密度的最大值相同且均出现在阴极表面, 然后从阴极到阳极逐渐降低. 同时与图5相对应, 施加纵向磁场后, 离子径向扩散减少导致等离子体被压缩, 使得轴线上离子数密度升高. 并且外施纵向磁场越大, 轴线上离子数密度越高. 在阴极斑点附近(z = 0—0.5 mm), 离子数密度的变化趋势受纵向磁场的影响不明显.

    图 6 不同外施纵磁条件下轴线上离子数密度变化\r\nFig. 6. Ion number density distributions along the axis under different external AMFs.
    图 6  不同外施纵磁条件下轴线上离子数密度变化
    Fig. 6.  Ion number density distributions along the axis under different external AMFs.

    图7展示了电弧电流为30 A时, 外施纵向磁场分别为0 mT和50 mT时, 离子沿x正方向的速度在x-z平面的相空间分布. 从图中的粒子分布可以看到, 在外施纵向磁场的作用下单阴极斑点等离子体射流沿径向的扩散减少, 说明外施纵向磁场对等离子体射流的扩散起到了很强的束缚作用, 这一现象和图5中的离子数密度分布相对应. 同时观察离子沿x方向的速度分布可以看出, 在外施纵向磁场的作用下, 离子沿径向速度的绝对值也因受到限制而变小.

    图 7 不同外施磁场条件下离子沿x方向速度在x-z平面的相空间分布\r\nFig. 7. Phase diagram of ion velocity along x-direction in x-z plane under different external AMFs.
    图 7  不同外施磁场条件下离子沿x方向速度在x-z平面的相空间分布
    Fig. 7.  Phase diagram of ion velocity along x-direction in x-z plane under different external AMFs.

    图8展示了外施纵向磁场为75 mT时不同电弧电流大小 (30, 60, 90, 120 A)下离子数密度在x-z平面的分布. 从图中可以看到, 当外施纵向磁场强度不变时, 随着电弧电流的逐渐增大, 等离子体射流的形状逐渐从圆柱形变成锥形, 离子在径向上的扩散逐渐增多. 这说明对于单阴极斑点等离子体射流来说, 增大电弧电流可以抑制外施纵向磁场对于等离子体射流的压缩效应.

    图 8 在Bz = 75 mT时不同电弧电流条件下离子数密度分布 (a) I = 30 A; (b) I = 60 A; (c) I = 90 A; (d) I = 120 A\r\nFig. 8. Ion number density distributions with different arc currents at Bz = 75 mT: (a) I = 30 A; (b) I = 60 A; (c) I = 90 A; (d) I = 120 A.
    图 8  Bz = 75 mT时不同电弧电流条件下离子数密度分布 (a) I = 30 A; (b) I = 60 A; (c) I = 90 A; (d) I = 120 A
    Fig. 8.  Ion number density distributions with different arc currents at Bz = 75 mT: (a) I = 30 A; (b) I = 60 A; (c) I = 90 A; (d) I = 120 A.

    可以用轴线上阳极处离子密度和阴极处离子密度的比值来衡量外施纵向磁场对等离子体射流的压缩效应. 不同情况下阳极处离子数密度与阴极处离子数密度的比值如图9所示. 当无外施纵向磁场时, 等离子体射流的扩散不受影响, 不同电弧电流所对应的比值几乎相等. 在同等电弧电流情况下, 随着纵向磁场强度的增加, 轴线上阳极处离子数密度与阴极处离子数密度的比值也随之变大. 由于电弧电流相同时阴极处离子数密度一致(图6), 这说明阳极处离子数密度随着外施纵向磁场的增大而逐渐升高. 由此可知纵向磁场对等离子体射流的压缩效应会随着纵向磁场强度的增加而越来越强. 在同等外施纵向磁场情况下, 随着电弧电流的逐渐增大, 比值随之变小. 这说明电弧电流的升高会抑制纵向磁场对等离子体射流的压缩效应.

    图 9 轴线上阳极处离子数密度与阴极处离子数密度的比值\r\nFig. 9. The ratios of the ion number density at the anode to that at the cathode on the axis.
    图 9  轴线上阳极处离子数密度与阴极处离子数密度的比值
    Fig. 9.  The ratios of the ion number density at the anode to that at the cathode on the axis.

    通过以上的仿真结果可以发现, 外施纵向磁场会压缩等离子射流, 抑制离子沿径向的扩散使得等离子体射流的扩散半径减小, 同时这种压缩效应不仅和外施纵向磁场的强度有关, 而且和电弧电流的自身大小有关. 该模拟结果和Wang[11]在实验中观察到的现象一致. 在本节中, 我们从能量守恒的角度对外施纵磁对等离子体射流的压缩效应进行理论分析.

    图3所示, 离子从阴极斑点喷射出之后呈圆锥形扩散. 其轴向速度使其从阴极向阳极运动, 其径向速度使其向四周扩散. 在单阴极斑点等离子体射流情况下, 由于电弧电流较小, 所以可以忽略自生环向磁场对其的影响. 等离子体和磁场的总能量由3部分组成: 电子能量、离子能量、磁场能量. 其中电子能量主要由电子温度决定(在真空电弧中电子的热速度远大于其漂移速度), 假设电子温度不变, 则电子能量保持恒定. 离子能量由离子沿径向速度的动能和沿轴向速度的动能组成, 又由于外施纵向磁场对轴向速度没有影响, 可以认为离子的轴向速度在整个过程中保持不变. 所以等离子体射流在外施纵向磁场中的运动过程可以简化为离子的径向动能和等离子体磁场能量之间的相互转换过程. 对于阴极斑点附近一个长为Δz的等离子体圆柱(Δz取值极小), 由能量守恒可得:

    12M0v2r0+B2z2μ0πr20Δz=12M0v2r+B2z2μ0πr2Δz,
    (30)
    M0=n0miπr20Δz,
    (31)

    式中, M0为该等离子体圆柱的质量, vr0为离子初始径向速度, r0为阴极斑点半径, vr为离子径向速度, r为等离子体射流半径. 从图7中可以发现, 在外施纵磁作用下, 离子的径向速度在阴极斑点附近最大. 随着等离子体射流的逐渐扩散, 等离子体射流半径逐渐增大, 离子的径向速度也逐渐变小. 所以当满足vr = 0时, 等离子体射流半径r取到最大值. 将vr = 0和(31)式代入(30)式可得:

    r2r20r20=μ0n0miv2r0B2z.
    (32)

    由于等离子体射流半径最大值r远大于阴极斑点半径r0, 所以r2r20r2. 由(32)式可得:

    r|vr0|r0|Bz|μ0n0mi.
    (33)

    将阴极斑点处离子密度n0的表达式(24)式代入(33)式可得:

    r|vr0|r0|Bz|μ0γj0vi0.
    (34)

    从(34)式中可以发现等离子体射流半径的最大值r和外施纵向磁感应强度Bz的取值成反比, 体现为纵向磁场对等离子体射流有压缩效应, 当外施纵向磁场的强度值增大时, 等离子体射流的扩散半径会减小. 同时从(34)式中还可以发现, 在外施纵向磁场强度不变Bz的情况下, 电流密度j0的增大可以减弱外施纵向磁场对等离子体射流的压缩效应, 该结论与仿真结果与图9相对应.

    本文基于三维混合等离子体仿真方法建立了外施纵向磁场下真空电弧等离子体射流仿真模型, 模型中将离子作为宏粒子考虑, 而电子作为无质量流体处理, 仿真计算了单阴极斑点等离子体射流在自生电磁场与外施纵向磁场共同作用下等离子体的分布运动状态. 仿真结果如下.

    1)单阴极斑点等离子体射流在离开阴极斑点后呈锥形扩散至极板间, 离子密度从阴极到阳极快速下降. 在本文的仿真条件下(30 A ≤ I ≤ 120 A), 等离子体射流产生的自生磁场不会对等离子体射流本身产生显著的影响.

    2)外施纵向磁场对等离子体射流的扩散有压缩效应. 在外施纵向磁场的作用下, 离子沿径向的运动受到抑制, 其径向速度变小导致等离子体射流的扩散半径变小, 同时使得沿轴线上的离子密度升高.

    3)外施纵向磁场对等离子体射流的压缩效应与外施纵向磁场的强度值及等离子体射流自身的电流大小有关. 在等离子体射流自身电流大小不变的情况下, 当纵向磁场强度值逐渐增大时, 压缩效应逐渐增强; 在外施纵向磁场强度值不变的情况下, 当等离子体射流自身的电流逐渐变大时, 压缩效应逐渐减弱.

    作者衷心感谢鞍山师范学院的王虹宇教授在计算方法方面的帮助.

    [1]

    Beilis I I 2001 IEEE Trans. Plasma Sci. 29 657Google Scholar

    [2]

    王建华, 耿英三, 刘志远, 闫静 2017 高压电器 53 1Google Scholar

    Wang J H, Geng Y S, Liu Z Y, Yan J 2017 High Volt. Appar. 53 1Google Scholar

    [3]

    Sanders D M, Anders A 2000 Surf. Coat. Technol. 133 78Google Scholar

    [4]

    Geng J Y, Chen Y C, Sun S R, Huang W D, Wang H X 2020 Plasma Sci. Technol. 22 094012Google Scholar

    [5]

    Keidar M, Zhuang T, Shashurin A, Teel G, Chiu D, Lukas J, Haque S, Brieda L 2014 Plasma Phys. Controlled Fusion 57 014005Google Scholar

    [6]

    王立军, 贾申利, 杨泽, 史宗谦 2017 高压电器 3 22Google Scholar

    Wang L J, Jia S L, Yang Z, Shi Z Q 2017 High Volt. Appar. 3 22Google Scholar

    [7]

    Rondeel W G J 1975 J. Phys. D: Appl. Phys. 8 934Google Scholar

    [8]

    Keidar M, Schulman M B 2001 IEEE Trans. Plasma Sci. 29 684Google Scholar

    [9]

    王立军, 贾申利, 史宗谦, 荣命哲 2005 中国电机工程学报 25 113Google Scholar

    Wang L J, Jia S L, Shi Z Q, Rong M Z 2005 Chin. Soc. for Elec. Eng. 25 113Google Scholar

    [10]

    Jia S L, Zhang L, Wang L J, Chen B, Shi Z Q, Sun W 2011 IEEE Trans. Plasma Sci. 39 3233Google Scholar

    [11]

    Wang C, Shi Z Q, Wu B Z, Gao Z P, Jia S L, Wang L J 2016 J. Phys. D: Appl. Phys. 49 135203Google Scholar

    [12]

    李晗蔚, 孙安邦, 张幸, 姚聪伟, 常正实, 张冠军 2018 物理学报 4 143Google Scholar

    Li H W, Sun A B, Zhang X, Yao C W, Chang Z S, Zhang G J 2018 Acta Phys. Sin. 4 143Google Scholar

    [13]

    Shmelev D L, Uimanov I V 2015 IEEE Trans. Plasma Sci. 43 2261Google Scholar

    [14]

    Shmelev D L, Oreshkin V I, Uimanov I V 2019 IEEE Trans. Plasma Sci. 47 3478Google Scholar

    [15]

    Li C, Ebert U, Hundsdorfer W 2010 J. Comput. Phys. 229 200Google Scholar

    [16]

    Arai K, Takahashi S, Morimiya O, Niwa Y 2003 IEEE Trans. Plasma Sci. 31 929Google Scholar

    [17]

    Kutzner J, Miller H C 1992 J. Phys. D: Appl. Phys. 25 686Google Scholar

    [18]

    Winske D, Omidi N 1991 Hybrid Codes: Methods and Applications (New Mexico: Los Alamos National Lab.) pp103−160

    [19]

    Beilis I I, Keidar M, Boxman R L, Goldsmith S 1998 J. Appl. Phys. 83 709Google Scholar

    [20]

    Tóth G 2000 J. Comput. Phys. 161 605Google Scholar

    [21]

    Harned D S 1982 J. Comput. Phys. 47 452Google Scholar

    [22]

    Schade E, Shmelev D L 2003 IEEE Trans. Plasma Sci. 31 890Google Scholar

    [23]

    Kutzner J, Miller H C 1989 IEEE Trans. Plasma Sci. 17 688Google Scholar

    期刊类型引用(5)

    1. 董华军,程靖洲,赵一鉴,库照宇,李东恒. 大电流真空电弧开断过程瞬态特性仿真分析. 电机与控制学报. 2024(01): 189-196 . 百度学术
    2. 蒋原,马速良,武建文,李擎,朱立颖. 中频真空电弧边缘击穿现象及仿真研究. 电工技术学报. 2024(09): 2887-2895 . 百度学术
    3. 姚芳,张鑫,申会员,彭美南,李冬,郭文谦. 墙壁插座直流燃弧特性测试技术研究. 电工电能新技术. 2022(02): 81-88 . 百度学术
    4. 黄小龙,孙韬,冉垠康,冯治国,赵莉华,王立军. 直流开断过程中真空电弧等离子体瞬态特性仿真研究. 广东电力. 2021(12): 28-37 . 百度学术
    5. 王宏达,王立军,施炜鑫,王凯,贾申利. 两种不同结构纵—横复合触头磁场计算及电弧实验研究. 高压电器. 2021(11): 10-17+26 . 百度学术

    其他类型引用(6)

  • 图 1  物理模型示意图

    Fig. 1.  The schematic of physical model.

    图 2  程序执行步骤

    Fig. 2.  Execution steps of the program.

    图 3  三维空间离子数密度分布

    Fig. 3.  The distribution of ion number density in 3D space.

    图 4  (a) 轴向电流密度分布; (b)自生磁感应强度分布 I = 30 A, Bz = 0 mT

    Fig. 4.  (a) The distribution of axial current density; (b) the distribution of self-generated azimuthal magnetic field I = 30 A, Bz = 0 mT

    图 5  I= 30 A时不同外施纵磁条件下离子数密度分布 (a) Bz = 0 mT; (b) Bz = 25 mT; (c) Bz = 50 mT; (d) Bz = 75 mT

    Fig. 5.  Ion number density distributions under different external AMFs at I= 30 A: (a) Bz = 0 mT; (b) Bz = 25 mT; (c) Bz = 50 mT; (d) Bz = 75 mT.

    图 6  不同外施纵磁条件下轴线上离子数密度变化

    Fig. 6.  Ion number density distributions along the axis under different external AMFs.

    图 7  不同外施磁场条件下离子沿x方向速度在x-z平面的相空间分布

    Fig. 7.  Phase diagram of ion velocity along x-direction in x-z plane under different external AMFs.

    图 8  Bz = 75 mT时不同电弧电流条件下离子数密度分布 (a) I = 30 A; (b) I = 60 A; (c) I = 90 A; (d) I = 120 A

    Fig. 8.  Ion number density distributions with different arc currents at Bz = 75 mT: (a) I = 30 A; (b) I = 60 A; (c) I = 90 A; (d) I = 120 A.

    图 9  轴线上阳极处离子数密度与阴极处离子数密度的比值

    Fig. 9.  The ratios of the ion number density at the anode to that at the cathode on the axis.

  • [1]

    Beilis I I 2001 IEEE Trans. Plasma Sci. 29 657Google Scholar

    [2]

    王建华, 耿英三, 刘志远, 闫静 2017 高压电器 53 1Google Scholar

    Wang J H, Geng Y S, Liu Z Y, Yan J 2017 High Volt. Appar. 53 1Google Scholar

    [3]

    Sanders D M, Anders A 2000 Surf. Coat. Technol. 133 78Google Scholar

    [4]

    Geng J Y, Chen Y C, Sun S R, Huang W D, Wang H X 2020 Plasma Sci. Technol. 22 094012Google Scholar

    [5]

    Keidar M, Zhuang T, Shashurin A, Teel G, Chiu D, Lukas J, Haque S, Brieda L 2014 Plasma Phys. Controlled Fusion 57 014005Google Scholar

    [6]

    王立军, 贾申利, 杨泽, 史宗谦 2017 高压电器 3 22Google Scholar

    Wang L J, Jia S L, Yang Z, Shi Z Q 2017 High Volt. Appar. 3 22Google Scholar

    [7]

    Rondeel W G J 1975 J. Phys. D: Appl. Phys. 8 934Google Scholar

    [8]

    Keidar M, Schulman M B 2001 IEEE Trans. Plasma Sci. 29 684Google Scholar

    [9]

    王立军, 贾申利, 史宗谦, 荣命哲 2005 中国电机工程学报 25 113Google Scholar

    Wang L J, Jia S L, Shi Z Q, Rong M Z 2005 Chin. Soc. for Elec. Eng. 25 113Google Scholar

    [10]

    Jia S L, Zhang L, Wang L J, Chen B, Shi Z Q, Sun W 2011 IEEE Trans. Plasma Sci. 39 3233Google Scholar

    [11]

    Wang C, Shi Z Q, Wu B Z, Gao Z P, Jia S L, Wang L J 2016 J. Phys. D: Appl. Phys. 49 135203Google Scholar

    [12]

    李晗蔚, 孙安邦, 张幸, 姚聪伟, 常正实, 张冠军 2018 物理学报 4 143Google Scholar

    Li H W, Sun A B, Zhang X, Yao C W, Chang Z S, Zhang G J 2018 Acta Phys. Sin. 4 143Google Scholar

    [13]

    Shmelev D L, Uimanov I V 2015 IEEE Trans. Plasma Sci. 43 2261Google Scholar

    [14]

    Shmelev D L, Oreshkin V I, Uimanov I V 2019 IEEE Trans. Plasma Sci. 47 3478Google Scholar

    [15]

    Li C, Ebert U, Hundsdorfer W 2010 J. Comput. Phys. 229 200Google Scholar

    [16]

    Arai K, Takahashi S, Morimiya O, Niwa Y 2003 IEEE Trans. Plasma Sci. 31 929Google Scholar

    [17]

    Kutzner J, Miller H C 1992 J. Phys. D: Appl. Phys. 25 686Google Scholar

    [18]

    Winske D, Omidi N 1991 Hybrid Codes: Methods and Applications (New Mexico: Los Alamos National Lab.) pp103−160

    [19]

    Beilis I I, Keidar M, Boxman R L, Goldsmith S 1998 J. Appl. Phys. 83 709Google Scholar

    [20]

    Tóth G 2000 J. Comput. Phys. 161 605Google Scholar

    [21]

    Harned D S 1982 J. Comput. Phys. 47 452Google Scholar

    [22]

    Schade E, Shmelev D L 2003 IEEE Trans. Plasma Sci. 31 890Google Scholar

    [23]

    Kutzner J, Miller H C 1989 IEEE Trans. Plasma Sci. 17 688Google Scholar

  • [1] 周雄峰, 陈彬, 刘坤. 氩气等离子体射流特性: 电压、气流、外磁场的综合影响. 物理学报, 2024, 73(22): 225201. doi: 10.7498/aps.73.20241166
    [2] 张雪雪, 贾鹏英, 冉俊霞, 李金懋, 孙换霞, 李雪辰. 辅助放电下刷状空气等离子体羽的放电特性和参数诊断. 物理学报, 2024, 73(8): 085201. doi: 10.7498/aps.73.20231946
    [3] 佟磊, 赵明亮, 张钰如, 宋远红, 王友年. 带有射频偏压源的感性耦合Ar/O2/Cl2等离子体放电的混合模拟研究. 物理学报, 2024, 73(4): 045201. doi: 10.7498/aps.73.20231369
    [4] 胡杨, 罗婧怡, 蔡雨烟, 卢新培. 外加磁场对螺旋等离子体的影响. 物理学报, 2023, 72(13): 130501. doi: 10.7498/aps.72.20222442
    [5] 牛中国, 许相辉, 王建锋, 蒋甲利, 梁华. 飞翼模型纵向气动特性等离子体流动控制试验. 物理学报, 2022, 71(2): 024702. doi: 10.7498/aps.71.20211425
    [6] 税敏, 席涛, 闫永宏, 于明海, 储根柏, 朱斌, 何卫华, 赵永强, 王少义, 范伟, 卢峰, 杨雷, 辛建婷, 周维民. 激光等离子体射流驱动亚毫米直径铝飞片及姿态诊断. 物理学报, 2022, 71(9): 095201. doi: 10.7498/aps.71.20212136
    [7] 牛中国, 许相辉, 王建峰, 蒋甲利, 梁华. 飞翼模型纵向气动特性等离子体流动控制试验研究. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211425
    [8] 杨丽君, 宋彩虹, 赵娜, 周帅, 武珈存, 贾鹏英. 大气压氩气刷形等离子体羽的特性研究. 物理学报, 2021, 70(15): 155201. doi: 10.7498/aps.70.20202091
    [9] 张亚容, 韩乾翰, 郭颖, 张菁, 石建军. 大气压脉冲放电等离子体射流特性及机理研究. 物理学报, 2021, 70(9): 095202. doi: 10.7498/aps.70.20202246
    [10] 张钰如, 高飞, 王友年. 低气压感性耦合等离子体源模拟研究进展. 物理学报, 2021, 70(9): 095206. doi: 10.7498/aps.70.20202247
    [11] 王鹏, 沈赤兵. 等离子体合成射流对超声速混合层的混合增强. 物理学报, 2019, 68(17): 174701. doi: 10.7498/aps.68.20190683
    [12] 郭恒, 张晓宁, 聂秋月, 李和平, 曾实, 李志辉. 亚大气压六相交流电弧放电等离子体射流特性数值模拟. 物理学报, 2018, 67(5): 055201. doi: 10.7498/aps.67.20172557
    [13] 郭恒, 苏运波, 李和平, 曾实, 聂秋月, 李占贤, 李志辉. 亚大气压六相交流电弧等离子体射流特性研究:实验测量. 物理学报, 2018, 67(4): 045201. doi: 10.7498/aps.67.20172556
    [14] 李刘合, 刘红涛, 罗辑, 许亿. 带状真空电弧磁过滤器等离子体分布特性及制备类金刚石膜研究. 物理学报, 2016, 65(6): 065202. doi: 10.7498/aps.65.065202
    [15] 杨成, 周昕. 液态水中的多种局域结构. 物理学报, 2016, 65(17): 176501. doi: 10.7498/aps.65.176501
    [16] 伍飞飞, 廖瑞金, 杨丽君, 刘兴华, 汪可, 周之. 棒-板电极直流负电晕放电特里切尔脉冲的微观过程分析. 物理学报, 2013, 62(11): 115201. doi: 10.7498/aps.62.115201
    [17] 廖瑞金, 伍飞飞, 刘兴华, 杨帆, 杨丽君, 周之, 翟蕾. 大气压直流正电晕放电暂态空间电荷分布仿真研究. 物理学报, 2012, 61(24): 245201. doi: 10.7498/aps.61.245201
    [18] 王晶, 马瑞玲, 王龙, 孟俊敏. 采用混合模型数值模拟从深海到浅海内波的传播. 物理学报, 2012, 61(6): 064701. doi: 10.7498/aps.61.064701
    [19] 吴 翊, 荣命哲, 杨 飞, 王小华, 马 强, 王伟宗. 引入6波段P-1辐射模型的三维空气电弧等离子体数值分析. 物理学报, 2008, 57(9): 5761-5767. doi: 10.7498/aps.57.5761
    [20] 方道腴. 在横向磁场中真空电弧后退速度和阴极温度的关系. 物理学报, 1983, 32(7): 838-844. doi: 10.7498/aps.32.838
  • 期刊类型引用(5)

    1. 董华军,程靖洲,赵一鉴,库照宇,李东恒. 大电流真空电弧开断过程瞬态特性仿真分析. 电机与控制学报. 2024(01): 189-196 . 百度学术
    2. 蒋原,马速良,武建文,李擎,朱立颖. 中频真空电弧边缘击穿现象及仿真研究. 电工技术学报. 2024(09): 2887-2895 . 百度学术
    3. 姚芳,张鑫,申会员,彭美南,李冬,郭文谦. 墙壁插座直流燃弧特性测试技术研究. 电工电能新技术. 2022(02): 81-88 . 百度学术
    4. 黄小龙,孙韬,冉垠康,冯治国,赵莉华,王立军. 直流开断过程中真空电弧等离子体瞬态特性仿真研究. 广东电力. 2021(12): 28-37 . 百度学术
    5. 王宏达,王立军,施炜鑫,王凯,贾申利. 两种不同结构纵—横复合触头磁场计算及电弧实验研究. 高压电器. 2021(11): 10-17+26 . 百度学术

    其他类型引用(6)

计量
  • 文章访问数:  6545
  • PDF下载量:  148
  • 被引次数: 11
出版历程
  • 收稿日期:  2020-10-14
  • 修回日期:  2020-11-10
  • 上网日期:  2021-02-24
  • 刊出日期:  2021-03-05

/

返回文章
返回