等离子体模拟器 返回
电磁学与等离子体

等离子体模拟器

实时模拟洛伦兹力 F = q(E + v×B) 驱动的带电粒子运动。可视化回旋运动、E×B漂移及等离子体参数空间。

场参数设置

统计信息

计算结果
1.00
ωc(电子)
1.00
拉莫尔半径
0.50
E×BVelocity
平均动能
Visualization
理论与主要公式

$$\omega_p = \sqrt{\frac{n_e e^2}{\varepsilon_0 m_e}}$$

電子プラズマ振動数:\(n_e\) 電子数密度 [m⁻³]、\(e\) 電荷、\(m_e\) 電子質量

$$\lambda_D = \sqrt{\frac{\varepsilon_0 k_B T_e}{n_e e^2}}$$

デバイ長(遮蔽長):電場が遮蔽される特性長さ [m]

$$\mathbf{F} = q(\mathbf{E} + \mathbf{v}\times\mathbf{B})$$

ローレンツ力:荷電粒子の運動方程式の右辺

什么是等离子体模拟器

🙋
这个模拟器里,那些飞来飞去的小光点是什么呀?它们好像被看不见的力推着走。
🎓
简单来说,那些是模拟的带电粒子,比如电子和离子。它们之所以会拐弯、加速,是因为受到了“洛伦兹力”。在实际工程中,比如在核聚变托卡马克装置里,就是用强大的磁场来约束这些高温粒子,不让它们撞到容器壁上。你试着拖动右边的“磁场强度B”滑块,看看粒子们的运动轨迹会发生什么变化?
🙋
诶,真的吗?我把磁场调大,它们转的圈圈立刻变小了!这是为什么?
🎓
这是因为磁场力只改变粒子的运动方向,不改变它的速度大小。粒子速度不变,但磁场越强,让它拐弯的力就越大,所以转的圈子就越小。这个圈子的半径叫做“拉莫尔半径”。工程现场常见的是,在粒子加速器里,我们就是通过精确控制磁场,让粒子沿着设计好的环形轨道高速运动的。你再把“粒子Velocityv”的滑块调小看看,会发现半径也会跟着变小。
🙋
哦!那旁边计算出来的“等离子体频率”又是什么意思?听起来好像它们自己在振动?
🎓
你的感觉很准!简单来说,如果把等离子体里的电子群想象成一块弹簧,等离子体频率就是这块弹簧自己振动的固有频率。比如在电离层(我们头顶的大气层),无线电波能不能反射回来,就和这个频率密切相关。改变参数后你会看到,当你把“电子密度n”滑块往上拉,这个频率会显著增大,这意味着等离子体对电磁场的响应更快了。

物理模型与关键公式

粒子运动的根本驱动力是洛伦兹力,它决定了带电粒子在电磁场中的每时每刻的加速度。

$$ \vec{F}= q(\vec{E}+ \vec{v}\times \vec{B}) $$

其中,$q$是粒子带电量,$\vec{v}$是粒子速度,$\vec{E}$是电场强度,$\vec{B}$是磁感应强度。$\vec{v}\times \vec{B}$是矢量的叉乘,结果是一个始终垂直于速度和磁场的力。

在均匀磁场中(忽略电场),粒子会做匀速圆周运动,其回旋半径(拉莫尔半径)和角频率由以下公式给出:

$$ r_L = \frac{m v_\perp}{|q| B}, \quad \omega_c = \frac{|q| B}{m} $$

$r_L$是拉莫尔半径,$m$是粒子质量,$v_\perp$是垂直于磁场的速度分量。$\omega_c$是回旋频率(角频率),它只取决于粒子的质荷比和磁场强度,与速度无关。

现实世界中的应用

磁约束核聚变(如托卡马克):利用强大环形磁场约束高达上亿摄氏度的等离子体,使氢同位素发生聚变反应。控制粒子的拉莫尔半径和漂移运动是实现稳定约束的关键。

粒子加速器与同步辐射光源:通过交变电场加速带电粒子,并用精确设计的磁场(二极磁铁、四极磁铁)引导和聚焦粒子束,使其在环形真空管道中高速运行。

空间等离子体与航天器防护:地球磁层、太阳风都是天然的等离子体。研究等离子体与磁场的相互作用,有助于预测太空天气,并设计保护卫星免受高能粒子轰击的措施。

等离子体加工与材料表面处理:在半导体工业中,利用低温等离子体进行刻蚀、镀膜和清洗。通过调节电场、磁场和气体密度,可以精确控制工艺效果。

常见误解与注意事项

首先,请明确本模拟器观测的是单粒子运动。实际工程中处理的等离子体是海量粒子的集合体,粒子间的碰撞与集体行为(如等离子体振荡)具有本质重要性。虽然通过此工具可以学习到E×B漂移“与粒子种类无关”,但实际等离子体中会产生由密度梯度和磁场曲率引起的其他漂移,这些漂移会显著增加约束难度。

其次,“实时”显示的局限性。模拟速度取决于计算负荷。例如当粒子质量极轻(如电子)或磁场较强时,回旋频率 $\omega_c$ 会变得极大,可能导致数值计算无法跟上而显示卡顿,或引发数值不稳定问题。在实际工程CAE中,如何稳定求解这类“刚性问题是关键所在。

最后,初始条件的重要性。本工具中粒子的初始位置和速度通常采用固定或随机设置。但在核聚变反应堆模拟等场景中,需要设定粒子按特定分布函数(如麦克斯韦分布)生成。切勿仅通过观察单个粒子轨迹就过早断言“等离子体会如此运动”。

使用指南

  1. 设置磁感应强度B(0-2T):调整Bslider或直接输入BvalNum数值,控制粒子的回旋半径r=mv/(qB)
  2. 配置电场分量:通过ExvalNum/Exslider设置Ex(-1000~1000 V/m),EyvalNum/Eyslider设置Ey,观察E×B漂移速度vd=(Ex×By)/(B²)
  3. 设定粒子数量N(1-500):选择NvalNum或Nslider确定仿真规模,点击运行按钮启动洛伦兹力驱动的动力学计算

具体计算示例

质子在强磁场中的运动:B=1.5T,Ex=500V/m,Ey=0V/m,粒子质量mp=1.67×10⁻²⁷kg,电荷q=1.6×10⁻¹⁹C。初速度v₀=10⁶m/s垂直于B,回旋频率fc=qB/(2πm)≈2.3×10⁷Hz,回旋半径r=mv₀/(qB)≈0.7cm,E×B漂移速度vd=Ex/B≈3.3×10⁵m/s。可视化显示螺旋轨迹叠加向y方向的系统漂移

实务注意事项

  1. 等离子体诊断中:B>0.5T时电子与离子轨迹分离明显,适合模拟托卡马克磁约束参数(典型值B=2-3T)
  2. E×B剪切:当Ex在500V/m量级变化时,漂移速度梯度产生剪切层,影响等离子体不稳定性增长率
  3. 粒子碰撞:N>300时计算资源占用显著增加,建议分段仿真多个粒子群而非单次运行大数量
  4. 边界条件:设定粒子反弹或周期性边界,防止轨迹逃逸导致数值不稳定