Axisymmetric reduction of the adjoint operator in ideal MHD equilibrium

技術要約：理想MHD平衡における随伴演算子の軸対称縮退

問題提起
理想磁気流体力学（MHD）平衡の計算は、磁場閉じ込め装置の設計や診断データからのプラズマ状態の再構成において極めて重要である。軸対称平衡（例：トカマク）はグラッド・シャフラノフ方程式によって支配されるが、三次元（3D）平衡（例：ステラレータ）には数値ソルバーが必要となる。両方の領域において、多数のパラメータに対する感度解析が必要となる最適化手法がますます利用されている。随伴法（Adjoint method）は、このような感度解析のための効率的な枠組みを提供する。しかし、完全な3D平衡問題の随伴定式化と、軸対称グラッド・シャフラノフ方程式の随伴との間の正確な数学的関係は、これまで明示的に示されてこなかった。具体的には、3D随伴演算子が軸対称性の仮定の下でどのように縮退するか、また、その結果得られる構造が、軸対称磁気エネルギーに固有の幾何学的重み因子を考慮した既存の軸対称随伴定式化と一致するかどうかが不明であった。

手法
著者らは、初期段階では対称性を仮定せずに、理想MHD平衡方程式の随伴を導出するために演算子論的アプローチを採用している。

1. 一般的な随伴定式化: MHD平衡は、状態変数として電流密度（$J$）、磁場（$B$）、および圧力（$p$）を持つ制約付き最適化問題として扱われる。平衡制約と目的関数を用いてラグランジアンが構成される。
2. 3D随伴の導出: ラグランジアンを状態変数に関して変分し、固定境界条件（$B \cdot n = 0$）の下でガウスの散度定理を適用することで、著者らは3D随伴方程式を導出する。物理的なラグランジュ変位とは異なる、静的平衡制約に対するラグランジュ乗数としての随伴変数（$\alpha, \beta, \gamma$）が導入される。
3. 軸対称縮退: 著者らは、導出された3D随伴系に対して軸対称性を課す。磁場および随伴場 $\alpha$ を、トロイダル幾何学におけるポロイダル磁束関数（$\psi$）とトロイダル角（$\phi$）を用いて表すことで、発散ゼロのベクトル場の標準的な表現を利用する。
4. 投影と簡略化: 縮退された系を、トロイダル方向に対して平行および垂直方向に投影し、スカラー成分を孤立させる。著者らは特に、磁束的な摂動と解釈されるスカラー量 $\eta = -\alpha \cdot \nabla \psi$ を分析する。
5. 演算子の比較: 得られたスカラー方程式を、摂動を受けたグラッド・シャフラノフ演算子と比較する。また、内積構造、特にポロイダル磁気エネルギーの変化から生じる幾何学的重み $1/R^2$ の役割についても検討する。

主な貢献と結果

• 明示的な縮退: 本論文は、完全な3D同次随伴系の軸対称縮退が、摂動を受けたグラッド・シャフラノフ演算子（Zakharov & Pletzer 1999）と正確に一致するスカラー演算を与えることを示している。この演算子は、平衡磁束関数の微小変化を支配する。
• 変数の対応: 本研究は、3D随伴変数と縮退されたスカラー変数の間の直接的な対応関係を確立している。縮退されたスカラー随伴変数は、力平衡方程式に関連する非圧縮随伴場 $\alpha$ を用いて $\eta = -\alpha \cdot \nabla \psi$ と同定される。
• 重み付き内積構造: 解析により、縮退されたスカラー演算子は、形式的に重み付き内積 $(u, v)_w = \int u v R^{-2} dV$ に関してのみ自己随伴であることが明らかになった。この重み因子 $1/R^2$ は、恣意的な課成ではなく、ポロイダル磁気エネルギーの軸対称表現から自然に生じるものであることが示された。
• 構造的一貫性: 本研究は、内積における幾何学的重みを欠いていた従来の定式化（例：Sanpei et al. 2021, 2023）との相違を解消する。著者らは、重みのない定式化は数値最適化においては実用可能であるが、3D随伴とその軸対称縮退との間の構造的互換性と自己随伴性を維持するためには、重み付きの定式化が必要であることを示している。

意義
本論文は、完全な3D随伴定式化を軸対称グラッド・シャフラノフ随伴問題へと結びつける、厳密な演算子レベルの一致性に関する記述を提供している。軸対称の随伴構造がいかに3Dフレームワークの中に埋め込まれているかを明らかにすることで、本研究は感度解析のための統一的な理論的基礎を提供する。これは以下の事項において特に重要である：

• 弱非軸対称構成: 3D随伴定式化は、軸対称状態の小さな摂動としての平衡を分析するための「親」フレームワークとして機能し、軸対称応答と3D補正を体系的に整理することを可能にする。
• 随伴ベースの再構成と最適化: 3D変数と軸対称変数の明示的な対応関係は、感度結果の解釈を助け、分岐解析や線形化平衡系における特異な応答の特性評価に寄つد可能性がある。

著者らは、この縮退が3Dおよび軸対称の理想MHD随伴問題の間の数学的関係を明確にし、適切な幾何学的重みが適用される場合に、次元を超えて構造的特性（自己随伴性など）が保持されることを保証すると結論付けている。