熱輸送

磁気閉じ込め核融合プラズマにおける熱輸送とは、加熱によって投入されたエネルギーがプラズマ中心から外周へ、そして壁へと逃げていく過程のことである。核融合炉の成立条件は「十分な高温を十分な長時間保つ」ことだが、その成否を実際に決めるのがこの熱輸送の大きさであり、閉じ込め磁場の強度や装置サイズの設計はすべて熱輸送の予測に依存している。熱輸送を論じる基本量は熱流束であり、温度勾配に比例する拡散型の関係式で記述されるのが普通だが、実験が示すところによれば、その比例係数は磁場配置・衝突頻度・温度勾配そのものに強く依存し、単純な拡散では説明できない現象も数多く観測されている。

磁気閉じ込めにおける熱輸送：新古典輸送と異常輸送

磁場中の荷電粒子は磁力線に沿って旋回しながら運動するため、垂直方向の輸送は原理的に極めて小さい。しかしトロイダル磁場では強度が不均一なため、磁力線の弱い側に捕らえられた「バナナ軌道」の粒子が上下にドリフトし、これが拡散の下限を与える。この粒子軌道論に基づく輸送が新古典輸送であり、その大きさは有効リップルと呼ばれる無次元量で磁場幾何の影響を特徴づけられる。ヘリカル装置の設計思想の中心はまさにここにあり、例えばW7-Xでは捕獲電子のバウンス平均半径方向ドリフトを抑制することで新古典輸送を低減し、体積平均ベータ約5%まで良好なMHD安定性を保つという最適化が図られている。さらに最大J特性（第二断熱不変量の極大を磁軸上に置く配置）は衝突無しの捕獲電子モードに対する耐性をもたらすと期待されている。

ところが、新古典理論が与える輸送は実測値の下限にすぎない。実際のプラズマでは、微視的乱流が引き起こす異常輸送が支配的となり、新古典レベルを一桁以上上回ることが珍しくない。W7-Xの初期高性能運転でも、低ベータ放電では電子エネルギー輸送が新古典レベルを有意に上回ることが測定されている。熱輸送の研究とは、突き詰めれば「どの領域で新古典輸送が支配的か、どこから異常輸送が支配的か、そして異常輸送を何が駆動しているのか」を問う営みである。

この対比を示すのが、W7-Xの「標準」配置と「高ミラー」配置で測定された電子・イオン温度分布と、パワーバランス解析から得られた電子熱流束である。電子熱流束は新古典予測値を大きく上回り、この差が乱流に起因する異常輸送の大きさを直接物語っている。

図の(a)には両配置の電子温度（丸印）とイオン温度（×印）の半径方向分布、およびECRH（電子サイクロトロン加熱）の正規化パワー堆積分布が示され、(b)にはパワーバランスから求めた電子熱流束と新古典熱流束が並べて描かれている。電子加熱が局所的に投入された領域で熱流束が急増し、それが新古典値をはるかに超えている様子が読み取れる。

乱流駆動機構：ITG、TEM、ETG、マイクロテアリングモード

異常輸送の主犯は、温度勾配や密度勾配を自由エネルギー源とする微視的不安定性、すなわちドリフト波乱流である。主要なモードは三つに大別される。イオン温度勾配（ITG）モードはイオン温度勾配が規格化された臨界値を超えると励起され、主にイオン熱輸送を駆動する。捕獲電子モード（TEM）は磁場不均一で捕獲された電子のバウンス運動と共鳴して成長し、電子熱輸送に寄与する。電子温度勾配（ETG）モードは電子スケール（イオン旋回半径より短い波長）で励起されるもので、波長は短いが電子熱輸送に直接効く。これらに加えて、有理面（磁力線が整数回の周回で元に戻る面）に局在する電磁的不安定性としてマイクロテアリングモードがあり、共鳴的な平行電流層を形成して磁気島を作り、入れ子状の磁面トポロジーを破る。高ベータ球状トカマクや通常トカマクのペデスタル領域で重要とされ、すべての分岐において電子温度勾配が不安定化の必要条件である。

どのモードが支配的になるかはプラズマ条件で切り替わる。例えばアルカトールC-ModのオーミックLモードでは、低衝突頻度のLOC領域では外側領域の乱流がTEM支配、高衝突頻度のSOC領域ではITG支配へと移行することが線形ジャイロキネティック計算と整合して示されている。またASDEX UpgradeのHモード実験では、イオン熱輸送が理論予測どおりITG不安定性によって駆動され、温度比（Te/Ti）や高速イオン集団が乱流安定性に与える影響がジャイロキネティック計算と一致することが確認された。ITGは電子熱流束にも寄与するが、電子加熱の全量を賄えない場合には別途電子系の不安定（ここではETG）が必要になることも示された。

