[論文レビュー] Temporal evolution of magnetic molecular shocks II. Analytics of the steady state and semi-analytical construction of intermediate ages

P. Lesaffre, J.-P. Chièze+2 more|ArXiv.org|Jul 28, 2004

Ionosphere and magnetosphere dynamics参考文献 9被引用数 21

ひとこと要約

本稿では、準定常状態近似を活用することで、時間に依存する磁気分子衝撃を定常状態コードのみでモデル化する半アナリティカル手法を開発する。磁気衝撃における密度、圧力、圧縮比についての解析的関係を導出し、時間に依存するシミュレーションと比較して検証した。その結果、安定な衝撃は実質的に切断された定常状態であることが示された。主な貢献は、低速衝撃（u < 20 km s⁻¹）のための中間年齢衝撃において、完全な時間に依存するシミュレーションを用いずにCJ型衝撃の挙動を正確に特定できる予測アルゴリズムの開発である。

ABSTRACT

In the first paper of this series (Paper I) we computed time dependent simulations of multifluid shocks with chemistry and a transverse magnetic field frozen in the ions, using an adaptive moving grid. In this paper, we present new analytical results on steady-state molecular shocks. Relationships between density and pressure in the neutral fluid are derived for the cold magnetic precursor, hot magnetic precursor, adiabatic shock front, and the following cooling layer. The compression ratio and temperature behind a fully dissociative adiabatic shock is also derived. To prove that these results may even hold for intermediate ages, we design a test to locally characterise the validity of the steady state equations in a time-dependent shock simulation. Applying this tool to the results of Paper I, we show that most of these shocks (all the stable ones) are indeed in a quasi-steady state at all times, i.e. : a given snapshot is composed of one or more truncated steady shock. Finally, we use this property to produce a construction method of any intermediate time of low velocity shocks (u &lt; 20 km/s) with only a steady-state code. In particular, this method allows one to predict the occurrence of steady CJ-type shocks more accurately than previously proposed criteria.

研究の動機と目的

• 横磁場を有する定常状態磁気分子衝撃における密度、圧力、圧縮比に関する解析的関係を導出すること。
• 時間に依存する磁気衝撃シミュレーションにおいて、定常状態近似の有効性を検証すること、特に中間年齢の状況において。
• 低速衝撃（u < 20 km s⁻¹）の時間に依存する衝撃構造を再構築するために、定常状態コードのみを用いる実用的手法を開発すること。
• 音速点形成に基づく既存の基準を再考することで、磁化された流れにおけるCJ型衝撃発生の予測をより信頼性高くすること。
• 定常状態モデルと時間に依存するモデルの間のギャップを埋め、観測された分子衝撃のより迅速かつ正確な診断を可能にすること。

提案手法

• 定常状態磁気衝撃における中性およびイオン化流体の力学的解析解を導出し、密度、圧力、磁気圧縮比との関係を含む。
• 時間に依存するシミュレーション内での定常状態方程式の有効性を評価するための局所的診断手法を導入し、時間微分項とフラックス勾配を比較する。
• Paper I の時間に依存する多流体シミュレーション結果にこの診断を適用し、安定な衝撃がその進化全般にわたり準定常状態に保たれていることを確認した。
• 定常状態モデルの速度および密度プロファイルを用いて、動的に整合性のある半径で定常状態解を切断することで、中間時間の衝撃プロファイルを構築した。
• 2段階のアルゴリズムを実装した：(1) 全定常状態プロファイルを計算し、(2) 時間に依存するシミュレーションにおける特定の年齢における局所的流れ条件と一致する半径で切断する。
• 化学ネットワークの柔軟性を確保するため、非動的支配種については事後処理により対応を拡張した。

実験結果

リサーチクエスチョン

• RQ1 時間に依存する磁気分子衝撃がその進化過程でどの程度準定常状態を維持するか。
• RQ2 磁場の存在下でも、中間年齢衝撃のダイナミクスを解析的定常状態関係がどれほど正確に記述できるか。
• RQ3 既存の基準よりも信頼性の高い方法で、磁化された流れにおけるCJ型衝撃の発生を予測するにはどうすればよいか。
• RQ4 準定常状態仮定のもとで、時間に依存する衝撃構造を定常状態コードのみで再構築できるか。
• RQ5 急速に反応する化学反応や強い不安定性を示す系では、準定常状態近似の限界は何か。

主な発見

• シミュレーションにおける安定な磁気分子衝撃は、すべての時点で準定常状態にあり、各時間ステップのスナップショットは1つ以上の切断された定常状態衝撃に対応している。
• 冷たい磁気前駆体、高温の磁気前駆体、断熱衝撃波面、冷却層における密度および圧力に関する解析的関係が導出され、検証された。
• 完全に解離する断熱衝撃の背後における圧縮比が解析的に導出され、衝撃モデルのベンチマークとしての役割を果たす。
• CJ型衝撃形成の基準は見直すべきである：磁気圧縮の増強およびピストン付近での速度再結合の影響により、従来の予想よりも低い速度でもCJ型衝撃が発生しうる。
• 本手法は、低速衝撃（u < 20 km s⁻¹）のための中間年齢衝撃構造を、定常状態コードのみを用いて成功裏に再構築でき、より迅速な診断が可能になった。
• 急激に変化する冷却剤や顕著な拡散効果を示す系、および強い不安定性を示す衝撃（例：弱く解離する衝撃）では、準定常状態仮定が失敗する。