図の左端に、薄い炭素の外殻(厚さ約40 μm)に包まれた直径3 mm級のDTペイロードからなる殻付きペレットの断面が示してあり、中央と右には、DIII-D放電の電子密度・電子温度プロファイルを初期条件として計算された粒子堆積プロファイルと、ペレット半径のプラズマ中での変化が並んでいる。ペレットがプラズマ内部へ進むにつれて半径が減少し、それに対応して燃料粒子が堆積していく様子が読み取れる。この種の予測計算は、ITERやCFETRのような燃焼プラズマ装置で必要とされる高密度化(グリーンウォルド密度に対する比が1近傍、場合によってはそれ以上)を実現するためのペレット燃料供給設計の基礎となっている。
磁化プラズマ中のペレットアブレーション:遮蔽、雲、ストリエーション
高温プラズマ中のペレットアブレーションの基本的な機構は次の通りである。ペレット表面は周囲の電子の熱流で加熱されて昇華し、生じた中性気体はさらに電離されて、ペレットを取り囲む「アブレーション雲」と呼ばれる低温・高密度のプラズマ層を形成する。この雲の密度は背景プラズマより典型的に3桁ほど高く、温度ははるかに低いため、電子の熱流は雲の中で放射冷却とともに強く減衰する。ペレット自身はこの雲に「遮蔽」されて直接の熱流から守られており、アブレーション率は遮蔽の程度によって自己制御される。遮蔽が弱ければ熱流が増えてアブレーションが加速し、雲が厚くなって再び遮蔽が強まる、という負帰還が働く。中性気体とプラズマの両段階の遮蔽を含むスケーリング則が理論的に導かれており、ペレットアブレーション率の予測の骨格となっている。
雲は磁場中にあるため、その振る舞いは磁場の勾配に強く支配される。雲内部の圧力凸部と磁場勾配の組み合わせが電場を作り、E×Bドリフトによって雲の物質が大半径方向へ移動する(∇Bドリフト)。このドリフトが堆積プロファイルをペレット軌道からずらすため、正確な燃料供給予測にはアブレーション相とその後の均質化相の両方を扱う必要がある。雲の形状は磁力線に沿って伸びた偏平な構造をとり、横方向の膨張は磁場に抑えられ、磁力線方向の膨張が支配的となる。
図はシミュレーションで得られたアブレーション雲の幾何形状を模式的に示したもので、磁力線方向(Z)に伸びた偏平な楕円体の雲と、ペレットの進行方向および∇Bドリフト方向の角度関係が描かれている。雲の粒子内容・寸法・密度・温度は、その後のドリフト相の初期条件として堆積予測の精度を直接左右する。
ペレットアブレーションで最も古くから知られ、いまだ機構が確定していない現象がストリエーション不安定性である。アブレーション雲の発光はほぼ常に強く変調されており、高速撮影された雲の像には磁力線に沿って伸びた明暗の帯が等間隔に並ぶ。これはアブレーション率そのものが周期的に変動していることを示唆し、堆積プロファイルが軌道に沿って時間的・空間的に周期的な構造を持つことになるため、燃料供給の予測精度に直結する。説明仮説としては、有理磁気面でのエネルギー枯渇、雲が「中空の管」状になりペレットが管壁を横切るたびに遮蔽が強弱するという模型、雲の静電ポテンシャル分布がもたらす高速回転に伴うレイリー・テイラー不安定性、そして∇Bドリフトで移動する電離雲の層が新たにアブレーションされた中性気体の膨張を拘束するという模型などが提案されてきた。有理面の効果は支配的でないことが示され、中空管模型は雲が磁力線に凍りつくという仮定の妥当性が疑われ、回転模型は∇Bドリフトを無視する点が課題とされ、3次元磁気流体シミュレーションは上記の中性気体拘束機構によるアブレーション率の振動を再現しない、といった具合に、各仮説は部分的な支持と反証を受けている。雲の直径を∇Bドリフト速度で割った時間は観測されるストリエーション周波数とほぼ一致することが知られており、ドリフト動力学が変調の時間スケールを設定することはほぼ確実だが、不安定性の本体は未解明のままである。
図はASDEX Upgradeでのペレット注入時にアブレーションモニターで記録された発光の時間変化であり、隣接する極大の間隔から読み取れるストリエーション周波数(61、91、118 kHzなど)が灰色の帯で示され、挿入図にはフーリエスペクトルが併せて示してある。ペレットがより高温のプラズマ内部へ進むにつれて変調周波数が上昇する(チャープする)傾向が観測されており、これは周波数が局所電子温度や加熱出力に依存することを示す指標として、各模型の検証に用いられている。
実験データの体系化も進んでいる。複数のトカマク(JET、TFTR、DIII-D、ASDEX Upgrade、T-10など)から水素・重水素ペレットの質量・速度、プラズマプロファイル、発光データを集めた国際ペレットアブレーションデータベースが整備され、アブレーション理論のスケーリング解析と検証の共通基盤となっている。また、ペレットが通過したと考えられる磁気面の位置から列状のフィラメント(細長い発光領域)が順次形成されることが観測されており、フィラメント位置がペレットの弾道運動と時間依存アブレーション模型から予測できること、すなわち定常模型ではなく時間依存の取り扱いが必要なことが示されている。
慣性核融合におけるレーザーアブレーション:アブレーション圧力と面のダイナミクス
慣性核融合では、数mm級の燃料カプセル表面に強力なレーザーを照射し、表面物質をアブレーションさせる。吹き飛ばされる物質の運動量の反作用としてアブレーション圧力が発生し、これがカプセル内へ衝撃波を駆動して燃料を圧縮・加熱する。この意味で、レーザーアブレーションの圧力と質量アブレーション率は慣性核融合の駆動効率そのものである。
レーザーアブレーションの構造は、レーザー吸収が起こる臨界面、そこからアブレーション面までの熱伝導帯、そして物質が実際に剥ぎ取られるアブレーション面、の三層で理解できる。レーザーエネルギーは臨界面付近で吸収され、電子の熱伝導によって熱が内側のアブレーション面へ運ばれ、そこで物質が加熱・加速されて外向きに噴出する。この熱伝導は古典的なスピッツァー・ヘルム模型の想定より平均自由行程が長く、熱流は局所的な勾配だけで決まらない非局所性を示す。実験データの説明にはしばしば熱流に上限(フラックスリミッター)を設ける近似が用いられてきたが、一定値や純粋に時間依存のリミッターではレーザーアブレーション全過程の時空間発展を同時に再現できないことが示されている。時空間に依存する非局所熱輸送模型を用いた放射流体シミュレーションによれば、非局所的な熱流の再分配により、臨界面では熱流がより強く抑制され、伝導帯では横方向輸送が抑えられ、アブレーション面の前方では予熱が起こる。これらが合わさってレーザー吸収効率の向上、伝導帯の伸長、質量アブレーション率と衝撃波速度の増大をもたらす。非局所性が単なる補正項ではなく、駆動効率を質的に変える要素であることを示す結果である。
アブレーション圧力の時間発展は、レーザー強度の時間整形(ランプ圧縮など)と組み合わせて、衝撃波の合体を制御する設計パラメータとして扱われる。指数則によるアブレーション圧力の近似がどの条件範囲で有効かの検証も行われており、圧力の正確な予測が圧縮設計の前提となっている。また、アブレーションプラズマ自体も研究対象であり、レーザーアブレーションで生成したプルームの大気中での膨張、磁場による閉じ込め、帯電粒子の最大放出フラックスなどは、標的照射の副次現象としてだけでなく、物質分析や核融合炉内のトリチウム・ダスト除去への応用も視野に入れられている。
アブレーション面の安定性とレイリー・テイラー問題
慣性核融合でアブレーションが本質的に関わるもう一つの問題が流体不安定性である。圧縮の際、低密度の駆動プラズマが高密度の燃料を押し込む構図は、重力場中で軽い流体が重い流体を支えるのと同じく、レイリー・テイラー不安定性の古典的な増幅条件を満たす。表面の微細な凹凸は加速の間に指数関数的に成長し、最悪の場合には燃料と周囲物質の混合を引き起こして点火を妨げる。
アブレーションがある場合、この不安定性は純粋なレイリー・テイラーとは異なる。アブレーション面から物質が絶えず流出し、その質量流が凸部を「削り」、さらに熱伝導による平滑化が働くため、成長率は古典値から下方修正される。波数の高い(波長の短い)摂動ほどアブレーションによる安定化が効き、成長率はある波数で最大をもち、それ以上では抑制される。この「アブレーション的に安定化されたレイリー・テイラー」の線形成長率と非線形段階の振る舞いを、熱伝導・放射・流体の自己無撞着な取り扱いで解析する理論が整備されており、高モード(高波数)摂動に対する標的感度の評価に用いられている。実際、高波数の表面粗さが最大の線形増幅率をもつモードで非線形段階に達するとアブレーション面が大きく変形することがシミュレーションで示され、これを遅らせる手段としてプリアブルス(前駆パルス)による整形が有効であることが確認されている。標的設計における安定性評価は、アブレーション面の物理と直接に結びついている。
なお、ペレットアブレーション雲の回転に伴うレイリー・テイラー不安定性がストリエーションの説明候補として提案されていることは前節で述べた通りで、アブレーションとレイリー・テイラー型不安定性の結びつきは磁化プラズマ側にも現れる。
道具としての、そして診断としてのアブレーション
アブレーションは核融合研究において能動的な手段としても使われる。第一の用途は燃料供給であり、ガス注入では達成できない中心密度の上昇をペレット注入が実現する。反応器設計では、ペレットの貫通深度と堆積位置をアブレーション模型で予測し、必要な注入頻度とペレットサイズを決定する。反応器級の高温ペデスタルでは深い貫通が難しいため、慣性核融合向けに開発されたペレット製造技術を応用し、耐熱性の固体殻でペイロードを覆って強度を高めるという案も検討されている。
第二の用途は急速停止(ディスラプション緩和)である。不純物(例えばアルゴン)の小型ペレットを注入すると、そのアブレーションと放射損失によって放電の温度崩壊と急速シャットダウンが誘起される。ここで興味深いのは、アブレーション率が診断として使える点である。熱プラズマの電子温度から予測される値を大きく超えるアブレーション率が観測された場合、その超過分は温度崩壊時に生成した非熱的脱走電子による追加加熱の証拠と解釈できる。実際、この超過アブレーション率から推定した脱走電子シード電流は1〜10 kA程度となり、雪崩理論に基づく推定とほぼ整合する一方、ドライサー式の解析推定は過小、ホットテール模型は過大となる。アブレーション率がプラズマ中の非熱的電子集団の定量診断として機能する例である。
第三に、アブレーション発光そのものが最も手軽なペレット診断である。アブレーション雲の発光強度は局所アブレーション率に比例するという仮定のもとで、発光の時間変化からペレットの貫通深度や堆積位置が推定されてきた。しかしこの比例関係の根拠は必ずしも確立しておらず、雲のスペクトルと発光モデルから合成スペクトル・合成像を作り実測と比較する検証が進められている。局所アブレーション率の時間発展と雲の輝線スペクトル(可視域の水素系列線や不純物線)を同時に計算し、フィルタ透過特性やカメラの伝達関数を畳み込んで実測像と突き合わせることで、発光強度とアブレーション率の関係、雲の密度・温度・粒子内容が定量的に検証されつつある。アブレーション雲の特性は局所プラズマ条件のみで決まり磁場配位には依存しないと考えられており、ヘリカル装置で得られた検証結果はトカマクにも適用できる。
このほか、レーザーアブレーションは炉内のトリチウムやダストの在庫管理・除去、再堆積炭素層の熱酸化による粒子放出評価、極端紫外レーザーを用いた飛行時間質量分析など、核融合炉の運転・保全に関わる応用にも広がっている。アブレーションは、燃料供給・圧縮駆動・急速停止・診断という、核融合プラズマ研究の複数の要所を貫く現象なのである。