静電乱流

静電乱流とは、磁場閉じ込めプラズマ中に生じる、密度・電位・電場の時間的・空間的な不規則変動のうち、磁場の擾乱を伴わず静電ポテンシャルの変動が支配的なものを指す。この乱流は、プラズマの圧力勾配や密度勾配が持つ自由エネルギーを駆動源として成長する微視的不安定の非線形に飽和した状態であり、磁気閉じ込め核融合装置では古典的・新古典的な輸送予測を大きく上回る「異常輸送」の主因となる。乱流渦に乗った粒子と熱のラジアル方向の動き、すなわち E×B ドリフトによる対流が、閉じ込め性能を直接左右するのである。

図は、内部輸送障壁を持つトカマクプラズマを対象としたジャイロキネティック・シミュレーションで得られた乱流構造を示す。(a)(b) は電位擾乱 δφ の二次元モード構造で、(a) が線形段階、(b) が非線形段階のものであり、(c)(d) はそれぞれに対応する時間平均された E×B せん断率 ω_{E×B} の半径方向分布である。線形段階では特定のモード構造が見えるのに対し、非線形段階では細かい渦が全面に広がり、それに伴って E×B せん断も半径方向に強く変動する様子が読み取れる。この「線形に成長した波動が非線形に乱流へ移行し、自らが生む流れのせん断と相互作用する」という描像が、静電乱流研究の中心にある。

静電乱流とは何か

磁場中の荷電粒子は磁力線の周りをラーマ回転するため、低周波の擾乱に対しては磁場と垂直な面内の運動が支配的になる。擾乱の周波数がイオンのラーマ周波数より十分低く、磁場擾乱が無視できる場合、粒子はほぼ静電ポテンシャル Φ の等高線に沿って E×B ドリフトで動く。このような条件下での乱流が静電乱流であり、電位の等高線が示す「渦」がそのまま流体渦のように振る舞う点が特徴である。擾乱の相対振幅はプラズマ中心部では数パーセント程度と小さいが、端部プラズマでは擾乱密度が背景密度に匹敵するほど大きくなる場合があり、この「強い乱流」領域では準線形理論では記述できない非線形効果が本質的になる。

静電乱流の研究は、トカマクやステラレータといった核融合装置だけでなく、単純磁化トーラスのような基礎実験装置、さらにオーロラ電離圏や宇宙空間のプラズマにも共通して見られる現象である。磁場の幾何学が異なっても「勾配を持つプラズマがドリフト波型の乱流を自発的に励起する」という構図は普遍的であり、装置の壁面近傍や磁力線が開いた領域では、乱流が「ブロブ」と呼ばれる高密度の塊として対流的に輸送される様子も観測されている。

駆動不安定：ドリフト波と交換モード

静電乱流の自由エネルギー源は、主に二種類の不安定に分けられる。

第一はドリフト波である。密度や温度の勾配があると、プラズマは勾配に垂直な方向（ダイヤモンド方向）に一様な反磁性ドリフトを持つ。この背景の上に波動が立つと、電子の応答が完全に断熱的であれば波動は安定化するが、磁力線の曲がりや粒子のトラップ、沿磁場方向の電子の有限な移動時間などが断熱応答を崩すと、勾配の自由エネルギーを吸収して波動が成長する。イオン温度勾配（ITG）モードやトラップ電子モード（TEM）はその代表であり、トカマク中心部の乱流輸送を担う主要な不安定として知られる。トラップ電子は、正の磁気せん断を持つ標準的な配位では ITG モードの成長率を高める方向に働く。これは、トラップ電子が断熱的な電子遮蔽効果を弱め、イオンの圧力摂動を打ち消す力が減るためである。一方、負の磁気せん断を持つ内部輸送障壁領域では、シャフラノフ・シフトや波動・粒子共鳴が絡み、挙動はより複雑になる。

