Development of a Reduced Multi-Fluid Equilibrium Model and Its Application to Proton-Boron Spherical Tokamaks

本論文は、陽子・ホウ素（p-$^{11}$B）核融合を目指す球状トカマクにおいて、回転による種分離や静電ポテンシャルといった多流体効果を考慮した低次元平衡モデルを開発し、実験・炉規模の構成における回転速度が平衡状態に与える影響を明らかにしたものです。

要約

タイトル： 「超高速回転する『究極のエネルギー工場』を設計するための、新しい計算ルール」

1. 背景：夢のエネルギー「プロトン・ボロン」

今、世界中の科学者が「究極のクリーンエネルギー」を探しています。その筆頭候補が**「プロトン（水素の仲間）」と「ボロン（ホウ素）」**をぶつけ合わせる核融合です。

これまでの主流（重水素・トリチウム）と違って、この方法は「放射能のゴミ（中性子）」がほとんど出ない、まさに「夢のエネルギー」です。しかし、これには**「ものすごい高温」と、それを閉じ込めるための「強力な磁石の檻（磁気閉じ込め）」**が必要になります。

2. 課題：回転しすぎて「中身がバラバラ」になる問題

このエネルギーを作る装置（球状トカマク）は、燃料を効率よく加熱するために、中身を猛烈なスピードで**「回転」**させます。

ここで問題が発生します。燃料には「軽いプロトン」と「重いボロン」の2種類が混ざっています。
これを、**「回転するメリーゴーラウンド」**に例えてみましょう。

• 軽いプロトンは、小さな子供のようなものです。回転しても、真ん中の方に留まろうとします。
• 重いボロンは、大きな荷物を持った大人です。遠心力で、外側の端っこ（メリーゴーラウンドの縁）にグイグイ押し流されてしまいます。

これまでの計算モデル（単一流体モデル）は、「中身はみんな一緒の液体だ」と決めつけて計算していました。しかし、実際には**「軽い奴は真ん中、重い奴は外側」と、燃料がバラバラに分かれてしまう**のです。これでは、燃料同士がぶつかり合えず、エネルギーがうまく作れません。

3. この研究のすごいところ：新しい「賢い計算モデル」の開発

研究チームは、この「バラバラ現象」を正確に予測できる、新しい計算方法（縮小多流体平衡モデル）を開発しました。

これまでの完璧すぎる計算モデルは、あまりに複雑すぎて、スパコンを使っても計算が止まってしまったり、時間がかかりすぎたりしていました（例えるなら、天気予報をするのに、地球上のすべての空気の粒子の動きを計算しようとして、計算が終わらないような状態です）。

今回のモデルは、**「大事なところ（回転による分離と、電気的なバランス）だけをピンポイントで計算する」**という賢いやり方です。これにより、エンジニアが設計図を作る際に、実用的かつ正確に「燃料がどう動くか」をシミュレーションできるようになりました。

4. 何がわかったのか？

研究チームが、実験用と将来の発電用（EHL-2とEHL-3B）の2つの装置モデルで実験したところ、驚きの事実がわかりました。

• 回転が速すぎると、燃料が「三日月形」に偏る： ボロンが外側に追い出され、中心部がスカスカになってしまうことが判明しました。
• 電気が発生する： 軽いものと重いものが分かれることで、装置の中に数千〜数万ボルトという、巨大な「電気の壁」が勝手にできてしまうことがわかりました。

5. まとめ：未来への一歩

この研究は、単に「計算ができました」というだけではありません。

「燃料がバラバラになって、エネルギーが作れなくなる！」という失敗を、実際に装置を作る前にコンピュータ上で予見できるようになったことが最大の成果です。

この「賢い計算ルール」を使えば、将来、人類がクリーンで無限のエネルギーを手に入れるための「設計図」を、より正確に、より早く描けるようになるのです。

技術概要

技術要約：プロトン・ホウ素球状トカマクにおける低減マルチ流体平衡モデルの開発

1. 背景と課題 (Problem)

プロトン・ホウ素（$p\text{-}^{11}\text{B}$）核融合は、中性子放出を抑えたクリーンなエネルギー源として期待されていますが、実現には極めて高いイオン温度（>100 keV）と強力な磁場閉じ込めが必要です。この条件を満たす候補として、高ベータ（$\beta$）特性を持つ球状トカマク（ST）が挙げられます。

