周辺乱流

磁気閉じ込め核融合装置において、周辺乱流とは、最後の閉じた磁気面（LCFS）近傍およびその外側のスクレイプオフ層（SOL）に生じる強いプラズマ密度・電位・温度ゆらぎのことである。閉じ込め領域の内部で支配的な微視的乱流と比べ、周辺乱流はゆらぎの相対振幅が桁違いに大きく、密度ゆらぎの相対振幅は典型的に δn/n ∼ 5–100 % に達する。周波数帯域はおよそ 10 kHz から 1 MHz、磁力線に垂直な構造の大きさは 0.1–10 cm 程度である一方、磁力線に沿っては数メートルにわたって伸びる。つまり周辺乱流の構造は、断面では小さくて長く伸びた「フィラメント（細糸）」であり、実質的に二次元的な性質を持つ。この乱流は粒子と熱の断面外輸送を直接担い、L モードから H モードへの遷移、SOL の熱流幅、そして密度限界といった、周辺物理の中心的な問題すべてに深く関わっている。

周辺乱流の性質と特徴

周辺乱流が何を物理的に見ているのかを示すため、まずダイバータ配位のトカマク周辺の全体像から見ていく。

図は、ポロイダル断面内の密度スナップショット（左）と重水素放射率分布（右）である。左図では、セパラトリックス（白破線）の内側に高い密度の閉じ込め領域が広がり、その外側の SOL に向かって密度が下がる様子、また X 点からダイバータに向かって磁力線に沿った構造が見て取れる。右図は中性ガスのイオン化に伴う放射が X 点近傍やダイバータ面上で強くなることを示し、周辺プラズマが中性粒子との相互作用と強く結合した領域であることを物語っている。

実験的な特徴付けは主にラングミュアプローブ、遠赤外レーザー散乱、反射計、そして高速可視カメラによって行われてきた。これらの測定が明らかにした普遍的な性質のひとつは、SOL に入るにつれて密度ゆらぎの歪度（skewness）が増大することである。これは密度の正のスパイク、すなわち「ブロブ」が SOL に存在することの足跡であり、複数の装置のデータを統合した国際周辺乱流データベースの解析によって、すべての装置で再現された。ゆらぎの相関長・相関時間は巨視的なドリフト波パラメータに対してべき乗則に従うが、その依存性は純粋なドリフト波乱流から期待されるものより弱く、結果として乱流拡散係数は磁場強度の逆数でスケールする。これはギャロ・ボーム型よりもボーム型に近いスケーリングであり、周辺乱流の輸送を予測する上で重要な指針となっている。また、閉じ込め領域端と SOL でほぼ同一のスケーリング則が得られたことは、閉じ込め領域内のドリフト波乱流と SOL のブロブとが強く結合していることを示している。

測定されたゆらぎは電位・電子温度にも及ぶ。球状トカマクでの高速撮像実験では、イオン温度ゆらぎが電子温度ゆらぎの約 2 倍になることが報告されており、これは平行方向の高速な電子熱伝導によって電子の温度ゆらぎが均される一方、イオン側のゆらぎが残るためと説明されている。このように周辺乱流は静電的な密度・電位ゆらぎにとどまらず、温度ゆらぎや場合によっては磁気ゆらぎを含む多次元的な現象である。

駆動機構：ドリフト波、抵抗性バルーニング、せん断流

周辺乱流の駆動源として最も基本的なのはドリフト波である。密度勾配と磁場の曲率が組み合わさることで、プラズマは磁力線に垂直な方向に伝播する波動を自発的に励起する。スラブ幾何での理論解析とトロイダル幾何での理論解析の双方から、ドリフト波の伝播が周辺乱流の源となることが示されてきた。ドリフト波乱流が「断熱的」であるかどうかは、断熱性パラメータ α = k∥v²ₜₕₑτₑᵢ/ω（ここで k∥ は磁力線平行方向の波数、τₑᵢ は電子・イオン衝突時間、ω はドリフト波の主要周波数）によって決まる。α が 1 を超えると電子はボルツマン応答に近づき乱流は弱まるが、α が 1 を下回ると流体極限へ移行し、非断熱的な抵抗性モードが支配的となる。

