衝撃波

衝撃波とは、流体やプラズマ中を伝わる圧縮性の不連続面であり、その面を境に密度・圧力・温度・流速が急峻に跳び、流れのエントロピーが増大する構造である。音速より速い圧縮擾乱が背景の流れに追いつき、エネルギーが散逸する形で集積されることで形成される。プラズマ物理学では、衝撃波は惑星圏や太陽風から超新星残骸、レーザー核融合や慣性閉じ込め核融合(I Cf)の標的内部に至るまで、ほぼすべてのスケールに現れる最も普遍的な現象の一つであり、その散逸機構と構造の理解は流体力学・統計力学・プラズマ物理が交差する中心課題である。

図は、衝撃波面の上流と下流における運動論的構造を示したものである。上段には電子密度の空間分布が、下段には位相空間におけるイオン分布が描かれ、上流域の高速度イオンがポテンシャル障壁で反射して「反射ビーム」を形成すること、下流域ではトラップされたイオンが「トラッピング渦」を構成することが見て取れる。衝撃波とは単なる流体の不連続面ではなく、反射・捕捉・加熱といった微視的な粒子挙動によって支えられているのである。

衝撃波とは何か：跳び、マッハ数、ランキン・ユゴニオ条件

衝撃波の本質は「突き上げ機構(overtaking)」にある。圧力擾乱の伝播速度は媒質の圧縮が強いほど大きい(あるいは流れ自体が超音速であれば背景に追いつく)ため、なめらかな圧縮波の波峰が後続部分に追い上げられ、次第に波面が急峻化し、最終的に有限幅の跳びとして安定する。波面の急峻化を抑えるのは散逸、すなわち衝突による粘性・熱伝導、あるいはプラズマでは集団的な波動散逸であり、跳びの幅と内部構造はこの散逸過程で決まる。

跳びの上流(upstream)と下流(downstream)の状態は、衝撃波面に固定した系で質量・運動量・エネルギーの流束をそれぞれ保存させることで結び付けられる。これがランキン・ユゴニオ条件であり、衝撃波の強さは入射流速の音速に対する比、すなわちマッハ数で特徴付けられる。マッハ数が大きいほど圧縮比・温度跳びは大きくなり、強い衝撃波では単原子気体の極限圧縮比 4(プラズマでは電子の自由度を含めて異なる値をとる)に漸近する。衝撃波はエントロピーを生成する散逸構造であり、したがって可逆的な音波とは根本的に異なる。

ただし、この古典的な描像が成立するには、上流と下流のそれぞれが熱力学的平衡状態として定義できることが前提である。衝突が支配的な媒体ではこれは満たされるが、後述のように衝突レスプラズマでは下流状態がマクスウェル分布から大きく外れ得るため、跳びの両側を熱力学量だけで特徴付ける記述は破綻する。生成過程の初期に現れる幾何学的構造(いわゆる「ムッセルシェル」構造)のように、衝突系と衝突レス系で共通する形成的特徴も存在するが、最終状態の普遍性という点では両者は著しく異なる。

衝突性衝撃波と衝突レス衝撃波

恒星内部や高密度のレーザー生成プラズマでは、粒子間衝突の平均自由行程が衝撃波幅より短く、粘性と熱伝導が散逸を担う古典的な衝突性衝撃波が実現する。一方、宇宙空間の大半のように希薄な媒質では、クーロン衝突の平均自由行程が系の寸法そのものに匹敵する。例として西暦1006年の超新星残骸では、下流の温度およそ 10⁷ K、密度 0.1 cm⁻³ という典型的な値から、イオンの衝突平均自由行程は 10¹⁹ cm と見積もられ、残骸の半径に相当し、衝撃波幅(後に述べるようにイオン慣性長程度、~10⁸ cm)より圧倒的に大きい。このような条件下では衝突による粘性散逸は機能せず、衝撃波の急峻化と散逸を支えるのは電磁的な変動やプラズマ波といった集団効果である。これが衝突レス衝撃波である。

衝突レス衝撃波は、散逸を担う媒体によって大きく二つに分けられる。一つは電磁衝撃波で、主に宇宙環境で重要であり、ウェーベル(フィラメンテーション)不安定性が自己生成磁場を介して散逸を供給する。もう一つは静電衝撃波で、強い静電ポテンシャル障壁が散逸と反射を担い、電子がポテンシャルに捕捉され、イオンは障壁で反射されることでエネルギーを失う。冒頭の図が示した構造、すなわち上流の反射イオンと下流のトラップ渦は、まさにこの静電衝撃波の典型的な内部構造である。

