ドリフト波乱流

ドリフト波乱流とは、磁場閉じ込めプラズマにおいて密度や温度の空間的不均一性を駆動力として励起されるドリフト波が、多数重ね合わされて非線形に相互作用する状態を指す。ドリフト波はプラズマの逆磁気方向に伝わる低周波の波であり、これが乱流化すると粒子や熱を磁場を横切って運ぶ異常輸送を引き起こす。核融合プラズマの閉じ込め性能を律速する主要因の一つと考えられており、トカマクやスタセレータの中心閉じ込め、ペデスタル、密度限界近傍の挙動まで、その理解は磁場閉じ込め研究の中心課題に位置づけられている。

ドリフト波とその自由エネルギー

ドリフト波の源は、密度勾配という自由エネルギーである。磁場に垂直な方向に密度がなだらかに減少するプラズマでは、電子とイオンが磁力線に沿って動きやすさを異にするため、粒子のラーモア運動や平行電子応答と組み合わさって、プラズマ全体が逆磁気方向に漂うような動きが生じる。これがドリフト波であり、その周波数はイオンのラーモア周波数よりずっと低く、波数スケールはイオンのラーモア半径あるいは音速で規格化したイオン音波半径 ρs 程度のオーダーにある。勾配が十分に急で、電子の平行応答が抵抗などによって完全でない場合、波はプラズマからエネルギーを汲み上げて成長する。これがドリフト波不安定であり、抵抗性ドリフト波不安定や、イオン温度勾配(ITG)不安定などがその代表例である。

乱流としてのドリフト波を理解する鍵は、この系が本質的に二次元的な非線形ダイナミクスを持つことにある。強い磁場に沿った運動は速く緩和されるため、磁場に垂直な断面内の電場 E×B ドリフトが主要な非線形過程を担う。E×B ドリフトは電場に垂直な方向へプラズマを運ぶため、ポテンシャルの空間構造がそのまま流体渦のように変形・輸送され、渦度方程式に従う二次元乱流に近い振る舞いを示す。ただしドリフト波系では、密度勾配に由来する波の伝播と、平行電子応答の緩和(断熱性)が非線形渦運動に重なるため、純粋な二次元乱流とは異なる固有の構造形成が起こる。

実際の乱流状態を時間発展で見ると、初期にほぼ一様な揺らぎから、渦の合体と引き伸ばしを経て、系のスケールに匹敵する大きな構造が形成されていく様子が数値計算で追える。次の図は、ポテンシャル・密度・渦度のスナップショットを時間順に並べたものであり、乱流が小スケールの揺らぎから渦構造へと自己組織化していく過程を示している。

図に見えるように、ポテンシャル場では時間の経過とともに細かい擾乱が合流して孤立した渦塊が現れ、密度場には双極子的な構造が、渦度場には強い渦の集中が形成される。このようなコヒーレント構造の生成は、二次元乱流に特有の逆カスケード(エネルギーが大きなスケールへ移る非線形輸送)と、ドリフト波の勾配駆動との競合によって実現するものであり、ドリフト波乱流が単なる雑然とした揺らぎではなく構造を持つ系であることを示している。

非線形乱流:渦、ストリーマー、クランプ

ドリフト波乱流の非線形構造は、その伸びる方向によって分類される。ラーモア半径から数えると、等方に近い渦(エディ)が基本だが、非線形相互作用が進むと、半径方向に細長く伸びた構造(ストリーマー)や、ポロイダル方向に広がる帯状構造(後述のゾーナルフロー)へとエネルギーが配分される。ストリーマーは半径方向輸送を直接担う構造であり、その振幅には波動運動エネルギーと勾配の関係から上限が議論される。また、乱流中で粒子を長時間閉じ込める位相空間の「クランプ」(塊)や、双極子渦の生成も報告されており、乱流が輸送を担う粒子を実空間でも位相空間でも局所的に捕獲しうることが示されている。

