含跨场流的弛豫磁流体力学（RxMHD）稳态可解性条件与环形几何中 m=1 至 m=2 磁岛转变

《Journal of Plasma Physics》：Relaxed magnetohydrodynamics with cross-field flow

　　

Dewar 等（Phys. Plasmas 27, 062507, 2020）的相空间拉格朗日模型（phase-space Lagrangian model）提供了一种将跨场流（cross-field flow）纳入弛豫平衡（relaxed equilibria）并同时保持理想磁流体力学（ideal magnetohydrodynamics, MHD）力平衡的模型。在此，研究人员刻画了稳态解空间，并识别出一个可解性条件（solvability condition），该条件通过度规张量（metric tensor）将规定的约束流（constrained flow）与几何耦合。利用该条件，研究人员在板（slab）、柱（cylindrical）和环形（toroidal）几何中构建了平衡。在环形几何中，跨场流与磁岛结构（magnetic-island structure）强烈相关：改变角频率会修正磁场径向分量的主傅里叶谐波（dominant Fourier harmonic），并能够驱动从主（m=1）岛向次（m=2）岛的转变。在板和柱几何中，流参数对岛宽影响较弱，但强烈修正平衡剖面（equilibrium profiles）。

研究背景：Kruskal & Kulsrud 的变分模型曾用于构建三维（3-D）磁流体力学（MHD）平衡，但理想 MHD 的冻结磁通约束使该模型在非轴对称几何中于有理面出现平行电流奇异性，且不允许磁岛形成。Taylor 弛豫、Finn & Antonsen、多区域弛豫 MHD（multiregion relaxed MHD, MRxMHD）以及阶梯压力平衡代码（stepped-pressure equilibrium code, SPEC）发展出可允许磁拓扑变化的弛豫模型。然而，早期 Taylor 弛豫与 MRxMHD 未包含等离子体流；Finn–Antonsen 通过环向角动量拉格朗日乘子引入流，但该约束仅在轴对称系统中有物理意义。Dewar 等提出基于相空间哈密顿作用量原理的 RxMHD 模型，可自然容纳跨场流，并以磁螺旋、交叉螺旋和熵为宏观约束，且其 Euler–Lagrange 方程可恢复理想 MHD 力平衡。但该理论的稳态问题仍欠定，计算上需要额外结构。研究人员因此针对二维板、柱和环形几何求解 RxMHD 稳态方程，识别可解性条件，并研究跨场流对平衡剖面与磁岛几何的影响。研究结论表明，跨场流可决定性影响磁岛结构，在环形几何中可由 m=1 主岛转变为 m=2 次岛；该工作为将更一般的流剖面纳入 SPEC 提供了依据。

关键技术方法：研究人员采用 Dewar 等的相空间拉格朗日 RxMHD 框架，在二维曲线坐标 {s,θ,ζ} 中令 ?_ζ=0，以度规张量编码板、柱和环形几何。稳态 Euler–Lagrange 方程与广义 Beltrami 方程、运动方程、连续性方程联立；引入纯环向约束流 v^ζ，利用可解性条件选择可容许流型。数值求解使用 Dedalus v3，tau-spectral 方法，s 方向 Chebyshev 展开、θ 方向 Fourier 展开，并进行力平衡残差收敛测试。无样本队列来源。

2 作用量原理与 Euler–Lagrange 方程：研究人员从 MHD 哈密顿密度出发，加入交叉螺旋、磁螺旋和熵三个宏观约束，构造相空间拉格朗日量并变分，得到 Euler–Lagrange 方程。结论是总流 u 为约束流 v 与弛豫场向流 u^Rx 之和，且由这些方程可导出理想 MHD 力平衡 ρu·?u=-?p+(?×B)×B；模型不强制理想 Ohm 定律，因而允许磁拓扑变化。

