[論文レビュー] A 3+1 Computational Scheme for Dynamic Spherically Symmetric Black Hole Spacetimes -- II: Time Evolution

Jonathan Thornburg|ArXiv.org|Jun 6, 1999

Black Holes and Theoretical Physics参考文献 86被引用数 7

ひとこと要約

本稿では、ブラックホール除去技術を用いてスカラー場と結合した球対称ブラックホール時空を安定かつ高精度に進化させる3+1数値スキームを提示する。自由進化を採用し、Eddington-Finkelsteinに類似した座標系を用い、空間および時間方向に4次精度の有限差分法を適用し、ルンゲ＝クッタ積分法を用いることで、t = 4000mを超える長時間にわたり安定した進化が達成され、計量成分および制約条件の誤差が10^-5未満となる。

ABSTRACT

This is the second in a series of papers describing a 3+1 computational scheme for the numerical simulation of dynamic black hole spacetimes. We discuss the numerical time-evolution of a given black-hole-containing initial data slice in spherical symmetry. We avoid singularities via the "black-hole exclusion" or "horizon boundary condition" technique, where the slices meet the black hole's singularity, but on each slice a spatial neighbourhood of the singularity is excluded from the domain of the numerical computations. After first discussing some of the key design choices which arise with the black hole exclusion technique, we then give a detailed description of our numerical evolution scheme for spherically symmetric scalar field evolution, assuming that a black hole is already present on the initial slice. We use a free evolution, with Eddington-Finkelstein-like coordinates and the inner boundary placed at a fixed coordinate radius well inside the horizon. Our numerical scheme is based on the method of lines (MOL), where spacetime PDEs are first finite differenced in space only, yielding a system of coupled ODEs for the time evolution of the field variables along the spatial-grid-point world lines. These ODEs are then time-integrated by standard methods. We use 4th order finite differencing in both space and time, with 5 and/or 6 point spatial molecules (off-centered near the grid boundaries), and a Runge-Kutta time integrator. The spatial grid is smoothly nonuniform, but not adaptive. We present numerical black hole + scalar field evolutions showing that this scheme is stable, can evolve "forever" (we have gone to t > 4000m), and is very accurate. At a resolution Delta_r/r = 3% near the horizon, typical errors in g_ij(K_ij) at t=100m are <= 1e-5(3e-7), and the energy constraint is < 3e-5.

研究の動機と目的

• 球対称性を有するスカラー場を伴う動的ブラックホール時空をシミュレートする、強固な3+1数値フレームワークの構築を目的とする。
• 計算領域からブラックホール特異点の近傍領域を除外することで、数値的特異点を回避する。
• 高次精度有限差分法と標準的なODE積分法を用いて、時間発展における高精度と長期間の安定性を達成する。
• 長期間にわたり収束性と制約保存性を示す数値的進化を用いて、スキームの妥当性を検証する。
• 制御可能な数値誤差を伴う球対称性下でのブラックホール形成および進化のシミュレーションの基盤を提供する。

提案手法

• 特異点近傍の空間的領域を計算領域から除外することで、ブラックホール除去技術を適用し、座標特異点を回避する。
• ホライズン付近での正則性を保証し、安定な進化を可能にするために、Eddington-Finkelsteinに類似した座標系を用いる。
• 時間方向に偏微分方程式（PDE）を離散化せず、空間方向にのみ離散化する「線法（MOL）」を採用し、世界線に沿った場の進化を記述する連立常微分方程式（ODE）系を導出する。
• 空間方向に5点または6点のステンシルを用い、4次精度の有限差分法を適用する。境界付近ではオフセンター化を施し、精度を維持する。
• 時間方向に4次精度のルンゲ＝クッタ積分法を用いて、ODE系を時間発展させる。
• 空間グリッドは滑らかに非一様であるが、適応的でない。ホライズン付近での解像度はΔr/r ≈ 3%である。

実験結果

リサーチクエスチョン

• RQ1 3+1数値スキームは、球対称性下でスカラー場と結合したブラックホール時空を、安定的かつ長期間にわたり進化させることができるか？
• RQ2 ブラックホール除去技術は、座標特異点を回避しつつ、数値的精度を維持するのにどの程度有効か？
• RQ3 この文脈において、4次精度有限差分法とルンゲ＝クッタ時間積分法を用いることで、どの程度の精度が達成可能か？
• RQ4 エネルギー制約のような制約条件は、長期間の進化においてどの程度保存されるか？
• RQ5 数値的発散や精度の低下が生じないよう、任意に遅い時刻までシステムを進化させることができるか？

主な発見

• t = 4000mを超える時間にわたり安定した進化が達成され、長期間にわたる数値的安定性が示された。
• ホライズン付近でΔr/r ≈ 3%の解像度を有する条件下で、t = 100mにおける計量成分g_ijの誤差は≤ 1×10^-5、外的曲率K_ijの誤差は≤ 3×10^-7に抑えられた。
• 同じ時間間隔にわたり、エネルギー制約の保存精度が3×10^-5未満であることが確認され、強い制約満たしの成立が示された。
• 4次精度の空間的および時間的離散化を用いることで、高い精度と収束性が確保された。
• 内境界をホライズンより内側の座標半径に固定した状態でも、スキームは安定的かつ正確に保たれた。
• 結果から、ブラックホール除去技術とEddington-Finkelsteinに類似した座標系の組み合わせが、座標特異点を回避する信頼性の高いブラックホール時空シミュレーションを可能にすることが確認された。