これらの構造形成を理解する枠組みとして、乱流を波数空間の「波の集団」とみなす波動運動論がある。ドリフト波のスペクトルを実空間と波数空間からなる位相空間上の流体とみなすと、非線形相互作用は波数空間での拡散として表され、乱流強度の局在や伝播を記述できる。この見方では、乱流は空間的に一様に分布せず、局所的に高強度の領域を形成して伝わる(乱流の広がり)ことが自然に導かれる。実際、測地音響モード(GAM)の励起を介して乱流が空間的に広がる現象や、乱流強度の間欠性を記述する非摂動的な統計理論が研究されており、乱流の統計性質そのものが研究対象となっている。低温プラズマ実験での統計的性質の測定は、こうした理論的描像を実験的に検証する場を提供している。

長谷川・三間(Hasegawa–Mima)と長谷川・脇谷(Hasegawa–Wakatani)モデル

ドリフト波乱流の理論研究を支えてきたのが、一対の簡約モデルである。長谷川・三間(Hasegawa–Mima)モデルは、電子が磁力線に沿って瞬時に(断熱的に)ボルツマン応答すると仮定し、イオンの分極を考慮した渦度方程式のみで系を閉じる。このモデルでは密度勾配は固定されており、ドリフト波は実質的に減衰する系として扱われるが、乱流から大規模構造へのエネルギー移動や、ゾーナルフロー生成の基本機構を最も明快に示す「パラダイム」として広く用いられてきた。ドリフト波の分散関係 ω = ω*/(1 + k⊥²ρs²)(ここで ω* は逆磁気ドリフト周波数)が基本スケールを与える。

一方、長谷川・脇谷(Hasegawa–Wakatani)モデルは、電子の平行応答を有限の抵抗を介して記述し、密度とポテンシャルの連成を陽に含む。これにより抵抗性ドリフト波不安定が自励的に駆動され、乱流が勾配からエネルギーを汲み上げながら飽和する自己無撞着な系が得られる。平行応答の強さを表す断熱性パラメータの値によって、系は断熱極限(電子がほぼ瞬時に応答し、乱流は減衰的なドリフト波の集団として振る舞う)から流体力学的極限(電子応答が遅く、不安定波が強く駆動される)まで連続的に変化する。この断熱性の違いは、乱流の性質だけでなく、後述するゾーナルフローとの空間的な関わり方までも変えることが、修正長谷川・脇谷(Hasegawa–Wakatani)モデルを用いた非線形シミュレーションと固有モード解析によって示されている。すなわち、断熱性が高い場合、乱流はゾーナルフローの速度が最大となる領域(「丘」)に局在するのに対し、断熱性が低い流体力学的な場合には、ゾーナルフローの谷に乱流が捕獲される(「谷」への捕獲)という分岐が起こる。どちらに局在するかは、ドリフト波のポロイダル位相速度とゾーナルフロー速度の大小という単純な規則で整理でき、密度限界近傍で衝突周波数が上がり断熱性が低下すると乱流・輸送の空間分布が質的に変化しうることを示唆する。

トカマクのトロイダル幾何を取り入れた三次元計算では、これらの簡約モデルで見られた渦運動と輸送が、実際の磁場配置の中でも定量的に追えることが示されており、簡約モデルと実機幾何をつなぐ橋渡しが進んでいる。また、多種イオンを含む系への拡張(例えば重水素・三重水素混合や不純物を含むプラズマ)でも、抵抗性ドリフト波乱流による異常輸送が調べられている。

ゾーナルフローによる自己制御と輸送抑制

ドリフト波乱流の最も重要な性質の一つは、乱流自身が自分を抑制する流れを生むことにある。乱流の E×B 揺らぎはレイノルズ応力を通じて、ポロイダル方向に一様で半径方向に変化する帯状の流れ、すなわちゾーナルフローを駆動する。ゾーナルフローは半径方向のシアーを持つため、乱流の渦を半径方向に引き伸ばして断片化し、勾配を横切る輸送を効率的に弱める。この「乱流がゾーナルフローを生み、ゾーナルフローが乱流を削る」という負帰還は、捕食者・被食者(prey–predator)型の力学としてモデル化でき、乱流強度とゾーナルフロー強度が時間的に振動しながら飽和する様子を定量的に記述する。

次の図は、この自己制御の効果を可視化したものである。平行電子応答の弱い(流体力学的な)場合には乱流の渦構造が空間に広く分布するのに対し、応答の強い(断熱的な)場合には乱流が半径方向の狭い帯に強く局在し、そこに鋭いゾーナルフローのシアーが形成されている。

