An $H^{-1}$ least-squares UnCut FEM on domains defined by a level set function

技術要約：レベルセット関数によって定義された領域における $H^{-1}$ 最小二乗 UnCut FEM

問題提起
本論文は、レベルセット関数 $\phi$（すなわち $\Omega = \{x : \phi(x) < 0\}$）によって領域が暗黙的に定義される、滑らかな境界を持つ $\mathbb{R}^d$ 内の領域 $\Omega$ 上でのポアソン方程式およびストークス方程式の数値解法を取り扱う。このような領域における主要な課題は、非適合有限要素法（unfitted FEM）に固有の「小さなカットセル」問題である。これは、境界が背景メッシュの要素を任意に横切ることで、極小または高度に異方的なサブ要素が生じる問題である。CutFEM や XFEM といった既存の手法は、ゴーストペナルティやセル集約を通じて安定性を確保しているが、これらは曲面カットセル上での正確な数値積分を必要とすることが多く、実装を複雑にする。さらに、最近提案された $\phi$-FEM のような従来の安定化された非適合手法は、安定性を確保するために十分な大きさの安定化パラメータを選択する必要があるが、実際にはその最小値が未知であることが多く、ユーザーによる調整を必要とする。

手法
著者らは、カットされたサブ領域上での数値積分を回避し、支配方程式の $H^{-1}$ 残差を最小化する最小二乗定式化を導入した、新しい UnCut 有限要素法を提案している。

1. 領域の拡張と関数空間: $\phi$-FEM と同様に、手法は計算領域を、物理領域と交差するすべてのフル背景要素からなる $\Omega_h$ へと拡張する。近似解は、有限要素関数とレベルセット関数 $\phi_h$ の積として構成され、これにより離散界面 $\Gamma_h = \{\phi_h = 0\}$ 上でディリクレ境界条件が満たされる。
2. 最小二乗定式化: 標準的なガラーキン定式化の代わりに、本手法は、拡張された作用素の離散 $H^{-1}(\Omega_h)$ ノルムの二乗を含む汎関数を最小化する。これは、離散リースの写像 $(I - \Delta_h)^{-1}$ を通じて定義される補助変数 $w_h$ を導入することによって達成される（ここで $\Delta_h$ は斉次ノイマン境界条件を持つ離散ラプラシアンである）。
   • ポアソン方程式の場合、スキームは、拡張された作用素の残差を補助変数の $H^1$ ノルムにおいて最小化するように設計されており、カット要素上の最小二乗項およびエッジ上のジャンプペナルティによって補強される。
   • ストークス方程式の場合、定式化は同様に運動量方程式の $H^{-1}$ 残差を最小化し、発散拘束およびカット要素上の最小二乗安定化項を含む。
3. 安定化: 手法は、安定性を確保するためにカット要素（$T \in \mathcal{T}_h^\Gamma$）およびエッジ（$F \in \mathcal{F}_h^\Gamma$）に対するペナルティ項を用いる。決定的なことに、最小二乗の構造により、これらの安定化項を固定された単位重みで用いても効果的に機能するため、ユーザーが調整すべき大きなペナルティパラメータを必要としない。

主な貢献

• パラメータフリーの堅牢性: 主な貢献は、安定化パラメータを「十分に大きく」する必要がない定式化を実現したことである。本手法は固定された単位重みで堅牢であり、最適なパラメータ値が事前に不明であるという $\phi$-FEM の実用的な限界に対処している。
• 曲面積分の回避: $\phi$-FEM と同様に、本手法は全背景要素上での体積積分を行うことにより、曲面カットセル上の複雑な数値積分を回避する。
• 理論的解析: 著者らは、以下の通り最適な次数の収束率を厳密に確立している：
  • ポアソン方程式については、最適な $H^1$ ノルム収束が証明されている。
  • ストークス方程式については、速度の最適な $H^1$ ノルム収束、および圧力の（定数を除いた）最適な $L^2$ ノルム収束が導出されている。
• 新しい安定化メカニズム: 本論文は、提案された UnCut 最小二乗 FEM の固有の安定化構造が、ストークス問題における圧力の安定性に悪影響を与えることなく、証明に十分な強さの安定化項を持てることを示している。これは、非適合離散化における一般的な困難である。

結果

• 理論的: レベルセット関数とメッシュの正則性に関する緩やかな仮定の下で、離散解の存在と一意性が証明されている。誤差評価は、メッシュサイズ $h$ およびレベルセット関数の近似誤差 $\|\phi - \phi_h\|_{H^1}$ に関して導出されている。
• 数値的: FEniCS を用いた 2D および 3D のポアソンおよびストークス問題の実験により、理論的な収束率が確認された。
  • 本手法は、ポアソン問題において $H^1$ および $L^2$ ノルムの両方で最適な収束を達成する。
  • ストークス問題では、速度は $H^1$ で、圧力は $L^2$ で最適に収束する。
  • 堅牢性の比較: 元の $\phi$-FEM との直接比較により、提案手法は幅広い安定化パラメータ $\sigma$ に対して精度を維持することが示された。対照的に、元の $\phi$-FEM は、特に圧力変数において、$\sigma$ が小さい場合または大きい場合に顕著な誤差の増大を示す。

意義と主張
本論文は、提案された $H^{-1}$ 最小二乗 UnCut FEM が、非適合有限要素法の堅牢性と簡便性を大幅に向上させることを主張している。$H^{-1}$ 残差を最小化するために最小二乗定式化を活用することで、本手法は、$\phi$-FEM の実装上の利点（カットセル上の積分を避けること）を維持しつつ、ユーザーが調整する大きな安定化パラメータへの依存を排除することに成功している。著者らは、このアプローチがパラメータフリーのスキームを提供し、メッシュ要素が任意に小さかったり不規則にカットされていたりする場合でも、カットの構成に対して堅牢であることを主張している。本研究は、ポアソンおよびストークス問題の両方に対して最適な次数の誤差評価を提供することで、これらの主張に対する厳密な理論的基礎を確立している。

著者らは、今後の研究方向として、最適な次数の $L^2$ 誤差評価の理論的導出、ノイマン／ロービン境界条件への拡張、および時間依存問題への適用を挙げている。