[論文レビュー] New numerical framework for the generalized Baumgarte-Shapiro-Shibata-Nakamura formulation: The vacuum case for spherical symmetry

M. A. Alcoforado, R. F. Aranha+2 more|arXiv (Cornell University)|May 19, 2021

Pulsars and Gravitational Waves Research参考文献 62被引用数 7

ひとこと要約

この論文では、球対称性における一般化BSSN形式を解くために、ガラーキン・コロケーション法に基づく新規のオープンソーススペクトルコードRIOを提示する。原点および空間無限遠における正則性を内蔵した基底関数を用いることで、明示的な正則化を施さずに制約の指数的収束を達成し、真空中の球対称テストにおいてハミルトニアンおよび運動量制約で準マシン精度に近い精度を示した。

ABSTRACT

Here we report a developed high performance and simplified version of the code denominated RIO, which can be easily extended, for the generalized BSSN formulation. We implement a code which is regular at the center of symmetry, without use a special procedure for regularization, as usual. We get exponential convergence for constraints. The numerical algorithm is based on the Galerkin-Collocation method developed successfully for diverse physical scenarios by the Numerical Relativity Group at UERJ. For the sake of clarity in presentation, we consider here the most simple case to display the most salient features of the procedure. Thus, we focus on the definite tests of the new numerical framework. The timing and performance of the code show that we can reach a better accuracy close to the machine precision, for the Hamiltonian and momentum constraints. RIO will be an open source code; currently is under continuous development to consider more general and realistic problems.

研究の動機と目的

• 球対称性における一般化BSSN形式のための高性能で拡張可能な数値フレームワークの開発。
• 原点および極軸における便宜的な正則化手順の必要性を排除するため、正則性をスペクトル基底関数に埋め込むこと。
• ガラーキン・コロケーションスペクトル法を用いて制約違反における指数的収束を達成すること。
• 将来の2次元および3次元問題への拡張、例えば重力崩壊やブリル波衝突のための基盤を築くこと。
• UERJにおける数値相対論研究のための統一的かつ保守可能なコードベース（RIO）を確立すること。

提案手法

• ガラーキン・コロケーション（GC）スペクトル法が用いられ、ガラーキン射影とガウス＝レジェンドル＝ロバッチョフ点におけるコロケーションを統合する。
• 基底関数はL²[0, ∞)ヒルバート空間から選択され、r=0および空間無限遠における内在的正則性を保証する。
• 偏微分方程式系は、これらの基底関数におけるスペクトル展開を経て、自己同様の力学系に変換される。
• 時間発展は、スペクトル係数に適用される固定ステップの4次ルンゲ＝クッタ積分法によって実行される。
• コロケーション点は特異性を避けるように配置され、境界条件は基底関数の選択によって自然に扱われる。
• 既知の初期データおよび制約を有する真空中の球対称ケースに対して、手法の妥当性が検証された。

実験結果

リサーチクエスチョン

• RQ1 スペクトルガラーキン・コロケーション法は、球対称性における一般化BSSN形式のハミルトニアンおよび運動量制約において、指数的収束を達成できるか？
• RQ2 基底関数の選択によって、原点および空間無限遠における正則性を内蔵的に実現でき、特別な正則化手順の必要性を排除できるか？
• RQ3 このフレームワークの性能および精度は、既存のBSSN実装と比較して、制約違反および収束率の観点でどのように異なるか？
• RQ4 このフレームワークは、重力崩壊や円柱座標系における2次元問題といったより複雑なシナリオに拡張可能か？
• RQ5 RIOコードは真空中の球対称テストにおいて、達成可能な精度と計算効率はどの程度か？

主な発見

• RIOコードはハミルトニアンおよび運動量制約において指数的収束を達成し、スペクトル精度を示した。
• 制約違反は準マシン精度にまで低減され、高い数値的精度を示した。
• 基底関数の選択によって、原点および空間無限遠における正則性が自然に強制され、特別な正則化技術の必要性が排除された。
• コードは安定しており、固定ステップの4次ルンゲ＝クッタ積分法を用いた時間発展において良好に動作した。
• フレームワークは拡張可能であり、現在も重力崩壊やブリル波衝突を含む2次元および3次元問題への将来的な応用に向けて積極的に開発が進められている。
• RIOコードはオープンソースであり、数値相対論研究における継続的開発および再利用を目的として設計されている。