マイクロテアリング乱流には特異な飽和機構がある。共鳴面に局在した磁気島が局所の電子温度勾配を平坦化し、それ自身の駆動源を弱めることで飽和するのである。乱流による熱流束は有理面近傍に強く局在するため、低次の有理面上で局所拡散率が十分大きくなると、そこだけ勾配が削り取られ、勾配は有理面から離れた場所に押しやられる。この効果は系のサイズに依存し、非ジャイロボーム型のスケーリングを生むため、マイクロテアリング乱流に支配された大型炉では予想より閉じ込めが劣化しうると議論されている。

臨界勾配、スティフネス、プロファイルの頑健性

乱流駆動熱輸送の最も重要な性質は、勾配に対する閾値の存在である。温度勾配が臨界値以下では乱流熱流束はほぼゼロだが、臨界値を超えると急激に増大する。この非線形な応答のため、加熱パワーを増やしても温度勾配はほとんど上がらず、熱流束だけが増える。これをスティフネス（剛性）と呼び、その帰結として温度プロファイルは加熱条件の変化に対して「頑健（レジリエント）」になる。トカマクで顕著に観測されるこの性質は、古典的ヘリカル装置では見られないとされてきた。

ASDEX Upgradeでは電子サイクロトロン加熱を用いて、エッジを固定したまま閉じ込め領域の電子熱流束だけを変化させる実験が行われ、∇Te/Teに明確な臨界値(∇Te/Te)臨界が存在すること、定常解析と変調解析の双方がそのモデルと一致することが示された。特に軸外加熱では、堆積点内側の領域がちょうど臨界値上に置かれるため、実測のTe分布から直接臨界値を読み取ることができ、その絶対値と半径依存性（捕獲電子割合と密度勾配で決まる）がTEM乱流の予測と一致し、ETG乱流の予測閾値は実験値より高く、この条件下ではETGは励起されないと結論された。

スティフネスの定量的な尺度として、熱パルスから得られる拡散率とパワーバランスから得られる拡散率の比が用いられる。W7-Xの三つの磁場配置（標準・高ミラー・高イオタ）での測定では、低衝突頻度側でこの比が1を大きく超え、衝突頻度の増加とともに1に近づく、すなわち衝突頻度が高いほど輸送が拡散的・局所的になる傾向が見出されている。また、捕獲電子がTEM乱流の駆動源であることから、最大J特性を深く捕獲粒子がサンプリングする配置ではTEM駆動の電子熱流束が小さくなると期待される。

熱パルス伝播と摂動輸送解析

熱輸送を実験的に分離する有力な手法が摂動解析である。プラズマに短時間の熱的な擾乱（パルス）を与え、その温度擾乱の伝播を測定することで、定常のパワーバランス解析とは独立に熱拡散率を求めることができる。擾乱の源としては、トカマクで自然に生じるソートゥース崩壊、ECRHの変調、レーザーアブレーションによる不純物注入（エッジ冷却パルス）などが使われる。

JETの初期の系統的研究では、ソートゥース起源の熱パルスを12チャンネルの電子サイクロトロン放射計で高時間分解能で追跡し、電子熱拡散率が加熱パワー全体（2〜13.5 MW）にわたってほぼ一定値（2.5±0.5 m²/s程度）をとることが示された。一方パワーバランスから得られる拡散率はオーミック加熱では熱パルス値の約2.5分の1であり、補助加熱の増加とともに両者は近づく。この差は、熱拡散率が温度勾配に対して急峻に増大する関数であること（スティフネス）で統一的に説明できる。勾配依存型の輸送モデルは、熱パルスの振る舞い、加熱増大に伴う閉じ込めの劣化、プロファイルの一貫性を同時に記述でき、しかも総入力パワーのような非局所パラメータを必要としない。この解釈に立てば、熱パルスのデータは閉じ込め時間がパワーの平方根に反比例するという経験則とは矛盾する。

摂動解析には注意点もある。熱パルス拡散率は一般にパワーバランス拡散率より大きく、両者の比がスティフネスの指標となるが、これは擾乱が背景プロファイルの勾配を「追加で」駆動するためであり、単純な拡散則からのずれの証拠でもある。W7-Xでの変調ECRH実験では、衝突頻度を掃引しながら熱パルス拡散率とパワーバランス拡散率を比較することで、TEM支配からETG/ITG支配への遷移や、準等ダイナミック配置によるTEM乱流の安定化が調べられている。次の図は、三つの磁場配置について衝突頻度に対して両者の拡散率を示したものである。