第二は交換（インターチェンジ）モードである。磁力線の曲率が悪い（曲率中心が圧力勾配と同じ側にある）領域では、重力場中の重い流体が軽い流体の下にあるレイリー・テイラー不安定と本質的に同じ機構で、プラズマ塊が外へ膨らむように動く。トカマクの端部や、単純磁化トーラスの実験では、この曲率駆動のモードが乱流の主要な駆動源であることが、ガスパフ撮像による発光変動の観測とラングミュア・プローブによる電位測定の同時解析から示されている。両者の構造の間には 90 度の位相差が見出され、これが曲率駆動モードの特徴的なシグネチャとなる。交換乱流が発達すると、密度の高い「ブロブ」が磁力線に沿った電流の閉回路を介して電荷分離を緩和しながら、その大きさと速度の間に普遍的なスケーリング則に従って対流的に外へ移動する。このブロブ輸送は、スクレイプオフ層での粒子・熱負荷を決める重要な過程である。

乱流輸送と異常閉じ込め損失

静電乱流が輸送をもたらす機構は、E×B ドリフトによる半径方向の対流である。擾乱電場 Ẽ と擾乱密度 ñ の位相関係がずれていれば、⟨ṽ_r ñ⟩ という形の粒子フラックス、⟨ṽ_r T̃⟩ という形の熱フラックスが生じる。古典拡散が衝突によるランダムウォークに基づくのに対し、乱流輸送は渦の相関時間と相関長で決まる拡散であり、その大きさは実験値を説明するのに必要なレベルに達する。このため、磁気閉じ込め装置のエネルギー閉じ込め時間は、衝突に基づく新古典理論の予測を大きく下回り、これが「異常輸送」と呼ばれてきた。

乱流輸送のスケーリングは、装置サイズが大きくなるにつれてボーム型からジャイロボーム型へ移行するとされる。ジャイロボーム・スケーリングでは輸送係数が温度勾配に比例するため、より高温のプラズマほど乱流の影響が相対的に増すことを意味し、反応炉設計上の重要な帰結である。乱流による輸送は粒子・熱だけでなく運動量にも及ぶ。乱流が運ぶ半径方向・極方向の運動量フラックス（レイノルズ応力）は、外部駆動なしにプラズマを回転させる「固有回転」の起源となりうる。さらに、擾乱振幅が大きい端部領域では、モード間結合や乱流の空間的広がり（スプレッディング）を表す非線形の運動量フラックスが、準線形のレイノルズ応力と同程度になりうることが理論的に示されており、強い乱流領域では準線形描像の限界が実際に現れる。

乱流輸送の定量的な記述には、輸送行列の概念が用いられる。粒子・熱・運動量の各フラックスを、密度勾配・温度勾配・回転勾配といった熱力学的な推進力に分解し、それぞれの係数（拡散係数、熱伝導率、ピンチ係数など）を乱流から求める枠組みである。これにより、例えば温度勾配が急になるほど熱フラックスが増えるという実験的事実を、乱流の駆動強度の変化として統一的に理解できる。内部輸送障壁（ITB）は、この輸送係数が局所的に新古典レベルまで低下した領域であり、そこでは温度・密度勾配が急峻になる。ITB の形成・維持・崩壊の理解は、静電乱流研究の中心的な応用課題である。

ゾーナルフローと乱流の自己制御

静電乱流の最も重要な非線形性の一つが、ゾーナルフローの生成である。ゾーナルフローとは、極方向にもトロイダル方向にも一様で、半径方向にのみ変化する E×B 流れのことである。ドリフト波乱流の渦は本来あらゆる波数を持つが、非線形相互作用の中で波数が互いに打ち消し合う成分、すなわち変調不安定を経て、波長の長い帯状の構造が選択的に増幅される。これは、乱流のレイノルズ応力が半径方向に不均一であることが駆動力となって起こる、乱流自身による流れの自発的生成と見なせる。

