磁场位形对稳定性的影响
在磁约束等离子体中,磁场位形不仅决定粒子轨道与约束性能,也直接决定宏观稳定性。磁流体力学(MHD)把等离子体视为导电流体,其稳定性可由能量原理判断:磁张力通常提供稳定作用,而磁压力、电流梯度与曲率驱动则可能提供自由能。若位形设计不当,微小的扰动会沿磁力线或径向发展成交换模、扭曲模、撕裂模或气球模,导致约束退化甚至破裂。
磁面、安全因子与磁剪切
理想磁约束位形通常由嵌套磁面描述。磁力线在磁面上的扭转程度用安全因子 (q) 表示:
- 在托卡马克中,(q = \frac{\text{环向圈数}}{\text{极向圈数}});
- 有理面满足 (q = m/n),其中 (m,n) 为整数;
- 磁剪切定义为 (s = \frac{r}{q}\frac{dq}{dr}),表示相邻磁面间磁力线的相对扭转。
磁剪切是稳定性的核心参数。强剪切使径向扰动在不同磁面上相位失配,从而抑制理想扭曲模和部分气球模。但剪切过强或电流剖面过峰,也可能在有理面附近驱动电阻撕裂模,形成磁岛。
曲率驱动:交换模与气球模
磁场曲率对稳定性有两面性。若曲率矢量与压力梯度方向不利,即所谓“坏曲率”,等离子体类似重流体位于轻流体之上,容易发生交换不稳定性。托卡马克外侧通常为坏曲率区,内侧为好曲率区。
稳定坏曲率的主要手段包括:
- 磁井:使磁场强度沿磁力线存在极小值,平均最小 (B) 位形可稳定交换模;
- 磁剪切:抑制扰动的径向扩展;
- 有限极向场:在托卡马克中通过极向场产生一定磁井。
气球模是压力梯度与坏曲率共同驱动的局域模,其稳定性阈值常用临界 (\beta) 表示。提高 (\beta) 可提升聚变功率密度,但会逼近气球模极限。因此,位形优化常在约束与稳定性之间折中。
磁剪切与电流剖面
电流剖面决定 (q) 剖面,进而影响扭曲模和撕裂模。最经典的例子是 Kruskal-Shafranov 极限:对于托卡马克,若边界安全因子 (q_a < 1),理想 (m=1,n=1) 扭曲模通常不稳定。因此,基本设计原则是保持 (q_a > 1)。
内部扭曲模与锯齿振荡密切相关。当芯部 (q) 穿过 1 时,(q=1) 有理面内可能发生周期性重联,导致芯部温度骤降。撕裂模则发生在 (q=2,3/2) 等有理面,电阻效应使磁力线重联并形成磁岛。新经典撕裂模由自举电流驱动,可显著降低约束,甚至引发破裂。
边界与三维位形
边界磁场位形同样关键。偏滤器位形中的 X 点可把开放磁力线引向靶板,控制杂质和热流。边缘局域模(ELM)与边界压力梯度和电流密度密切相关。共振磁扰动(RMP)通过施加非轴对称场,在边界形成小磁岛,可主动抑制 ELM,但需避免过度降低约束。
仿星器通过三维磁场设计,可产生磁井、准对称和优化粒子轨道。与托卡马克相比,仿星器无需大等离子体电流,因而避免电流驱动扭曲模,但三维线圈制造和输运优化更复杂。场反位形、紧凑环等则利用不同曲率和剪切组合探索替代路径。
建模要点
在数值建模中,通常先求平衡,再做稳定性分析:
- 用 Grad-Shafranov 方程求轴对称平衡,或用 VMEC 等代码求三维平衡;
- 提取 (q) 剖面、压力剖面、磁剪切和曲率;
- 用理想/电阻 MHD 本征值代码计算增长率与模结构;
- 扫描 (\beta)、(q_a)、电流剖面和边界场,寻找稳定窗口。
例如,在圆柱近似下取 (q(r)=q_0+(q_a-q_0)(r/a)^2)。若 (q_a<1),(m=1,n=1) 扭曲模可能增长;若 (q_a>1),理想扭曲模被稳定,但 (q=2) 有理面仍可能发生撕裂。增加磁剪切通常提高气球模 (\beta) 阈值,但需同时控制有理面位置。
小结
磁场位形通过曲率、磁剪切、磁井、电流剖面和边界三维场共同影响稳定性。设计原则可概括为:避免大范围坏曲率,保持足够 (q_a),合理分布电流与磁剪切,利用磁井和三维优化,并控制边界共振。理解这些几何因素,是等离子体稳定性分析与磁约束装置建模的基础。