密度が高く衝突性が増した環境では、抵抗性バルーニングモード（RBM）が主要な駆動源となる。これは圧力勾配と磁場曲率によるバルーニング駆動に、周辺プラズマの有限抵抗が絡む不安定であり、ドリフト波乱流とは異なる非線形発展をたどる。理論モデル（RDZ モデルやその拡張であるセポス（SepOS）モデル）は、密度の増加に伴い周辺乱流が静電的なドリフト波乱流から電磁的な RBM 乱流へと遷移することを予言し、これは複数のトカマク実験および 3 次元シミュレーションによって支持されている。電磁効果のうち磁気フラッター（摂動磁場に沿ったプラズマ流）は、輸送に二重の影響を与える。高ベータ環境では磁気面を壊してラジアル方向の電子熱輸送を直接増大させる一方、低ベータ環境ではむしろ静電的な E×B 輸送を間接的に抑制する効果が知られている。フルサイズの周辺乱流シミュレーションでは、フラッターを含めることで乱流輸送が大幅に減少することが示されており、予測的シミュレーションにおいて電磁項を落とすことの危険性が指摘されている。

図は、磁気フラッターを含まないモデル（左）と含むモデル（右）で、瞬間的な密度分布（上段）と電位ゆらぎ（下段）を比較したものである。フラッターを含む場合、閉じ込め領域端の乱流が弱まり、SOL に伝播するフィラメントの数が減少し、電位ゆらぎの強度が下がって渦のサイズが大きくなる様子が見て取れる。

乱流の飽和には E×B せん断流が本質的な役割を果たす。周辺では、平均流（径向電場に由来する）とゾーナルフロー（乱流自身がレイノルズ応力を通じて駆動する帯状流）の両方が存在し、これらのせん断が渦を引き裂くことで乱流の成長が抑えられる。シミュレーションでは、フィラメントの源が自発的に形成された E×B せん断層の近く（閉じ込め領域端から約 2 cm 以内）にあることが示されており、乱流の発生位置とせん断層の位置が密接に結びついていることがわかる。

フィラメント、ブロブ、間欠性

周辺乱流の最も顕著な現れ方が、フィラメント（ブロブ）と呼ばれる局在構造である。高速撮像実験では、局所的に強い光放射を示す「ブロブ」がポロイダル方向とラジアル方向の両方に移動する様子が捉えられており、その典型的な速度は約 10⁵ cm/s である。ブロブは SOL における間欠的な大きなバーストとして観測され、条件付き平均解析によれば、すべての装置で外向きの粒子輸送に寄与する一方、負の密度スパイクである「ホール」には一様な伝播パターンが見られない。

フィラメントのモード数やサイズはプラズマパラメータに依存する。L モードでの自動画像解析によれば、モード数は密度の増加とともに増え、エッジ安全因子 q₉₅ の減少とともに減る。フィラメント幅は逆の依存性を示す。3 次元二流体シミュレーションは観測されたフィラメントと同程度のサイズとラジアル速度を再現しており、数値モデルが周辺乱流の実体を捉えつつあることを示している。

統計的には、周辺乱流は強い間欠性を示す。ウェーブレット解析や高次統計量の解析から、ゆらぎの分布はガウス分布から大きく外れ、稀で大きなイベントが輸送の大部分を担うことが明らかになっている。この間欠性は、SOL での粒子・熱輸送を平均的な拡散係数だけでは記述できないことを意味し、周辺輸送モデルの構築上の重要な制約となっている。

乱流輸送：L-H 遷移、拡がり、熱流幅

周辺乱流輸送が最も劇的に変化するのが L モードから H モードへの遷移である。L-H 遷移では E×B せん断流の急増によって乱流が抑制され、実験でも周辺乱流レベルの顕著な減少が観測される。H モードでは SOL の E×B せん断が乱流をほぼ消し去るほど強く、その結果として新古典的な磁気ドリフトに由来する極めて狭い熱流幅が現れる。ヒューリスティック・ドリフトモデルは熱流幅を λq ∼ V_Dτ∥ ∼ ερ_θ,i（ε は逆アスペクト比、ρ_θ,i はポロイダルイオン・ジャイロ半径）と予測し、これは H モードでの広範な測定と一致するが、あまりに狭いため工学的に受け入れがたい。

一方、オーミック加熱や L モードのプラズマでは、熱流幅は SOL 内の局所的なラジアル輸送と磁力線に沿った熱流の積で決まると伝統的に理解されてきた。すなわち λq ∼ (χ⊥τ∥)^(1/2) であり、χ⊥ は乱流熱拡散係数、τ∥ ∼ Rq/Vₜₕ,ᵢ は平行通過時間である。しかし近年、SOL の乱流が必ずしも SOL 内で局所的に生成されるとは限らず、閉じ込め領域端（ペデスタル）で生成された乱流エネルギーがラジアル方向に伝播する「乱流拡がり（turbulence spreading）」の寄与が重要であることが認識されてきた。生産比 Rₐ（LCFS を横切る乱流エネルギーフラックスの、SOL 内での積算生産に対する比）を用いた解析では、Rₐ がしばしば 1 を超え、Rₐ が大きいほど熱流幅が大きくなるという相関が見出された。これは多くの場合、SOL の乱流レベルと熱流幅の決定において、エッジからの拡がりが局所生産と同等以上に重要であることを示している。