図の上段はポテンシャル揺らぎのスナップショット、下段はゾーナルフローの半径方向プロファイルであり、乱流の局在とゾーナルフローのシアーが対応して現れることが分かる。

ゾーナルフローによる制御は、磁場閉じ込めの高性能化モード(LモードからHモードへの遷移)の中心的なパラダイムとされている。Hモードでは強いシアーを持つ大規模流が乱流輸送を抑制し、エッジに勾配の急な障壁が形成される。遷移には外部入力電力の閾値があり、その閾値は密度やプラズマ組成に依存するが、シアーの発達と乱流抑制の間の直接的な因果関係を実験で確定することは容易ではなく、電場の井戸の深さと乱流レベルの相関などが測定されている。近年は、時間依存の確率密度関数と情報理論を用いて、乱流・ゾーナルフロー系が生成する統計状態の多様さを定量化する解析も行われている。

ゾーナルフローとの相互作用には、トロイダル性に由来する側面もある。トロイダル系ではゾーナルフローは測地曲率を通じて有限周波の測地音響モード(GAM)と結合する。ドリフト波乱流が駆動するGAMの半径方向伝播方向は、乱流スペクトルのピーク波数に依存して定在波、外向き伝播、内向き伝播と変化することが波動運動論に基づく解析で示されており、乱流と固有モードの非線形結合がモードの空間構造そのものを決定しうる。また、ゾーナルフローによる乱流の捕獲・局在化は、前節で述べた断熱性依存の分岐を通じて、密度限界近傍の輸送変調とも結びついている。

さらに、高速イオン(追加加熱で生じるエネルギー粒子)がこの制御系に与える影響も注目されている。高速イオンが熱イオンを「希釈」すると、ドリフト波の分散関係が変化してゾーナルフロー生成を妨げる周波数のミスマッチが減り、ゾーナルフロー活動がむしろ強まる。その結果、飽和乱流強度は希釈率と密度勾配の関数として強く抑制され、輸送が新古典輸送に近い水準まで下がりうることが、波動運動論に基づく自己無撞着な計算で示されている。これは、高速イオンが主導する内部輸送障壁(ITB)を持つ高性能閉じ込めモードの理解につながる描像である。

異常輸送と閉じ込めへの意義

ドリフト波乱流がなぜ核融合研究の中心に置かれるかといえば、その答えは異常輸送にある。古典的・新古典的な衝突に基づく輸送予測を大きく上回る粒子・熱の漏出が磁場閉じ込めプラズマで観測され、その主因が乱流による E×B 輸送であることが、数値シミュレーションと実験の両面から確立されてきた。ドリフト波乱流は密度・温度勾配を下流へ平坦化する方向に輸送を生むため、閉じ込め性能は「駆動される勾配」と「乱流が許す勾配」の均衡で決まる。混合長的な見積もりでは、輸送係数は乱流の相関長(イオン音波半径程度)と揺らぎ振幅の積で与えられ、乱流強度の抑制がそのまま閉じ込めの改善に直結する。

この観点から、ドリフト波乱流の研究は三つの方向で閉じ込めに接続している。第一に、ゾーナルフローによる乱流抑制は L–H 遷移や内部輸送障壁の物理の中核であり、輸送の分岐(ビファケーション)として理解される。第二に、高速イオンによる希釈効果のように、乱流の駆動と制御のバランスを変える機構の解明は、新たな高性能運転モードの設計指針となる。第三に、密度限界近傍では衝突性の増加が断熱性を下げ、乱流・ゾーナルフロー系の空間構造と輸送の変調をもたらすため、密度限界の理解に乱流輸送の視点が不可欠になっている。

また、ドリフト波乱流は他の揺動とも結合する。乱流はアルフベン固有モードを散乱・減衰させることが示されており、高速イオン物理と熱スケール乱流の跨ぐスケール結合の経路の一つとなっている。さらに、乱流中の重い不純物粒子の運動をラグランジュ的に解析する研究では、乱流の多スケール構造が不純物輸送に直接影響することが調べられている。このようにドリフト波乱流は、勾配駆動という単純な自由エネルギー源から出発しながら、渦構造の形成、ゾーナルフローによる自己制御、そして閉じ込め性能の決定までを貫く、磁場閉じ込めプラズマ物理の中心的な現象なのである。