[論文レビュー] A new 3D transport and radiation code for galactic cosmic rays

M. W. Werner, R. Kissmann+2 more|arXiv (Cornell University)|Aug 13, 2013

Astrophysics and Cosmic Phenomena被引用数 5

ひとこと要約

本論文は、従来の2次元シミュレーションの制限を克服するため、完全な3次元、非等方的拡散および並列計算を可能にする新しい3次元宇宙線（CR）輸送・放射線コードを提示する。高度な数値スolverを内蔵し、収束制御機能を備えており、特に高エネルギー電子や非軸対称な源分布に対して、GALPROPより優れた精度と効率を示している。

ABSTRACT

We show the necessity for a new approach towards comprehensive and consistent simulations of the propagation of galactic cosmic rays. Our developments are optimised for addressing the spatially 3-dimensional inhomogeneous diffusion problem and utilise contemporary numerical methods. We aim to address the transport problem in a full 3-dimensional environment. For that, we test the transition from 2D to 3D simulation results within an existing propagation code. We present sub-kpc scale simulations that allow the investigation of small- scale structures regarding different model conditions such as variety regarding non-axisymmetric cosmic ray source distributions. These results are discussed critically and motivate our development of a new transport code for galactic cosmic rays. The capabilities of this code are outlined.

研究の動機と目的

• 現在の2次元宇宙線（CR）伝搬コードが、等方的拡散や軸対称な源分布といった簡略化仮定に依存しているという制限を解消する。
• GALPROPのような既存コードが収束制御や数値フィードバックを欠いていること、ユーザーが手動で収束テストを実施する必要があるという点を克服する。
• 空間的に変化する拡散テンソルや非軸対称な源分布（例：スパイラルアームのパターン）を伴う、高精度な3次元CR伝搬シミュレーションを可能にする。
• 一般的なオペレータ分割ではなく、物理学に特化したスolverを用いた最適化されたコードアーキテクチャを構築し、定常状態および時間依存解の両方を扱えるようにする。
• MPIを活用して現代の高性能コンputングを活用し、CRモデリング分野全体での一貫性を高める。

提案手法

• 宇宙線輸送方程式の各物理項に高度なスolverを適用し、従来のCrank-Nicolsonスキームに代わる、現代の数値法を用いた新しい3次元輸送コードを実装する。
• 次元別オペレータ分割を、対流、空間的拡散、運動量空間のエネルギー損失に特化したスolverに置き換え、精度と安定性を向上させる。
• スパース行列スolverを用いて、完全な3次元非等方的拡散に起因する非三重対角系を処理し、定常状態輸送方程式の効率的解法を実現する。
• 定常状態解を初期条件として用いる時間依存スイッチャーを統合し、一時的源や時間変化するCR分布のモデル化を可能にする。
• 収束モニタリング機能をコードに内蔵し、数値的安定性を確保するためにユーザーがtimestep_repeatなどのパラメータを手動で調整する必要がなくなる。
• メッセージパッシングインターフェース（MPI）標準を活用し、高性能コンputングクラスタ上でスケーラブルかつ並列計算を実現する。

実験結果

リサーチクエスチョン

• RQ1 2次元と3次元の宇宙線シミュレーションは、特に高エネルギー領域において、陽子および電子分布の予測にどの程度の差を示すか？
• RQ2 スパイラルアームのパターンを模倣する非軸対称な宇宙線源分布は、2次元モデルと比較して3次元シミュレーションにおけるCRの空間的形状にどのように影響を与えるか？
• RQ3 新しい数値スolverは、GALPROPのような既存コードに見られるユーザー依存の収束問題を解消できるか？
• RQ4 3次元CR輸送シミュレーションにおいて、一般化されたオペレータ分割と比較して、物理学に最適化されたスolverを用いることで得られる性能および精度の向上はどの程度か？
• RQ5 新しいコードは、現代の並列コンputングアーキテクチャ上で、サブkpcスケールのシミュレーションに対してどの程度効率的にスケーリングできるか？

主な発見

• 100 MeVの陽子では、GALPROPにおける2次元と3次元のシミュレーションが1％未満の偏差で一致しており、低エネルギー領域では一貫性が確認され、3次元実装の妥当性が裏付けられた。
• 1.2 TeVの電子では、2次元と3次元のシミュレーション間に顕著な差異（最大約10％）が生じ、2次元モデルが高エネルギー電子伝搬を正確に捉えていないことを示している。
• 新しいコードのスolverは、ユーザー定義の基準が不要で収束を達成しており、timestep_repeatなどのパラメータの手動チューニングの必要性が解消された。
• 時間積分を用いないことで、定常状態解が効率的に計算され、時間離散化に起因する数値誤差が低減され、精度が向上した。
• コードは完全な3次元非等方的拡散および非軸対称な源分布をサポートしており、構造的銀河環境におけるCR伝搬の新たな研究を可能にした。
• MPIの活用により、高性能コンputングシステム上でスケーラブルなパフォーマンスが実現され、サブkpcスケールの3次元シミュレーションが実現可能かつ効率的になった。