衝突レス衝撃波で特筆すべきは、跳びの両側の状態がランキン・ユゴニオ条件だけでは決まらないことである。下流状態は「衝突レスプラズマ不安定性に対して安定」という広い制約の中から選ばれ、必ずしもマクスウェル分布とは限らない。系の全体進化を追わなければ最終状態が特定できないため、衝突性衝撃波と比べて普遍性が失われ、場合によっては全体流と衝撃波遷移部のスケール分離そのものが崩れる。二つの相貫流プラズマが衝突して衝撃波を形成する過程では、形成途中の下流・上流域でイオン音波不安定性が励起され、透過イオンと反射イオンのそれぞれを熱化させる。ただし流入速度が十分大きいと、重なり層内で斜めイオン音波不安定性が強く発達し、かえって衝撃波の形成を阻害し得る。このように衝突レス衝撃波の形成と散逸は、運動論的不安定性の非線形発展と密接に結び付いている。

静電衝撃波・イオン音響衝撃波の構造と運動論的効果

電子・イオンプラズマを低周波で見ると、電子はボルツマン分布に従って速く応答し、イオンの慣性と電子の圧力が音波(イオン音波)を支える。その非線形発展として、圧縮の跳びであるイオン音響衝撃波が現れる。散逸が弱ければ跳びは孤立波(ソリトン)の列へ分解され、散逸が十分に強ければ単一の単調な衝撃波となる。散逸の源泉は、荷電粒子間の衝突のほか、衝突レス環境では電子のランダウ減衰やイオン捕捉(位相空間の渦)などの運動論的効果であり、理論上はこの辺りで見られるように、超熱的(コッパ分布)電子や二温度イオン、ダスト粒子、電子・陽電子対の存在など、媒質の組成を変えるだけで衝撃波の構造や分散性・散逸性のバランスが大きく変わる。

運動論的な内部構造が決定的な意味を持つのは、跳びの幅と反射率である。跳びは数デバイ長あるいは数イオン慣性長の幅を持ち、その内部の電位障壁の高さが入射イオンの運動エネルギーを上回ればイオンは反射される。反射イオンは上流にビームとして残り、背景イオンとの間でイオン音波不安定性を駆動し、その非線形発展が上流の加熱と衝撃波面の緩和を進める。冒頭の図の位相空間に見える捕捉渦と反射ビームは、まさにこの自己組織化した運動論的構造である。したがって静電衝撃波の下流状態は「どれだけのイオンが反射し、どう熱化したか」という粒子分布の問題であり、流体記述では押さえきれない。

実験室での衝撃波：レーザー駆動衝撃波と慣性核融合

強レーザーを固体標的やプラズマに照射すると、焼灼圧力(アブレーション圧力)が超音速の圧縮波を駆動し、数マイクロメートル幅の衝撃波が標的内を進行する。現代の高強度レーザーでは、エネルギー密度が超新星爆発周辺の衝撃波に匹敵する速度の衝撃波を生成でき、しかも形成過程を高時空間分解能で計測できるため、宇宙衝撃波の理解に貢献する実験室天体物理学の舞台となっている。ただし注意が必要なのは、レーザー生成の非磁化衝撃波には集団的な静電力が介在し粒子衝突が無視できる場合があり、一方で大規模な宇宙衝撃波を記述する流体モデルは衝突による局所熱平衡を仮定しているため、レーザー生成衝撃波が常に宇宙衝撃波の相似モデルになるわけではない点である。それでも、衝突系と衝突レス系の双方に共通して働く物理過程――例えば衝撃波面の屈曲を増幅する非線形薄膜不安定性――は実験室で分離・検証でき、相似性が成り立つ範囲の理解を深めている。

慣性核融合(ICF)では、衝撃波の制御は圧縮の成否を直接左右する。間接駆動ではホルラーム内のX線がアブレーターを加熱して複数段の衝撃波を駆動し、その到達時刻の同期が燃料の圧縮効率を決める。衝撃波速度と放射温度の間のスケーリング則を逆に使えば、アブレーター内部を貫く衝撃波の計測からホルラーム内のピーク放射温度を推定でき、ソフトX線分光計やX線検出器の視野・観測方向による系統誤差を避ける診断手法として発展している。

ホルラーム内部では、壁から膨張する高Zプラズマと充填ガスや標的コロナの相互作用が衝撃波を駆動する。ホルラームは軸対称であるため、複数の衝撃波は軸上で相互作用する。二つの衝撃波が対向して衝突すると、両者の電磁場は互いに通り抜け、衝突後には高密度のプラトーが形成される。この過程は放射流体シミュレーションで通常は流体として扱われるが、低密度低Z領域では二峰性のイオン分布が観測されるなど運動論的効果が顕在化しており、実験とシミュレーションの不一致の原因の一つに挙げられている。

