ジャイロ運動論

ジャイロ運動論（gyrokinetic theory）は、強く磁化されたプラズマ中の荷電粒子の運動を記述するための漸近的理論体系であり、磁気核融合研究における乱流・輸送解析の事実上の標準手法となっている。その基本思想は、ラーモア回転という最も速い周期運動を解析的に平均化してしまい、粒子を「荷電リング」あるいはガイドセンター（旋回中心）の分布関数として扱う点にある。これにより、時間ステップや格子幅をサイクロトロン周期とラーモア半径に縛られることなく、乱流の成長に支配的な低周波・長波長ダイナミクスだけを数値的に追跡できる。核融合プラズマではラーモア半径 ρ は平衡場の特性長 L に比べて十分小さく（ε = ρ/L ≪ 1）、この小パラメータに基づく摂動展開が理論の出発点となる。

ジャイロ運動論のオーダリングとジャイロ平均方程式

ジャイロ運動論は、ε = ρ/L ≪ 1 を小ささの尺度とする系統的なオーダリングから出発する。分布関数は平衡マクスウェル分布 F₀ と一次の摂動部分に分けられ、摂動は電位 φ に対するボルツマン応答と、非ボルツマン成分 h とに分解される。h はガイドセンター座標で表すとジャイロ角に依存しない。これが理論の核心で、磁場に垂直な速い回転運動は完全に消去され、残るのは平行運動と、磁場曲率・勾配に起因するドリフト（磁ドリフト v_D）および E×B ドリフト v_E である。

支配方程式であるジャイロキネティック方程式は、ガイドセンター分布 h に対する輸送方程式の形をとる。位相空間の流れは平行運動、ミラー力、ドリフト速度から構成され、電磁場はジャイロ平均 ⟨·⟩ を施された形で効く。ジャイロ平均とは、ガイドセンターを固定してジャイロ角について 2π にわたり積分する操作であり、これにより場の空間変化はラーモア半径程度のスケールで平滑化される。この平滑化が、例えばイオン温度勾配モードのような波長が ρ_s 程度のドリフト波乱流を、イオンのラーモア運動を解かずに記述することを可能にする。電荷と電流の自己無撞着性は、ジャイロ平均された電荷・電流から準中性条件とアンペール則（並進不変量としての平行アンペール則）を通じて求められる。

電磁的な拡張では、摂動磁場の平行成分 δB_∥ と、磁力線に沿ったベクトルポテンシャル成分 A_∥ を場の変数に加える。平衡電流を含む完全な定式化では、線形長波長極限で理想磁気流体力学（MHD）理論が正しく再現され、バルーニングモード、キンクモード、シェアアルフヴェン波がジャイロキネティック方程式系から導かれる。これは、ジャイロ運動論が単なる乱流理論ではなく、圧力勾配駆動・電流駆動の両方の不安定性を統一的に扱う運動論的枠組みであることを示している。一方でオーダリングの帰結として、ジャイロ運動論は原理的に高周波・短波長の物理（サイクロトロン共鳴そのもの、電子の慣性スケールを超える効果の一部）を含まず、適用範囲は ε ≪ 1 が成り立つ領域に限られる。

数値手法：デルタf・フルf の粒子法とオイラー法

ジャイロキネティック方程式は 5 次元位相空間（3 次元空間＋平行速度＋磁気モーメント）上の非線形輸送方程式であり、6 次元の完全運動論に比べれば格段に扱いやすいとはいえ、依然として大規模数値計算を要する。主要な離散化手法は二つに大別される。

第一は粒子インセル（PIC）法である。ジャイロキネティック方程式を特性曲線に沿ったガイドセンターの運動方程式として解き、分布関数を多数のマクロ粒子で表現する。摂動が平衡に比べて小さい場合にはデルタf 法が用いられ、平衡部分を厳密に差し引いた摂動分布 δf だけを粒子で運ぶことで統計ノイズを大幅に抑える。一方、平衡プロファイル自体が乱流輸送によって時間発展する状況や、定常状態からの大きな偏りを扱う場合には全分布を運ぶフルf（total-f）法が使われる。電磁場は格子点上のポアソン方程式・アンペール方程式で解き、粒子と格子の間で電荷・電流を分配する。PIC 法では、平衡電荷と摂動電荷の大きな桁の差を正確に差し引く「キャンセレーション問題」が電磁計算特有の難点として知られ、これを回避するための定式化や離散化の工夫が発展してきた。