3 二维平衡与纯环向约束流 v
3.1 坐标系统：研究人员定义曲线坐标 {s,θ,ζ}，二维假设使所有场量不依赖 ζ，并由边界线性插值构造位置矢量。通过场线积分获得 Poincaré 图。结论是二维假设限制扰动为环向模数 n=0，唯一可能共振为 ι=0，高阶 m 谐波只改变共振层结构而不产生新岛链。
3.2 广义 Beltrami 方程：研究人员用矢势和规范条件处理广义 Beltrami 方程，边界条件给出总的极向和环向磁通 ψ_θ、ψ_ζ，并推导出 μ_BA_ζ+νu_ζ-B_ζ=C 以及平行 Alfvén Mach 数 M_∥^Rx。结论是该式需与 ζ 分量方程联立求解。
3.3 运动方程：基于新经典考虑，研究人员设约束流仅有环向分量 v=v^ζe_ζ，并分解总流的平行与跨场分量。结论是磁岛 O 点处 B^s=B^θ=0，跨场流局部消失；同时推导出可解性条件 ?v^ζ×?u_ζ^Rx+v^ζ?v^ζ×?g₃₃=0，该条件限制可容许的 v^ζ。
3.4 数值方法与收敛图：研究人员使用 Dedalus v3、tau-spectral 方法，在 s 方向用 Chebyshev 展开、θ 方向用 Fourier 展开求解。收敛测试表明力平衡残差先指数下降，后平台至约 10^-12，说明数值解可靠；同时指出 Dedalus v3 不支持自定义坐标，且依赖直接矩阵求逆，不能处理欠定系统。
4 HKT 板几何与柱几何中的示例解
4.1 HKT 板几何：研究人员构造 Hahm–Kulsrud–Taylor（HKT）板几何，其中 g₃₃=1，可解性条件给出 v^ζ=F(νu_ζ^Rx)，并在线性近似下取 v^ζ=αu_ζ^Rx+ω。数值解显示旋转变换在 s≈0.3 穿过 ι=0 并形成岛，磁场和流范数在岛位置最大，跨场流在 O 点消失。
4.2 柱几何：研究人员构造空心柱几何，仍有 g₃₃=1，因此方程形式类似。结果表明旋转变换呈凸剖面，岛移向 s≈0.17；柱向曲率使 s=0 附近磁场和流范数升高，跨场流比降低，O 点跨场流仍消失。
4.3 α 对跨场流比、平衡剖面和磁岛的影响：研究人员令 ω=0 并改变 α。结论是 α=0 时跨场流为零，类似带流的 MRxMHD；α 增大时跨场流比增大，压力、流和磁场剖面的变化增强，且 α→∞ 时 u_⊥/u_∥→|B_p/B_ζ|。旋转变换变化较小，岛宽对 α 不敏感，主要受体边界几何决定；在 α=-1、ω=0 时环向流可消失，流变为纯极向，这与未强制 ?×(u×B)=0 有关。
5 环形几何与 Finn–Antonsen 平衡
研究人员构造同心环几何，g₃₃=R²，取 v^ζ=ω 以满足可解性条件，得到广义 Bernoulli 方程，并发现其等价于 Finn–Antonsen 平衡，因此 Finn–Antonsen 平衡是本文模型的一个子类。该结果还与 Maschke & Perrin 的旋转平衡理论联系，密度和压力出现离心增强，不再在磁面上为常数；但本文模型不假设全局磁面叶理，因此允许磁岛。
5.1 角频率 ω 对跨场流比和磁岛几何的影响：研究人员扫描 ω 从 -3 到 3。结论是 ω=-3 时存在大主岛，位置约 s≈0.35；增大 ω 后主岛缩小，随后出现次岛，ω=-0.65 时主岛消失而仅余次岛，ω=-0.5 时次岛合并并在 θ=0 附近重建主岛；继续增大 ω 后主岛增长。B^s 的一阶谐波 b_m=1^s 随 ω 增大而减小，二阶谐波 b_m=2^s 近似不变，并在 -1?ω?-0.4 区间主导。跨场流比与主岛宽度相关。近似模型将 b_m=1^s 化为平移绝对值函数，其平移与 ψ_θ 相关；ψ_θ=0 时的小差异被归因于大纵横比和一阶展开。
6 结论与讨论
讨论部分总结：MHD 平衡问题，尤其含流时，通常欠定，需要额外结构；本文通过预设约束流形式提供该结构。该模型允许显著环向流而不完全抑制极向分量，并重现磁岛 O 点跨场速度消失的实验一致行为。模型本身不直接包含新经典或微湍流效应。未来方向包括全三维推广、沿准对称方向指定约束流、发展多区域 RxMHD 并与 SPEC 集成，以及与 Voigt 正则化 MHD 或电阻撕裂模饱和比较。

研究结论（译自原文）：在这项工作中，研究人员展示了求解 Dewar 等提出的“半弛豫”平衡模型的若干方法，探索了不同几何中的解空间，并表明等离子体流与几何内在耦合。特别地，研究人员证明等离子体流，尤其当前弛豫 MHD 模型所允许的跨场流，能够决定平衡剖面和 ι=0 磁岛几何。在板与柱几何中，增大 v^ζ 与 u_ζ^Rx 之比会增强跨场流并陡化平衡梯度，但 ω 不显式出现在控制方程中，因而不影响磁压剖面。在环形几何中，ω 的影响显著：在大多数可达范围内，增大 |ω| 导致主岛加宽；但在小负 ω 时趋势逆转，主岛缩小并最终分裂为两个次岛；继续增大 |ω| 时次岛合并，重新形成主岛。该非单调磁岛动力学可由 B^s 的一阶傅里叶谐波解释。与理想 MHD 不同，RxMHD 允许无需显式电阻的磁拓扑变化。