生成されたゾーナルフローは、その半径方向のせん断によって乱流渦を引き裂き、位相相関を断ち切る。E×B せん断率が線形成長率の最大値を上回れば、乱流は実効的に抑制される。この「乱流がゾーナルフローを生み、ゾーナルフローが乱流を抑える」という負帰還機構により、乱流は系のパラメータで決まる準定常レベルに自己制御される。輸送のスケーリングがボーム型からジャイロボーム型へ変わるのも、この自己制御の働きによる部分が大きい。L モードから H モードへの遷移（L–H 遷移）の際に観測されるリミットサイクル振動は、乱流とゾーナル電場の間の捕食者・被食者関係の実験的現れと解釈されており、複数のトカマクで確認されている。

図は、同じ時刻の二次元乱流構造を、ゾーナルフローを含む場合 (a) と含まない場合 (b) で比較したものである。ゾーナルフローを除いた (b) では大規模な渦構造がはっきりと残るのに対し、含む (a) では渦がせん断によって細分化され、構造が大きく崩れている。乱流の帯状構造の半径方向への断片化と極方向への長波長構造の生成という、ゾーナルフロー生成の描像を可視化した例といえる。

なお、ゾーナルフローを駆動するのはレイノルズ応力だけではない。近年の理論・シミュレーション研究は、乱流によるエネルギーフラックスやトロイダル・レイノルズ応力もゾーナル電場の駆動に寄与しうることを示しており、イオン圧力勾配の振動が実験で観測される L–H 遷移のリミットサイクルと結びつくという観測事実を説明する枠組みとして注目されている。また、乱流スペクトルの非対称性がなくても非線形運動量フラックスは生じうるため、極方向回転の説明にはレイノルズ応力だけでは不十分な場合がある。

各装置での実験と数値シミュレーション

静電乱流の研究は、装置の規模も目的も異なる多様な実験・理論・計算の組み合わせで進められている。

実験面では、トカマク端部や反転磁場ピンチ（RFP）の端部で、ガスパフ撮像による発光変動とラングミュア・プローブによる電位・密度変動の同時測定が標準的な手法となっている。条件付き平均解析により乱流構造の形と位相関係を抽出し、駆動機構を同定する。端部乱流は多重フラクタル性を持つことが統計解析から示されており、単純な拡散過程ではない非自己相似的な時空間構造を持つことが分かっている。確率的磁場中での乱流輸送、円筒プラズマでの自己組織化、単純磁化トーラスでのブロブ生成とその電極バイアスによる制御、高速イオンの乱流による散乱など、基礎装置での研究は核融合装置では分離しにくい素過程を切り出す役割を果たしている。宇宙プラズマでは、オーロラ領域の密度・電場の不均一が静電乱流の源となり、広帯域の静電擾乱が荷電粒子を散乱する過程が研究されている。

数値シミュレーションでは、第一原理に基づくジャイロキネティック・コードが主流である。これは、ラーマ運動を摂動として扱い、磁場と垂直な面内の速い回転運動を平均化した上で、分布関数の緩慢な発展を解く手法である。実験で得られた平衡配位（磁気せん断、シャフラノフ・シフト、トラップ粒子の割合など）をそのまま用いて、線形の分散関係・モード構造から非線形の飽和・輸送までを追うことができる。トカマクの ITB プラズマやステラレータの準等磁場配位など、実際の装置パラメータでのシミュレーションが行われ、トラップ電子の非断熱効果や磁気せん断の符号、全磁束面上の効果が乱流の振る舞いにどう影響するかが調べられている。一方で、数値コードの検証・妥当化（ベリフィケーション・バリデーション）も活発な課題であり、数値リートランスフォームのような手法の精度確認が進められている。実験・理論・シミュレーションの三者が同じ物理量を比較できるようになって初めて、乱流輸送の定量的予測が可能になるというのが、この分野の現在の到達点である。