Axisymmetric reduction of the adjoint operator in ideal MHD equilibrium

技术摘要：理想 MHD 平衡中伴随算符的轴对称归约

问题陈述
理想磁流体动力学（MHD）平衡计算是磁约束装置设计以及从诊断数据重建等离子体状态的基础。虽然轴对称平衡（如托卡马克）由 Grad–Shafranov 方程控制，但三维（3D）平衡（如仿星器）需要数值求解器。两种机制都越来越多地采用基于优化的方法，其中需要针对众多参数进行灵敏度分析。伴随方法为这种灵敏度分析提供了一个高效的框架。然而，三维全量平衡问题的伴随形式与轴对称 Grad–Shafranov 方程的伴随形式之间，尚未建立起明确的数学关系。具体而言，目前尚不清楚在假设轴对称的情况下，三维伴随算符如何归约，以及其产生的结构是否与已建立的轴对称伴随形式一致，特别是关于轴对称磁能中固有的几何权重因子问题。

方法论
作者采用算符理论方法，在最初不假设对称性的情况下，推导了理想 MHD 平衡方程的伴随。

1. 通用伴随形式： 将 MHD 平衡视为一个受约束的优化问题，其状态变量为电流密度 ($J$)、磁场 ($B$) 和压力 ($p$)。利用平衡约束和目标泛函构建拉格朗日量。
2. 三维伴随推导： 通过对拉格朗日量关于状态变量求变分，并在固定边界条件 ($B \cdot n = 0$) 下应用散度定理，作者推导出了三维伴随方程。引入伴随变量 ($\alpha, \beta, \gamma$) 作为对静态平衡约束的拉格朗日乘子，这些变量与物理拉格朗日位移不同。
3. 轴对称归约： 作者将轴对称性施加于推导出的三维伴随系统中。他们利用环形几何中无散矢量场的标准表示法，将磁场和伴随场 $\alpha$ 表示为极向磁通函数 ($\psi$) 和环向角 ($\phi$) 的函数。
4. 投影与简化： 系统沿环向和横跨环向进行投影，以隔离标量分量。作者专门分析了标量量 $\eta = -\alpha \cdot \nabla \psi$，将其解释为一种类通量扰动。
5. 算符比较： 将得到的标量方程与扰动的 Grad–Shafranov 算符进行比较。作者还检查了内积结构，特别是源于极向磁能变化的几何权重 $1/R^2$ 的作用。

主要贡献与结果

• 显式归约： 本文证明了全量三维齐次伴随系统的轴对称归约产生了一个标量算符，该算符与扰动的 Grad–Shafranov 算符（Zakharov & Pletzer 1999）完全一致。该算符控制着平衡通量函数的无穷小变化。
• 变量对应关系： 研究建立了三维伴随变量与归约标量变量之间的直接对应关系。归约标量伴随变量被识别为 $\eta = -\alpha \cdot \nabla \psi$，其中 $\alpha$ 是与力平衡方程相关的不可压缩伴随场。
• 加权内积结构： 分析阐明了归约标量算符仅在加权内积 $(u, v)_w = \int u v R^{-2} dV$ 下是形式自伴随的。该权重因子 $1/R^2$ 被证明自然地源于极向磁能的轴对称表达式，而非人为强加的。
• 结构一致性： 该工作解决了先前文献（如 Sanpei 等人 2021, 2023）中遗漏内积中几何权重的问题。作者表明，虽然未加权的公式在数值优化中是切实可行的，但加权公式对于保持三维伴随与其轴对称归约之间的结构兼容性和自伴随性是必要的。

意义
本文提供了一个严谨的算符级一致性陈述，将全量三维伴随形式与轴对称 Grad–Shafranov 伴随问题联系起来。通过澄清轴对称伴随结构是如何嵌入三维框架中的，这项工作为灵敏度分析提供了统一的理论基础。这对于以下方面尤为重要：

• 弱非轴对称构型： 三维伴随形式可以作为分析轴对称状态微扰的“父”框架，从而允许系统地组织轴对称响应和三维修正。
• 基于伴随的重建与优化： 三维与轴对称变量之间的显式对应有助于解释灵敏度结果，并可能有助于分叉分析或表征线性化平衡系统中的奇异响应。

作者得出结论，这种归约澄清了三维与轴对称理想 MHD 伴随问题之间的数学关系，确保了在应用适当几何权重时，其结构特性（如自伴随性）能够在不同维度间得以保持。