Analysis of a finite element method for second order uniformly elliptic PDEs in non-divergence form

Weifeng Qiu による論文「非発散型の 2 階一様楕円型偏微分方程式に対する有限要素法の解析」の詳細な技術的概要を以下に示す。

1. 問題の定義

本論文は、リプシッツ多面体領域 $\Omega \subset \mathbb{R}^d$ ($d=2, 3$) 上で定義された 2 階偏微分方程式 (PDE) の 2 種類のクラスに対する数値解析を取り扱っている：

1. 非発散型の線形楕円型 PDE:
   $$A : D^2u + b \cdot \nabla u + cu = f \quad \text{in } \Omega, \quad u = 0 \text{ on } \partial\Omega$$
   ここで、$A$ は対称かつ一様正定値行列であり、$b, c$ は有界な係数である。

2. 一様楕円型ハミルトン・ヤコビ・ベルマン (HJB) 方程式:
   $$\sup_{\alpha \in \Lambda} [A_\alpha : D^2u + b_\alpha \cdot \nabla u + c_\alpha u - f_\alpha] = 0 \quad \text{in } \Omega, \quad u = 0 \text{ on } \partial\Omega$$
   ここで、$\Lambda$ はインデックス集合であり、係数は一様楕円性と有界性の条件を満たす。

主要な課題:

• 非発散型: 標準的な弱定式化（部分積分に基づく）は、$A$ または領域の高い正則性を仮定しない限り、2 階項 $A:D^2u$ を試験関数へ容易に移行できないため、直接適用できない。
• 解空間: 粘性解ではなく、ソボレフ空間 $W^{2,p}(\Omega)$ における強解に解析の焦点が当てられている。
• 領域の幾何学: 既存の文献でしばしば要求される $C^{1,1}$ 領域や凸領域よりも一般的である、2 次元の非凸多角形を含むリプシッツ多面体領域の扱いを目的としている。
• 係数の正則性: 最適収束率を維持しつつ、HJB 方程式に対しては特に不連続な係数にも対応できる必要がある。

2. 手法

著者は、線形 PDE と HJB 方程式の両方に対して、統合された有限要素法 (FEM) を提案している。

A. 離散化と演算子

• メッシュ: 準一様、単体、適合的な三角分割 $\mathcal{T}_h$。
• 有限要素空間: 多項式次数 $r \geq 2$ に対して $V_h = H^1_0(\Omega) \cap P_r(\mathcal{T}_h)$。
• 離散勾配 ($\nabla_h$): この手法は、メッシュ要素上での部分積分と内部面をまたぐジャンプ項を通じて定義される特定の離散勾配演算子を利用する。この演算子により、基底関数の $C^1$ 連続性を必要とせずにヘッシアン $D^2u$ の近似を可能にする。
  $$(\nabla_h v)_i = \partial_{h, x_i} v$$
• 重み関数: 非発散型を扱うために、方程式は重み関数 $\gamma$（HJB の場合は $\gamma_\alpha$）を乗算する。
  • 線形の場合: $\gamma(x) = \frac{\text{Tr} A(x)}{|A(x)|^2}$。
  • これにより、方程式は離散演算子に適した形 $\tilde{A} : D^2 u + \dots = \tilde{f}$ に変換される。

B. 数値スキーム

• 線形の場合 (式 14): すべての $v_h \in V_h$ に対して、$u_h \in V_h$ を見つける：
  $$\int_\Omega (\tilde{A} : \nabla_h(\nabla u_h) + \tilde{b} \cdot \nabla u_h + \tilde{c} u_h) v_h \, dx = \int_\Omega \tilde{f} v_h \, dx$$
• HJB の場合 (式 16): すべての $v_h \in V_h$ に対して、$u_h \in V_h$ を見つける：
  $$\left( \sup_{\alpha \in \Lambda} \gamma_\alpha [A_\alpha : \nabla_h(\nabla u_h) + b_\alpha \cdot \nabla u_h + c_\alpha u_h - f_\alpha], v_h \right)_\Omega = 0$$
  注: $\Lambda$ が単一要素を含む場合、HJB スキームは線形スキームと完全に一致する。

C. 解析的枠組み

解析は、以下のように定義される離散 $W^{2,p}$ ノルムに依存している：
$$\|v\|_{W^{2,p}_h(\Omega)} = \|D^2_h v\|_{L^p(\Omega)} + \left( \sum_{F \in \mathcal{F}_h^I} h_F^{1-p} \|\llbracket \nabla v \rrbracket\|_{L^p(F)}^p \right)^{1/p}$$
主要な解析ツールには以下が含まれる：