図は、対向する二つの静電衝撃波の相互作用を粒子シミュレーションで追ったものである。上段と中段は相互作用の前・最中・後の電場とイオン位相空間で、衝突時に電場のピークが一時的に合体して 1+2′ のような構造が現れ、その後ふたたび分離していく様子が分かる。下段は密度と温度のプロファイルで、相互作用後に密度プラトーが形成されることが示されている。

さらに、ホルラーム充填ガスの密度は真空から 1.6 mg cc⁻¹ まで変わり得るため、衝撃波を常に衝突レスと見なすのは適切でない。衝突は弱くても長時間スケールでは累積的に衝撃波構造を変化させ、イオンの反射率にも影響する。イオンを粒子法とモンテカルロ衝突モデルで、電子を質量なし流体で扱うハイブリッド手法は、この弱衝突性・長時間の領域を扱う有力な枠組みとなっている。衝撃波相互作用の理解は、ホルラーム内のプラズマ状態評価を通じて核融合設計に直結する問題なのである。

宇宙・天体物理における衝撃波：粒子加速と観測

宇宙では衝撃波は遍在する。太陽風中の「局所的」な衝撃波(惑星弓状衝撃波、コロナ質量放出の先行衝撃波)から、銀河団スケールの宇宙論的衝撃波まで、マッハ数と長さスケールは多種多様である。中でも超新星残骸は最も強力な衝突レス衝撃波の源であり、Ia型超新星で 10⁵¹ エルグ、コア崩壊型では 10⁵³ エルグ のエネルギー放出に伴い、最大数万 km s⁻¹ に達する放出物が星間物質を押し拡げる。周囲の物質を圧縮・加熱する順衝撃波と、放出物自身を減速・加熱する逆衝撃波の二重構造が基本像である。

衝突レス衝撃波の観測的特徴は、X線で見れば、衝撃波面直後の温度(典型例で 10⁷ K 程度)まで加熱された希薄プラズマの熱放射と、衝撃波面に局在する非熱的なシンクロトロン放射の共存である。興味深いのは、加熱と圧縮は衝撃波面という極めて薄い領域(イオン慣性長や反射イオンのジャイロ半径程度、例えば 10⁸ cm のオーダー)で起こるのに対し、その後の平衡への緩和は桁違いに遅いことである。電離平衡への緩和時間は約 10¹²/n 秒(n は下流密度)で、典型的な残骸では残骸の年齢を超えるため、X線スペクトルは一般に「過電離でない」プラズマの署名を示す。電子と陽子の温度平衡に至る時間も同様に長く、衝撃波直後のプラズマは種ごとに異なる温度を保ったまま緩やかに緩和する。この非平衡性そのものが、衝撃波面での加熱が衝突ではなく集団過程で起きている証拠である。

衝撃波面はまた、宇宙で最も効率的な粒子加速器である。荷電粒子は衝撃波面を前後するたびに収束する磁場乱流中で散乱され、往復のたびに平均してエネルギーを得る。これが一次フェルミ加速、すなわち拡散型衝撃波加速であり、超新星残骸の殻では GeV 電子からの電波シンクロトロン放射が系統的に観測され、若い残骸の衝撃波面では 10¹³ eV に達する電子からのX線シンクロトロン放射も観測されている。少数の事例では高エネルギー陽子の証拠も得られており、宇宙線の起源としての衝撃波加速の役割を裏付けている。さらに、加速粒子が占めるエネルギーが大きくなると、加速自身が衝撃波構造に逆作用し、圧縮や温度跳びを修正する。この効率的加速による衝撃波の変調は、観測された放射の空間分布から読み取ることができる。

図は超新星残骸 SN 1006 の合成画像である。左の全景(寸法 3×10¹⁹ cm)では、ソフトX線(赤)が衝撃加熱された放出物のリップル状の構造を、硬X線(青)が北東・南西の縁におけるシンクロトロン放射を示す。右の拡大図は、北西の衝撃波面を鋭く画定するバルマー Hα フィラメント(幅 3×10¹⁶ cm)であり、宇宙の衝撃波面が実際に観測された姿である。

なお、衝撃波現象は完全電離プラズマに限らない。ダストを含むプラズマでの塵音響衝撃波、量子プラズマやスピン効果を含む理論的展開、さらには変圧器油中のパルス放電のような非プラズマ系まで、衝撃波の概念は流体の圧縮性と散逸が共存するあらゆる場に拡張されている。宇宙から実験室、工業応用に至るこの普遍性こそ、衝撃波物理がプラズマ科学の横断的中心的主題であり続ける理由である。