[論文レビュー] Discontinuous Galerkin methods for the $p$--biharmonic equation from a discrete variational perspective

Tristan Pryer|CentAUR (University of Reading)|Sep 18, 2012

Advanced Numerical Methods in Computational Mathematics参考文献 26被引用数 18

ひとこと要約

本稿では、有限要素ヘシィアンに基づく離散変分式を用いた、$p$–双調和方程式に対する不連続ガレルキン（dG）有限要素法を提案する。この手法は、離散解およびヘシィアンが連続解に収束することを保証し、$p=2$ では最適収束率が理論的に証明され、$p>2$ では数値的に観察されており、$p \in (1,\infty)$ において堅牢性と精度を示している。

ABSTRACT

We study discontinuous Galerkin approximations of the $p$--biharmonic equation from a variational perspective. We propose a discrete variational formulation of the problem based on a appropriate definition of a finite element Hessian and study convergence of the method (without rates) using a weak lower semicontinuity argument. We present numerical experiments aimed at testing the robustness of the method. We also note a superconvergence effect for some values of $p$.

研究の動機と目的

• 第四階非線形偏微分方程式である $p$–双調和方程式に対して、安定かつ収束する不連続ガレルキン有限要素法を構築すること。既存の有限要素法の定式化が存在しない。
• 非一致型で局所的に滑らかな有限要素空間におけるヘシィアンの近似という課題に取り組み、第四階問題の変分的離散化にとって不可欠である。
• 弱下半連続性および離散関数解析を用いて、離散解およびヘシィアンが連続解に収束することを確立すること。
• 異なる $p$ 値および多項式次数において、手法の堅牢性と収束率を数値的に検証すること。
• 曲率に基づく偏微分方程式など、より複雑な第二階変分問題への拡張のための枠組みを構築すること。

提案手法

• 補助変数として有限要素ヘシィアンを導入することで、変分的枠組み内でのヘシィアンの弱い強制を可能にする離散変分式を構築する。
• ヘシィアンを有限要素空間で近似する混合定式化を採用し、可能な限り低次の表現を用いることで計算コストを低減する。
• 離散問題は変分原理から導出され、適切な数値フラックスおよび安定化項を用いることで、一貫性と安定性を保証する。
• 収束性は [DPE10] の枠組みを用い、離散ソボレフ空間における弱下半連続性およびコンパクトネスの議論に依拠して証明する。
• ヘシィアンは、局所多項式上での離散勾配作用素によって定義され、$C^0$-適合だが不連続なヘシィアン行列の近似を可能にする。
• FEniCS/DOLFIN を用いた数値実験では、メッシュの適応的 refinement を行い、自由度を削減するための補助ヘシィアン変数を用いた混合系を解く。

実験結果

リサーチクエスチョン

• RQ1 離散変分的アプローチを用いて、$p$–双調和方程式に対する不連続ガレルキン法を定式化できるか？
• RQ2 非一致型設定において、局所滑らかな有限要素関数のヘシィアンを一貫的かつ安定的に近似する方法は何か？
• RQ3 離散解およびヘシィアン近似の理論的・数値的収束率は、それぞれどのように確立できるか？
• RQ4 $p \in (1,\infty)$ の異なる値、特に $p=2$ および $p>2$ に対しても、この手法は堅牢か？
• RQ5 この枠組みは、曲率や高階微分作用素を含むより複雑な変分問題へ拡張可能か？

主な発見

• $p=2$ の場合、最適収束率が達成される：$k=2$ のとき $\|u - u_h\|_{L^2} = \mathcal{O}(h^2)$ であり、$k>2$ のとき $\mathcal{O}(h^{k+1})$ となる。ここで $k$ は多項式次数を表す。
• 離散ヘシィアンの収束速度は $\|\Delta u - \mathscr{D}[u_h]\|_{L^p} = \mathcal{O}(h^{k-1})$ であり、理論的期待と整合的である。
• 数値実験により、$p=2$、$p=2.1$、$p=10$ において理論的収束率が確認され、さまざまな非線形性において堅牢性が示された。
• $p$ が 2 から大きく離れた場合でも、手法は安定かつ収束を示し、非線形第四階問題への広範な適用可能性を示している。
• ヘシィアンに $\mathbb{W} = \mathbb{P}^{k-1}$ の低次の有限要素空間を用いることで、$\mathbb{V}^{d\times d}$ 全体の空間を用いる場合と同等の収束率が得られ、計算効率が向上する。
• ヘシィアンに補助変数を導入した実装により、保存領域および計算コストが削減されつつ、精度や収束性に影響を与えない。