第二はオイラー的手法で、位相空間を格子で直接離散化し、方程式を偏微分方程式として解く。ノイズが本質的にない代わりに 5 次元格子の記憶容量と計算量が課題となり、局所（フラックスチューブ）近似で次元を落として用いられることが多い。どちらの系統でも、コード間の比較（ベンチマーク）によって線形成長率や非線形飽和レベルの一致が確認されており、手法間の差異は物理ではなく離散化の問題として管理されている。

衝突の取り扱いも数値手法の重要な一部である。ジャイロキネティック方程式自体は衝突なしでも定式化できるが、乱流の定常状態を得るには散逸項が必要であり、弱衝突性プラズマであっても衝突は飽和機構として本質的な役割を果たす。実装される衝突演算子は通常、線形化されたフォッカー・プランク演算子を簡略化したものであり、運動量とエネルギーの保存を保証する近似が施される。しかしこの簡略化は衝突効果が強い場合に不正確になることが知られており、高衝突性の領域や高電荷の不純物イオン（衝突周波数は電荷の 4 乗に比例する）を扱う場合には、より忠実な演算子の実装が進められている。ヘリカル系では最適化された装置（例えばウェンデルシュタイン 7-X）で新古典輸送が強く抑制され、乱流輸送が新古典輸送と同程度になりうるため、衝突を含むジャイロキネティック乱流計算の重要性が相対的に高まっている。

微小乱流、輸送、ゾーナルフロー

ジャイロキネティック理論の最大の応用分野は、温度・密度勾配によって駆動される微小乱流による異常輸送の理解と予測である。イオン温度勾配（ITG）モード、電子温度勾配（ETG）モード、捕捉電子モードなどが代表的な駆動機構であり、これらはドリフト波型の微不安定性として、勾配が臨界値を超えると成長する。非線形段階では乱流渦が径方向に熱・粒子を輸送し、その輸送係数は古典・新古典拡散を桁違いに上回りうる。

非線形飽和を支配する鍵概念がゾーナルフローである。これはフラックス面平均された帯状の E×B 流で、乱流自身が生成する。生成機構は乱流渦同士の非線形相互作用（レイノルズ応力）であり、生成されたゾーナルフローのせん断が乱流エディを径方向に引き裂いて相関を失わせることで、乱流レベルと輸送が自己規制される。このフィードバック機構により、勾配を上げても輸送があまり増えない「スティフネス」がコアプラズマの特徴として現れる。ゾーナルフローのダイナミクスが乱流飽和を決定する様子は、ジャイロキネティック計算で電位の時間発展を追うと直接観察できる。乱流エディが発達すると帯状構造の流れが成長し、エディがせん断によって崩壊していく過程が可視化される。

図はイオン温度勾配乱流における静電ポテンシャルの時間発展を三つの時刻で示したもので、左から右へ乱流エディが発達し、やがて帯状の構造（ゾーナルフロー）がエディを崩壊させていく過程が見て取れる。ゾーナルフローのダイナミクスが電磁 ITG 乱流の飽和を決定している。

乱流計算には局所近似と大域近似の二つの流儀がある。局所計算は一つのフラックス面上の周期箱で乱流を解き、輸送係数を効率よく求める。大域（global）計算はトーラス全体の径方向構造を解き、プロファイルの緩和、平衡シアや磁気配位の効果、モード間の大域的結合を扱える。磁気配位の効果は特にヘリカル系で重要であり、例えばウェンデルシュタイン 7-X の磁気配置空間における ITG 乱流の大域計算や、ヘリカル系の多イオン種プラズマでの乱流輸送計算が行われている。また、トライアングラリティ（断面の三角形度）のような断面形状パラメータが温度勾配駆動乱流に与える影響を系統的に調べる研究も進んでおり、負のトライアングラリティが乱流を安定化する方向に働くことと合わせ、形状制御による輸送制御の指針が得られている。実験との比較では、位相コントラストイメージングのような乱流計測とシミュレーションの直接照合や、多イオン種プラズマでの密度ピーク形成の再現などが行われている。

多スケール電磁物理：高速粒子、アルフヴェンモードと MHD 結合

燃焼プラズマでは、核融合生成のアルファ粒子や中性粒子入射・高周波加熱に由来する高速イオンが大量に存在する。これらは熱プラズマの乱流とは全く異なるスケールのマクロモード、すなわちアルフヴェン固有モード（AE）やエネルギー粒子モードを駆動しうる。ジャイロ運動論は、熱プラズマの微小乱流（周波数・波長とも小さい）と、高速粒子が共鳴して駆動するマクロなアルフヴェンモード（周波数はイオンサイクロトロン周波数よりはるかに低いが乱流より高い）を、同一の方程式系で扱える数少ない枠組みである。このため、微視的乱流と巨視的モードの非線形結合を調べる多スケール計算の基盤となっている。