図の(a)がパワーバランス解析による電子熱拡散率、(b)が熱パルス伝播解析によるそれで、いずれも規格化衝突頻度に対して「高ミラー」（×印）、「高イオタ」（星印）、「標準」（丸印）の三配置が重ねて描かれている。パワーバランス値が配置間でほぼ重なるのに対し、熱パルス値は低衝突頻度側で配置差が顕著に現れる。

熱パルス拡散率とパワーバランス拡散率の比、すなわちスティフネスの指標そのものを衝突頻度に対して示したのが次の図である。

図の(a)はスティフネス比 χHP/χPB、(b)は熱パルスの振幅から得た拡散率と位相から得た拡散率の比（減衰の指標）であり、同じく三配置のデータが色分けされている。低衝突頻度側でスティフネス比が1を超えて増大し、衝突頻度の増加とともに1へ収束する傾向が読み取れる。

非局所輸送とアバランシェ型熱輸送

局所拡散モデルの前提は、ある点の熱流束がその点の温度勾配だけで決まるというものである。しかし実験はこの前提が破れる現象を複数報告してきた。代表例が非局所効果である。エッジを冷却すると、擾乱が到達するはずのない中心部の電子温度がむしろ急上昇するというもので、アルカトールC-ModのオーミックLモードでレーザー不純物注入によるエッジ冷却を用いて観測された。低衝突頻度のLOC領域ではこの現象が明瞭に現れ、中心拡散率の急激な低下か見かけの熱ピンチの発生かのいずれかで説明される。一方、高衝突頻度のSOC領域では輸送は「局所的」に戻り、中心温度は閉じ込め時間のスケールで低下する。興味深いことに、非局所性が現れる転移密度は回転反転の臨界密度と一致し、熱輸送と運動量輸送、そしてTEMからITGへの乱流支配の転移が結びついている可能性を示唆する。繰り返し冷パルスを注入してフーリエ解析した結果でも、LOCでは電子熱輸送が純拡散でないのに対し、SOCでは拡散的になることが示された。

非局所性を説明する主要な概念がアバランシェ型輸送である。これは、局所的な臨界励起が隣接領域を順次励起しながら長距離を伝播する間欠的な輸送イベントであり、砂山モデルに代表される自己組織化臨界性（SOC）の典型的な挙動である。砂山に砂を少しずつ加えると、勾配が臨界角を超えた場所から崩壊が連鎖し、個々の崩壊は局所的でも系全体としては長距離の砂の輸送が起こる。磁化プラズマでは、この「崩壊」が熱や粒子のバリスティックな輸送として現れる。フラックス駆動型シミュレーションでは、アバランシェがレジリエントな圧力プロファイルやE×B階段構造の形成に支配的役割を果たすことが報告されており、プロファイルの頑健性が臨界勾配型の拡散輸送ではなく非拡散的なアバランシェによって維持されうることが示唆されている。

実験的には、JT-60Uでアバランシェ型電子熱輸送の定量的評価が行われた。反射計で測定した密度変動強度の特徴的なバーストが電子温度擾乱と同期することを利用し、条件付き平均によって電子温度擾乱を抽出したところ、その擾乱が直径磁気ドリフト速度の数分の一の速度で、プラズマ半径の半分以上を一挙に伝播する様子が捉えられた。擾乱の時間スケールは1/f冪乗則に従い、SOCの特徴と整合する。重要なのはその輸送量であり、アバランシェによる熱流束は平均加熱パワーが高いほど急増し、高加熱極限では増分加熱パワーの80%に達する。これは、勾配を固定してもプロファイルが保たれる「レジリエント・プロファイル」の実体が、まさにこの長距離アバランシェ輸送であることを示すものである。次の図は、11 MW放電で条件付き平均により抽出された電子温度擾乱の伝播を示したものである。

図の(a)は密度変動の包絡線（−20 kHz〜−250 kHz帯）で、(b)〜(i)は参照半径 ρ = 0.49 と 0.54 から内外へ条件付き平均された電子温度変動 ΔTe/Te である。青い丸印と矢印は相互相関解析で求めた「山」と「谷」の時間遅れを示し、擾乱の塊が参照点から半径方向へ順次伝わっていく様子が見て取れる。

このように熱輸送の研究は、新古典論による下限の理解から出発し、微視的乱流の臨界勾配とスティフネス、そして局所拡散を超える非局所・アバランシェ的輸送へと広がってきた。炉設計の観点では、異常輸送をどこまで抑制できるか、そして非拡散的輸送がプロファイル制御にどの程度の自由度を残すかが、依然として中心的な未解決問題である。