[論文レビュー] A strongly hyperbolic first-order CCZ4 formulation of the Einstein equations and its solution with discontinuous Galerkin schemes

Michael Dumbser, Federico Maria Guercilena+3 more|arXiv (Cornell University)|Jul 31, 2017

Numerical methods for differential equations被引用数 8

ひとこと要約

本稿では、制約減衰を備えたCCZ4系から導かれた、強く双曲型の一次形式であるFO-CCZ4を提示する。これにより、高次離散ガラーキン（DG）法を用いた安定なシミュレーションが可能となる。本研究では、一次形式を用いて、初めての3次元高次DG法による移動穿孔を伴う二重ブラックホール系のシミュレーションを達成した。これは、多次元真空中のテストを通じて検証され、メッシュの適応的細分化とサブセルリミッターによって強化されている。

ABSTRACT

We present a strongly hyperbolic first-order formulation of the Einstein equations based on the conformal and covariant Z4 system (CCZ4) with constraint-violation damping, which we refer to as FO-CCZ4. As CCZ4, this formulation combines the advantages of a conformal and traceless formulation, with the suppression of constraint violations given by the damping terms, but being first order in time and space, it is particularly suited for a discontinuous Galerkin (DG) implementation. The strongly hyperbolic first-order formulation has been obtained by making careful use of first and second-order ordering constraints. A proof of strong hyperbolicity is given for a selected choice of gauges via an analytical computation of the entire eigenstructure of the FO-CCZ4 system. The resulting governing partial differential equations system is written in non-conservative form and requires the evolution of 58 unknowns. A key feature of our formulation is that the first-order CCZ4 system decouples into a set of pure ordinary differential equations and a reduced hyperbolic system of partial differential equations that contains only linearly degenerate fields. We implement FO-CCZ4 in a high-order path-conservative arbitrary-high-order-method-using-derivatives (ADER)-DG scheme with adaptive mesh refinement and local time-stepping, supplemented with a third-order ADER-WENO subcell finite-volume limiter in order to deal with singularities arising with black holes. We validate the correctness of the formulation through a series of standard tests in vacuum, performed in one, two and three spatial dimensions, and also present preliminary results on the evolution of binary black-hole systems. To the best of our knowledge, these are the first successful three-dimensional simulations of moving punctures carried out with high-order DG schemes using a first-order formulation of the Einstein equations.

研究の動機と目的

• 高次不連続ガラーキン（DG）法に適した一次形式で強く双曲型のアインシュタイン方程式の形式を構築すること。
• 共形的かつトレースレス形式の利点を、制約破壊の減衰と組み合わせることで、数値的安定性を向上させること。
• 特に二重ブラックホール系を対象とした、ブラックホール時空の高次で安定なDG法によるシミュレーションを可能とすること。
• ブラックホール進化における特異性といった数値的課題を、高度なリミッターと適応的メッシュ細分化を用いて解決すること。

提案手法

• 一次および二次の順序制約を慎重に用いて、CCZ4系の一次形式を導出することで、強い双曲型性を保証する。
• 58個の進化変数を有する非保存型双曲型PDE系を構築し、純粋な常微分方程式と線形的に退化した場を有する低次元双曲型系に分解する。
• 局所時間ステッピングと適応的メッシュ細分化を備えた、任意の高次法（ADER）-DGスキームにこの系を実装する。
• ブラックホールシミュレーションで生じる特異性に対処するため、三階ADDER-WENOサブセル有限体積リミッターを統合する。
• 一次元、二次元、三次元の標準的真空中テストを通じて、形式の妥当性を検証する。
• 二重ブラックホール系の進化にこのスキームを適用し、一次形式を用いた高次DG法による初めての3次元移動穿孔シミュレーションを達成する。

実験結果

リサーチクエスチョン

• RQ1 高次DG離散化に対応できる一次形式で強く双曲型なアインシュタイン方程式の形式を構築できるか？
• RQ2 FO-CCZ4形式は標準的なゲージ選択のもとでどのように強く双曲型性を保証するのか？また、その固有構造はいかなるものか？
• RQ3 適応的メッシュ細分化と局所時間ステッピングを併用した高次DG法が、移動穿孔を伴う二重ブラックホール系を効果的に進化させられるか？
• RQ4 ADER-WENOサブセルリミッターは、ブラックホール特異性付近の数値的不安定性をどれほど効果的に抑制できるか？
• RQ5 FO-CCZ4形式は多次元真空中テストおよび二重ブラックホール進化において、性能と精度の両面でどの程度の効果を示すか？

主な発見

• 解析的計算により全固有構造を確認した結果、特定のゲージ選択のもとでFO-CCZ4形式は強く双曲型であることが示された。
• この系は、純粋な常微分方程式と線形的に退化した場のみを含む低次元双曲型系に分解され、数値的取り扱いが簡素化される。
• 本形式により、高次不連続ガラーキン法を用いた、初めての移動穿孔を伴う二重ブラックホール系の3次元シミュレーションが達成された。
• 一次元、二次元、三次元の真空中テストにおいて、本スキームは高い精度と安定性を示した。
• 適応的メッシュ細分化、局所時間ステッピング、およびADER-WENOサブセルリミッターの組み合わせは、ブラックホール特異性付近の数値的不安定性を効果的に管理した。
• 二重ブラックホール進化に関する予備的結果から、一次形式を用いた高次DG法が、複雑な時空力学を扱う上で実用的であることが示された。