Augmented Lagrangian methods produce cutting-edge magnetic coils for stellarator fusion reactors

技術要約：ステラーレータコイル最適化のための拡張ラグランジュ法

問題提起
ステラーレータ核融合炉用磁気コイルの設計は、制御された核融合を実現するための決定的なボトルネックとなっている。ステラーレータは（トカマクとは異なり）誘導プラズマ電流を必要とせずに定常運転が可能であるが、その複雑な三次元磁場は、入り組んだ外部コイルシステムを必要とする。これらのコイルの最適化は、不良設定（ill-posed）かつ高度に非凸な逆問題である。従来の設計ワークフローは、プラズマ境界の最適化（ステージI）とコイル形状の最適化（ステージII）を分離しているため、この複雑さに対応することが困難である。

現在の手法は通常、各目的関数や制約条件（磁場精度、コイルの曲率、コイル間距離など）に対して手動で重みを割り当てることで、制約付き最適化問題を非制約の多目的問題へと変換することに依存している。このアプローチには以下のようないくつかの限界がある：

1. 手動チューニング： 重みには物理的な意味がなく、実現可能な解を見つけるために、膨大な計算コストを要する手動のチューニングが必要となる。
2. 局所解： 非凸なランドスケープにより、最適化アルゴリズムが劣劣な構成に陥りやすい。
3. 工学的実行不可能性： 生成されたコイルセットの多くは、過度な曲率、相互連結、あるいは材料の許容範囲を超える力などの理由から実用的ではなく、Wendelstein 7-Xのようなプロジェクトにおける建設の遅延や、NCSXのキャンセルを招いている。

手法
著者らは、これらの課題に対処するために、SIMSOPT最適化スイート内に実装された拡張ラグランジュ（AL）法を提案している。従来の加重和アプローチとは異なり、AL法は制約を動的に扱う。

• 定式化： 最適化は、ラグランジュ関数 $L_A(x, \lambda, \mu) = f(x) - \lambda^\top c(x) + \frac{1}{2} \|\sqrt{\mu} \circ c(x)\|_2^2$ を最小化する。ここで、$f(x)$ は目的関数、$c(x)$ は制約、$\lambda$ はラグランジュ乗数、$\mu$ はペナルティパラメータである。
• 動的更新： 固定されたユーザー定義の重みの代わりに、ラグランジュ乗数 ($\lambda$) とペナルティパラメータ ($\mu$) が、制約の違反に基づいて最適化プロセス全体を通じて自動的に更新される。これにより、手動の重みチューニングが不要になる。
• 制約処理： 二乗磁束誤差（磁場偏差の尺度）を、主要な目的関数ではなく制約として扱う。閾値（例：$f_{SF} \lesssim 10^{-6}$）を設定することで、アルゴリズムは、理論的に完璧だが実用的に不可能なゼロ誤差の解を追い求めるのではなく、磁場精度が十分になった段階で工学的制約を満たすことを優先する。
• 工学的制約： 本フレームワークには、コイルの複雑さに関する標準的な指標（最小コイル間距離およびコイル・表面間距離、総コイル長、曲率制限、平均二乗曲率、連結数（インターリンキング防止）、および電磁力）が含まれている。
• 最適化アルゴリズム： 内側の最小化ステップには、L-BFGS-B クイック・ニュートン法を利用する。勾配は、解析的な微分とベクトル・ヤコビアン積による自動微分を組み合わせて計算される。

主な貢献と結果
論文では、理論的概念から完全に建設された装置に至るまで、5つの異なるステラーレータ構成に対してコイルセットを生成することにより、AL法の有効性を実証している。すべてのケースにおいて、本手法は手動の重みチューニングなしに、以前に発表されたセットと同等またはそれを上回るパレート最適解を生成した。

1. 準軸対称（QA）ステラーレータ：

   • Landreman & Paul の QA 構成に対して、コイルの力と磁場精度の間のパレートフロントをマッピングした。
   • 従来の4コイル解（Wechterm et al.）よりも、コイル数が減ることでリップルが増大するという予想に反して、磁場精度とコイル間距離の両面で優れた1/2フィールド周期あたり3コイルの解を生成した。
   • 総長制約を緩和することで新たなパレート境界を辿り、以前のスキャンよりも一桁優れた磁場精度を達成した。

2. 準ヘリカル対称（QH）ステラーレータ：

   • Landreman & Paul の QH 構成に適用し、物理性能において参照用の5コイル解（Wiedman et al.）に匹敵しながら、工学的指標（例：コイル長の短縮、間隔の拡大）を改善した4コイル解を生成した。
   • 純粋に工学的特性を最適化した第3の構成は、物理的忠実度を犠牲にする代わりに、コイル・表面間距離を38%増加させ、HTS材料コストを10%削減した。これは、製造可能性のために物理性能をトレードオフできる能力を示している。

3. Stellaris (SQuID) 構成：

   • Proxima Fusion の Stellaris コンセプトに対し、1/2フィールド周期あたり5および6コイルの代替案を見出した。
   • これらの解は、オリジナルの設計と比較して、磁場誤差を25%改善し、コイル間間隔を13 cm広げると同時に、HTSターンへのピーク力を29%減少させた。
   • 工学的特性を最適化した構成は、推定HTSコストを6%削減し、幾何学形状を大幅に単純化した。

4. Wendelstein 7-X (W7-X)：

   • 既存の W7-X 原子炉に対する代替コイルセットを生成した。
   • 構成 #1： コイル数を50個から40個（1/2フィールド周期あたり5個から4個）に減らしつつ、磁場精度を向上させ、総コイル長を短縮した。
   • 構成 #2： 総コイル長をわずか5%増加させるだけで、電磁力を20%削減した。
   • 構成 #3： 長さや力にほとんど変化を与えることなく、磁場誤差と曲率を劇的に改善した。
   • すべての代替案は、オリジナル設計よりも低い曲率と滑らかな幾何学形状を示した。

5. Helically Symmetric eXperiment (HSX)：

   • 既存の HSX デバイスに対し、オリジナルの6コイル設計と比較して、総コイル長を25%、曲率を45%削減しつつ、妥当な磁場精度を維持する4コイル解を生成した。これは、研究規模の装置における製造の簡素化の可能性を強調している。

意義と主張
本論文は、拡張ラグランジュ法が以下の点において、ステラーレータコイル設計の最先端を根本的に進化させると主張している：

• 重みチューニングの自動化： 現在のワークフローにおける非効率性と劣劣な結果の主な原因である、ヒューリスティックな手動の重み調整への依存を排除した。
• 堅牢性： 多様な磁気構成（QA, QH, QI）において、非凸な最適化ランドスケープをナビゲートし、実現可能なパレート最適解を見出すことに成功した。
• 汎用性： 理論的概念から既存の製造済み原子炉（W7-X, HSX）まで、幅広い適用性を示し、確立された工学的ソリューションを改善できることを証明した。
• 効率性： 以前の確率的または加重和を用いた手法で必要とされた大規模なパラメータスキャン（例：8,500回の実行）よりも、大幅に少ない計算オーバーヘッドでこれらの結果を達成した。

著者らは、この手法が、プラズマ物理学の要件と、原子炉規模の統合に必要な厳格な工学的制約とのバランスを取るための、計算効率が高く堅牢な経路を提供すると結論付けている。