VEQ: a fast parametric Grad--Shafranov solver for fixed-boundary tokamak equilibria with flexible source profiles

本論文は、柔軟なソースプロファイルと調和展開を利用することで、輸送幾何学結合ワークフローにおける繰り返しのクエリに適した低遅延かつ高精度な平衡再構成を実現する、固定境界トカマク平衡のための高速なパラメトリックソルバーであるVEQを導入する。

要約

トカマク（ドーナツ型の核融合炉）を、超高温のプラズマで満たされた巨大で見えない風船だと想像してみてください。この風船が破裂したり崩壊したりするのを防ぐために、科学者たちは、それをつなぎとめている磁気の「皮膚」がどのような形をしているか、そして内部の力がどのようにバランスしているかを正確に知る必要があります。これは「平衡（エクリブリウム）」を見つけると呼ばれます。

通常、この平衡を計算することは、触れるたびに形が変わる巨大な3Dジグソーパズルを解くようなものです。それは正確ですが、非常に時間がかかります。もし、原子炉を制御するために1秒間に数百回も形状をチェックする必要があるなら、その時間は長すぎます。

この論文は、これらの計算の「早送りボタン」となるように設計された新しいツール、VEQ（Veloce EQuilibrium）を紹介しています。その仕組みを、簡単な比喩を使って説明します。

1. 「シェイプシフター」対「ピクセル絵師」

従来の計算手法は、ピクセル絵師のようなものです。彼らは、数百万もの微細なグリッド点（ピクセル）における値を計算することで、磁場全体をマッピングしようとします。それは詳細ですが、重くて時間がかかります。

VEQは、より**シェイプシフター（変身体）**に近いものです。すべてのピクセルを塗る代わりに、いくつかの主要なつまみ（パラメータと呼ばれるもの）で構成された、柔軟な数学的「骨格」を用いてプラズマの形状を記述します。

• プラズマの形状を粘土の塊だと考えてください。
• 従来のソルバーは、粘土の上のあらゆる小さな凹凸を彫刻しようとします。
• VEQは、あらかじめ定義された「引き伸ばし」や「曲げ」のツール（数学的な高調波や多項式）を使用します。あなたは単に、粘土を正しい形にするために、いくつかのダイヤルを回すだけです。ピクセルを塗るよりも回すべきダイヤルの数が圧倒的に少ないため、この計算は驚異的に速くなります。

2. 「ユニバーサル・トランスレーター（万能翻訳機）」

核融合研究における最大の悩みの一つは、異なるコンピュータプログラムが異なる言語を話していることです。あるプログラムは「圧力」を、別のプログラムは「電流」を、また別のプログラムは「安全係数」（安定性の指標）を出力するかもしれません。

VEQは、ユニバーサル・トランスレーターとして機能します。論文では、VEQには6つの異なる「入力ルート」（USBポートのようなもの）があることが示されています。これら異なるソースからのデータを入力すれば、VEQはそれらをすべて独自の内部言語へと翻訳し、問題を解決します。

• 主張： 著者らは、同じ完璧なデータを6つのポートすべてに投入してテストを行いました。その結果、どのポートを使用しても、VEOは全く同じ結果を生成しました。これは、異なる入力ケーブルを使用したからといって、翻訳機がエラーを引き起こしたりしないことを証明しています。

3. 「速度と精度のトレードオフ」

この論文は、単に「速い」と言っているだけではありません。「スケッチ」と「写真」の選択のように、速度と精度の間のトレードオフについても探求しています。

• スケッチ（低パラメータ）： ダイヤルを非常に少なく使います。一瞬（ミリ秒単位）で済み、素早い確認には十分ですが、細かいディテールを見逃す可能性があります。
• 写真（高パラメータ）： 多くのダイヤルを使います。少し時間はかかりますが（それでも非常に速く、複雑な形状でも約19ミリ秒）、高い精度で形状を捉えます。
• 結果： 著者らは、3種類の異なるプラズマ形状（標準的な「D型」、高性能な「Hモード」、そして複雑な形状の「X点」を持つもの）でテストを行いました。その結果、少ないダイヤル数であっても、VEQは複雑な形状を、目に見えないほど小さな誤差（原子炉のサイズの0.2%未満）で再現できることがわかりました。

4. 「ストレス・テスト（限界の検証）」

著者らは、限界についても正直に述べています。力のバランス（プラズマを保持する張力）が完璧な場所と、そうでない場所をチェックしました。

