[論文レビュー] Building Space-Time Meshes over Arbitrary Spatial Domains

Jeff Erickson, Damrong Guoy+2 more|ArXiv.org|Jun 1, 2002

Computational Geometry and Mesh Generation参考文献 4被引用数 9

ひとこと要約

本稿では、任意の次元における任意の単体的領域に対して、局所的な空間メッシュ品質および特徴サイズに応じて要素の持続時間を適応的に調整する一般化された「テントピッチャー」アプローチを用いて、段階的な時空メッシュを生成する新規なアルゴリズムを提示する。この手法は、1要因あたり定数時間で高品質で非循環的な時空メッシュを構築し、コーン制約を満たし、グローバルな時間ステップを回避することで、時空不連続ガラーキン有限要素問題を線形時間で効率的に解くことを可能にする。

ABSTRACT

We present an algorithm to construct meshes suitable for space-time discontinuous Galerkin finite-element methods. Our method generalizes and improves the `Tent Pitcher' algorithm of Üngör and Sheffer. Given an arbitrary simplicially meshed domain M of any dimension and a time interval [0,T], our algorithm builds a simplicial mesh of the space-time domain Mx[0,T], in constant time per element. Our algorithm avoids the limitations of previous methods by carefully adapting the durations of space-time elements to the local quality and feature size of the underlying space mesh.

研究の動機と目的

• 不連続ガラーキン（DG）有限要素法に適した、任意の空間領域に対してスケーラブルで汎用的な時空メッシュ生成アルゴリズムの開発。
• 従来の手法が固定されたグローバル時間ステップに依存するか、メッシュに制限的な幾何的制約を課えるという限界を克服すること。
• 時空要素がコーン制約を満たすようにすることにより、結合された要素依存関係を回避し、線形時間での解法戦略を可能にすること。
• 波速の変動と局所的な空間メッシュ品質に対応できるように、テントピッチャーアルゴリズムを一般化し、段階的な時間方向のメッシュ細分化を可能にすること。
• 誤差推定値に応じて、局所的な再メッシュ、細分化、粗化などの動的メッシュ操作をサポートすること。

提案手法

• アルゴリズムは、空間時空を通過する前進フロント戦略を用い、局所的な空間メッシュ品質および特徴サイズに従って、段階的に四面体要素を構築する。
• 各フロント頂点が十分に時間的に進むことを保証する修正された進行制約を導入し、各ステップでの最小進行量を制御する調整可能なパラメータ ε を用いる。
• 要素の持続時間は、各点における局所的な空間メッシュ品質および波速に応じて適応的に選択され、固定時間ステップを回避する。
• コーン制約は、隣接する要素間の境界面が影響と依存性のコーンを分離するような形で、逆方向の情報伝播を許可しないようにすることで強制される。
• 元のテントピッチャーを、非一様な時間ステッピングと任意次元の単体的メッシュを扱えるように一般化した。
• メッシュの非循環性を維持し、依存関係グラフが部分順序を形成することを保証することで、要因ごとの効率的解法を可能にする。

実験結果

リサーチクエスチョン

• RQ1 固定されたグローバル時間ステップに依存せず、局所的な空間メッシュ品質に応じて要素の持続時間を適応的に調整する時空メッシュ生成アルゴリズムを設計できるか？
• RQ2 コーン制約をどのように強制すれば、非循環的依存関係を保証し、時空DG問題の線形時間解法を可能にできるか？
• RQ3 テントピッチャーアルゴリズムを、平面的領域に限らない任意の次元における任意の単体的空間領域に一般化できるか？
• RQ4 進行制約パラメータ ε の値の変化が、メッシュ品質、要素数、および時間方向の段階的細分化に与える影響は何か？
• RQ5 誤差推定値に応じて、局所的再メッシュ、細分化、粗化などの動的メッシュ操作をサポートできるように、この手法を拡張できるか？

主な発見

• アルゴリズムは1要因あたり定数時間で時空メッシュを構築し、大規模なシミュレーションにおいて非常に効率的である。
• 進行制約パラメータ ε の異なる値においても、メッシュ品質および要素数は数パーセントの範囲で変動し、より小さい ε 値ではわずかに少ない要素数が得られた。
• 主に鋭角な三角形からなる高品質な空間メッシュでは、進行制約の影響は孤立した領域を除き最小限に抑えられ、メッシュ品質の変動に対して高いロバスト性を示した。
• コーン制約を満たす段階的な時空メッシュを効果的に生成でき、非循環的依存関係を保証し、DG系の線形時間解法を可能にした。
• 波速の変動と局所的な空間メッシュ品質に応じた可変時間ステッピングを可能にし、任意次元の単体的領域を扱えるように一般化したことで、従来の領域幾何や波速挙動に関する制限を克服した。
• 実験結果から、複雑な空間領域に対しても時間方向に段階的なメッシュが生成され、結果として得られる時空要素の持続時間に明確な時間方向の段階的細分化が確認された。