1. 大域的 $W^{2,p}$ 評価: リプシッツ多面体上の連続問題に対する強解の適切性を証明する。
2. 離散コンパクト性: 離散 $W^{2,p}_h$ ノルムで有界な数列が $W^{1,p}$ において収束部分列を持つことを確立する。
3. 一意性の議論: 以前の作品で使用された特殊な双対性議論を回避し、強解の一意性を用いて背理法により安定性を証明する。
4. 上限の可測性: HJB 方程式について、インデックス集合 $\Lambda$ 上の上限が可測であり、可算部分集合 $\tilde{\Lambda}$（仮定 6）によって近似可能であることを証明する。

3. 主要な貢献

A. 統合された枠組み

本論文は、単一の FEM 定式化が線形非発散型 PDE と完全非線形 HJB 方程式の両方を効果的に解くことを示している。インデックス集合が単元集合である場合、HJB スキームは自然に線形スキームに退化する。

B. 領域および係数の仮定の緩和

• 領域: 最適収束結果を $C^{1,1}$ 領域または厳密に凸な領域からリプシッツ多面体領域へ拡張する。
  • $d=3$ の場合、$\Omega$ は凸でなければならない。
  • $d=2$ の場合、$\Omega$ は非凸となり得る（$p$ に関する制限付き）。
• 係数: 本手法は、可算な稠密部分集合を含む特定の近似条件（仮定 6）を満たす限り、HJB 方程式において不連続な係数を許容する。これは、多くの既存の DG 法が要求する連続性の仮定（$A \in C^0$）を緩和するものである。

C. $W^{2,p}$ における最適収束

本論文は、広範な $p$ に対して離散 $W^{2,p}_h$ ノルムにおける最適収束率を証明している：

• 凸多面体 ($d=2,3$): $p \in (1, 2]$。
• 非凸多角形 ($d=2$): $p \in (4/3 - \epsilon_1, 4/3 + \epsilon_2)$。
  これは、$L^2$ または $H^2$（$p=2$）での収束のみを保証する手法に対する重要な改善である。

D. 新規な安定性証明

連続問題の適切性に依存する特殊な双対性議論に依存する以前の作品（Feng らなど）とは異なり、本論文は離散コンパクト性と強解の一意性を用いて安定性を確立している。このアプローチは、非凸領域や不連続な係数に対してより頑健である。

4. 主要な結果

1. 適切性 (定理 3.3, 4.1):

   • 指定された $A, b, c$ および領域の条件の下で、線形 PDE と HJB 方程式の両方に対して、$W^{2,p}(\Omega) \cap W^{1,p}_0(\Omega)$ における強解の存在と一意性が証明されている。
   • HJB 方程式については、可算部分集合の仮定を通じて不連続な係数を許容するよう、Cordes 条件が一般化されている。

2. 安定性 (定理 3.2, 4.3):

   • 離散解 $u_h$ は安定性評価を満たす：$\|u_h\|_{W^{2,p}_h} \leq C \|f\|_{L^p}$。
   • この安定性は、領域の幾何学に依存する特定の $p$ の範囲で成り立つ。

3. 最適収束 (定理 3.3, 4.3):

   • 誤差評価が導出される：
     $$\|u - u_h\|_{W^{1,p}} + \|u - u_h\|_{W^{2,p}_h} \leq C \left( \inf_{\chi_h \in V_h} (\|u - \chi_h\|_{W^{1,p}} + \|u - \chi_h\|_{W^{2,p}_h}) + \|P_h(D^2u) - D^2u\|_{L^p} \right)$$
   • これは、厳密解が十分に滑らかな場合、手法が離散 $W^{2,p}$ ノルムにおいて最適収束率を達成することを意味する。

5. 意義

• 理論と実践の架け橋: 本論文は、離散化が極めて困難である非発散型 PDE を解くために、標準的な有限要素空間（$C^0$ 連続性を持つ）を使用するための厳密な理論的基盤を提供している。
• 非凸性の扱い: 非凸な 2 次元多角形への解析を拡張することにより、この手法は、領域が完全に凸であることが稀である、より広範な実用的な工学および物理的問題に適用可能となる。
• HJB への応用: 係数の連続性の仮定の緩和は、係数が不連続となり得る（例えば、スイッチング・レジームなど）確率制御問題に適した手法を可能にする。これは、以前の $C^0$ 適合法が失敗したシナリオである。
• 統合されたアプローチ: 1 つの数値スキームが線形および完全非線形の両方のケースを処理することを示すことは、複雑なシステムの実装と理論的解析を簡素化する。

要約すると、Weifeng Qiu の研究は、緩和された正則性要件を持つ一般的な多面体領域上で強ノルムにおいて最適収束を達成する、堅牢で統合された有限要素枠組みを確立することにより、非発散型方程式に対する数値 PDE の最先端を進展させた。