The Tempered Finite Element Method

技術要約：テンパード有限要素法 (TFEM)

問題提起
有限要素法 (FEM) は、非構造格子による複雑な形状の扱いに長けていることから、計算物理学および工学において支配的な数値手法である。しかし、標準的な FEM 定式化は、特定のメッシュ構成、特に「キャップ」（角度が180°に近づく三角形）がバンド状に配置されている場合に、脆弱性を示す。孤立した退化した要素（細長い要素や孤立したキャップ）は収束を妨げないが、キャップのバンドは「ロッキング」現象を引き起こす。このような構成では、数値解が真の解とは無関係に、そのバンドに沿って線形になるよう強制される。その結果、有効なメッシュサイズは、バンドの厚さではなくバンドの長さによって決定されてしまう。これは、標準的な FEM における収束性の喪失を招く。歴史的に、「最大角条件」（すべての角度が180°から離れていること）は収束に不可欠であると考えられてきたが、近年の研究では、孤立した要素についてはそれが厳密には必要ではないことが示されている。課題は、界面追跡手法（X-Mesh など、メッシュのリレーイングによって界面を追跡する手法）において発生しうる、ゼロ測度またはほぼ退化した要素を含むメッシュをどのように扱うかである。

手法
著者らは、キャップのバンドを含むメッシュ上での収束を可能にする、古典的な FEM の修正版であるテンパード有限要素法 (TFEM) を提案している。核心となるアイデアは、冶金における「テンパリング（焼き戻し）」から着想を得ている。焼き戻しが硬化鋼の内部応力を緩和して柔軟性を高めるのと同様に、TFEM は退化したバンド内の材料特性を局所的に「軟化」させることで、数値解が連続的な問題に適合できるようにする。

ポアソン方程式のような楕円型問題の場合、標準的な剛性行列 $K_{ij}$ は、要素写像のヤコビアン行列式 $J$ の逆数を含む。退化した要素では、$J \to 0$ となり、ゼロ除算や無限の固有値を引き起こす。TFEM は、パラメータ $J_{min}$ を用いてヤコビアンをゼロから遠ざけることで、局所的な剛性行列を修正する：
$$\tilde{K}_{ij} = \int_{\hat{E}} \frac{1}{\max(J, J_{min})} [\dots] d\xi$$
これは、弱定式化において可変拡散係数 $d(x)$ を導入することと同等である：
$$d(x) = \begin{cases} 1 & x \in \Omega \setminus B \\ D = \frac{J}{J_{min}} & x \in B \end{cases}$$
ここで $B$ は退化した要素のバンドである。$J_{min}$ は、拡散係数の減少（ロッキングを回避するため）と近似特性の保持とのバランスをとるように、メッシュサイズ $h$ に依存して（具体的には $J_{min} \sim h^{d+1}$、ここで $d$ は空間次元）選択される。

要素の測度が厳密にゼロ（$J=0$）の場合、著者らは TFEM スキームの極限が、界面における解のジャンプを罰則化するモルター法 (mortaring method) に対応することを実証している。

主な貢献

1. 収束の理論的証明： 本論文は、キャップのバンドに対する TFEM の厳密な誤差解析を提供している。最適なテンパリング・パラメータの選択 $D \sim \bar{h}/h^2$ （ここで $\bar{h}$ はバンド幅）を用いることで、本手法が規則的なメッシュ上の標準的な FEM と同様の $O(h)$ の $H^1$ セミノルムにおける収束率を達成することを証明している。
2. モルター法との等価性： 幅が厳密にゼロのバンドに対して、TFEM がジャンプ罰則を伴う非適合モルター法と等価であることを証明し、退化したメッシュの取り扱いと界面結合技術の間の理論的な架け橋を築いた。
3. 複雑な物理への拡張： 本手法は、スカラー・ラプラシアンを超えて以下に拡張された：
   • 線形弾性： 界面を横断する応力解析において、標準的な FEM と同等の精度を示す。
   • 非適合メッシュのモルター法： TFEM が、遷移要素を用いることなく異なる解像度のメッシュをシームレスに接続できることを示す。
   • 高次要素： $J_{min}$ が要素次数に応じて適切にスケールされる限り、多項式次数が 4 に達しても最適な収束率が維持されることを検証した。
   • 移流： 流線アップウィンド/ペトロフ・Galerkin (SUPG) 安定化を用いた移流拡散問題に TFEM を適用し、テンパリングは安定化項に対してのみ必要であることを指摘した。
4. ロバスト性解析： 数値実験により、関係式 $J_{min} = C h^{d+1}$ における定数 $C$ に対する本手法のロバスト性が確認された。誤差は広い範囲の $C$ 値（ほぼ 1 桁の範囲）にわたって最適に近い状態を維持しており、標準的な前処理器を用いれば、線形システムの条件数は管理可能な範囲に留まる。

結果

• 数値実験： 製造解を用いた 2D および 3D ポアソン問題において、キャップのバンドを持つメッシュ上の TFEM は、規則的なメッシュ上の標準的な FEM の収束率を回復する。具体的には、2D では $J_{min} \sim h^3$、3D では $J_{min} \sim h^4$ と選択することで、最適な $L^2$ および $H^1$ 誤差収束が得られる。
• 誤差の感度： $H^1$ 誤差は $L^2$ 誤差よりも $J_{min}$ の選択に対する感度が低く、広いパラメータ範囲でプラトー（停滞）を示す。
• 条件数： バンド幅が減少するにつれて剛性行列の条件数は増大するが、不完全コレスキー前処理を使用することで、反復解法の収束はバンドの厚さに依存しない状態が保たれる。
• ランダムテスト： 2D における 100 回のランダムテストにより、単純な選択（$C=1$）であっても、TFEM の誤差は規則的なメッシュ上の標準的な FEM の誤差とほぼ同一であることが確認された。

意義および主張
本論文は、TFEM が退化した要素に関連する厳格なメッシュ品質要件を緩和することにより、有限要素解析の柔軟性を大幅に向上させると主張している。これにより、以前は動作不能と考えられていた特定のメッシュトポロジー（キャップのバンド）を使用するための、ゼロ測度またはほぼ退化した要素（動的な界面追跡によって生成されるものなど）の使用が可能になる。

著者らは、本手法が既存の FEM コードへの実装が極めて容易であり、ゼロ除算を防ぐためのヤコビアン評価の修正のみを必要とすることを強調している。TFEM は、高品質なメッシュに対する標準的な FEM の代替品としてではなく、これまで動作不能とされていた特定のメッシュ構成（キャップのバンド）の使用を可能にする堅牢な拡張として位置づけられている。本手法は、整合的なメッシュの維持が困難または計算コストが高い X-Mesh などの界面追跡法において特に価値が高い。結論として、パラメータ $C$ の理論的な最適値は導出されているものの、本手法の実用的なロバスト性により、幅広い許容範囲が存在するため、即座に適用可能であるとしている。