[論文レビュー] Spectral methods in general relativity -- toward the simulation of 3D-gravitational collapse of neutron stars

S. Bonazzola, J. Frieben+2 more|ArXiv.org|Apr 12, 1996

Pulsars and Gravitational Waves Research参考文献 30被引用数 11

ひとこと要約

本稿では、一般相対性理論におけるアインシュタイン方程式を解くためのスペクトル法を、一般化された球対称に似た座標系で提示している。この手法により、中性子星の崩壊を極めて高精度にシミュレートすることが可能になる。このアプローチは、有限差分法と比較して数個の桁の高い精度を達成しており、球対称な崩壊および磁化され回転する中性子星への応用に成功しており、テンソルに基づくスキームを用いた3次元非軸対称シミュレーションの基盤を築いている。

ABSTRACT

Several applications of spectral methods to problems related to the relativistic astrophysics of compact objects are presented. Based on a proper definition of the analytical properties of regular tensorial functions we have developed a spectral method in a general sphericallike coordinate system. The applications include the investigation of spherically symmetric neutron star collapse as well as the solution of the coupled 2D-Einstein-Maxwell equations for magnetized, rapidly rotating neutron stars. In both cases the resulting codes are efficient and give results typically several orders of magnitude more accurate than equivalent codes based on finite difference schemes. We further report the current status of a 3D-code aiming at the simulation of non-axisymmetric neutron star collapse where we have chosen a tensor based numerical scheme.

研究の動機と目的

• コンパクトな相対論的対象物に適用可能な、一般相対性理論におけるアインシュタイン方程式を解く高精度なスペクトル法の開発を目的とする。
• 有限差分法よりも精度を向上させた、3次元重力的崩壊を伴う中性子星のシミュレーションにおける数値的課題に対処することを目的とする。
• スペクトル技法を、磁場および回転効果を含む非軸対称な3次元中性子星配置に拡張することを目的とする。
• 2次元および3次元における結合されたアインシュタイン＝マクスウェル方程式を解くための、テンソル的スペクトル法に基づく堅牢な数値フレームワークを構築することを目的とする。
• 中性子星の崩壊や重力波放射といった複雑な天体物理学的現象をシミュレートするのに適した、計算効率と精度を達成することを目的とする。

提案手法

• 本手法は、一般化された球対称に似た座標系においてスペクトルコロケーション技法を用いて、アインシュタイン場方程式を離散化する。
• 曲がった時空における正則なテンソル関数の厳密な数学的取り扱いに基づいており、安定性と収束性を保証する。
• スペクトル的手法は、空間的微分を高次精度で表現するためのグローバル基底関数（例えば、チェビシェフ多項式やルジャンドル多項式）を用いる。
• 非軸対称な崩壊のための3次元方程式系を処理するため、テンソルに基づく数値スキームが実装されている。
• 磁化され高速に回転する中性子星の2次元アインシュタイン＝マクスウェル方程式を、角度方向および半径方向のスペクトル分解を用いて解いている。
• コードは球対称なケースで検証され、適応的グリッド解像度を用いた擬似スペクトル的手法を用いて3次元に拡張されている。

実験結果

リサーチクエスチョン

• RQ1 スペクトル法は、相対論的中性子星崩壊のシミュレーションにおいて、有限差分スキームと比較して著しく高い精度を達成できるか？
• RQ2 スペクトル法は、中性子星に関連する非自明な球対称に似た座標系におけるテンソル場を扱うために一般化可能か？
• RQ3 スペクトル法を3次元非軸対称な重力的崩壊に拡張する可能性と性能は何か？
• RQ4 スペクトル法は、磁化され回転する中性子星の結合されたアインシュタイン＝マクスウェル方程式をどの程度の精度で解けるか？
• RQ5 強い重力場および磁場を伴う複雑で非線形な系、例えば3次元中性子星崩壊において、スペクトル法は安定性と収束性を維持できるか？

主な発見

• 球対称な中性子星崩壊のシミュレーションにおいて、スペクトル法は同等の有限差分コードと比較して、数個の桁の高い精度を達成している。
• 磁化され高速に回転する中性子星の2次元アインシュタイン＝マクスウェル系は、高い精度と安定性で正しく解かれている。
• 3次元コードフレームワークは開発中であり、テンソルに基づくスペクトル的手法を用いた非軸対称崩壊のシミュレーションに有望である。
• 一般化された球対称に似た座標系におけるスペクトル法の使用により、曲がった時空におけるテンソル場の正確な表現が可能になった。
• スペクトル法は、有限差分法と比較して、優れた収束速度と低い数値的拡散を示している。
• 滑らかな解を伴う問題、例えば平衡状態および準平衡状態の中性子星構成に対して、本手法は計算的に効率的である。