[論文レビュー] A Multiple-Grid-Patch Evolution Scheme for 3-D Black Hole Excision

Jonathan Thornburg|ArXiv.org|Dec 4, 2000

Pulsars and Gravitational Waves Research被引用数 4

ひとこと要約

本稿では、3次元ブラックホール除去における標準的なxyzグリッドの代わりに、r×(角度座標)トポロジーを用いた複数グリッドパッチ方式を提案する。これにより、滑らかで階段状でない除去境界が得られ、数値的不安定性が低減される。4次精度の有限差分法とADM方程式を用い、コードはカー初期データを約60Mまで時間発展させ、パッチの隅では制約誤差が約1%に達する。安定性を向上させるためにBSSN形式の実装が進行中である。

ABSTRACT

When using black hole excision to numerically evolve a fully generic black hole spacetime, most 3-D 3+1 codes use an $xyz$-topology (spatial) grid. In such a grid, an $r = \constant$ excision surface must be approximated by an irregular and non-smooth "staircase-shaped" excision grid boundary, which may introduce numerical instabilities into the evolution. In this paper I describe an alternate scheme, which uses multiple grid patches, each with topology $\{r imes ({ m angular coordinates})\}$, to cover the slice outside the $r = \constant$ excision surface. The excision grid boundary is now smooth, so the evolution should be less prone to instabilities. With 4th order finite differencing, this code evolves Kerr initial data to ${\sim} 60M$ using the ADM equations; I'm currently implementing the BSSN equations in it in the hope that this will improve the stability.

研究の動機と目的

• 標準的な3次元$3+1$有限差分コードにおける不規則で階段状の除去境界が引き起こす数値的不安定性を是正すること。
• 一般のブラックホール時空を安定かつ高精度に時間発展させるためのマルチパッチグリッドシステムを用いたスケーラブルなスキームの開発。
• 標準的なxyzグリッドトポロジーを、r×(角度座標)パッチに置き換えることで滑らかな除去表面を実現すること。
• パッチ境界の重なり領域において内部データのみを用いた1次元補間を可能にすることで、数値的安定性を向上させること。
• 新規スキームにBSSN形式を実装・テストし、現在のADM形式に比べて長期的な安定性を向上させること。

提案手法

• 本スキームは、$r = \text{constant}$の除去表面の外側を覆う複数のグリッドパッチを用いる。各パッチはトポロジー$ r \times (\text{角度座標}) $を持つ。
• 各パッチは角度方向のゴーストゾーンを有し、これらの変数は隣接するパッチからの1次元テンソル補間によって設定される。
• 隣接するパッチは$r$および1つの角度座標を共有しており、補間の一貫性とグリッドの連続性が保証される。
• グリッドは径方向に滑らかに非一様であり、空間方向には4次精度の有限差分法が適用される。
• 時間積分には、先行する球対称実装と整合する4次精度ルンゲ・クッタ法の方法が用いられる。
• パッチを複数のグリッド間隔分重ねることで、補間にパッチ内部のデータのみが使用されるようにし、安定性が向上する。

実験結果

リサーチクエスチョン

• RQ1 r×(角度座標)トポロジーを持つマルチパッチグリッドシステムは、標準的なxyzトポロジーのグリッドと比較して、3次元ブラックホール除去における数値的不安定性を低減できるか？
• RQ2 重なり領域における1次元補間の使用が、有限差分時間発展の安定性および精度に与える影響は何か？
• RQ3 ADM形式は、複数のパッチが接する部分において、この新規スキームでどの程度の性能を示すか？
• RQ4 BSSN形式は、ADM形式と比較して、このマルチパッチフレームワークにおいて安定性をさらに向上させられるか？
• RQ5 滑らかな除去境界は、階段状の境界に起因する不安定性をどの程度解消できるか？

主な発見

• 複数グリッドパッチ方式により、標準的なxyzグリッドコードが示す不規則で階段状の除去境界が、滑らかで$r = \text{constant}$の境界に置き換えられた。
• 4次精度の有限差分法とADM方程式を用い、コードは座標時間で約$60M$までカー初期データを時間発展させた。
• 3つのパッチがすべて集まる隅では、有限差分法の不安定性が顕在化し、$t \sim 60M$で制約違反が約1%に達する。
• パッチ接合部での不安定性は、このマルチパッチ構成下での現在のADM形式の限界を示唆している。
• 著者は、長期的な数値的安定性の向上を調査するために、コード内にBSSN形式の実装を進行中である。
• 内部データのみを用いた重なり領域のゴーストゾーンを用いることで、補間が安定し、境界アーチファクトが回避される。