図は、熱流幅がどのように決まるかを説明する概念図である。エッジ乱流と乱流輸送、E×B せん断、SOL 内の局所生産と LCFS を横切る乱流拡がりのエネルギー生産比 Rₐ、SOL のゆらぎと乱流混合、新古典ドリフト、垂直輸送と平行通過といった要素がどう連鎖して最終的な熱流幅 λq を形作るかが示されている。

乱流拡がりは密度限界近傍でも重要な役割を果たす。グリーンウォルド密度限界に近づくと、乱流拡がりの指標となる乱流内部エネルギーのラジアルフラックスが増大し、エッジ冷却と同時に観測される。この場合、乱流粒子フラックス自体は顕著に変化しないこともあり、乱流拡がりはエッジ冷却の指標として粒子輸送よりも優れている可能性がある。間欠的な対流イベント（ブロブとホール）が乱流拡がりの運び手として機能していることが示唆されており、エッジと SOL をつなぐ物理として注目されている。

周辺乱流と密度限界

核融合炉の経済性は高密度運転を要求するが、実験的には密度に上限（密度限界）が存在する。トカマクではこの限界は経験的にグリーンウォルド密度限界 nGW = Ip/πa² で記述され、線平均密度がこれに近づくと閉じ込めの劣化が進み、超過するとディスラプションに至る。密度限界の物理は一般にエッジプラズマの冷却と結びつけられており、冷却が MARFE を生み、電流分布の収縮を経て大規模な MHD 不安定を誘発すると考えられている。しかし、その冷却の引き金となるのが周辺乱流の増強である可能性が高い。

密度を上げていくと、まず周辺乱流がドリフト波乱流から抵抗性バルーニング乱流へと遷移する。この遷移はグリーンウォルド係数の比較的低い段階（複数の実験で fGW ∼ 0.2 前後）で観測される。続いてラジアル磁気ゆらぎレベルの強い増大が見られ、さらにプラズマゆらぎの飽和が起こるが、これらのイベントはいずれも実際のディスラプションからはまだ遠い。このように密度限界に至るまでの過程で複数の乱流レジーム変化が段階的に起こることが、近年の実験とシミュレーションで明らかになりつつある。

もうひとつの有力な機構がエッジせん断流の崩壊である。密度の増加とともにゾーナルフローのエネルギーが減衰し、E×B せん断による乱流抑制が弱まることで乱流レベルと外向き輸送が急増する。この機構は、長距離相関（ゾーナルせん断流の指標）の減少、レイノルズパワーの低下、乱流粒子輸送の増大として複数の装置で観測されている。理論的には、衝突性の増大が断熱性を下げ、ゾーナルフローの生成を減らすことでせん断流崩壊が起こると説明される。

図は、高密度プラズマで考えられるフィードバック機構をまとめた概念図である。電子密度の上昇が衝突性の増大とエッジプロファイルの平坦化の両方を引き起こし、それぞれが乱流遷移・増強と E×B せん断流（平均流）の低下につながり、両者が合わさって外向き輸送をさらに強めるという正のフィードバックループが描かれている。

ただし、乱流遷移やせん断流崩壊は密度限界の十分条件ではない可能性がある。せん断層の崩壊や静電的から電磁的 RBM への遷移が、実際の密度限界よりかなり低い密度で観測された例もあり、密度限界に至るまでの過程には複数の段階があると考えられている。また、中性ガスの役割も無視できない。リサイクリング（壁での中性化と再イオン化）は密度・温度プロファイルの形成に本質的であり、シミュレーションに中性ガスモデルを組み込むと、イオン化源の非一様な分布がポロイダル方向の密度非対称性を生み、乱流の駆動モードがバルーニング型からイオン温度勾配型へと変わることが示された。中性粒子による希釈と放射冷却はダイバータに向かう電子流体を冷やし、衝突性の増大を通じて平行熱伝導を弱め、その結果として垂直乱流輸送がラジアル温度プロファイルを平坦化する。ダイバータ中性密度がさらに上がると遠方 SOL の密度が上昇し、「密度ショルダー」が形成される。このように周辺乱流は、中性粒子、放射、せん断流、電磁効果が絡み合う複合系の一部として密度限界物理を支配しており、その完全な理解と予測が次世代融合炉の排熱・排粒子設計の鍵となっている。