• 内部： プラズマの内部において、VEQは極めて優秀です。力はほぼ完璧にバランスしています。
• エッジ（縁）： 最外殻の境界付近では、誤差がわずかに大きくなります。これは、VEQが滑らかで柔軟な骨格を使用しているのに対し、実際のプラズマの境界はギザギザしていたり、鋭い角（X点のようなもの）を持っていたりするためです。
• 教訓： VEQは原子炉の「中央部」には最適です。もし境界の極めて正確な状態を知る必要があるなら、より遅い、より詳細なツールでダブルチェックする必要があるかもしれません。しかし、ほとんどの制御タスクにおいては、VEQは十分に高速かつ正確です。

5. 「輸送（トランスポート）テスト」

最後に、VEQによる形状の小さな誤差が、熱がプラズマを通じてどのように移動するかを予測しようとする際（「輸送」テスト）、大きな問題を引き起こすかどうかをテストしました。

• 結果： 誤差は極めて小さく（1%未満）、実用上の問題はありませんでした。これは、少し曲がった定規で部屋を測ったとしても、その測定ミスが「どのサイズのラグを買うべきか」という決定に影響を与えないようなものです。

まとめ

VEQは、核融合炉のための新しい、超高速な計算機です。すべての点をマッピングする代わりに、柔軟な数学的骨格を使用してプラズマの形状を記述します。

• 高速： 複雑な形状をミリ秒単位で解決できます。
• 柔軟： 多様なソースからのデータを受け入れます。
• 信頼性： 原子炉の内部に対してよく機能し、最も外側のエッジに注意を払えば、ほとんどの制御システムにおいて十分に正確です。

著者らは、VEQが「今、プラズマはどのような形をしているのか？」という問いを何度も繰り返す必要があるシステム、例えばリアルタイムの原子炉制御や、最適な動作条件を見つけるための何千回ものシミュレーション実行に最適であると結論づけています。

技術概要

技術要約：VEQ – 固定境界トカマク平衡のための高速パラメトリック Grad–Shafranov ソルバー

問題提起
トカマクの統合モデリング・ワークフロー（IMAS、RAPTOR、SuperCodeなど）では、平衡幾何学、計量係数、および診断プロファイルの繰り返し的なクエリが必要となる。これらのワークフローは、一度限りの高忠実度再構成とは異なる。それらが求めているのは、パラメータ・スキャンや制御指向のイテレーションに適したコンパクトな表現でありながら、単純な解析的閉鎖を超えて有用となる、十分な2次元成形性と力学的バランスの情報を保持できる表現である。既存のツールは、多くの場合、二極化した領域に位置している。低次の工学的閉鎖（Millerモデルなど）は高速だが残留診断機能に欠け、一方で完全な平衡ソルバー（CHEASE、EFIT、DESCなど）は、高忠実度のネイティブ出力を提供するものの、内ループでの使用には計算コストが高すぎるか、あるいはインフラストラクチャが重すぎる。したがって、連続的な固定境界表現でありながら、有用な2次元成形性と残留診断を保持しつつ、再生コストを低く抑えられる手法へのギャップが存在する。

手法
本論文では、この中間領域を埋めるために設計されたコンパクトなパラメトリック・フレームワークである VEQ（Veloce EQuilibrium）を紹介する。ベンチマークされた具体的な実装である VEQPy は、軸対称の固定境界 Grad–Shafranov ソルバーである。

• パラメトリック表現: 平衡は、パラメトリックなフラックス面写像 $(\rho, \theta) \to [R(\rho, \theta), Z(\rho, \theta)]$ によって定義される。幾何学は、MXH型フラックス面調和関数（歪んだ極角）とシフトされたチェビシェフ多項式による径方向プロファイルを用いてパラメータ化されている。能動的な未知数は、グリッドベースのポロイダル磁束マップの節点値ではなく、これらの調和関数およびプロファイルの係数である。
• 変分定式化: ソルバーは、変分的に誘導された投影 Grad–Shafranov 残差を強制する。強形式の 방정식을 各点で解く代わりに、この手法は、能動的な係数によって生成される許容される変分に制限を課す。これにより、形状ブロックが変分構造を持つ投影（ペトロフ・ガルツキン条件）となり、プロファイル・ブロックがモーメント投影によって閉じられる、有限次元系が得られる。
• 柔軟な入力ルート: 様々な入力（解析的研究、再構成ファイル、制御モデル）との互換性を確保するため、VEQは6つの入力ルート（PF, PP, PI, PJ1, PJ2, PQ）をサポートしている。これらのルートは、圧力勾配、トロイダル磁場関数、ポロイダル磁束勾配、囲まれたトロイダル電流、電流密度、または安全因子（q）プロファイルを受け取る。ルート固有の閉鎖処理により、これらのヘテロジニアスな入力を共通の正規化されたソース表現へとマッピングし、すべてのルートが同一の有限次元残差演算子へと供給されるようにする。
• 実装: ソルバーは、非線形ソルバー（SciPyのPowellハイブリッド法）と加速されたバックエンド（Numba）を使用している。セットアップ（グリッド依存の前計算、演算子データ）と解法フェーズを分離することで、低レイテンシの繰り返しクエリを可能にしている。

