[論文レビュー] Incorporation of cosmic ray transport into the ZEUS MHD code. Application for studies of Parker instability in the ISM

M. Hanasz, H. Lesch|ArXiv.org|Sep 24, 2003

Solar and Space Plasma Dynamics参考文献 31被引用数 40

ひとこと要約

本稿では、宇宙線輸送を拡散・拡散方程式に組み込むための数値的手法を提示しており、平行および垂直方向の拡散係数を用いて磁場に沿ったおよび垂直方向の異方性拡散を可能にする。この手法は、現実的な銀河系拡散係数（例：$3 \times 10^{28}\ \text{cm}^2\text{s}^{-1}$）の範囲で数値的に安定であり、有限の拡散がパーカー不安定性の成長率を低下させることを示しており、線形近似に反するものであり、宇宙線駆動型銀河系空間力学における非線形効果の重要性を浮き彫りにしている。

ABSTRACT

We present a numerical algorithm for the incorporation of the active cosmic ray transport, into the ZEUS-3D magnetohydrodynamical code. The cosmic ray transport is described by the diffusion-advection equation. The applied form of the diffusion tensor allows for anisotropic diffusion of cosmic rays along and across the magnetic field direction, which is controlled by two parameters: the parallel and perpendicular diffusion coefficients. The implemented numerical algorithm is tested by comparison of the diffusive transport of cosmic rays to analytical solutions of the diffusion equation. Our method is numerically stable for a wide range of diffusion coefficients, including the realistic values inferred from the observational data for the Milky Way of about $ 6 imes 10^{28} \cm^2 \s^{-1}$. The presented algorithm is applied for exemplary simulations of the the Parker instability triggered by cosmic rays injected by a single SN remnant.

研究の動機と目的

• ZEUS-3D MHDコードに、活性的な宇宙線輸送を組み込むための数値的に安定な手法を開発すること。
• 別々の平行方向（$K_\parallel$）および垂直方向（$K_\perp$）の拡散係数を用いて、磁場に沿ったおよび垂直方向の異方性輸送をモデル化すること。
• アルゴリズムの妥当性を確認するため、拡散方程式の解析的解との比較による検証。
• 拡張されたコードを用いて、1つの超新星残骸から生じる宇宙線によって引き起こされるパーカー不安定性を研究すること。
• 有限の拡散係数が、パーカー不安定性の非線形的進化に与える影響を調査し、線形解析からの仮定に疑問を呈すること。

提案手法

• ステガードメッシュ上で、オペレータ分割と有限差分法を用いて、ZEUS-3D MHDコードを宇宙線拡散・拡散方程式を解くように変更する。
• 別々の平行方向（$K_\parallel$）および垂直方向（$K_\perp$）の拡散係数を有するテンソルを用いて異方性拡散を実装する。
• 宇宙線輸送シミュレーション中に磁場の発散がゼロとなる条件を保つため、制約付き伝搬形式を用いる。
• 磁場が存在する重力的層化した銀河系空間媒体（ISM）への、1つの超新星残骸からの宇宙線注入をシミュレートする。
• 幅広い$K_\parallel$値の範囲で、シミュレートされた宇宙線拡散を拡散方程式の解析的解と比較することで、アルゴリズムの妥当性を検証する。
• さまざまな$K_\parallel$値（現実的な値：$3 \times 10^{28}\ \text{cm}^2\text{s}^{-1}$ と低い値：$3 \times 10^{27}\ \text{cm}^2\text{s}^{-1}$）を用いた3次元パーカー不安定性のシミュレーションを実施する。

実験結果

リサーチクエスチョン

• RQ1 現実的な$K_\parallel \approx 3 \times 10^{28}\ \text{cm}^2\text{s}^{-1}$を有する有限の宇宙線拡散は、パーカー不安定性の成長率にどのように影響を与えるか？
• RQ2 宇宙線拡散が制限されている場合、パーカー不安定性の進化において非線形効果がどれほど支配的になるか？
• RQ3 異なる拡散係数を用いた場合、時間経過に伴い宇宙線圧力の空間的分布はどのように変化するか？
• RQ4 観測データから推定される値を含む広範な拡散係数範囲において、数値的アルゴリズムが宇宙線輸送を正確にシミュレートできるか？
• RQ5 有限の拡散を組み込むことで、パーカー不安定性のモード構造（例えば、交差モードと波状モードの相対的寄与）がどのように変化するか？

主な発見

• 数値的アルゴリズムは、現実的な銀河系の値（$3 \times 10^{28}\ \text{cm}^2\text{s}^{-1}$）を含む広範な拡散係数範囲で、高い数値的安定性を示して宇宙線輸送をシミュレートできた。
• 現実的な拡散係数（$K_\parallel = 3 \times 10^{28}\ \text{cm}^2\text{s}^{-1}$）の場合、宇宙線はより均一に拡散され、ピークの宇宙線圧力が低下し、低拡散ケースと比較して初期の浮力力が弱まった。
• 現実的な拡散係数の場合、成長率が低い拡散ケース（$K_\parallel = 3 \times 10^{27}\ \text{cm}^2\text{s}^{-1}$）よりも遅くパーカー不安定性が成長したが、これは線形解析からの予想に反する。
• 現実的な拡散係数の場合、時間経過とともに波状モードが優勢になる（約$t = 150$付近）、一方で初期には交差モードが支配的であり、非線形的発展が遅延していることが示された。
• 現実的な拡散係数の場合、$t = 100$における最大の垂直速度は低拡散ケースよりも低かったが、後に波状モードの出現に伴い不安定性が加速した。
• 結果から、低拡散係数の場合、線形近似の適用性に疑問が呈され、強い非線形効果（例：垂直速度が音速に近づく）が顕著に現れるため、完全な非線形シミュレーションが不可欠であることが示唆された。