Wiki / 現象 レイリー・テイラー不安定性 （Rayleigh-Taylor instability、Rayleigh-Taylor）は、核融合・プラズマ物理学の現象。密度勾配と加速度が逆向きの界面で生じる流体力学的不安定性で、慣性核融合のアブレーション面やZピンチの圧縮過程で成長する。これにより界面が混合・乱流化して核融合性能を低下させるため、その抑制が重要である。この用語集が追う論文のうち、タイトルにこの主題が現れるものは107本ある。比重が最も大きいのは直近の2025–2029年で、その期間に発表された論文の0.37%を占める。
目次
1 解説 2 推移 3 研究の変遷 4 この数字の作り方 解説 レイリー・テイラー不安定性（Rayleigh–Taylor instability、略してRTI）は、密度の異なる二つの流体の界面が、重い流体を軽い流体が押し上げる方向に加速度を受ける（あるいは重い流体が軽い流体に支えられる）ときに生じる流体力学的な不安定性である。界面にわずかな凹凸があると、重い流体の突出部（スパイク）が下方へ、軽い流体の突出部（バブル）が上方へ成長し、界面は次第に混ざり合う構造へと変形する。天体現象から慣性核融合 、Zピンチ・プラズマ に至るまで、加速される界面が存在するあらゆる場で現れる普遍的な不安定性であり、特に慣性核融合（ICF）ではターゲットの圧縮を破壊しうる最大の脅威の一つとされる。図は、レーザー加熱されるプラズマ中でこの不安定性を扱う典型的な計算設定を示したものである。
fig.1 Enhancement of self-similar growth in multimode ablative Rayleigh–Taylor instability by self-generated magnetic fields (2026) · CC BY 4.0
(a) The schematic diagram of the simulation setup. (b) The initial density and temperature profiles along the z axis.
図の(a)は、加速方向に密度の階段状の界面（約5 g/cm³の高密度領域と低密度領域の境界）を置き、そこへ熱流束を入射させるシミュレーション構成の模式図、(b)はその初期の密度と温度のプロファイルであり、密度が急峻に落ちるアブレーション 前面と、温度がなだらかに下がる高温コロナの共存、すなわちRTIが問題となる空間構造が読み取れる。
古典的なレイリー・テイラー不安定性と線形成長
最も単純な設定は、一様重力下で重い流体（密度 ρ_h）が軽い流体（ρ_l）の上に載っている静的な界面である。界面を波数 k のじょう乱で撹乱すると、じょう乱は指数関数的に成長し、その成長率は古典理論により
γ = √(A g k)
と与えられる。ここで g は加速度、A = (ρ_h − ρ_l)/(ρ_h + ρ_l) はアトウッド数である。重要なのは成長率が √k に比例すること、すなわち波長が短いじょう乱ほど速く成長することである。これは古典RTIでは波長の下限が存在しないことを意味し、実際には粘性や表面張力、あるいは流体の圧縮性が短波長側の成長を抑えるが、原理的には任意の短波長じょう乱が増幅される。物理的には、界面の「谷」では重い流体が軽い流体に置き換わるため系の位置エネルギーが下がり、その解放エネルギーがじょう乱を駆動する。バロクリニックな渦生成、すなわち密度こう配と圧力こう配が平行でなくなるときに生じる渦度が、界面の変形を自己増殖させるのが本質的な機構である。
この古典的な枠組みは、加速される界面一般にそのまま翻訳できる。慣性系で見れば、加速度 g は等価な重力として働き、重い流体が「下」、軽い流体が「上」に来る配置であれば界面は必ず不安定になる。
レイリー・テイラー不安定性
種別 現象
論文数 107 本
初出 1970 –1974 年（2本）
最盛期 2025 –2029 年 0.37%
ほかの綴り Rayleigh-Taylor 索引にもどる
アブレーション安定化と慣性核融合におけるRTI ICF では、レーザー光やX線がターゲット表面を加熱して激しいアブレーション （表面の蒸発・吹き出し）を起こし、その反作用でシェルを加速する。この質量アブレーションは単なる駆動機構ではなく、RTIに対する強力な安定化機構として働く。アブレーション前面の参照系で見ると、プラズマは界面を横切って連続的に流れ去っており、界面のじょう乱はこの流れによって運ばれ、成長が遅らされる。渦度がアブレーション流によって界面から対流で持ち去られること（渦度の対流）と、アブレーションによる動的な過剰圧力がじょう乱を押し戻すことが、主要な安定化機構として理論的に同定されている。
その結果、線形成長率は古典的な √(gk) から修正され、しばしばタカベの公式と呼ばれる経験式
γ = α √(gk / (1 + k L_m)) − β k V_a
でよく近似される。ここで L_m は密度こう配の最小スケール長、V_a はアブレーション速度、α と β はアブレータ材料に依存する無次元係数である。第2項はアブレーションによる安定化を表し、波数 k が大きい（波長が短い）ほど効く。そのため成長率はある波数で最大となり、それより短い波長では成長率がゼロを横切る。このカットオフ波長より短いじょう乱は線形に完全に安定化される。古典RTIでは成長率が k とともに単調に増えたのに対し、アブレーション性RTI（ARTI）では短波長じょう乱が抑制される点が本質的な違いである。
この安定化の程度は熱伝導の詳細に強く依存する。電子の平均自由行程が温度こう配スケール長の1%程度を超えると、古典的なシュピッツァー・ヘルムホルツ（局所）熱伝導は破綻し、熱流は非局所的になる。非局所輸送はアブレーション前面の内側への予熱を通じて V_a と L_m を増大させ、ARTIの成長をさらに抑える方向に働く。数値シミュレーションでは、非局所効果を取り入れると実測成長率が理論予測より下がる実験的事実が説明され、特に短波長じょう乱に対する抑制が強まりカットオフ波長が長波長側へ移ることが示されている。一方で、単に V_a と L_m を増やすだけでは非局所輸送による安定化を完全には説明できず、非局所熱流束そのものがじょう乱を直接減衰させる効果を含めた線形理論が必要とされる。ARTIの成長率予測と実験観測の間には依然として溝が残っており、非局所輸送の影響は現在も活発に研究されているテーマである。
Zピンチおよび磁化プラズマにおける磁気レイリー・テイラー不安定性 Zピンチ では、軸方向に流す大電流が作る方位角磁場のローレンツ力（j×B力）がプラズマライナーを内側へ加速する。この加速界面、すなわちインプロージョンするプラズマの外表面は、軽い磁場領域が重いプラズマを押し込む構造になっており、磁気レイリー・テイラー不安定性（MRTI）が不可避的に発生する。物理的には古典RTIと同じくらせ状の界面じょう乱が成長するが、磁場の張力がじょう乱を引き戻す方向に働くため、磁場の向きとじょう乱波数の関係によって成長率が修正される。特に軸方向（縦）磁場を外部から印加すると、磁力線の張力が短波長じょう乱を安定化させ、成長を大幅に遅らせることが理論・実験の両面で示されている。
MRTIはZピンチを用いた慣性核融合 、X線放射源、中性子源 など、高速Zピンチのあらゆる応用にとって支配的なリスクである。特にZピンチ・ダイナミック・ホールラウムでは、MRTIの発達が放射温度を直接左右するため、その制御が成否を分ける。線形理論では単一モードの成長はインプロージョン距離のみで決まるが、高速Zピンチでは不安定性が速やかに強非線形段階へ進むため、線形・弱非線形理論ではインプロージョン全過程を追えない。実験的には、電流立ち上がり時間が長く初期半径が大きいインプロージョンほどMRTIが完全に発達しやすいこと、すなわち非線形じょう乱振幅がインプロージョン速度を固定すれば立ち上がり時間（または初期半径）にほぼ線形に比例して増大することが、シミュレーションと理論から示され、ワイヤアレイ実験によっても検証されている。図は、ワイヤアレイZピンチ実験で観測された停滞前のX線自己発光像と、二次元放射磁気流体シミュレーションで再現された像を並べたものである。
図の左と中の列が異なるショットでの実験データ、右の列が対応するシミュレーション結果であり、上から下へインプロージョンが進むにつれ（各パネルの数値は停滞までの残り時間をnsで示す）、ライナー表面に軸方向に伸びるMRTIのじょう乱構造が発達していく様子が、実験・シミュレーションともに一貫して見て取れる。
磁化プラズマ 一般では、磁場とRTIの相互作用は多様である。有限ベータプラズマでは二流体効果や有限ラーモア半径 効果が成長率に影響し、磁気剪断や粘性もモードの安定性を変える。部分電離プラズマやダストプラズマ 、強く相関のある量子プラズマ といった極限状態でもRTIの線形理論が展開されており、磁場がじょう乱を安定化するか逆に増幅するかは、熱伝導の異方性や磁力線の曲がりを介した熱流の変形など、複数の経路に依存して決まる。
非線形発展：バブル、スパイク、自己相似混合
線形段階を過ぎると、RTIの界面は特徴的な非線形構造へ発達する。重い流体は指状の「スパイク」となって軽い流体中を落下し、軽い流体は丸い「バブル」となって重い流体中を上昇する。単一モードでは、バブルの貫通速度はポテンシャル流理論が予測する一定の終端速度に一旦飽和するが、その後再加速される現象が古典・アブレーション 性の両RTIで報告されている。この再加速は、バブル内部への渦度の蓄積によって説明される。アブレーション性RTIでは、アブレーションがスパイク近傍の質量流出を減らしバブル幅を狭めることでバブル内部の渦度を強め、バブルの再加速を促進する（渦加速機構）。すなわちアブレーションは線形段階では安定化に働くが、非線形段階ではバブル貫通をむしろ速めるという逆説的な役割を果たす。
多モードの状況、すなわち現実の実験に近い多数のじょう乱モードが共存する状況では、バブル同士の競合と合体が支配的となる。バブル貫通距離 h_b は、自己相似段階で
h_b = α_b A S
というスケーリング則に従う。ここで S は時間積分 ∫√g dt の二乗で定義される長さ次元の量、A はアトウッド数、α_b は自己相似成長係数である。α_b はICF ターゲット設計上の最重要パラメータの一つであり、シェルの完全性を保つために許容される最大の飛行中アスペクト比を決めるからである。α_b の初期じょう乱振幅 h_0 への依存性に応じて、古典RTIには二つのレジームがある。h_0 が小さい「バブル競合」レジームでは α_b が h_0 に敏感で、浮力・抗力モデルから導かれた式で近似される。h_0 が臨界値を超えると「バブル合体」レジームへ移行し、α_b は h_0 に依らなくなる。
アブレーションが存在する多モードARTIでは、この様相が変わる。アブレーションは小さい h_0 では線形安定化を通じて α_b を古典値より下げるが、大きい h_0 ではアブレーション駆動の渦がバブル速度を加速し、バブル競合からバブル合体への移行を妨げるため、十分大きな h_0 では α_b が古典値を上回る。つまり多モードARTIは事実上常にバブル競合レジームに留まる。さらに、レーザープラズマ 中ではビアマン電池機構によってRTI自体が磁場を自己生成することが知られており、プロトンラジオグラフィ実験はメガガウス級の自己生成磁場を観測している。この磁場は電子の熱流を磁化させることでアブレーションを変調し、スパイク近傍のアブレーションを弱めバブル幅を狭めて渦度を強め、短波長モードの非線形バブル貫通速度を高める。多モードシミュレーションでは、自己生成磁場がアブレーション安定化を弱めて α_b を増大させ、混合層の発達を加速することが示されている。図は、同じ多モードARTIを磁場のフィードバックあり・なしで比較した密度コンターである。
図の上段(a)(b)は磁場フィードバックなし、下段(c)(d)はありの場合で、左列が非線形初期（h_b = 4.8 μm）、右列が深非線形段階（h_b = 16.5 μm）に対応する。磁場ありではバブルが細く鋭く突出し、スパイク構造がより顕著になっているのが分かる。
緩和と制御の戦略
RTIの成長を抑えることは、ICF やZピンチ の性能を左右する中心的な課題であり、複数の戦略が発展してきた。第一は駆動プロファイルの最適化である。Zピンチでは、電流立ち上がり時間や負荷半径を短く・小さく抑えることで、MRTIが完全に発達する前に停滞に至らせることができる。実際、同じインプロージョン速度でも立ち上がり時間が50 ns、100 ns、150 nsと変わるだけで非線形じょう乱振幅が系統的に増大することがワイヤアレイ実験で確認されており、パルス形状の設計がそのまま不安定性制御の手段となる。ガスパフ Zピンチでは初期密度プロファイルの工夫によってMRTIの成長が緩和された例もある。
第二は磁場の利用である。軸方向磁場の印加はZピンチのライナー・ガスパフ両実験でMRTIの成長を実際に低減することが示されており、磁気化ライナー慣性核融合（マグリフ、MagLIF）では構造材料の強度や小半径と相まって不安定性の発達を抑えている。ただしICF文脈では磁場の効果は一方向的ではない。外部磁場は電子熱伝導を抑制して質量アブレーション を減らし、その結果アブレーション安定化を弱めてARTI成長をかえって増幅するレジームがシミュレーション・実験で示されている。一方で磁場印加はホットスポットの磁気閉じ込めを通じて収縮性能を向上させる効果もあり、NIF での間接駆動実験やマグリフ、磁化直接駆動実験で融合収率・ホットスポット温度の向上が確認されている。すなわち磁場は収縮性能には有利だが特定の条件下ではRTIには不利であり、そのトレードオフの定量的評価が設計上の課題となっている。
第三はターゲット・負荷構造そのものによる制御である。Zピンチでは多層シェル負荷や誘電体コーティングがMRTIの抑制に用いられてきた。ICFでは、アブレーション安定化を最大限に活かすため短波長じょう乱をカットオフ以下に抑える表面仕上げ、熱伝導の異方性や非局所輸送を利用した安定化、そして多モード発展を支配する α_b を小さく保つ駆動設計が組み合わされる。熱伝導の異方性がARTIに安定化効果を持つこと、流れの曲率が磁化プラズマ の線形RTI成長率に影響することなど、制御の観点からの効果も理論的に整理されつつある。RTIは原理的に排除できない不安定性である以上、その成長をどこまで許容し、どの機構で抑えるかという設計の最適化こそが、この分野の研究の核心であり続けている。
レーザー核融合とレーザープラズマの文脈で、粘性、磁気シア、密度勾配、自己生成磁場がレイリー・テイラー不安定性に及ぼす影響が調べられている。直接駆動レーザー核融合や多層構造標的における不安定性の研究が特徴的である。
レイリー・テイラー不安定性（Rayleigh-Taylor instability）の研究史 — FusionPapers
![Comparison of the density contours in case iii without ((a) and (b))/with ((c) and (d)) magnetic feedback. (a) and (c) are in the early nonlinear stage [h_b = 4.8~\mu \mathrm{m}]. (b) and (d) are in the deeply nonlinear stage [h_b = 16.5~\mu \mathrm{m}].](https://fusionpapers-figures.s3.ap-northeast-1.amazonaws.com/nfae315cf2.jpg)