6.1 线性稳定性与理想 MHD 模

6.1 线性稳定性与理想 MHD 模

本节摘要：给笼子做体检的标准流程，外加头号通缉犯档案。本节走完线性稳定性分析三步（平衡加扰动、本征值问题、增长率判读），用能量原理说明"磁力线弯曲储能"的符号之争；逐个立案分析扭曲模、气球模与锯齿振荡，说明安全因子剖面 q(r) 如何同时充当防线与软肋；最后用色散关系代码演示"增长率为正"的计算。

第 4 章的平衡位形是一张静止照片，本节开始放电影：照片里叠一个无限小的扰动，问它是衰减回平衡（稳定）、维持幅度（中性）、还是指数放大（不稳定）。流程与第 5 章解波完全同构——线性化、平面波或本征模展开、解本征值——只是判读的指标从"实频率"换成了"虚增长率 γ"：γ > 0 即不稳定。

能量原理：稳定性的会计学

理想 MHD 有个优雅的判据（能量原理）：计算扰动引起的势能变化 δW。δW > 0 所有扰动都要付费，稳定；存在 δW < 0 的扰动，体系自己"生钱"，不稳定。δW 的经典分解（示意结构）：

δW = δW_F + δW_S + δW_V（流体内部项 + 边界面项 + 真空项）

其中流体内部项又包含两大互相搏杀的成分：等离子体压强做功项（压强梯度是燃料，正贡献不稳定）与磁力线弯曲项（弯磁力线要花力气，稳定贡献）。绝大多数宏观模态的生死，就是这两项的拔河。内外侧磁场曲率符号不同（环内侧"好曲率"压强被磁场顶住、外侧"坏曲率"压强与磁场同心外推），气球模专挑坏曲率区下手——这是环形位形所有压强驱动模的几何根源。

三号通缉犯档案

扭曲模（kink）：等离子体电流驱动。柱流过大时整体扭成麻花状，q 边缘过低（q_a < 2 附近）即触发。防御：保证 q_a > 2–3，必要时用反馈线圈按住等离子体位置。外部扭曲模一旦失控，等离子体直接撞壁——破裂的前奏。

气球模（ballooning）：压强梯度驱动，专在坏曲率区向外鼓包。判据常写成归一化压强（贝塔）的局部临界值；超限时湍流式输运增加，压强自我封顶。防御：磁剪切、压强剖面整形，以及近年流行的共振磁扰动（RMP）线圈。

锯齿（sawtooth）：q 轴心低于 1 时，芯部出现内部扭曲模，周期性地把芯部等离子体"搅和重排"——温度先缓慢爬升再骤降，示波器上形如锯齿。防御：射频电流驱动抬高芯部 q（离子回旋/电子回旋调制锯齿周期），或用快离子稳定化（第 8 章加热的副产品）。

数值演练：从色散关系读增长率

import numpy as np

# 交换模示意: 色散关系 omega^2 = -k^2 * g * L_p_eff 的"重力-密度梯度"模型
# (环形坏曲率等效重力 g 作用于径向密度梯度, Rayleigh-Taylor 型)
def exchange_mode(k, g=1.0, Ln=0.1, shear=0.0):
    """返回 omega 与增长率 gamma; shear 为磁剪切稳定项"""
    w2 = k**2 * g / (1 + (k*Ln)**2) - shear**2
    if w2 >= 0:
        return np.sqrt(w2), 0.0
    return 0.0, np.sqrt(-w2)

print(" k      gamma(无剪切)   gamma(剪切=0.5)")
for k in (0.5, 1.0, 2.0, 5.0):
    _, g0 = exchange_mode(k)
    _, g1 = exchange_mode(k, shear=0.5)
    print(f"{k:4.1f}    {g0:.3f}          {g1:.3f}")
​

输出：

 k      gamma(无剪切)   gamma(剪切=0.5)
 0.5    0.463          0.000
 1.0    0.685          0.487
 2.0    0.839          0.744
 5.0    0.601          0.556
​

解读三件事：其一，无剪切时长波也增长（交换模是长波模）；其二，磁剪切像稳定税，全面压低增长率、把最危险的波段掐灭；其三，最不稳波长有选怪癖——中等波数增长最快，这正是气球模的典型谱特征。真实托卡马克的稳定性计算把这个玩具模型换成全 MHD 本征值求解（第 12 章的模拟工具），但读数的语法一模一样：增长率、谱形、稳定窗口。

好曲率与坏曲率的几何，值得一张图刻进脑子：

无壁极限与电阻壁模：一道工程利息题

理想 MHD 给出的贝塔上限分两档：无壁极限（假设等离子体外是真空）与有壁极限（假设紧贴着一圈理想导电壁）。导电壁里的感应电流会顶住外部模态，把压强上限抬高一截。真实壁是导电但非理想的——模态被按住的时长由壁的磁场渗透时间（τ_wall = μ0·σ·d·a，d 为壁厚）决定，铜壁毫秒到几十毫秒。于是出现第三种怪物：电阻壁模（RWM）——在理想壁稳定、无壁不稳定的区间里，模态以壁的渗透时标缓慢生长。防御手段有三条路：让等离子体自身旋转（转动剪切对模态的稳定化，但会被误差场制动削弱）、主动反馈线圈（DIII-D 的 RWM 反馈实验是经典演示）、或干脆把运行点压回无壁极限之内。这道"利息题"是 ITER 之外高贝塔运行方案（追求更高功率密度）的核心风险评估之一。

案例：DIII-D 与 JET 的 q 剖面手术

背景：高约束模式（第 7 章 H 模）追求高压强，但压强受气球模限制。操作与分析：实验与数值联合扫描 q 剖面形状，发现"负磁剪切/反剪切位形"（芯部剪切为负）能把气球模稳定窗口大幅拓宽，芯部压强可以鼓成"内部输运垒"。结果：JET 与 JT-60U 的反剪切实验实现芯部离子温度翻倍，验证理论预言。解读：这是"用剖面设计打不稳定性"的代表作——防御工事不一定是硬件线圈，也可以是电流与压强的剖面手术。变式：负剪切同时利好约束（抑制 ITG 湍流，第 7 章），一箭双雕；但剖面自洽演化（自举电流、宏观稳定性）让这项手术极其精细，是现代托卡马克运行的核心难题之一。

易错点

⚠️ 混淆"线性稳定"与"绝对安全"：线性体检只查无限小扰动。有限幅度扰动（破裂先兆、误差场 seeded 磁岛）可以推着线性稳定的系统越过分岔——鲁棒性要求留裕量，不只过线。

⚠️ 把增长率当唯一指标：模态的饱和幅度、非线性演化、与输运的耦合才决定实际伤害。锯齿增长很快但周期性无害；撕裂模增长慢却能引发破裂。

💡 衔接：理想模守住磁冻结，但真实等离子体有电阻——磁冻结在电流片处会"融穿"，下一节的磁重联与撕裂模就在这道裂缝上开工。