实时模拟洛伦兹力 F = q(E + v×B) 驱动的带电粒子运动。可视化回旋运动、E×B漂移及等离子体参数空间。
粒子运动的根本驱动力是洛伦兹力,它决定了带电粒子在电磁场中的每时每刻的加速度。
$$ \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中,如何稳定求解这类“刚性问题是关键所在。
最后,初始条件的重要性。本工具中粒子的初始位置和速度通常采用固定或随机设置。但在核聚变反应堆模拟等场景中,需要设定粒子按特定分布函数(如麦克斯韦分布)生成。切勿仅通过观察单个粒子轨迹就过早断言“等离子体会如此运动”。
质子在强磁场中的运动: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方向的系统漂移