[論文レビュー] A space-time FEM for PDEs on evolving surfaces

Joerg Grande, Maxim A. Olshanskii+1 more|arXiv (Cornell University)|Mar 2, 2014

Fluid Dynamics and Thin Films参考文献 13被引用数 13

ひとこと要約

本稿では、表面パラメータ化やボリューム領域への拡張を必要とせず、変化する表面におけるPDEの解法を目的とした空間時間有限要素法（FEM）を提示する。この手法は、固定されたEuler型グリッド上での連続的区分線形空間的および不連続的区分線形時間的有限要素関数のトレースを用い、空間および時間で2次精度を達成するとともに、ドロップレットの衝突などのトポロジー変化を最小限の質量損失で安定して扱える。大きな時間ステップでも同様に有効である。

ABSTRACT

The paper studies a finite element method for computing transport and diffusion along evolving surfaces. The method does not require a parametrization of a surface or an extension of a PDE from a surface into a bulk outer domain. The surface and its evolution may be given implicitly, e.g., as the solution of a level set equation. This approach naturally allows a surface to undergo topological changes and experience local geometric singularities. The numerical method uses space-time finite elements and is provably second order accurate. The paper reviews the method, error estimates and shows results for computing the diffusion of a surfactant on surfaces of two colliding droplets.

研究の動機と目的

• 表面パラメータ化やボリューム領域への拡張を必要とせず、変化する表面におけるPDEの解法を対象とした、強固で精度が高く効率的な有限要素法の開発。
• ドロップレットの衝突などのトポロジー変化を伴う表面における輸送および拡散現象のシミュレーションを可能にする。
• 空間時間有限要素法を用いて、空間および時間で2次収束を確保する。
• レベルセットに基づく表面表現に起因する幾何的近似誤差が存在する中でも、弱い意味での質量保存を維持する。
• 固定されたEuler型グリッドと大きな時間ステップを用いて、幾何的特異点を伴う問題に対しても本手法の強健性を示す。

提案手法

• 本手法は、空間時間内に固定されたEuler型メッシュ上で定義された標準的な体積有限要素関数のトレースに基づく空間時間有限要素法を採用する。
• 空間的有限要素は連続的区分線形、時間的有限要素は不連続的区分線形であり、自然な時間積分手順を可能にする。
• 表面はレベルセット関数を用いて暗黙的に表現され、リメッシュなしにトポロジー変化や特異点を扱える。
• 空間時間多様体上の積分を近似するための局所的数値積分手順が採用され、四面体を細分化し、レベルセット関数の線形補間のゼロ等高線を計算する。
• 物質微分および表面微分作用素（接線勾配、Laplace-Beltrami作用素）を含む弱形式を用い、トレース空間上でGalerkin離散化を実施する。
• 離散系は空間時間ピラミッド上での要素単位のアセンブリにより構築され、表面の変化に伴い局所的にメッシュを細分化することで幾何的精度を向上させる。

実験結果

リサーチクエスチョン

• RQ1 表面パラメータ化やボリューム領域への拡張を回避しつつ、変化する表面におけるPDEに対して2次精度を維持できる空間時間FEMを開発できるか？
• RQ2 表面が合体や衝突などのトポロジー変化を経験する場合、本手法はどのように性能を示すか？
• RQ3 レベルセット表現に起因する幾何的近似誤差が、数値解における質量保存性にどの程度影響を及えるか？
• RQ4 複雑な表面変化の過程で、固定されたEuler型グリッドと大きな時間ステップを用いても、本手法は精度と安定性を維持できるか？
• RQ5 幾何的特異点や不連続性を伴う表面変化に対しても、本手法は十分に強健であると言えるか？

主な発見

• 理論的誤差推定と数値実験により、本手法は空間および時間で2次収束を達成していることが確認された。
• 衝突するドロップレット上での界面活性剤拡散問題において、Δt = 1/128 の細かいメッシュで相対的質量損失は約 0.09%（9×10⁻⁴）にまで低下した。
• 粗いメッシュ細分化（l=4, Δt=1/64）では、相対的質量損失は約 0.065 に増加し、細かいメッシュの約 3.6 倍に上昇した。
• 大きな時間ステップ Δt = 1/4（t=1 までに4ステップ）でも、安定かつ精度の高い解が得られ、さらにメッシュの細分化に伴い質量損失は減少した。
• ドロップレットの融合といったトポロジー的特異点に対しても、本手法は固定されたEuler型グリッド上で精度と安定性を維持する強健性を示した。
• 数値結果から、幾何的誤差が質量損失の主な要因であることが確認されたが、メッシュの細分化に伴いその影響は急速に減少した。