[論文レビュー] Generating Initial Data in General Relativity using Adaptive Finite Element Methods

Burak Aksoylu, David Bernstein+2 more|ArXiv.org|Jan 21, 2008

Advanced Numerical Methods in Computational Mathematics参考文献 104被引用数 5

ひとこと要約

本稿では、一般相対性理論における初期データ生成のための適応型有限要素法（AFEM）を提示している。アインシュタイン制約方程式の共形形式化を解くことにより、弱形式、残差に基づく後験的誤差指標、およびFETKソフトウェアパッケージを用いて、最適収束性を達成し、複雑な二重コンパクト物体の幾何構造に対しても高い精度と効率性を発揮する。

ABSTRACT

The conformal formulation of the Einstein constraint equations is first reviewed, and we then consider the design, analysis, and implementation of adaptive multilevel finite element-type numerical methods for the resulting coupled nonlinear elliptic system. We derive weak formulations of the coupled constraints, and review some new developments in the solution theory for the constraints in the cases of constant mean extrinsic curvature (CMC) data, near-CMC data, and arbitrarily prescribed mean extrinsic curvature data. We then outline some recent results on a priori and a posteriori error estimates for a broad class of Galerkin-type approximation methods for this system which includes techniques such as finite element, wavelet, and spectral methods. We then use these estimates to construct an adaptive finite element method (AFEM) for solving this system numerically, and outline some new convergence and optimality results. We then describe in some detail an implementation of the methods using the FETK software package, which is an adaptive multilevel finite element code designed to solve nonlinear elliptic and parabolic systems on Riemannian manifolds. We finish by describing a simplex mesh generation algorithm for compact binary objects, and then look at a detailed example showing the use of FETK for numerical solution of the constraints.

研究の動機と目的

• 一般相対性理論におけるアインシュタイン制約方程式に起因する結合非線形楕円型系を解くための堅牢な数値フレームワークの開発。
• 特に複雑な幾何構造を示すコンパクトな二重物体系のシミュレーションに適した正確な初期データの構築という課題に取り組む。
• 残差に基づく誤差指標を用いた適応的メッシュ・リファインメントを実装し、計算コストを最小限に抑えながら解の精度を向上させる。
• 適応的でマルチブロック、四面体メッシュ上でのADM質量や共形キリングベクトルといった物理的量の高精度計算を可能にする。
• 低レベルの有限差分離散化を回避するため、弱形式の指定のみを要する柔軟で使いやすい有限要素実装（MC/FETK経由）を提供する。

提案手法

• ヨークの共形分解を用いて、結合されたハミルトニアンおよび運動量制約方程式の弱形式を導出し、ガラーキン型有限要素離散化を可能にする。
• 適応的メッシュリファインメントを導くために、事前および事後誤差推定を適用し、局所的メッシュリファインメントの意思決定に残差に基づく事後誤差指標を用いる。
• コンパクトな二重物体を表すマルチブロック・スフィア型領域における四面体メッシュリファインメントのための適応的単体二等分戦略を実装する。
• FETK有限要素ソフトウェアパッケージを用いて、複雑な多様体トポロジー（マルチチャートアトラスおよび境界条件を含む）を処理する。
• 非線形代数系を効率的に解くために、代数的多重グリッドソルバー（HBMG）を用いたニュートンベースの続次法を採用する。
• 弱形式およびその線形化をC言語で直接コード化し、ユーザー作業を最小限に抑え、相対性理論の応用において高い抽象化と再利用性を実現する。

実験結果

リサーチクエスチョン

• RQ1 適応型有限要素法は、一般相対性理論における非線形的で結合されたアインシュタイン制約方程式系にどのように効果的に適用可能か？
• RQ2 共形形式化された制約の文脈において、メッシュリファインメントを誘導するのに最も効果的な後験的誤差指標は何か？
• RQ3 高次の精度と最適収束性は、複雑なマルチブロック領域上での適応型有限要素法を用いて、ハミルトニアンおよび運動量制約に対して達成可能か？
• RQ4 適応的で非構造的四面体メッシュ上でのADM質量および共形キリングベクトルは、コンパクトな二重物体系に対してどのように高精度に計算可能か？
• RQ5 従来の有限差分法やスペクトル法と比較して、FETKベースの実装はどの程度実装負荷を低減するか？

主な発見

• 適応型有限要素法は、理論的事前および事後誤差推定の支援を受けて、アインシュタイン制約方程式に対して最適収束率および最適計算複雑性を達成する。
• 残差に基づく事後誤差指標は、高誤差領域を的確に特定し、追加コストを最小限に抑えながら解の精度を向上させる有効なメッシュリファインメントを可能にする。
• MC/FETKフレームワーク上での実装により、ユーザーは弱形式およびその線形化のみを指定すればよく、従来の有限差分スキームと比較して著しくコード作成の負担が軽減される。
• R_{\text{out}} = 30、R_{\text{med}} = 7、a_c = 1.5 を有する13パッチのマルチブロック領域上での数値結果は、N = 32, 36, 40 の解像度でブリル波初期データに対して収束を示し、N = 16, 24, 36 を用いたx軸上での点ごとの収束が確認された。
• 本手法は、適応的四面体メッシュ上でのADM質量および共形キリングベクトルの計算に成功し、二重コンパクト物体系の初期データの物理的解釈を高精度で可能にする。