主な貢献

1. 有限次元パラメトリック表現: コンパクトなMXH–チェビシェフ・フラックス面表現の定義、および（フラックス面の移動に伴う対流変分を含む）変分的に誘導された投影 Grad–Shafranov 残差の導出。
2. ルート固有のソース・クロージャ: ヘテロジニアスな1次元入力（圧力、電流、安全因子など）を共通の残差系へとマッピングするフレームワークの導入。これにより、入力ルートの一貫性と幾何学的近似誤差を区別している。
3. デュアル診断アプローチ: 最小化された投影残差ノルム（ソルバーによって最小化されるもの）と、サンプリングされた点ごとの強形式の Grad–Shafranov 残差の両方を報告する。これにより、簡約系の収束と、点ごとの力学的バランス・スクリーニングを区別している。
4. ベンチマーク: Solov'ev、CHEASE、および EFIT 由来の G-EQDSK ケースを用いた包括的なテストを行い、入力ルートの一貫性、フラックス面再構成誤差、精度とコストのトレードオフ、および繰り返し解法のレイテンシを定量化した。

結果

• ルートの一貫性: 共通のリファレンス平衡から生成された、滑らかで相互に互換性のある入力を用いた制御されたテストにより、すべての6つの入力ルートが、数値的なリカバリパス（積分/微分ステップなど）に起因する差異はあるものの、同一の固定境界パラメトリック平衡に収束することが示された。
• 精度とレイテンシ: ソルバーは、3つの代表的なケース（D型 Solov'ev、Hモード CHEASE、および X点 EFIT）でテストされた。
  • D型: 9個の能動パラメータを用いて、マイナー半径正規化形状誤差 $1.4 \times 10^{-3}$ を達成し、ソルブのみのメディアン時間は 1.6 ms であった。
  • Hモード: 65個のパラメータを用いて、誤差 $1.1 \times 10^{-3}$ を達成し、メディアン時間は 19 ms であった。
  • X点: 94個のパラメータを用いて、誤差 $1.9 \times 10^{-3}$ を達成し、メディアン時間は 15 ms であった。
• 残留診断: 点ごとの強形式の診断により、能動的な表現を豊かにすることは、主に内部の力学的バランスを改善することが明らかになった。HモードおよびX点のケースでは、グローバルなRMSおよび最大残差値は依然として境界付近の寄与に支配されており、これは固定境界の仮定が、内部の形状変化のみでは境界適合の不一致や特異構造（X点など）を完全に補償できないことを示している。
• 輸送結合: 孤立した1次元輸送・幾何学結合テストにおいて、温度プロファイルの応答は、VEQ の幾何学誤差に対して約1%未満に留まった。

意義と主張
本論文は、境界の精緻化、局所的な補正、またはより高精度なソルバーを必要とするケースをスクリーニングするために、点ごとの診断を保持することを条件として、VEQ が統合モデリングおよび制御ワークフローにおける繰り返しの平衡幾何学クエリに適していると主張している。

著者らは、VEQ を高精度な再構成ソルバーの普遍的な代替品ではなく、補完的なものとして位置づけている。その主な価値は、低コストで再生可能な、明示的な投影残差とサンプリングされた診断を備えた、連続的かつコンパクトな表現を提供することにある。論文では、示されたタイミングの結果はセットアップ後の「ソルブのみ」のレイテンシを反映したものであり、前処理、I/O、および離散化を含むフルパイプラインのタイミングの置き換えとして解釈すべきではないと明記している。本手法は、解析的閉鎖とフルソルバー・ネイティブ・パイプラインとの間のギャップを埋めるためのものであり、幾何学的な豊かさと計算効率のバランスを提供するものである。