[論文レビュー] Adaptive Finite Elements and Colliding Black Holes

Douglas N. Arnold, Arup Mukherjee+1 more|ArXiv.org|Sep 15, 1997

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

ひとこと要約

本稿では、一般相対性理論における衝突するブラックホールシミュレーションの初期データ問題を解くために、非線形楕円型境界値問題として定式化された非有界3次元領域における、適応的非構造的四面体メッシュを用いた有限要素法を提示する。この手法は、メッシュの適応的 refinement、多重グリッドソルバ、および特化した変分式を組み合わせることで、LIGOによる検出を支援する重力波源の高精度なモデル化を実現する。

ABSTRACT

According to the theory of general relativity, the relative acceleration of masses generates gravitational radiation. Although gravitational radiation has not yet been detected, it is believed that extremely violent cosmic events, such as the collision of black holes, should generate gravity waves of sufficient amplitude to detect on earth. The massive Laser Interferometer Gravitational-wave Observatory, or LIGO, is now being constructed to detect gravity waves. Consequently there is great interest in the computer simulation of black hole collisions and similar events, based on the numerical solution of the Einstein field equations. In this note we introduce the scientific, mathematical, and computational problems and discuss the development of a computer code to solve the initial data problem for colliding black holes, a nonlinear elliptic boundary value problem posed in an unbounded three dimensional domain which is a key step in solving the full field equations. The code is based on finite elements, adaptive meshes, and a multigrid solution process. Here we will particularly emphasize the mathematical and algorithmic issues arising in the generation of adaptive tetrahedral meshes.

研究の動機と目的

• 一般相対性理論における衝突するブラックホールの初期データ問題を解く際の計算的課題に対処すること。
• 2つのブラックホール系に生じる非線形的かつ非有界な3次元楕円型境界値問題を扱える数値フレームワークの開発。
• LIGOの検出作業を支援するため、重力波源の正確なシミュレーションを可能にすること。
• ブラックホール特異点付近での解像度を向上させるために、適応的メッシュの細分化を有限要素法と統合すること。
• 非構造的四面体要素上の多重グリッド解法と変分式を用いて、安定性と効率性を確保すること。

提案手法

• 初期データ問題を非有界3次元領域における非線形楕円型境界値問題として定式化する。
• 非構造的四面体メッシュ上での連続的で一次の基底関数を用いた有限要素法を用いる。
• 事後誤差推定器に基づく適応的メッシュの細分化を実装し、ブラックホール特異点付近に解像度を集約する。
• 有限要素離散化から生じる大規模な線形系を効率的に処理するため、多重グリッド反復ソルバを適用する。
• アインシュタイン制約方程式の物理的制約を尊重する変分式を構築する。
• 計算領域の非有界性を扱うために、ドメイン分割法を採用する。

実験結果

リサーチクエスチョン

• RQ1 一般相対性理論における衝突するブラックホールの初期データ問題に、適応的有限要素法をどのように効果的に適用できるか。
• RQ2 ブラックホールホライズン付近の特異的挙動を捉えるために、最も効果的な誤差推定とメッシュの細分化戦略は何か。
• RQ3 非有界領域における非線形楕円型偏微分方程式の有限要素離散化から生じる大規模かつスパースな線形系に対して、多重グリッドソルバをどのように最適化できるか。
• RQ4 適応的細分化を可能にするために、アインシュタイン制約方程式と整合性を持つ変分式は何か。
• RQ5 メッシュの適応的細分化が、2つのブラックホール初期データのシミュレーションにおける精度と計算コストに与える影響は何か。

主な発見

• 適応的有限要素法は、均一な細分化と比較して、ブラックホール特異点付近での解像度を著しく向上させるとともに、計算コストを削減する。
• 事後誤差推定と局所的メッシュ細分化の組み合わせにより、初期データを支配する非線形楕円型系の高精度な近似が可能になった。
• 多重グリッドソルバは、有限要素離散化から生じる大規模な線形系に対して、効率的な収束を達成した。
• 適切な漸近的境界条件とドメイン切断を用いることで、非有界領域を効果的に処理できた。
• このフレームワークは、重力波源のシミュレーションに堅牢な基盤を提供し、LIGOのデータ解析への応用が可能である。
• 非構造的四面体要素の使用により、複雑なブラックホール配置を幾何学的に柔軟に解像できるようになった。