Towards $C^0$ finite element methods for fourth-order elliptic equation. Part I: general boundary conditions

技術要約：高次楕円方程式に対する $C^0$ 有限要素法に向けて、第I部：一般境界条件

問題設定
本論文は、多角形領域 $\Omega \subset \mathbb{R}^2$ における高次楕円方程式、具体的には双調和方程式 $\Delta^2 u = f$ の数値近似を扱う。境界 $\partial \Omega$ は、互いに素なディリクレ（$\Gamma_D$）およびノイマン（$\Gamma_N$）セグメントに分割される。主な課題は、この問題に対して効果的な $C^0$ 有限要素法を開発することにある。ナヴィエ境界条件（$\Gamma_D = \partial \Omega$）を伴う双調和方程式に対する標準的なアプローチは、しばしば問題をポアソン方程式の系へとデカップリング（分離）することに依存しているが、この「素朴な（naive）」混合定式化は、非凸領域や特定の混合境界条件の設定において、真の解へ収束しない。この失敗は、**サポンギャン・パラドックス（Sapongyan paradox）**に起因しており、これは、デカップリングされた混合定式化の解空間が、元の双調和問題の解空間と一致しないために起こる。これは、再入角（re-entrant corners）から生じる低正則性項によるものである。

手法
著者らは、解が正しいソボレフ空間（$H^2(\Omega)$）内に留まることを保証するために、双調和問題をポアソン方程式の系へと分解し、さらに追加の中間ポアソン問題を付加した修正混合定式化を提案している。

1. 解空間の解析： 本論文では、$H^2(\Omega)$（適切な境界条件を伴う）から $L^2(\Omega)$ へのラプラス作用素 $\Delta$ の値域を解析している。これにより、多角形領域において、大きな内部角を持つ場合、値域 $S = \Delta(V^2)$ は $L^2(\Omega)$ の真の部分閉部分空間であり、有限次元の直交補空間 $S^\perp$ を持つことが特定される。
2. 基底の構成： 著者らは、$S^\perp$ のための基底として、最大の内部角 $\omega$ の近傍で定義される関数 $\xi_m$ を構成している。これらの関数は、特異項 $s_m(r, \theta)$（ラプラス作用素の角度成分の固有値から導出）と、補助的なポアソン問題を解くことで得られる正則な補正項 $\zeta_m$ によって構成される。直交補空間の次元 $d_\perp$ は、最大の内部角の大きさと、隣接するエッジにおける具体的な境界条件（ディリクレ／ノイマン）の組み合わせに依存する。
3. 修正アルゴリズム： 提案されている $C^0$ 有限要素アルゴリズム（アルゴリズム 3.1）は、以下の手順で進行する：
   • 補助関数 $w$ のための主要なポアソン問題を解く。
   • 直交補空間のための基底関数 $\xi_{m,h}$ を計算するために、補助的なポアソン問題を解く。
   • $w$ を部分空間 $S$ へ射影する（グラム行列システムを用いて係数を計算することで、$w$ から $S^\perp$ 内の成分を取り除く）。
   • 射影されたソース項 $w_S$ を用いて第2のポアソン問題を解き、最終的な近似 $u_h$ を得る。
4. ノイマン条件への拡張： 本手法は、適合条件 $\int_\Omega f = 0$ が必要となる純ノイマンケース（$\Gamma_N = \partial \Omega$）にも拡張されている。また、再入角において直交補空間の次元が 1 になることが示されている。

主な貢献

• サポンギャン・パラドックスの解決： 本論文は、一般的な多角形領域および混合境界条件において、双調和方程式の素朴なデカップリングが失敗することを実証する厳密な理論的枠組みを提供している。著者らは、真の解への収束を保証する修正定式化を提案している。
• 一般化された分解： ナヴィエ条件に限定されていた従来の作業とは異なり、本手法は任意のディリクレ条件とノイマン条件の組み合わせを扱うことができる。追加で必要となるポアソン方程式の数（0、1、または2）は、最大の内部角と局所的な境界条件のタイプによって明示的に決定される。
• 誤差解析： 提案された $C^0$ 線形有限要素法の厳密な誤差評価を確立している。解 $u_h$ の $H^1$ ノルムにおける収束率が $O(h^{\min\{1, 2\alpha\}})$ であることを証明している。ここで $\alpha$ は内部角と境界条件に依存する。具体的には、混合境界条件で角度 $\omega > \pi$ の場合、収束率は特異性に制限されるが、ナヴィエ／ノイマンの場合には、しばしば最適なレートが回復される。
• 数値検証： 様々な多角形領域（L字型や再入角を持つ幾何学的形状を含む）を用いた広範な数値実験により、理論的予測を検証している。これらの結果は、直交補空間が非自明である場合に素朴な手法が誤った解に収束する一方で、提案された修正手法が参照解（$C^0$-IPDGにより計算されたもの）に収束することを裏付けている。

結果と意義
本論文は、提案された修正混合定式化が、一般的な多角形領域における双調和方程式を解くために適格（well-posed）かつ効果的であることを示している。数値実験は、提案手法が予測された収束率を達成し、素朴な混合定式化で見られる不一致を正常に解決することを証明している。

本研究の意義は、高次の問題（これまで $C^1$ 要素の複雑さや単純なデカップリング戦略の収束失敗によって困難であったもの）に対して、$C^0$ 有限要素法（標準的な連続区分的多項式を利用するもの）の計算効率を拡張できる能力にある。境界条件と領域の幾何学が解空間に与える具体的な影響を特徴付けることで、著者らは、実装の単純さを損なうことなく数学的正当性を保証する堅牢なアルゴリズムの枠組みを提供している。本論文は一連の研究の第I部であり、後続の著作においてより複雑な低次項を伴う高次方程式を扱うための理論的基礎を築くものである。