アルフヴェン固有モードの側から見ると、トーラス固有モード（TAE）、逆シア固有モード（RSAE）、キネティックバルーニングモード（KBM）などが線形スペクトルとして解析され、非線形飽和機構が理論・計算の両面から研究されている。例えば RSAE は安全係数 q プロファイルの極小値近傍で励起され、電子非線形ダイナミクスによる飽和や、ITG 乱流との波－波相互作用によるエネルギー交換が調べられている。ドリフト波乱流はアルフヴェンモードを短波長のキネティックアルフヴェン波へ散乱して電子ランダウ減衰で消散させる経路を提供しうる一方、アルフヴェン固有モードが駆動するゾーナルフローがイオンスケール乱流を抑制する効果も見いだされている。高速粒子と乱流の共鳴相互作用による乱流抑制は、JET の重水素－三重水素実験を対象とした大域シミュレーションで定量的に評価されており、波－粒子共鳴による熱プラズマの乱流輸送の大幅な低減、電磁ベータ安定化による追加の抑制、そして不安定な TAE が駆動するゾーナルフロー・ゾーナルカレントによる深部コアでの輸送低減が分離・定量化されている。

図は大域電磁ジャイロキネティック計算で得られたトーラス方向モード数 n = 8 の静電ポテンシャルの小断面であり、コア領域に局在したモード構造が最外磁気面の内側に描かれている。大域計算ではこのように平衡配位の実形状の中でモードの空間構造が直接得られる。

さらに電磁ジャイロキネティック模型は、平衡電流を含む定式化により、理想 MHD の長波長極限（バルーニング・キンクモード、シェアアルフヴェン波）を回復する。これにより、衝突なしテアリングモードのような MHD 摂動と電磁乱流の共存・相互作用、圧力勾配駆動と電流駆動の不安定性の非線形結合を、一つの枠組みで計算できる。トカマクとヘリカル系の双方で、電磁乱流の大域粒子シミュレーションの実行可能性が示され、アルフヴェン固有モードと乱流が共存する系での周波数変化や、乱流が AE を調節して高速粒子輸送を変化させる効果などが調べられている。一方、エッジの H モードペデスタルや ELM（エッジ局在モード）物理では、ジャイロキネティック・ドリフトキネティック模型が MHD 模型（ペリング・バルーニングモードの線形解析、非線形 MHD 計算）や新古典輸送模型と並び、ペデスタルの微不安定性と輸送を記述する主要な道具となっており、MHD・ジャイロ流体・ジャイロキネティック各模型の統合が予測モデル構築の方向として意識されている。

検証・妥当性確認と予測モデリング

ジャイロキネティック・シミュレーションが設計根拠として信頼されるためには、検証（ベリフィケーション（verification）：方程式を正しく解けているか）と妥当性確認（バリデーション（validation）：解いた結果が実験と一致するか）が不可欠である。検証の側では、線形固有値問題の収束試験、コード間ベンチマーク、衝突演算子の精度試験（例えばスピッツァー問題の輸送係数の再現）などが行われる。妥当性確認の側では、トカマク実験で計測された乱流揺らぎの振幅・スペクトル・相関長と、シミュレーションを診断装置の応答に通した「合成診断」と比較する手法が確立されている。DIII-D での L モード検証研究では、コア中部（正規化半径 0.4～0.75 程度）では輸送と揺らぎレベルの予測が実験とよく一致する一方、外側領域ではエッジとの結合や電子・イオンスケールの同時解像が必要になることが示された。Alcator C-Mod のオーミック・I モードプラズマ、ASDEX Upgrade の負トライアングラリティ放電など、多様な条件下で同様の比較が行われている。

もう一つの実務的な流れは、非線形計算の結果を縮約した準線形輸送モデルである。非線形ジャイロキネティック計算は高価であるため、線形成長率と非線形飽和の経験則を組み合わせて輸送係数を高速に与える準線形モデルが構築され、粒子輸送や角運動量輸送のモデリング、プロファイル予測（与えられた加熱・燃料供給条件下での温度・密度分布の予測と実験検証）に用いられている。ペデスタル構造の予測モデル（幅・勾配・高さの三要素で特徴づける枠組み）も、MHD 安定化解析と微不安定性・新古典輸送のジャイロキネティック的知見を統合する形で発展している。これらの取り組みは、ITER のような燃焼プラズマでの第一原理に基づく予測能力の確立を共通の目標としており、ジャイロ運動論はその理論的・数値的中核であり続けている。