しかし、STでは中性粒子入射（NBI）による強力なトロイダル回転が誘発されます。$p\text{-}^{11}\text{B}$プラズマでは、陽子（$m_p$）とホウ素イオン（$\sim 11m_p$）の間で質量差が非常に大きいため、遠心力が重いホウ素イオンを低磁場側（LFS）へ押し出す**種分離（Species Separation）**が発生します。これにより電荷の偏極が生じ、静電ポテンシャル（$\Phi$）が発生します。
従来の単一流体MHDモデルでは、この種分離やそれに伴うポテンシャル構造を記述できず、一方で、極めて詳細なマルチ流体モデルは数値的な不安定性（数学的な硬さ）や計算負荷が大きく、実用的な設計ツールとしては不十分であるという課題がありました。

2. 研究手法 (Methodology)

本研究では、物理的な忠実度と計算の堅牢性のバランスを取った**「低減マルチ流体平衡モデル（Reduced Multi-fluid Equilibrium Model）」**を開発しました。

• モデルの簡略化戦略:
  • 採用: 強力なトロイダル回転による慣性と、準中性条件を維持するための自己整合的な静電ポテンシャル（$\Phi$）の結合。
  • 省略: 数値的な特異点（音速共鳴）を避けるため、ポロイダル流の慣性と圧力異方性を無視。これにより、支配方程式を楕円型（Grad-Shafranov方程式の拡張版）として維持。
• 数学的定式化:
  • 各種のベルヌーイ関係式（遠心力ポテンシャルを含む）と、準中性条件（$\sum q_s n_s = 0$）を組み合わせた、一般化されたGrad-Shafranov方程式を構築。
• 数値アルゴリズム (NFEQ):
  • 有限差分法（FDM）を用い、ポロイダル磁束（$\psi$）を解く外側のピカール反復ループと、静電ポテンシャル（$\Phi$）を解く内側のニュートン・ラフソン法を組み合わせた入れ子状反復スキームを開発。
  • 全電流（$I_p$）を一定に保つためのフィードバック制御機構を実装。

3. 主な貢献 (Key Contributions)

1. 実用的な理論基盤の確立: 複雑なマルチ流体物理を、エンジニアリング的な設計（パラメータスキャン等）に耐えうる計算コストで扱えるモデルとして定式化。
2. 数値コード「NFEQ」の開発: 高回転・高質量差の条件下でも数値的に安定して収束する、自己整合的な平衡計算アルゴリズムの提供。
3. $p\text{-}^{11}\text{B}$特有の物理現象の解明: 回転がプラズマの種分布とポテンシャルに与える影響を定量的に示した。

4. 結果 (Results)

ENNグループが設計した2つのデバイス（実験機EHL-2および炉型EHL-3B）への適用により、以下の知見が得られました。

• マッハ数による挙動の変化: イオンのマッハ数（$M$）が閾値（$M \approx 0.5$）を超えるとマルチ流体効果が顕著になる。
• 種分離の発生: 高回転領域（$M > 2$）では、重いホウ素イオンが低磁場側（LFS）に強く蓄積し、三日月型の密度分布を形成する。これにより、コア部での燃料（陽子とホウ素）の密度ミスマッチが生じ、核融合出力密度に影響を与える可能性がある。
• 大規模な静電ポテンシャル: 炉型EHL-3Bでは、種分離に伴い約10～15 kVという極めて大きな静電ポテンシャルが生成される。これは、マイクロ不安定性の抑制に関わる$E \times B$シアーなどの輸送現象に重要な影響を与える。

5. 意義 (Significance)

本研究は、$p\text{-}^{11}\text{B}$核融合炉の設計において、単一流体モデルでは見落とされる「燃料の空間的偏り」や「大規模な電位構造」を予測するための不可欠なツールを提供しました。
また、このモデルは$p\text{-}^{11}\text{B}$に限らず、標準的なD-T核融合における重金属不純物（タングステン等）の遠心力による偏析や、同位体分離、高度トカマクモード（ATモード）における輸送障壁の解析など、幅広い核融合研究に応用可能な汎用性を備えています。