[論文レビュー] Obtaining higher-order Galerkin accuracy when the boundary is polygonally approximated

Todd F. Dupont, Johnny Guzmán+1 more|arXiv (Cornell University)|Jan 9, 2020

Advanced Numerical Methods in Computational Mathematics参考文献 9被引用数 4

ひとこと要約

本稿では、多角形で近似された領域上でポアソン方程式を解くための対称的高次ガラーキン法を提示している。標準的な有限要素法では幾何的誤差が収束速度を低下させるが、真の境界からの符号付き距離に基づく境界層補正を変分形式に組み込むことで、この問題を是正する。この手法により、例えば2次要素では$O(h^3)$の最適収束速度を達成しながらも、対称性を保ったまま、ブラムブル・デュポン・トーメー法のような従来の非対称的手法の限界を克服する。

ABSTRACT

We study two techniques for correcting the geometrical error associated with domain approximation by a polygon. The first was introduced some time ago \cite{bramble1972projection} and leads to a nonsymmetric formulation for Poisson's equation. We introduce a new technique that yields a symmetric formulation and has similar performance. We compare both methods on a simple test problem.

研究の動機と目的

• 滑らかな領域を多角形で近似した場合に生じる標準的ガラーキン有限要素法の収束が最適でない問題に対処すること。
• 高次多項式近似（例えば2次および3次要素）において最適収束速度を回復する対称的変分形式を開発すること。
• ブラムブル・デュポン・トーメー法のような既存の高次補正手法における非対称な双線形形式の問題を克服すること。
• 領域近似に起因する幾何的誤差が、符号付き距離関数および法線微分を含む境界積分項によって是正可能であることを示すこと。

提案手法

• 符号付き距離$\delta$および法線微分を含む境界項を追加することで、標準ガラーキン形式を修正する対称的双線形形式$M_h(u,v)$を導入する。
• 整合性と対称性を確保するため、$\int_{\partial\Omega_h} \left(\gamma \frac{\delta}{h} - 1\right) \left( \delta \frac{\partial u}{\partial n} \frac{\partial v}{\partial n} + \frac{\partial u}{\partial n} v + \frac{\partial v}{\partial n} u \right) ds$ のような補正項を含む。
• 境界付近における解のテイラー展開を用いて、多角形境界と真の境界とのずれを補償する補正項を導出する。
• ディリクレ境界条件を伴うポアソン方程式にこの手法を適用し、体積積分および境界積分を含む変分形式を採用する。
• 安定性と最適収束を保証するために、安定化パrameter$\gamma$を用いたロビン型境界条件を変分形式に組み込む。
• 曲線境界を多角形で近似した滑らかな領域に対して、数値実験を実施。標準法および非対称補正法と比較する。

実験結果

リサーチクエスチョン

• RQ1 領域が多角形で近似された場合に、最適収束速度を回復する対称的変分形式を構築可能か？
• RQ2 対称性を保ちつつ高次精度を達成するための境界補正項の最適構造は何か？
• RQ3 対称的高次補正手法を用いた場合、標準法や非対称法と比較して誤差分布はどのように変化するか？
• RQ4 符号付き距離関数$\delta$を用いることで、標準的なニツェ型手法と比較して収束性はどの程度向上するか？

主な発見

• 提案された対称的手法$M_h$は、2次有限要素を用いる場合に$H^1$ノルムで最適な$O(h^3)$収束を達成し、多角形境界近似に対しても最適次数を回復する。
• 数値実験の結果、対称的手法は非対称ブラムブル・デュポン・トーメー法と同等の収束速度を示し、$k=2$では$L^2$誤差が$O(h^3)$、$k=3$では$O(h^4)$で減少する。
• 対称手法の誤差分布は境界付近に集中しており、内部の誤差はゼロに収束するため、汚染効果が効果的に抑制されていることが示唆される。
• 中空領域のテストケースでは、$h_{\text{max}} = 0.036$、$k=2$のとき$L^2$誤差が$1.54 \times 10^{-5}$に達し、最適収束が確認された。
• $\delta^\ell \partial^n u / \partial n^n$項を用いて高次補正を組み込むことで、安定性と対称性を維持しつつ、高次精度を実現する。
• 非対称系に起因する計算上の不利要因を回避するため、反復ソルバーや大規模シミュレーションに適した対称形式を採用している。