Wiki / 現象 プラズマ乱流 （Plasma turbulence）は、核融合・プラズマ物理学の現象。磁化プラズマ中で生じる非線形な揺動現象で、密度や電位などの空間・時間的な乱流的変動を指す。核融合装置では辺縁部からコア部まで広く現れ、乱流による異常輸送が閉じ込め性能を左右する。計測やジャイロ運動論的シミュレーションによりその性質が調べられている。この用語集が追う論文のうち、タイトルにこの主題が現れるものは187本ある。比重が最も大きかったのは1975–1979年で、その期間の論文の0.47%を占めた。直近の2025–2029年では0.19%。
目次
1 解説 2 推移 3 研究の変遷 4 よく登場する装置 5 この数字の作り方 解説 プラズマ乱流とは、磁場で閉じ込められたプラズマ中に自然に発生する、広い空間・時間スケールをもつ密度、温度、電位の激しい揺らぎと、それが引き起こす非線形な輸送現象のことである。熱核融合の文脈では、この乱流が磁面を横切る粒子と熱の拡散を引き起こし、閉じ込め性能の上限を決める中心的な要因となる。乱流そのものは流体乱流と似た統計性を持つが、荷電粒子が磁力線に沿って動き、密度勾配や温度勾配がドリフト波 型の不安定性を駆動するという、流体にはない固有の機構に支配されている。
プラズマ乱流の性質と駆動源
磁場による閉じ込めでは、粒子は磁力線の周りをらせん運動しながら、系に存在する圧力勾配や磁場の曲がりをエネルギー源として、摂動が自発的に成長する。トカマクやステラレータ の高温プラズマで最も重要な駆動源は、イオン温度勾配 に起因する不安定性と、電子温度勾配 や捕獲電子の効果による不安定性であり、いずれも勾配が臨界値を超えると成長する。乱流が発達すると、渦(エディ)やフィラメントと呼ばれる磁力線に沿った構造が形成され、それらが径方向に動くことで粒子と熱が閉じ込め領域から外へ運ばれる。このため乱流輸送 の理解は、単に揺らぎを観測することではなく、どの不安定性がどのスケールでエネルギーを供給し、それがどのように小さいスケールへ移って散逸するかという、エネルギーカスケードの問題として扱われる。
乱流の様子を具体的に描くには、数値シミュレーションで得られる三次元スナップショットが役に立つ。閉じ込め領域の最外側磁面と、そこに沿って伸びる一本の磁力線を重ねた密度の三次元分布を見ると、乱流構造が磁面に沿ってどのように分布し、磁力線の向きとどう関係しているかが直接読み取れる。
fig.2 Global fluid simulation of plasma turbulence in a stellarator with an island divertor (2022) · CC BY 4.0
Three-dimensional snapshot of plasma density, representing the geometry of the simulation considered here. The surface depicted corresponds to the LCFS and a field line (black line) is superimposed.
プラズマ乱流
種別 現象
論文数 187 本
初出 1970 –1974 年（1本）
最盛期 1975 –1979 年 0.47%
直近 2025 –2029 年 0.19% 索引にもどる
図は、閉じ込め領域の境界をなす最外閉磁面の三次元像に、一本の磁力線(黒線)を重ねたもので、面上の色はプラズマ密度の揺らぎを表す。密度の高い構造と低い構造が磁面の周囲に不均一に分布し、それらが磁力線に沿って伸びている様子がわかる。このように乱流構造は等方的ではなく、平行方向には長く、磁面横断方向には短いという強い異方性を持つ。これがプラズマ乱流の幾何学的な基本像である。
さらに、乱流は単一のスケールで閉じているわけではない。イオン スケールの乱流と電子スケールの乱流が互いに結合し、ゾーナルフローと呼ばれる帯状の電位構造を介してエネルギーをやり取りすることで、熱損失が単一スケールの理論から予測される量を上回り得ることが、イオンと電子の両スケールを同時に解像するシミュレーションによって示されている。このような跨スケール結合は、実験で観測される電子熱輸送 の説明に不可欠と考えられている。
核での乱流輸送と閉じ込め
核プラズマでは、乱流による横方向の熱輸送 がエネルギー閉じ込め 時間を直接制限する。重要な性質の一つは「剛性」である。温度勾配を臨界値を少し超えても大きく上げると、熱流が急激に増えるため、温度分布を自由に選べず、勾配が臨界値近くに張り付くように振る舞う。このため核の温度プロファイルは、加熱分布よりも乱流の性質によって強く規定される。電子熱損失が理論的予測を上回る問題は長年「大きな未解決問題」と呼ばれてきたが、跨スケールのジャイロ運動論 的シミュレーションが実験のイオン・電子両方の熱流を再現できるようになり、電子スケール乱流との結合が有力な説明として浮かびつつある。
輸送の理論的な扱いでは、乱流による粒子の動きを速度空間の拡散として定式化するフォッカー・プランク 的な記述や、リウヴィル変換を用いた輸送理論などが発展してきた。また、乱流がいつ、どのようにして衝突レスな温度勾配レスな温度勾配駆動不安定性の閾値を越えて遷移するかを分岐理論として扱う研究もある。これらはいずれも、乱流輸送 係数を経験則ではなく第一原理から決めようとする試みである。
境界とスクレイプオフ層の乱流: ブロブと排熱
閉じ込め領域の外側、磁力線が開いて壁に接続されるスクレイプオフ層 (SOL)では、乱流の振る舞いが質的に変わる。ここでは分離線近傍で生まれた高密度の構造(ブロブ)が、SOLを横切って壁方向へ弾道的に飛翔する。ブロブは径方向速度が毎秒数キロメートルに達することもあり、SOLでの粒子・熱輸送 の主要な担い手である。ブロブの生成位置は分離線内側の圧力勾配が最大となる領域と対応し、ペデスタル の構造がブロブの放出量を決めるため、境界乱流は核の閉じ込め状態とも結び付いている。
SOL乱流の最も重要な帰結が、ダイバータ (ダイバータ板)への熱負荷 の幅を決めることである。熱流の径方向減衰長 λq が狭いほど、同じ排熱をより小さな面積に集中させなければならず、壁材料の限界との兼ね合いで炉設計の制約になる。シミュレーションと実験の比較から、現在の装置では磁場ドリフトが主に熱流幅を決める「ドリフト支配」の領域にあるのに対し、ITER のような大型炉では機械が大きくなるほどドリフトによる径方向輸送が減り、逆に乱流による熱拡散が増えるため、「乱流支配」の領域へ移行することが示されている。つまり、現在の装置の熱流幅の経験則をそのまま大型炉へ外挿することはできず、境界乱流の予測精度が将来炉のダイバータ設計を直接左右する。
装置の種類によっても様相が異なる。ステラレータ では、W7-X の島ダイバータ 構成における観測で、フィラメントの径方向速度が毎秒百メートル程度とトカマクのブロブより一桁以上遅く、構造がほぼ磁面に張り付いて留まることがわかった。揺らぎの統計がほぼ正規分布に従うこととも合わせ、ステラレータ境界の乱流は外から運ばれた構造ではなく局所で生まれるものが支配的であると解釈されている。一方で、1〜10 kHz の広帯域スペクトルや、10〜50 kHz の準コヒーレントモードの存在など、トカマクと共通する特徴も多い。
抑制と制御: せん断流、断面形状、磁場摂動
乱流を抑制して閉じ込めを改善する、あるいは境界での熱負荷 を分散させるという二つの目的から、乱流の能動制御が研究されてきた。
第一の手段はせん断流である。径方向に変化する電場が作る E×B ドリフトのせん断は、乱流渦を引き伸ばして破壊する働きがあり、せん断率が十分大きいと乱流レベルと乱流粒子束が抑制される。実験では、プラズマ中に電極を挿入してせん断率を放電ごとに変え、その抑制のスケーリングを測定することができる。抑制はせん断率の二次関数として表され、乱流と粒子束の間の位相差(クロスフェーズ)が粒子束の抑制に鍵を握ること、つまり揺らぎの大きさだけでなく揺らぎと輸送の位相関係が制御されていることがわかった。この機構は、L モード から H モード への遷移で境界に局在した強せん断層が形成され乱流が一斉に静まるという、閉じ込め改善の中心的な描像の基礎を与える。
第二の手段はプラズマ断面の制御である。トカマクの断面の三角形度 (三角形状)は乱流に系統的な影響を与える。負の三角形度(δ<0)ではイオン温度勾配 モードや捕獲電子モードが安定化され、L モードのまま H モードに匹敵する閉じ込めが得られることが、TCV をはじめとする装置の実験とシミュレーションで示されてきた。境界では、負の三角形度が SOL 乱流を安定化し、熱流減衰長 λq を狭める方向に働く。これは閉じ込め改善という利点と、排熱面積を減らすという欠点が同じ乱流抑制から生じることを意味し、将来炉での負三角形度 運用を考える上で境界乱流の精密な理解が不可欠である。実際、負三角形度と正三角形度の磁気平衡を用いた非線形シミュレーションでは、両者の間で境界乱流の構造と強度に明瞭な差が現れる。
図は、負三角形度(左)と正三角形度(右)のプラズマについて、平衡電子圧力で規格化した電子圧力揺らぎのスナップショットを断面で示したものである。分離線(破線)の内側と外側で揺らぎの空間構造がどう分布し、三角形度の向きによって乱流の模様がどう変化するかが読み取れる。断面形状という平衡のパラメータ一つで乱流の様相が変わることは、乱流が「固定された性質」ではなく制御対象であることを示す好例である。
第三の手段は三次元磁場摂動 である。共振磁場摂動(RMP )は、コイルで外部から小さな非軸対称磁場を加え、境界にカオス的な磁力線領域を作ることでエッジ局在モード(ELM )を抑制する手法である。同時に RMP は境界乱流にも影響する。シミュレーションでは、RMP 印加によって密度揺らぎや乱流輸送 が減少する例が複数のモデルで報告されており、その機構は平衡の密度・電場プロファイルが変わることによる安定化と説明される。一方で、RMP は磁力線に沿った径方向輸送経路も作るため、ダイバータターゲット への熱負荷分布への影響は単純でない。ダイバータ 構成でのグローバルな二流体シミュレーションでは、RMP 印加がトロイダル平均した熱流を減らす方向に働くが、その効果はまだ控えめであり、コイル最適化の余地が大きい。いずれにせよ、ELM 制御 という本来の目的と乱流・排熱への副作用を同時に扱う視点が不可欠である。
シミュレーションと診断の方法
プラズマ乱流の研究は、数値シミュレーションと実験診断の両輪で進められている。
シミュレーションでは、対象とする領域によって用いるモデルが変わる。核では、渦運動を高速回転に平均化したジャイロ運動論 が標準であり、イオン・電子スケールを同時に解くマルチスケール計算が行われている。境界では衝突性が高いため、ドリフト近似したブラギンスキー流体方程式を用いるのが主流で、GBS のようなグローバル・フラックス駆動コードは、核からの粒子・熱の流入、乱流輸送 、壁での損失の三者の釣り合いを時間発展として解く。この種のコードは、単純化された線形装置や小型トーラスでの検証(バリデーション)を経てトカマク境界へ拡張され、さらにステラレータ の非軸対称磁場へも適用が広がっている。ステラレータの島ダイバータ シミュレーションでは、トカマクとは異なり低い極方向モード数のコヒーレントなモードが輸送を支配し、しかもその輸送が高磁場側で強くなるという、単純な曲率駆動の直観に反する結果が得られ、幾何学を取り込んだ線形理論によって説明が与えられている。
境界シミュレーション特有の難しさは幾何学にある。プラズマのダイナミクスは磁力線方向に強く異方性であるため、計算格子を磁力線や磁面に沿わせるのが望ましいが、その一方で壁は磁力線に沿っていないため、複雑な第一壁 形状を境界条件として正確に扱うことが難しい。磁場に沿った座標系、磁場に依存しない局所整合格子、ペナルティ法による埋め込み境界条件など複数の戦略が競合しており、近年はスプラインに基づく境界適合格子を用いて、バッフル(ガス排気を高めるダイバータ の囲い込み)を含む現実の壁形状をそのまま計算に取り込む手法が開発されている。ダイバータの閉じ込みは中性圧力を上げて脱離着(dettachment)への移行を容易にし、ターゲット熱負荷 を下げる有力な手段とされており、壁形状を忠実に扱える乱流シミュレーションはその設計を支える道具となる。
実験診断の側では、乱流を空間的に分解して測る手法が中心である。ビームエミッション分光(BES)による二次元測定から乱流の相関性を求める手法、上部混交共鳴(UHR)を使った後方散乱で小型乱流を調べる手法、集団光散乱による境界乱流の分析など、散乱・撮像系の診断が乱流の波数スペクトルと伝播方向を直接与える。また、条件付き平均やウェーブレット解析、ラグランジュ的な条件付き統計といった統計手法により、ブロブのような個別の構造を抽出し、その速度・形状・輸送への寄与を定量化する解析が定着している。これらの測定は、乱流がランダムな雑音ではなく、条件付きで捉えられる組織的な構造(条件付きエディ)の集まりであるという描像を実験的に裏付けてきた。
乱流が測定そのものに与える影響も無視できない。集団トムソン散乱 のような診断では、乱流揺らぎが測定信号に混入して精度を制限し得るため、乱流の性質を理解することが診断設計の側にも必要となる。逆に、乱流が強い状態ではテラヘルツ帯の強い電磁放射が発生し得るなど、乱流は観測対象であると同時に装置動作への影響源でもある。境界乱流の統計的性質を多装置データベースと突き合わせてスケーリング則として確立する試み、機械学習による低次元モデルで大規模シミュレーションを高速化する試みなど、予測精度を上げるための方法論の整備が現在も続いている。
PHYSICAL REVIEW LETTERS 1977
プラズマ乱流（Plasma turbulence）の研究史 — FusionPapers
