Parameter Scan of Multi-Fluid Equilibria in Rotating p-11B Plasmas: Effects on Fusion Power and Bremsstrahlung Losses

技術要約：回転するp-11Bプラズマにおけるマルチフルイド平衡のパラメータスキャン

問題提起
プロトン・ホウ素11（p-11B）核融合反応は、中性子が少ないという性質から魅力的であるが、クーロン障壁を克服するために必要なイオン温度（数百keV）が非常に高いという大きな課題に直面している。これらの温度域では、相対論的制動放射が支配的なエネルギー損失チャネルとなり、従来の熱的シナリオでは核融合出力が制動放射を上回ることが多い。さらに、これらの損失を軽減するために非平衡プラザ状態を維持するには、運動量交換に伴う電力コストが発生する。p-11Bプラズマにおいては、プロトンとホウ素イオンの質量差（$m_B \approx 11m_p$）が大きいため、単一フルイドモデルでは不十分である。既存のマルチフルイド平衡モデルの多くは、種依存的なトロイダル回転を含む複雑なパラメータ空間を効率的にスキャンする能力に欠けている。本論文は、回転する球状トカマク（ST）構成において、規定された種依存的回転が密度再分布、制動放射損失、および核融合出力にどのように影響するかを分析するための、堅牢かつ高速なスペクトル・フレームワークの必要性に対処するものである。

手法
著者らは、規定された種依存的トロイダル回転を持つ二次元軸対称マルチフルイド平衡のための高速スペクトル・パラメータスキャン・フレームワークであるVEQ-MFを開発した。この手法は、以下の3つのコアコンポーネントを統合している：

1. 支配方程式： モデルは、一般化されたボルツマン密度応答（遠心力を考慮）、準中性静電分極、および一般化されたグラッド・シャフラノフ（GS）方程式を結合させる。ここでは、トロイダル回転が支配的であり、磁力面内は等温であり、電子は質量ゼロの極限を仮定し、トランジション・シングラリティを避けるためにポロイダル流の慣性を無視している。
2. 数値ソルバー： VEQ-MFは、簡略化されたパラメータGS表現およびVEQスペクトルソルバーを拡張したものである。これは二層の非線形構造を採用している。内ループではベクトル化されたニュートン・ラフソン法を用いて局所的な準中性を解き、外ループではシフトされたチェビシェフ多項式とハイブリッド・ブロイデン法を用いてグローバルなGS演算子を解く。回転周波数を段階的に増加させるために、ホモトピー継続法が用いられる。
3. 出力評価： ポストプロセッシングでは、相対的なトロイダル流を考慮したドリフト・マックスウェル型反応率係数（$\langle \sigma v \rangle_{DM}$）を用いて核融合出力（$P_{fus}$）を計算し、相対論的な解析的放射モデルを用いて制動放射出力（$P_{brem}$）を計算する。主要な指標は、制動放射のみによる出力比 $R_{fb} = P_{fus}/P_{brem}$ である。
4. 検証： ソルバーは高解像度の有限差分法（FDM）ソルバーと比較検証され、スペクトルパラメータおよび格子細分化に関する収束チェックが行われた。

主な貢献

• VEQ-MFフレームワーク： 規定された種依存的トロイダル回転を持つ結合マルチ種回転平衡を、高い計算効率（MATLABでの非ゼロケースで1点あたり0.032–0.050秒）で扱うことができるスペクトルソルバーの導入。
• 診断的分解： 密度再分布（幾何学的重なり）とドリフト誘起の反応率係数の両方が全核融合出力に与える寄与を分離する新しい手法。
• パラメータスキャン： 2つのSTジオメトリ（EHL-2: $R_0=1.1$ m および EHL-3B: $R_0=3.2$ m）について$(\Omega_p, \Omega_B)$平面の包括的なスキャンを行い、微分回転が核融合と制動放射のバランスを変化させる領域を特定した。

結果

• 等回転（$\Omega_p = \Omega_B$）： 等回転下では、重いホウ素イオンは遠心力により低磁場側（LFS）に強く蓄積する。これにより静電分極が誘起され、電子が外側にシフトすることで、体積積分された $n_e^2$ が増加する。その結果、制動放射損失は核融合出力よりも速く増大し、$R_{fb}$ は回転数の増加とともに著しく減少する（例：EHL-2では約0.40から0.10未満へ低下）。
• 微分回転（$\Omega_p \neq \Omega_B$）：
  • 制動放射の低減： プロトンに対するホウ素の回転周波数（$\Omega_B$）を減少させると、遠心力による分極が弱まる。これにより電子密度プロファイルが平坦化され、$n_e^2$ の積分値が減少し、$P_{brem}$ が低下する。
  • 核融合出力の回復： より大きなEHL-3Bジオメトリでは、相対的なトロイダル流が顕著なドリフト運動エネルギー（$E_d$）を生成する。これが相対速度分布をシフトさせ、ドリフト・マックスウェル型反応率係数 $\langle \sigma v \rangle_{DM}$ を増大させる。
  • 競合する効果： 診断的分解により、微分回転が増大するにつれて2つの競合する傾向があることが明らかになった：(1) プロトンとホウ素の密度重なり積分が単調に減少し（核融合ポテンシャルを減少させる）、(2) 密度重み付きの有効反応率係数が着実に上昇する。EHL-3Bでは、反応率係数の上昇が最終的に密度の損失を相殺し、$P_{fus}$ が極小値から回復することを可能にする。
• 出力比の交差： EHL-3Bにおいて、制動放射の低減と核融合反応性の向上の相乗効果により、制動放射のみによる出力比 $R_{fb}$ は、最大の相対トロイダル流（$\gamma_{Bp} \approx 0$）において0.18（等回転時）から約1.3まで上昇する。EHL-2では、$\langle \sigma v \rangle_{DM}$ を有意に変えるほどのドリフトエネルギーが不足しているため、効果は制動放射の低減に限定される。
• 数値的安定性： ベンチマークの結果、貯蔵エネルギー、制動放射、および核融合出力において、VEQ-MFとFDMの間で1.7–3.4%の全域的な差が見られた。収束チェックは、スペクトルおよび格子の細分化に対する感度が1%未満であることを示している。

意義および主張
本論文は、その目的がp-11B核融合の原子炉レベルの実現可能性を評価することではないことを明示している。むしろ、パラメータスキャン可能なマルチフルイド平衡計算をテストするための、要求水準の高いモデルケースを提供することにある。その意義は、回転するマルチフルイドプラズマにおける密度重なりと反応率係数の寄与を、正常に分離できるベンチマーク済みの平衡スキャン・ワークフローを提供したことにある。

著者らは、報告された $R_{fb} = 1$ の交差が、規定されたプロファイル下での「減少した制動放射のみの出力比」の結果であることを強調している。これには、補助動力、輸送損失、および（最大のケースで推定される約5.2 GWの運動量交換に必要な）大幅な動力コストは含まれていない。したがって、これらの結果は、規定された微分回転によって生じる平衡の結果を特徴づけるものであり、そのような状態を維持するための輸送制限的なプロセスを記述するものではない。本研究は、これらの平衡計算を、プロファイル進化とエネルギーコストの輸送レベルの解析と結合させる将来の研究のための基礎を築くものである。