Error estimates for the interpolation and approximation of gradients and vector fields on protected Delaunay meshes in $\mathbb{R}^d$

技術要約：保護されたデローネメッシュにおける補間および近似の誤差推定

問題提起
本論文は、次元 $d > 2$ における単体メッシュ上の補間および近似に関する数学的理解における決定的な空白を扱っている。標準的なデローネメッシュは、二次元（$\mathbb{R}^2$）においては広く利用され、十分に理解されているが、高次元におけるその挙動はそれほど好ましくない。具体的には、$\mathbb{R}^d$ ($d>2$) における標準的なデローネ・三角形分割には、「スリバー（sliver）」と呼ばれる要素が含まれることが多い。これらは、最小高度と最大辺長の比が極めて小さい単体である。これらのスリバーは、厚みがゼロに限りなく近づくことがあり、それが勾配やベクトル場の補間および近似の品質を低下させる原因となる。

既存の文献では、標準的なデローネメッシュ上での区分線形補間（粗さの最小化や特定の汎関数など）に関する堅牢な誤差推定を提供している。しかし、高次の区分多項式補間に関する結果は乏しく、それらは暗黙的なメッシュ依存の定数に依存しているか、あるいはデローネの設定に明示的に適合していないことが多い。さらに、標準的なデローネメッシュは特定の汎関数を最小化するが、高次元において要素の厚さの下限を保証するものではないため、誤差推定の不安定性を招く可能性がある。

手法
著者らは、**保護されたデローネメッシュ（protected Delaunay meshes）**における勾配およびベクトル場の補間と近似に関する明示的な誤差推定を開発している。保護されたデローネメッシュは、各単体の外接球に半径 $r \ge \delta$ の球面バッファ領域を付加した修正された空球条件を満たすものである（ここで $\delta$ は「保護」パラメータである）。この保護により、病理的なスリバーの形成が防止される。

手法は以下の通りである：

1. 粗さとノルムの等価性： 著者らは、「粗さ（roughness）」の概念（もともと $\mathbb{R}^2$ において Rippa によって定義されたもの）を $\mathbb{R}^d$ へと拡張している。彼らは、エッジに沿った方向微分平方和に基づく「粗さ汎関数」$\Psi(\nabla v)$ を定義する。そして、この汎関数が勾配の標準的な $L^2$ ノルムと等価であることを証明する。このとき、等価定数はメッシュの幾何学的性質（具体的には最小の厚さ $C_\Xi$ と最大辺長汎関数 $C_\Upsilon$）に依存する。
2. 補間誤差の導出： 古典的な補間理論（ラグランジュ補間）と上述の等価性を用いて、勾配補間の誤差の上界を導出する。彼らは、誤差を関数のヘッセ行列、メッシュ間隔 $h$、および特定のメッシュ依存定数に関連付けている。
3. 楕円型問題に対する近似誤差： 楕円型問題（具体的には $-\nabla^2 u = f$）に対して、セアの補題（Cea's Lemma）を利用して、勾配近似誤差 ($\nabla(u - u_h)$) を抑え込む。彼らは、この誤差をメッシュの正則性パラメータ $\sigma(K)$ および最大最小包含球半径 $R_{\max}$ に結びつけている。
4. ベクトル場への一般化： この解析をスカラー勾配から一般的なベクトル場へと一般化し、同様の誤差境界が適用されることを示す。

主要な貢献と結果
本論文は、高次の多項式補間および近似に関する、幾何学的に依存した明示的な誤差推定を提供している。主な結果は以下の通りである：

• 明示的な誤差境界： 著者らは以下の形式の不等式を導出している：
  $$\|\nabla v - I_h(\nabla v)\|_{L^2(\Omega)} \le C \tilde{C}(\mathcal{T}_h) h^r \|\nabla(\nabla v)\|_{L^{2\varrho}(\Omega)}$$
  ここで、$\tilde{C}(\mathcal{T}_h)$ はメッシュの特性（Rajanの汎関数 $\Theta$、最小の厚さ $C_\Xi$、直径比 $C_\Delta$、および最大最小包含半径 $R_{\max}$）によって明示的に定義される定数である。
• 保護の役割： 解析は、補間の品質が単体の最小の厚さ（$C_\Xi$）に大きく依存することを実証している。
  • 標準的なデローネメッシュ ($d>2$) では、$C_\Xi$ がゼロに近づく可能性があり、誤差定数を増大させる。
  • 保護されたデローネメッシュでは、$C_\Xi$ は保護パラメータ $\delta$ の関数によって下限が抑えられる（具体的には $C_\Xi \ge \delta^2 / (8d\varepsilon^2)$）。
  • したがって、保護されたメッシュ上では、メッシュの正則性パラメータ $C_\sigma$（厚さに反比例するもの）が上限に抑えられる。
• メッシュパラメータの最適化： 本論文は、標準的なデローネメッシュがRajanの汎関数 $\Theta$（平方辺長の加重和）および最大最小包含半径 $R_{\max}$ を最小化することを確認している。しかし、保護されたデローネメッシュは、保護パラメータ $\delta$ を通じて最小の厚さ $C_\Xi$ の下限を保証するために不可欠である（標準的なメッシュではこれを達成できない）。この組み合わせにより、誤差推定における重要な定数が制御された状態に保たれる。

意義と主張
本論文の主要な貢献は、任意の次元における保護されたデローネメッシュ上での高次補間および近似に関する、明示的で解釈可能な誤差推定を提供することであると主張している。

• 空白の解消： 区分線形補間に焦点を当てたり、暗黙的な定数を用いたりした従来の成果とは異なり、本研究は高次の誤差境界をメッシュの幾何学的特性（厚さと正則性）に明示的に結びつけている。
• 保護の正当化： これらの結果は、高次元における保護されたデローネメッシュの使用を数学的に正当化している。非ゼロの保護 $\delta$ を確保することで、最小の厚さがゼロから離れた値に抑えられることが保証され、それによって誤差定数が任意に大きくなることを防いでいる。
• 実用的な意味合い： 著者らは、これらの推定値が「数値誤差を低減するための数学的に裏付けられたアイデア」を提供すると述べている。有限要素法を用いて勾配の補間やベクトル場の近似を行う実務者は、特に高次の多項式を使用する場合、安定性と精度を確保するために、保護されたデローネメッシュの採用を真剣に検討すべきであると彼らは主張している。

本論文の範囲は控えめであり、これらの結果は十分に滑らかな関数（区分 $W^{2,p}$ または $W^{1,m}$）に適用されるものであり、より不規則な解に対しては他の戦略が必要になる可能性があると述べている。本研究は新しいメッシュ生成アルゴリズムを提案するものではなく、既存の保護メッシュ生成技術を検証するための理論的な誤差解析を提供するものである。