$L^p$ Estimates for Numerical Approximation of Convex Hamilton-Jacobi Equations

技術要約：凸なハミルトン–ヤコビ方程式の数値近似における $L^p$ 誤差評価

問題設定
本論文は、$d$ 次元の平坦なトーラス $\mathbb{T}^d$ 上における凸なハミルトン–ヤコビ（HJ）方程式の数値近似を取り扱う。対象とする方程式は以下の通りである：
$$\begin{cases} \partial_t u + H(Du) = 0 & \text{in } Q := \mathbb{T}^d \times (0, T), \\ u(x, 0) = g(x) & \text{on } \mathbb{T}^d, \end{cases}$$
ここで、ハミルトニアン $H$ は滑らかで凸かつ強圧的であり、初期データ $g$ はリプシッツ連続かつ半凹（semiconcave）である。著者らは、単調数値スキームの2つのクラス、すなわち有限差分法とセミ・ラグランジュ法に焦点を当てている。古典的な結果（Crandall–Lions や Barles–Souganidis など）は、凸性の仮定を置かない一般的な単調スキームに対して $1/2$ 次の $L^\infty$ 収束率を確立しているが、本研究では、ハミルトニアンの凸性と解の半凹性を利用することで、より鋭い $L^p$ ノルム（$1 \le p < +\infty$）における誤差評価を確立することを目的としている。

手法
核となる手法は、随伴法（adjoint method）と半凹性評価の組み合わせに基づいている。このアプローチは、標準的な粘性解の技法（通常は変数の二重化法を用いる）とは異なり、スカラ保存則のために開発され、HJ方程式へと適応された双対性の議論を利用している。

1. 随伴フレームワーク： 著者らは、非線形HJ方程式を線形化された随伴方程式と結合させる。有限差分スキームの場合、誤差 $z = u^\varepsilon - u$ は、切断誤差を表すソース項を持つ輸送型の方程式を満たす。この方程式に、終端データ $\text{sgn}(z)$ を持つ後退随伴方程式の解 $\rho$ を乗じることで、$L^1$ ノルムにおける誤差は、切断誤差と随伴解の積分として表現される。
2. 半凹性の保存： 本解析の重要な構成要素は、数値スキームが初期データの半凹性を保存することを示すことである。
   • 有限差分スキームの場合： 著者らは、連続空間・離散時間型の補助スキームを構築する。彼らは、単調更新演算子が連続的な意味での半凹性を保存することを示し、これにより離散的な半凹性の再構成を必要とせずに済むようにしている。
   • セミ・ラグランジュ法の場合： 解析は、離散値関数の動的計画法の表現を利用する。証明は、値関数のリプシッツ連続性と半凹性の性質、および（最適制御の一意性に結びついた）ほとんど至る所での微分可能性を利用している。
3. 誤差分解： 切断誤差は以下の2つの部分に分解される：
   • 時間補間による誤差（離散的な時間ステップ間）。
   • 空間的一貫性誤差（離散勾配と連続勾配の差）。
     保存された半凹性を用いることで、数値解の二階微分に関する $L^1$ ノルムを抑えることができ、これにより空間的一貫性誤差を制御することが可能となる。
4. 補間： $L^1$ 誤差評価が 1 次のオーダーで確立された後、著者らは $L^1$ と古典的な $L^\infty$ 評価（次数 $1/2$）の間の補間不等式を用いて、すべての $1 \le p < +\infty$ に対する $L^p$ 誤差評価を導出する。

主要な貢献および結果

• 有限差分スキーム：

  • 著者らは、一貫性のある凸な数値ハミルトニアンを用いた単調有限差分スキームに対し、$O(h + \Delta t)$ の $L^1$ 誤差評価を証明した。
  • 標準的な $L^\infty$ 評価 $O((h+\Delta t)^{1/2})$ との補間により、以下の $L^p$ 誤差評価を導出している：
    $$\|u^\varepsilon_n - u(\cdot, n\Delta t)\|_{L^p} \le C (h + \Delta t)^{\frac{1}{2} + \frac{1}{2p}}.$$
  • この結果は、輸送–射影構造に依存することなく、多次元設定およびより広範なクラスの凸な数値ハミルトニアンへの適用を可能にし、先行研究（Cagnetti–Gomes–Tran など）を拡張するものである。

• セミ・ラグランジュ法：

  • $H(x, p) = H_0(p) - V(x)$ の形式のハミルトニアンを持つセミ・ラグランジュ法について、時間離散スキームに対する $L^1$ 誤差評価 $O(\Delta t)$ を確立している。
  • 完全離散スキーム（空間グリッド補間を含む）については、$L^1$ 誤差評価は $O(\Delta t + h/\Delta t)$ である。
  • 対応する $L^p$ 収束率は以下のように導かれる：
    $$\| \tilde{u}^\varepsilon_n - u(\cdot, n\Delta t) \|_{L^p} \le C \left( \Delta t + \frac{h}{\Delta t} \right)^{\frac{1}{p}} \left( \Delta t^{\frac{1}{2}} + \frac{h}{\Delta t} \right)^{1-\frac{1}{p}}.$$

意義および主張
本論文は、凸なHJ方程式に対する $L^p$ 誤差評価の体系的な導出を提供し、一般的な単調スキームで知られている標準的な $L^\infty$ 速度を改善することを主張している。その意義は以下の点にある：

1. 収束速度の鋭鋭化： 凸性と半凹性の仮定の下で、$L^1$（したがって $p>1$ における $L^p$）における収束速度が 1 次であることを示し、これは $L^\infty$ で見られる典型的な $1/2$ 次よりも優れている。
2. 手法の適応： 保存則や消散粘性近似に伝統的に用いられてきた随伴法を、局所的な単調スキーム（Crandall–Lions 型）やセミ・ラグランジュ・スキームに直接適用して分析することに成功した。
3. 一般化： 特定の 1 次元または輸送–射影構造に限定されることなく、多次元トーラスおよび広範なクラスの凸な数値ハミルトニアンへの適用性を拡張した。

著者らは、自らの結果が標準的な仮定（$H$ の滑らかさ、凸性、強圧性；$g$ のリプシッツ連続性と半凹性）に依拠していることを述べ、控えめなトーンを維持しており、$L^1$ メカニズムは半凹性–双対性の原理の直接的な適用であり、新規性は検討された特定の数値スキームに対するこれらの仮定の検証にあるとしている。