CN110278648B - 磁场位形对等离子体减压效应的影响研究方法 - Google Patents

磁场位形对等离子体减压效应的影响研究方法 Download PDF

Info

Publication number
CN110278648B
CN110278648B CN201910559066.7A CN201910559066A CN110278648B CN 110278648 B CN110278648 B CN 110278648B CN 201910559066 A CN201910559066 A CN 201910559066A CN 110278648 B CN110278648 B CN 110278648B
Authority
CN
China
Prior art keywords
magnetic field
particles
plasma
motion
particle
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
Application number
CN201910559066.7A
Other languages
English (en)
Other versions
CN110278648A (zh
Inventor
毛保全
白向华
李程
杨雨迎
宋瑞亮
朱锐
张天意
赵其进
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Academy of Armored Forces of PLA
Original Assignee
Academy of Armored Forces of PLA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Academy of Armored Forces of PLA filed Critical Academy of Armored Forces of PLA
Priority to CN201910559066.7A priority Critical patent/CN110278648B/zh
Publication of CN110278648A publication Critical patent/CN110278648A/zh
Application granted granted Critical
Publication of CN110278648B publication Critical patent/CN110278648B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • H—ELECTRICITY
    • H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05H—PLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
    • H05H1/00—Generating plasma; Handling plasma
    • H05H1/0006—Investigating plasma, e.g. measuring the degree of ionisation or the electron temperature
    • H—ELECTRICITY
    • H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05H—PLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
    • H05H1/00—Generating plasma; Handling plasma
    • H05H1/02—Arrangements for confining plasma by electric or magnetic fields; Arrangements for heating plasma
    • H05H1/10—Arrangements for confining plasma by electric or magnetic fields; Arrangements for heating plasma using externally-applied magnetic fields only, e.g. Q-machines, Yin-Yang, base-ball
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/10—Nuclear fusion reactors

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Mathematical Analysis (AREA)
  • Plasma & Fusion (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Optics & Photonics (AREA)
  • Plasma Technology (AREA)

Abstract

本发明公开了磁场位形对等离子体减压效应的影响研究方法，包括以下步骤：S1：均匀磁场减压研究，带电粒子在圆管中均匀磁场区域中运动时，等离子体中单个带电粒子在外加磁场作用下的运动，其他粒子对它的作用被完全忽略，这是等离子体的单粒子运动模型；S2：非均匀磁场减压研究；S3：模型验证，减压效应是圆管结构中等离子体受到磁场约束作用时的集体行为。本发明的研究方法更加的科学合理，粒子在垂直于磁场平面上的运动轨迹是由有限半径的拉莫尔回旋运动和磁场不均匀引起的漂移运动的叠加，回旋运动使带电粒子在垂直于磁场方向的运动受到抑制，当一个离子和中性原子碰撞后，离子会向不同的方向移动。

Description

磁场位形对等离子体减压效应的影响研究方法
技术领域
本发明涉及磁场领域，具体为磁场位形对等离子体减压效应的影响研究方法。
背景技术
经过多年的发展，等离子体技术已在材料、微电子、化工、受控核聚变、磁流体发电、机械及环保等众多学科领域中得到较广泛地应用，并已初步形成等离子体工业体系，在军事领域，等离子体隐身、通信和探测技、等离子体拦截等技术近年来的发展十分迅速，在火炮方面，国内外对等离子体点火以及火药燃烧产生等离子体已有一些深入的研究。
理想磁场位形划分为均匀磁场和非均匀磁场，均匀磁场是指圆管内部轴向和径向的磁感应强度是均匀的磁场，非均匀磁场主要包括梯度、弯曲等几种典型的磁场类型，带电粒子回旋中心垂直于磁力线方向的运动称为漂移，等离子体在这些不均匀磁场中的漂移现象也会对径向受力产生影响。
发明内容
本发明的目的在于提供磁场位形对等离子体减压效应的影响研究方法，以解决上述背景技术中提出的问题。
为实现上述目的，本发明提供如下技术方案：磁场位形对等离子体减压效应的影响研究方法，包括以下步骤：
S1：均匀磁场减压研究，带电粒子在圆管中均匀磁场区域中运动时，等离子体中单个带电粒子在外加磁场作用下的运动，其他粒子对它的作用被完全忽略，这是等离子体的单粒子运动模型，粒子在外加磁场中的运动服从牛顿(Newton)方程：
Figure GDA0003598236580000011
假设圆管中施加的外磁场沿Z轴方向且不随时间发生变化，可近似认为B＝Bez，E＝f＝0，粒子的运动方程为：
Figure GDA0003598236580000021
解上式得：
Figure GDA0003598236580000022
Figure GDA0003598236580000023
z＝vPt+z0其中rL＝v⊥/|ωc|＝mv⊥/|q|B，vP、v⊥分别为平行磁场方向和垂直磁场方向的速度分量，则由上式可得：
(x-x0)2+(y-y0)2＝rL 2
由上式可知，在均匀磁场中，粒子垂直于感应强度B的速度分量v⊥，使粒子在垂直于B的平面以半径rL＝mv⊥/qB做匀速圆周运动，平行于磁感应强度B的分量vP不会受到磁场的作用，使粒子保持直线运动状态，通常情况下粒子的回旋半径rL与圆管的半径相比非常小，所以粒子在垂直磁场方向上受到磁场的约束；
等离子体中通常含有一种以上的正离子，当正离子和电子没有达到平衡，则正离子和电子应作为两种不同的粒子体系考虑，采用简化双流体模型的方式建立等离子体双流体方程组：
Figure GDA0003598236580000024
式中
Figure GDA0003598236580000031
--为运流导数；
由于等离子体满足准电中性，因此认为离子电荷数与电子电荷数相等，为了减小计算难度简化模型，不考虑等离子体本身的电荷、电流产生的场，假设：
ρ＝nimi+neme≈n(mi+me)
Figure GDA0003598236580000032
Figure GDA0003598236580000033
Figure GDA0003598236580000034
Figure GDA0003598236580000035
上式方程组包括等离子体的质量守恒、动量守恒、能量守恒方程，将上式与Maxwell方程组联立，加上欧姆定律以及局部热平衡状态条件，最后得到等离子体的磁流体方程组：
Figure GDA0003598236580000036
式中σc为电导率，由于不考虑电场的作用，且外加磁场是恒定的，因此上式中忽略了电场项和位移电流项；
根据上述等离子体的磁流体方程组，在comsol软件中建立二维轴对称模型，r＝0为对称轴，流体区域长200mm，宽15mm，流体从下端进入圆管中，从上端流出，模型中忽略重力的影响，模型中线圈用正方体表示，调用COMSOL软件中的AC/DC模块中的磁场模块，设置线圈导体模型为均匀多匝，每个的线圈匝数为100匝，线圈导线截面积为10-6m2，线圈材料设置为铜，电导率为6×107S/m，求解域设置为空气，等离子体流体模型采用等离子体模块和和单相层流模块共同构建，单相流体的层流模块用来分析作为流体的等离子体的运动，两个物理场模块是通过等离子体所受到的洛伦兹力耦合，磁源与等离子体、流体通过多物理场中的等离子体电导耦合场进行耦合计算；
圆管中部区域0.05m至0.15m处的的磁通密度变化较小，0.05m和0.15m处的磁通密度均为0.0092T，0.1m处的磁通密度为0.0095T，变化幅度为3.1％，圆管中轴线和靠近线圈处的圆管内壁磁感应强度相同，表明圆管内部的磁场分布是均匀分布；
获得均匀磁场后，对等离子体流体进行瞬态分析，层流模块中设置流体为马赫数小于0.3的可压缩流，此时式中流体的运动方程变为：
Figure GDA0003598236580000041
式中F为流体单元受到磁场作用的洛伦兹力，磁场由通电线圈和等离子体流体运动共同产生，设置流体流质量流量SCCM为10mL/min，等离子体的初始电子密度为1020m-3，初始电子温度10eV，环境压力为1[Torr]，因为稳态磁场已经完成计算，因此直接计算瞬态流体流动，得到不同时刻流体的电子密度分布；
S2：非均匀磁场减压研究，带电粒子在磁场中的运动，通常会因磁场的变化而产生的漂移，例如在缓变磁场和不均匀磁场中，带电粒子都有一些独特的性质，等离子体中的带电粒子在梯度磁场中运动时，由于磁场的不均匀性，运动方程为：
Figure GDA0003598236580000051
此时磁场是关于回旋半径的函数，假设圆管中非均匀磁场为缓变磁场，即粒子回旋半径、螺旋轨道的螺距远小于非均匀性的特征长度，满足磁场不随时间变化时的缓变条件
Figure GDA0003598236580000052
则可将粒子在非均匀磁场中的运动轨迹看成是无数个均匀磁场中的漂移运动的叠加，当磁场的梯度方向垂直于自身方向时，粒子沿粒子
Figure GDA0003598236580000053
的运动，漂移速度为：
Figure GDA0003598236580000054
式中W⊥为带电粒子垂直动能；μ为磁矩；q为带电粒子电荷量，说明梯度漂移与粒子的电荷量和符号有关，正离子与电子的梯度漂移方向相反，
当磁力线弯曲时，等离子体的回旋中心沿弯曲磁力线运动时会产生指向曲率中心的曲率漂移，记磁力线的曲率半径为
Figure GDA0003598236580000055
粒子沿
Figure GDA0003598236580000056
方向漂移，漂移速度为：
Figure GDA0003598236580000057
由上式可知，曲率漂移同样与粒子的电荷符号相关，正离子与电子的曲率漂移方向相反，
通常情况下，梯度漂移和曲率漂移是同时存在的，带电粒子在非均匀磁场中的漂移总速度为
Figure GDA0003598236580000058
磁场关于Z轴对称，而且沿Z轴方向缓慢汇聚，磁场强度逐渐增大，由于洛伦兹力始终垂直于粒子运动速度，因而粒子动能不变，粒子动能W可以表示为平行磁场方向的动能WP与垂直磁场方向的动能W⊥之和，当带电粒子在缓变磁场中由弱磁场区域向强磁场区域运动时，随着磁场B的不断增强，W⊥也不断增大，由于磁场力与粒子的运动速度始终垂直，洛伦兹力对粒子不做功，粒子总能量W＝W⊥+WP保持守恒，故随着W⊥增大，WP将逐渐减小，此时由
Figure GDA0003598236580000061
可知粒子的横向运动vP减小，粒子回旋运动的半径及角速度为
Figure GDA0003598236580000062
与
Figure GDA0003598236580000063
当带电粒子运动到强磁场区域时，r变小而ω变大，带电粒子在非均匀磁场中的运动是回旋半径与轴向速度逐渐减小的螺旋运动；
取柱坐标系，
Figure GDA0003598236580000064
粒子速度
Figure GDA0003598236580000065
则轴向运动方程为：
Figure GDA0003598236580000066
对于缓变磁场，在回旋半径粒子的rc范围内，可得：
Figure GDA0003598236580000067
对于正离子，在柱坐标系中vθ＝-v⊥，和负离子vθ＝v⊥，则可得
Figure GDA0003598236580000068
式中
Figure GDA0003598236580000069
为平行磁场方向的梯度，上式表明带电粒子在纵向不均匀磁场中运动时，将受到从强磁场区域指向弱磁场区域的作用力；
磁场由两个平行的载流圆环线圈产生，磁场沿z轴方向先减小后增大，当粒子运动到磁场足够强的区域时，如果有W＝W⊥，即径向动能等于总动能，此时带电粒子的平行向动能WP＝0，粒子就不能继续沿轴线方向前进，但粒子仍受力FP的作用，因此会出现带电粒子运动“反射”的现象，于是粒子在两端强磁场区域之间来回反射；
并非所有的带电粒子都能被反射，当粒子运动到圆环线圈截面处，此时若速度vz没有减小为零，则带电粒子从磁镜中穿过，而不会被反射；
单粒子轨道描述法可以绘制出带电粒子在电磁场中的运动轨迹，对于温度和压力较低的等离子体，由于粒子密度比较低，在研究粒子在磁场中的运动时，可以忽略粒子间的相互作用，当然，单粒子轨道描述并不能给出等离子体的集体行为，但可以用于研究粒子在非均匀磁场中的运动规律；
利用COMSOL软件中的带电粒子追踪模块和AC/DC模块耦合建立带电粒子在磁场中的运动模型，圆管结构半径15mm，内壁面边界条件设置为反弹，带电粒子在圆管的一端释放，粒子以圆管轴线为中心轴线向圆管内部运动，在磁场中将圆环结构设置为线圈，线圈类型为圆形，匝数20匝，通入电流10A，在圆管端部释放20000个带电粒子，粒子质量设为电子质mp，电荷数设置为-1，粒子以圆柱轴线为中心轴线向通道内部运动，法向速度设为800m/s，粒子进入磁场区后，受到洛伦兹力的作用，对于速度方向不与轴线平行的带电粒子，磁场通过拉莫尔回旋使粒子沿轴线方向运动；
圆管内部不同区域的磁通密度分布曲线，图中0mm、5mm、10mm、15mm表示距圆管轴线等距离平行线，线圈截面上的磁感应强度最大，相邻两线圈的中点截面上的磁感应强度最小，而轴线处的磁通密度较为均匀，平均值为0.055T；在圆管的入口处和出口处，磁场衰减迅速，第一个线圈与第二个线圈间磁感应强度变化比较大，距离轴线的越远的直线上，磁通密度幅值变化越大，从图中可以看出，圆管的径向和轴向都形成了梯度磁场；
粒子进入磁场区后，受到洛伦兹力的作用，对于速度方向不与轴线平行的粒子，磁场通过拉莫尔回旋使粒子沿轴线方向运动；
S3：模型验证，减压效应是圆管结构中等离子体受到磁场约束作用时的集体行为，宏观上表现为壁面处的气体压强降低，当等离子体处于高温高压状态时，圆管壁面的压强数据较难采集，涉及到在圆管上开孔以及气密性等问题，而当等离子体处于低温低密度时，通常较为容易观测到等离子体在磁场中的运动变化，为了测试磁场对圆管中等离子体气体压强的影响，采用低气压低温等离子体试验平台进行等离子体在磁场中的压强变化试验；
试验平台包括设备台架，石英圆管，等离子体电源、真空计，压力计等，台架用于安装各种设备和仪表，等离子体电源可以在大气压下、低气压下和各种气氛环境中产生稳定的电弧放电、介质阻挡放电、辉光放电，石英圆管由真空泵抽至低压状态，两端电极通过辉光放电产生等离子体，辉光放电的视觉特征呈现均匀的雾状放电。
优选的，所述S1中角标α--表示不同的正离子或电子；
Figure GDA0003598236580000081
——为压强张量；
Figure GDA0003598236580000082
——为粘性应力张量；
Figure GDA0003598236580000083
——为摩擦阻力；Q——为不同类型离子碰撞引起的热能交换。
优选的，所述S1中ρ——约化密度；
Figure GDA0003598236580000084
——为约华速度；
Figure GDA0003598236580000085
——为电流密度。
优选的，所述S1中当各粒子的温度各向同性但不相等时，
Figure GDA0003598236580000086
式中mi为第i种粒子的质量，k为玻尔兹曼常数，qi为第i种粒子的电荷量，
Figure GDA0003598236580000087
优选的，所述S1中当各粒子温度相等，即热力学平衡状态时，
Figure GDA0003598236580000088
优选的，所述S2中分析不同线圈间距粒子运动的影响，改变线圈间距为30mm、20mm，使线圈截面处的磁感应强度不变。
优选的，所述S2中非均匀磁场中，带电粒子沿轴线方向幅值呈周期性变化，在线圈界面处的强磁场区域时，粒子受到的洛伦兹力较大，而在两线圈之间的弱磁场区域，粒子受到的洛伦兹力较小，由此可推断出等离子体在多圆环线圈组产生的磁场中，在线圈截面附近区域减压效应较好，而两线圈之间区域减压效应不强。
优选的，所述S3中实验步骤为调节石英管中的气压，采用进压或抽压的方式实现气体压强的调节，调节范围0-3000Pa；两端电极辉光放电产生等离子体，等离子体电流和功率通过控制台调节；将永磁铁固定于一端空真计附近，通过控制距离调节磁感应强度；待石英管中的气压稳定后，记录控制台中的压强数据，重复上述步骤。
与现有技术相比，本发明的有益效果是：本发明的研究方法更加的科学合理，当带电粒子在磁场中运动时，由于受到洛伦兹力，粒子在垂直于磁场平面上的运动轨迹是由有限半径的拉莫尔回旋运动和磁场不均匀引起的漂移运动的叠加，回旋运动使带电粒子在垂直于磁场方向的运动受到抑制，但由于中性粒子的存在，当一个离子和中性原子碰撞后，离子会向不同的方向移动，但它围绕磁场的旋转运动不会改变，因此大量中性粒子与离子、电子碰撞后也会间接受到磁场的影响，中性粒子对圆管壁面的碰撞减少，在宏观上就表现为对圆管内壁面压力的降低。
附图说明
图1为本发明的研究方法的流程图；
图2为本发明的正离子与电子的回旋运动图；
图3为本发明的圆管中粒子的回旋运动图；
图4为本发明的轴线上磁感应强度变化曲线图；
图5为本发明的圆管轴向压强变化曲线图；
图6为本发明的圆管径向压强变化曲线图；
图7为本发明的粒子在非均匀磁场中的运动图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一
请参阅图1，本发明提供一种技术方案：磁场位形对等离子体减压效应的影响研究方法，包括以下步骤：
S1：均匀磁场减压研究，带电粒子在圆管中均匀磁场区域中运动时，等离子体中单个带电粒子在外加磁场作用下的运动，其他粒子对它的作用被完全忽略，这是等离子体的单粒子运动模型，粒子在外加磁场中的运动服从牛顿(Newton)方程：
Figure GDA0003598236580000101
假设圆管中施加的外磁场沿Z轴方向且不随时间发生变化，可近似认为B＝Bez，E＝f＝0，粒子的运动方程为：
Figure GDA0003598236580000102
解上式得：
Figure GDA0003598236580000103
Figure GDA0003598236580000104
z＝vPt+z0其中rL＝v⊥/|ωc|＝mv⊥/|q|B，vP、v⊥分别为平行磁场方向和垂直磁场方向的速度分量，则由上式可得：
(x-x0)2+(y-y0)2＝rL 2
由上式可知，在均匀磁场中，粒子垂直于感应强度B的速度分量v⊥，使粒子在垂直于B的平面以半径rL＝mv⊥/qB做匀速圆周运动，平行于磁感应强度B的分量vP不会受到磁场的作用，使粒子保持直线运动状态，通常情况下粒子的回旋半径rL与圆管的半径相比非常小，所以粒子在垂直磁场方向上受到磁场的约束；
等离子体中通常含有一种以上的正离子，请参阅图2，当正离子和电子没有达到平衡，则正离子和电子应作为两种不同的粒子体系考虑，采用简化双流体模型的方式建立等离子体双流体方程组：
Figure GDA0003598236580000111
式中
Figure GDA0003598236580000112
--为运流导数；
由于等离子体满足准电中性，因此认为离子电荷数与电子电荷数相等；为了减小计算难度简化模型，不考虑等离子体本身的电荷、电流产生的场，假设：
ρ＝nimi+neme≈n(mi+me)
Figure GDA0003598236580000113
Figure GDA0003598236580000114
Figure GDA0003598236580000115
Figure GDA0003598236580000121
上式方程组包括等离子体的质量守恒、动量守恒、能量守恒方程，将上式与Maxwell方程组联立，加上欧姆定律以及局部热平衡状态条件，最后得到等离子体的磁流体方程组：
Figure GDA0003598236580000122
式中σc为电导率，由于不考虑电场的作用，且外加磁场是恒定的，因此上式中忽略了电场项和位移电流项；
根据上述等离子体的磁流体方程组，在comsol软件中建立二维轴对称模型，r＝0为对称轴，流体区域长200mm，宽15mm，流体从下端进入圆管中，从上端流出，模型中忽略重力的影响，模型中线圈用正方体表示，调用COMSOL软件中的AC/DC模块中的磁场模块，设置线圈导体模型为均匀多匝，每个的线圈匝数为100匝，线圈导线截面积为10-6m2，线圈材料设置为铜，电导率为6×107S/m，求解域设置为空气，等离子体流体模型采用等离子体模块和和单相层流模块共同构建，单相流体的层流模块用来分析作为流体的等离子体的运动，两个物理场模块是通过等离子体所受到的洛伦兹力耦合，磁源与等离子体、流体通过多物理场中的等离子体电导耦合场进行耦合计算；
圆管中部区域0.05m至0.15m处的的磁通密度变化较小，0.05m和0.15m处的磁通密度均为0.0092T，0.1m处的磁通密度为0.0095T，变化幅度为3.1％，圆管中轴线和靠近线圈处的圆管内壁磁感应强度相同，表明圆管内部的磁场分布是均匀分布；
获得均匀磁场后，对等离子体流体进行瞬态分析，层流模块中设置流体为马赫数小于0.3的可压缩流，此时式中流体的运动方程变为：
Figure GDA0003598236580000131
式中F为流体单元受到磁场作用的洛伦兹力，磁场由通电线圈和等离子体流体运动共同产生，设置流体流质量流量SCCM为10mL/min，等离子体的初始电子密度为1020m-3，初始电子温度10eV，环境压力为1[Torr]，因为稳态磁场已经完成计算，因此直接计算瞬态流体流动，得到不同时刻流体的电子密度分布，角标α--表示不同的正离子或电子；
Figure GDA0003598236580000132
——为压强张量；
Figure GDA0003598236580000133
——为粘性应力张量；
Figure GDA0003598236580000134
——为摩擦阻力；Q——为不同类型离子碰撞引起的热能交换，ρ——约化密度；
Figure GDA0003598236580000135
——为约华速度；
Figure GDA0003598236580000136
——为电流密度；
当各粒子的温度各向同性但不相等时，
Figure GDA0003598236580000137
式中mi为第i种粒子的质量，k为玻尔兹曼常数，qi为第i种粒子的电荷量，
Figure GDA0003598236580000138
当各粒子温度相等，即热力学平衡状态时，
Figure GDA0003598236580000139
S2：非均匀磁场减压研究，带电粒子在磁场中的运动，通常会因磁场的变化而产生的漂移，例如在缓变磁场和不均匀磁场中，带电粒子都有一些独特的性质，等离子体中的带电粒子在梯度磁场中运动时，由于磁场的不均匀性，运动方程为：
Figure GDA0003598236580000141
此时磁场是关于回旋半径的函数，假设圆管中非均匀磁场为缓变磁场，即粒子回旋半径、螺旋轨道的螺距远小于非均匀性的特征长度，满足磁场不随时间变化时的缓变条件
Figure GDA0003598236580000142
则可将粒子在非均匀磁场中的运动轨迹看成是无数个均匀磁场中的漂移运动的叠加，当磁场的梯度方向垂直于自身方向时，粒子沿粒子
Figure GDA0003598236580000143
的运动，漂移速度为：
Figure GDA0003598236580000144
式中W⊥为带电粒子垂直动能；μ为磁矩；q为带电粒子电荷量，说明梯度漂移与粒子的电荷量和符号有关，正离子与电子的梯度漂移方向相反；
当磁力线弯曲时，等离子体的回旋中心沿弯曲磁力线运动时会产生指向曲率中心的曲率漂移，记磁力线的曲率半径为
Figure GDA0003598236580000145
粒子沿
Figure GDA0003598236580000146
方向漂移，漂移速度为：
Figure GDA0003598236580000147
由上式可知，曲率漂移同样与粒子的电荷符号相关，正离子与电子的曲率漂移方向相反；
通常情况下，梯度漂移和曲率漂移是同时存在的，带电粒子在非均匀磁场中的漂移总速度为
Figure GDA0003598236580000148
磁场关于Z轴对称，而且沿Z轴方向缓慢汇聚，磁场强度逐渐增大，由于洛伦兹力始终垂直于粒子运动速度，因而粒子动能不变，粒子动能W可以表示为平行磁场方向的动能WP与垂直磁场方向的动能W⊥之和，当带电粒子在缓变磁场中由弱磁场区域向强磁场区域运动时，随着磁场B的不断增强，W⊥也不断增大，由于磁场力与粒子的运动速度始终垂直，洛伦兹力对粒子不做功，粒子总能量W＝W⊥+WP保持守恒，故随着W⊥增大，WP将逐渐减小，此时由
Figure GDA0003598236580000151
可知粒子的横向运动vP减小，粒子回旋运动的半径及角速度为
Figure GDA0003598236580000152
与
Figure GDA0003598236580000153
当带电粒子运动到强磁场区域时，r变小而ω变大，带电粒子在非均匀磁场中的运动是回旋半径与轴向速度逐渐减小的螺旋运动；
取柱坐标系，
Figure GDA0003598236580000154
粒子速度
Figure GDA0003598236580000155
则轴向运动方程为：
Figure GDA0003598236580000156
对于缓变磁场，在回旋半径粒子的rc范围内，可得：
Figure GDA0003598236580000157
对于正离子，在柱坐标系中vθ＝-v⊥，和负离子vθ＝v⊥，则可得：
Figure GDA0003598236580000158
式中
Figure GDA0003598236580000159
为平行磁场方向的梯度，上式表明带电粒子在纵向不均匀磁场中运动时，将受到从强磁场区域指向弱磁场区域的作用力；
磁场由两个平行的载流圆环线圈产生，磁场沿z轴方向先减小后增大，当粒子运动到磁场足够强的区域时，如果有W＝W⊥，即径向动能等于总动能，此时带电粒子的平行向动能WP＝0，粒子就不能继续沿轴线方向前进，但粒子仍受力FP的作用，因此会出现带电粒子运动“反射”的现象，于是粒子在两端强磁场区域之间来回反射；
并非所有的带电粒子都能被反射，当粒子运动到圆环线圈截面处，此时若速度vz没有减小为零，则带电粒子从磁镜中穿过，而不会被反射；
单粒子轨道描述法可以绘制出带电粒子在电磁场中的运动轨迹，对于温度和压力较低的等离子体，由于粒子密度比较低，在研究粒子在磁场中的运动时，可以忽略粒子间的相互作用，当然，单粒子轨道描述并不能给出等离子体的集体行为，但可以用于研究粒子在非均匀磁场中的运动规律；
利用COMSOL软件中的带电粒子追踪模块和AC/DC模块耦合建立带电粒子在磁场中的运动模型，请参阅图3，圆管结构半径15mm，内壁面边界条件设置为反弹，带电粒子在圆管的一端释放，粒子以圆管轴线为中心轴线向圆管内部运动，在磁场中将圆环结构设置为线圈，线圈类型为圆形，匝数20匝，通入电流10A，在圆管端部释放20000个带电粒子，粒子质量设为电子质mp，电荷数设置为-1，粒子以圆柱轴线为中心轴线向通道内部运动，法向速度设为800m/s，粒子进入磁场区后，受到洛伦兹力的作用，对于速度方向不与轴线平行的带电粒子，磁场通过拉莫尔回旋使粒子沿轴线方向运动；
圆管内部不同区域的磁通密度分布曲线，请参阅图4-6，图中0mm、5mm、10mm、15mm表示距圆管轴线等距离平行线，由图可知，线圈截面上的磁感应强度最大，相邻两线圈的中点截面上的磁感应强度最小，而轴线处的磁通密度较为均匀，平均值为0.055T；在圆管的入口处和出口处，磁场衰减迅速，第一个线圈与第二个线圈间磁感应强度变化比较大，距离轴线的越远的直线上，磁通密度幅值变化越大，请参阅图7，从图中可以看出，圆管的径向和轴向都形成了梯度磁场，分析不同线圈间距粒子运动的影响，改变线圈间距为30mm、20mm，使线圈截面处的磁感应强度不变；
粒子进入磁场区后，受到洛伦兹力的作用，对于速度方向不与轴线平行的粒子，磁场通过拉莫尔回旋使粒子沿轴线方向运动，非均匀磁场中，带电粒子沿轴线方向幅值呈周期性变化，在线圈界面处的强磁场区域时，粒子受到的洛伦兹力较大，而在两线圈之间的弱磁场区域，粒子受到的洛伦兹力较小，由此可推断出等离子体在多圆环线圈组产生的磁场中，在线圈截面附近区域减压效应较好，而两线圈之间区域减压效应不强；
S3：模型验证，减压效应是圆管结构中等离子体受到磁场约束作用时的集体行为，宏观上表现为壁面处的气体压强降低，当等离子体处于高温高压状态时，圆管壁面的压强数据较难采集，涉及到在圆管上开孔以及气密性等问题，而当等离子体处于低温低密度时，通常较为容易观测到等离子体在磁场中的运动变化，为了测试磁场对圆管中等离子体气体压强的影响，采用低气压低温等离子体试验平台进行等离子体在磁场中的压强变化试验；
试验平台包括设备台架，石英圆管，等离子体电源、真空计，压力计等，台架用于安装各种设备和仪表，等离子体电源可以在大气压下、低气压下和各种气氛环境中产生稳定的电弧放电、介质阻挡放电、辉光放电，石英圆管由真空泵抽至低压状态，两端电极通过辉光放电产生等离子体，辉光放电的视觉特征呈现均匀的雾状放电，实验步骤为调节石英管中的气压，采用进压或抽压的方式实现气体压强的调节，调节范围0-3000Pa；两端电极辉光放电产生等离子体，等离子体电流和功率通过控制台调节；将永磁铁固定于一端空真计附近，通过控制距离调节磁感应强度；待石英管中的气压稳定后，记录控制台中的压强数据，重复上述步骤。
通过对上述实施例进行研究，本发明的研究方法更加的科学合理，当带电粒子在磁场中运动时，由于受到洛伦兹力，粒子在垂直于磁场平面上的运动轨迹是由有限半径的拉莫尔回旋运动和磁场不均匀引起的漂移运动的叠加，回旋运动使带电粒子在垂直于磁场方向的运动受到抑制，但由于中性粒子的存在，当一个离子和中性原子碰撞后，离子会向不同的方向移动，但它围绕磁场的旋转运动不会改变，因此大量中性粒子与离子、电子碰撞后也会间接受到磁场的影响，中性粒子对圆管壁面的碰撞减少，在宏观上就表现为对圆管内壁面压力的降低。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。

Claims (1)

1.磁场位形对等离子体减压效应的影响研究方法，其特征在于，包括以下步骤：
S1：均匀磁场减压研究，带电粒子在圆管中均匀磁场区域中运动时，等离子体中单个带电粒子在外加磁场作用下的运动，其他粒子对它的作用被完全忽略，这是等离子体的单粒子运动模型，粒子在外加磁场中的运动服从牛顿（Newton）方程：
Figure 57395DEST_PATH_IMAGE001
，
假设圆管中施加的外磁场沿Z轴方向且不随时间发生变化，可近似认为
Figure 473333DEST_PATH_IMAGE002
，
Figure 740367DEST_PATH_IMAGE003
，粒子的运动方程为：
Figure 939267DEST_PATH_IMAGE004
解上式得：
Figure 608145DEST_PATH_IMAGE005
Figure 601509DEST_PATH_IMAGE006
Figure 355839DEST_PATH_IMAGE007
其中
Figure 561692DEST_PATH_IMAGE008
，
Figure 85077DEST_PATH_IMAGE009
、
Figure 514922DEST_PATH_IMAGE010
分别为平行磁场方向和垂直磁场方向的速度分量，则由上式可得：
Figure 756547DEST_PATH_IMAGE011
由上式可知，在均匀磁场中，粒子垂直于感应强度B的速度分量
Figure 31671DEST_PATH_IMAGE012
，使粒子在垂直于B的平面以半径
Figure 268617DEST_PATH_IMAGE013
做匀速圆周运动，平行于磁感应强度B的分量
Figure 603783DEST_PATH_IMAGE009
不会受到磁场的作用，使粒子保持直线运动状态，通常情况下粒子的回旋半径
Figure 332705DEST_PATH_IMAGE014
与圆管的半径相比非常小，所以粒子在垂直磁场方向上受到磁场的约束；
等离子体中通常含有一种以上的正离子，当正离子和电子没有达到平衡，则正离子和电子应作为两种不同的粒子体系考虑，采用简化双流体模型的方式建立等离子体双流体方程组：
Figure 411519DEST_PATH_IMAGE015
式中
Figure 909497DEST_PATH_IMAGE016
--为运流导数；由于等离子体满足准电中性，因此认为离子电荷数与电子电荷数相等，为了减小计算难度简化模型，不考虑等离子体本身的电荷、电流产生的场，假设：
Figure 681144DEST_PATH_IMAGE017
Figure 100624DEST_PATH_IMAGE018
上式方程组包括等离子体的质量守恒、动量守恒、能量守恒方程，将上式与Maxwell方程组联立，加上欧姆定律以及局部热平衡状态条件，最后得到等离子体的磁流体方程组：
Figure 983129DEST_PATH_IMAGE019
式中
Figure 335613DEST_PATH_IMAGE020
为电导率，由于不考虑电场的作用，且外加磁场是恒定的，因此上式中忽略了电场项和位移电流项；
根据上述等离子体的磁流体方程组，在comsol软件中建立二维轴对称模型，r=0为对称轴，流体区域长200mm，宽15mm，流体从下端进入圆管中，从上端流出，模型中忽略重力的影响，模型中线圈用正方体表示，调用COMSOL软件中的AC/DC模块中的磁场模块（mf），设置线圈导体模型为均匀多匝，每个的线圈匝数为100匝，线圈导线截面积为10-6m2，线圈材料设置为铜，电导率为6×107S/m，求解域设置为空气，等离子体流体模型采用等离子体模块（plas）和和单相层流（spf）模块共同构建，单相流体的层流模块用来分析作为流体的等离子体的运动，两个物理场模块是通过等离子体所受到的洛伦兹力耦合，磁源与等离子体、流体通过多物理场中的等离子体电导耦合场进行耦合计算；
圆管中部区域0.05m至0.15m处的的磁通密度变化较小，0.05m和0.15m处的磁通密度均为0.0092T，0.1m处的磁通密度为0.0095T，变化幅度为3.1%，圆管中轴线和靠近线圈处的圆管内壁磁感应强度相同，表明圆管内部的磁场分布是均匀分布；
获得均匀磁场后，对等离子体流体进行瞬态分析，层流（spf）模块中设置流体为马赫数小于0.3的可压缩流，此时式中流体的运动方程变为：
Figure 278161DEST_PATH_IMAGE021
式中，
Figure 716096DEST_PATH_IMAGE022
为粘性系数;
Figure 136713DEST_PATH_IMAGE023
为剪切应力;F为流体单元受到磁场作用的洛伦兹力，磁场由通电线圈和等离子体流体运动共同产生，设置流体流质量流量SCCM为10mL/min，等离子体的初始电子密度为1020m3，初始电子温度10eV，环境压力为1Torr，因为稳态磁场已经完成计算，因此直接计算瞬态流体流动，得到不同时刻流体的电子密度分布；
S2：非均匀磁场减压研究，带电粒子在磁场中的运动，通常会因磁场的变化而产生的漂移，在缓变磁场和不均匀磁场中，带电粒子都有一些独特的性质，等离子体中的带电粒子在梯度磁场中运动时，由于磁场的不均匀性，运动方程为：
Figure 202758DEST_PATH_IMAGE024
此时磁场是关于回旋半径的函数，假设圆管中非均匀磁场为缓变磁场，即粒子回旋半径、螺旋轨道的螺距远小于非均匀性的特征长度，满足磁场不随时间变化时的缓变条件
Figure 316207DEST_PATH_IMAGE025
，则可将粒子在非均匀磁场中的运动轨迹看成是无数个均匀磁场中的漂移运动的叠加，当磁场的梯度方向垂直于自身方向时，粒子沿粒子
Figure 241438DEST_PATH_IMAGE026
的运动，漂移速度为：
Figure 465746DEST_PATH_IMAGE027
式中
Figure 527243DEST_PATH_IMAGE028
为带电粒子垂直动能；
Figure 811594DEST_PATH_IMAGE029
为磁矩；q为带电粒子电荷量，说明梯度漂移与粒子的电荷量和符号有关，正离子与电子的梯度漂移方向相反，
当磁力线弯曲时，等离子体的回旋中心沿弯曲磁力线运动时会产生指向曲率中心的曲率漂移，记磁力线的曲率半径为
Figure 692962DEST_PATH_IMAGE030
，粒子沿
Figure 720961DEST_PATH_IMAGE031
方向漂移，漂移速度为：
Figure 902543DEST_PATH_IMAGE032
由上式可知，曲率漂移同样与粒子的电荷符号相关，正离子与电子的曲率漂移方向相反，
通常情况下，梯度漂移和曲率漂移是同时存在的，带电粒子在非均匀
磁场中的漂移总速度为
Figure 92216DEST_PATH_IMAGE033
，磁场关于Z轴对称，而且沿Z轴方向缓慢汇聚，磁场强度逐渐增大，由于洛伦兹力始终垂直于粒子运动速度，因而粒子动能不变，粒子动能
Figure 992039DEST_PATH_IMAGE034
可以表示为平行磁场方向的动能
Figure 151625DEST_PATH_IMAGE035
与垂直磁场方向的动能
Figure 453294DEST_PATH_IMAGE036
之和，当带电粒子在缓变磁场中由弱磁场区域向强磁场区域运动时，随着磁场B的不断增强，
Figure 813868DEST_PATH_IMAGE037
也不断增大，由于磁场力与粒子的运动速度始终垂直，洛伦兹力对粒子不做功，粒子总能量
Figure 200987DEST_PATH_IMAGE038
保持守恒，故随着
Figure 570788DEST_PATH_IMAGE036
增大，
Figure 461384DEST_PATH_IMAGE039
将逐渐减小，此时由
Figure 727280DEST_PATH_IMAGE040
可知粒子的横向运动
Figure 336116DEST_PATH_IMAGE041
减小，粒子回旋运动的半径及角速度为
Figure 509608DEST_PATH_IMAGE042
与
Figure 520290DEST_PATH_IMAGE043
，
当带电粒子运动到强磁场区域时，
Figure 488246DEST_PATH_IMAGE044
变小而
Figure 584378DEST_PATH_IMAGE045
变大，带电粒子在非均匀磁场中的运动是回旋半径与轴向速度逐渐减小的螺旋运动；
取柱坐标系，
Figure 889457DEST_PATH_IMAGE046
，粒子速度
Figure 754645DEST_PATH_IMAGE047
，则轴向运动方程为：
Figure 893502DEST_PATH_IMAGE048
对于缓变磁场，在回旋半径粒子的
Figure 476930DEST_PATH_IMAGE049
范围内，可得：
Figure 726646DEST_PATH_IMAGE050
对于正离子，在柱坐标系中
Figure 711919DEST_PATH_IMAGE051
，和负离子
Figure 490520DEST_PATH_IMAGE052
，则可得
Figure 561244DEST_PATH_IMAGE053
式中
Figure 614650DEST_PATH_IMAGE054
为平行磁场方向的梯度，上式表明带电粒子在纵向不均匀磁场中运动时，将受到从强磁场区域指向弱磁场区域的作用力；
磁场由两个平行的载流圆环线圈产生，磁场沿z轴方向先减小后增大，当粒子运动到磁场足够强的区域时，如果有
Figure 454431DEST_PATH_IMAGE055
，即径向动能等于总动能，此时带电粒子的平行向动能
Figure 935090DEST_PATH_IMAGE039
=0，粒子就不能继续沿轴线方向前进，但粒子仍受力
Figure 821007DEST_PATH_IMAGE056
的作用，因此会出现带电粒子运动“反射”的现象，于是粒子在两端强磁场区域之间来回反射；
并非所有的带电粒子都能被反射，当粒子运动到圆环线圈截面处，此时若速度
Figure 678104DEST_PATH_IMAGE057
没有减小为零，则带电粒子从磁镜中穿过，而不会被反射；
单粒子轨道描述法可以绘制出带电粒子在电磁场中的运动轨迹，对于温度和压力较低的等离子体，由于粒子密度比较低，在研究粒子在磁场中的运动时，可以忽略粒子间的相互作用，当然，单粒子轨道描述并不能给出等离子体的集体行为，但可以用于研究粒子在非均匀磁场中的运动规律；
利用COMSOL软件中的带电粒子追踪模块（cpt）和AC/DC模块耦合建立带电粒子在磁场中的运动模型，圆管结构半径15mm，内壁面边界条件设置为反弹，带电粒子在圆管的一端释放，粒子以圆管轴线为中心轴线向圆管内部运动，在磁场（mf）中将圆环结构设置为线圈（coil），线圈类型为圆形，匝数20匝，通入电流10A，在圆管端部释放20000个带电粒子，粒子质量设为电子质
Figure 372391DEST_PATH_IMAGE058
，电荷数设置为-1，粒子以圆柱轴线为中心轴线向通道内部运动，法向速度设为800m/s，粒子进入磁场区后，受到洛伦兹力的作用，对于速度方向不与轴线平行的带电粒子，磁场通过拉莫尔回旋使粒子沿轴线方向运动；
粒子进入磁场区后，受到洛伦兹力的作用，对于速度方向不与轴线平行的粒子，磁场通过拉莫尔回旋使粒子沿轴线方向运动；
S3：模型验证，减压效应是圆管结构中等离子体受到磁场约束作用时的集体行为，宏观上表现为壁面处的气体压强降低，当等离子体处于高温高压状态时，圆管壁面的压强数据较难采集，涉及到在圆管上开孔以及气密性等问题，而当等离子体处于低温低密度时，通常较为容易观测到等离子体在磁场中的运动变化，为了测试磁场对圆管中等离子体气体压强的影响，采用低气压低温等离子体试验平台进行等离子体在磁场中的压强变化试验；
试验平台包括设备台架，石英圆管，等离子体电源、真空计，压力计等，台架用于安装各种设备和仪表，等离子体电源可以在大气压下、低气压下和各种气氛环境中产生稳定的电弧放电、介质阻挡放电、辉光放电，石英圆管由真空泵抽至低压状态，两端电极通过辉光放电产生等离子体，辉光放电的视觉特征呈现均匀的雾状放电；
所述S1中角标
Figure 23952DEST_PATH_IMAGE059
--表示不同的正离子或电子；
Figure 69269DEST_PATH_IMAGE060
——为压强张量；
Figure 464478DEST_PATH_IMAGE061
——为粘性应力张量；
Figure 482112DEST_PATH_IMAGE062
——为摩擦阻力；
Figure 570154DEST_PATH_IMAGE063
——为不同类型离子碰撞引起的热能交换；所述S1中
Figure 837187DEST_PATH_IMAGE064
——约化密度；
Figure 36088DEST_PATH_IMAGE065
——为约华速度；
Figure 439387DEST_PATH_IMAGE066
——为电流密度；所述S1中当各粒子的温度各向同性但不相等时，
Figure 698330DEST_PATH_IMAGE067
；式中
Figure 46135DEST_PATH_IMAGE068
为第
Figure 783147DEST_PATH_IMAGE069
种粒子的质量，
Figure 306532DEST_PATH_IMAGE070
为玻尔兹曼常数，
Figure 736376DEST_PATH_IMAGE071
为第
Figure 712422DEST_PATH_IMAGE072
种粒子的电荷量，
Figure 253125DEST_PATH_IMAGE073
；所述S1中当各粒子温度相等，即热力学平衡状态时，
Figure 365438DEST_PATH_IMAGE074
；
所述S2中分析不同线圈间距粒子运动的影响，改变线圈间距为30mm、20mm，使线圈截面处的磁感应强度不变；所述S2中非均匀磁场中，带电粒子沿轴线方向幅值呈周期性变化，在线圈界面处的强磁场区域时，粒子受到的洛伦兹力较大，而在两线圈之间的弱磁场区域，粒子受到的洛伦兹力较小，由此可推断出等离子体在多圆环线圈组产生的磁场中，在线圈截面附近区域减压效应较好，而两线圈之间区域减压效应不强；
所述S3中实验步骤为调节石英管中的气压，采用进压或抽压的方式实现气体压强的调节，调节范围0-3000Pa；两端电极辉光放电产生等离子体，等离子体电流和功率通过控制台调节；将永磁铁固定于一端空真计附近，通过控制距离调节磁感应强度；待石英管中的气压稳定后，记录控制台中的压强数据，重复上述步骤。
CN201910559066.7A 2019-06-26 2019-06-26 磁场位形对等离子体减压效应的影响研究方法 Active CN110278648B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910559066.7A CN110278648B (zh) 2019-06-26 2019-06-26 磁场位形对等离子体减压效应的影响研究方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910559066.7A CN110278648B (zh) 2019-06-26 2019-06-26 磁场位形对等离子体减压效应的影响研究方法

Publications (2)

Publication Number Publication Date
CN110278648A CN110278648A (zh) 2019-09-24
CN110278648B true CN110278648B (zh) 2022-07-26

Family

ID=67963270

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910559066.7A Active CN110278648B (zh) 2019-06-26 2019-06-26 磁场位形对等离子体减压效应的影响研究方法

Country Status (1)

Country Link
CN (1) CN110278648B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112651194B (zh) * 2020-12-29 2023-03-21 广东电网有限责任公司电力科学研究院 一种等离子体的特性参数空间分布调控方法及装置
CN113643680B (zh) * 2021-08-12 2024-02-09 西北工业大学 多级环-腔耦合模型、声学超表面材料及反射式隐身结构
CN114018326B (zh) * 2021-11-03 2024-04-16 国网湖南省电力有限公司 基于微系统传感器阵列的低压台区环境多参量检测方法
CN114184578B (zh) * 2021-12-03 2023-06-23 西安交通大学 一种可实现瞬态及稳态测量的介质阻挡放电实验平台
CN115124106B (zh) * 2022-06-14 2024-01-19 西安电子科技大学 一种等离子体活化水荷电喷雾装置
CN115100830B (zh) * 2022-07-01 2023-05-16 中国人民解放军国防科技大学 基于等离子体及光敏二极管的高功率微波探测告警平台

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150097487A1 (en) * 2013-10-03 2015-04-09 Daniel Prater Plasma Confinement Device
CN103995923A (zh) * 2014-04-25 2014-08-20 中国人民解放军空军工程大学 用于磁化等离子体的增强型移位算子时域有限差分方法
CN104697397B (zh) * 2015-03-26 2016-06-15 中国人民解放军装甲兵工程学院 一种磁化等离子体火炮
CN106040554A (zh) * 2016-05-25 2016-10-26 东营威玛石油钻具有限公司 一种涂层感应加热方法
CN106057252A (zh) * 2016-06-21 2016-10-26 哈尔滨工业大学 用于多极磁阱磁约束装置的等离子体环向力平衡控制方法
CN106980710A (zh) * 2017-02-28 2017-07-25 上海电力学院 Esp工作电压贡献率与磁感应强度曲线拟合方法
CN109270374B (zh) * 2018-09-14 2019-12-31 电子科技大学 一种监测回旋行波安全工作的方法
CN109858192A (zh) * 2019-03-15 2019-06-07 中国人民解放军陆军装甲兵学院 基于fluent mhd模块的数值仿真计算方法
CN109933911B (zh) * 2019-03-15 2023-04-07 中国人民解放军陆军装甲兵学院 密绕螺线管内金属圆筒电磁场有限元分析方法

Also Published As

Publication number Publication date
CN110278648A (zh) 2019-09-24

Similar Documents

Publication Publication Date Title
CN110278648B (zh) 磁场位形对等离子体减压效应的影响研究方法
Hutchinson Ion collection by a sphere in a flowing plasma: I. Quasineutral
CN110287600B (zh) 磁控等离子体在圆筒内的流动及压力分布研究方法
Olsen et al. Investigation of plasma detachment from a magnetic nozzle in the plume of the VX-200 magnetoplasma thruster
Sukhinin et al. Plasma anisotropy around a dust particle placed in an external electric field
Shang et al. Plasma dynamics for aerospace engineering
Chai et al. Extreme ultra-violet burst, particle heating, and whistler wave emission in fast magnetic reconnection induced by kink-driven Rayleigh-Taylor instability
CN110377963A (zh) 均匀磁场对等离子体减压效应的影响仿真分析方法
Sekerak et al. Hall thruster plume measurements from high-speed dual langmuir probes with ion saturation reference
Imai et al. Deflections of dynamic momentum flux and electron diamagnetic thrust in a magnetically steered rf plasma thruster
Chen et al. Diagnostics of a microhollow cathode discharge at atmospheric pressure
Maiorov et al. On the interaction of microparticles with ion flux in gas discharge plasma
Sung et al. Studies on the optimum condition for the formation of a neutral loop discharge plasma
Shigeta et al. Numerical simulation of a potassium-seeded turbulent RF inductively coupled plasma with particles
Brandt et al. Particle-in-cell simulation of a down-scaled HEMP thruster
Benilov et al. Ion saturation currents to spherical and cylindrical electrostatic probes in collisional plasmas
Hong et al. Particle-in-cell simulation for effect of anode temperature on discharge characteristics of a Hall effect thruster
Haixin et al. Numerical study of the influence of dielectric tube on propagation of atmospheric pressure plasma jet based on coplanar dielectric barrier discharge
Stark et al. Elliptical dust growth in astrophysical plasmas
Urusov et al. Numerical simulation of a helical shape electric arc in the external axial magnetic field
Fu et al. Kinetic simulations of low-pressure inductively coupled plasma: an implicit electromagnetic PIC/MCC model with the ADI-FDTD method
Hu et al. Effects of combined coils on the uniformity of RF ion source
Urusov et al. Mechanisms of the formation of an electric arc of a helical shape in an external axial magnetic field
Fomin et al. Study of MHD-Interactions in Hypersonic Streams
Pogorelov et al. Comment on “On the interaction of the solar wind with the interstellar medium: Field aligned MHD flow” by R. Ratkiewicz and GM Webb

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant