A Residual-Certified Product-Integration Method for Nonlinear Caputo–Volterra Equations

技術要約：非線形Caputo–Volterra方程式に対する残差保証付き積分法

問題設定
本論文は、非線形Caputo初期値問題の数値解法を扱っており、これはAbel型の核を持つ弱特異なVolterra積分方程式として再定式化できる。具体的には、以下の式に焦点を当てる：
$$C D^\alpha_{0+} u(t) = F(t, u(t)), \quad u(0) = u_0, \quad 0 < \alpha < 1$$
ここで、$F$は局所リプシッツ条件を満たす。主な課題は、離散化から生じる非線形代数系の求解である。特に、グローバルなリプシッツ定数が大きい場合（標準的なバナッハの不動点定理が必要とする閾値を超える場合）、および解が初期特異性を示す場合に困難が生じる。本研究は、高次の積分則を因果的に実装しつつ、不正確な非線形解に対して計算可能かつ厳密な境界を提供できるかどうかを検討している。

手法
著者らは、非線型項 $F(t, u(t))$ の非一様メッシュ上での区分線形補間を用いた暗黙的積分スキームを提案している。

1. 離散化: Abel核を線形補間関数に対して正確に積分する。これにより、現在のノード $U_i$ が暗黙的に現れる下三角非線形系が得られる。
2. 可解条件: 純粋な因果的左端点則とは異なり、現ノードの項は対角係数 $\kappa_i L_i$ を導入する。本手法は、メッシュ局所的な可解性テストに依存している。もし $\mu_h = \max_i (\kappa_i L_i) < 1$ であれば、離散系は一意の解を持つ。この条件は、標準的なグローバル条件 $S = L^* T^\alpha / \Gamma(\alpha+1) < 1$ よりも緩い。
3. 反復ソルバー: グローバルなPicard反復と、ノードごとの（逐次的な）Picard反復の両方が、グローバルな条件 $S > 1$ であっても、局所的な条件 $\mu_h < 1$ のもとで収束することが証明されている。
4. 残差保証: 主要な構成要素は、成分ごとの代数誤差界の導出である。任意の不正確な離散ベクトル $W$ に対して、誤差 $|W - U_h|$ は $(I - C_h)^{-1} |W - T_h W|$ によって抑えられる。ここで、$C_h$ はリプシッツ定数と求積重みから導かれる非負の下三角行列である。この境界は、単一の前方代入によって計算され、非線形反復の不完全さに起因する誤差と離散化誤差を分離する。
5. 収束解析:
   • 滑らかな解に対して、本手法は十分細かいメッシュ上で2次収束を達成する。
   • べき乗型の初期特異性（$|g''(t)| \le C t^{\sigma-2}$）を持つ解の場合、勾配メッシュ（$t_j = T(j/M)^r$）上での誤差は $O(M^{-\min\{r\sigma, 2\}})$ のオーダーとなる。

主要な貢献
本論文の新規性は、求積則そのもの（これは古典的である）や勾配メッシュの使用（あるいは平滑化変換）にあるのではなく、以下の要素の組み合わせにあると主張している。

• 正確な対角テスト: グローバルなリプシッツ定数が小さくなくても機能する、明示的なメッシュ局所的条件。
• 重み付き縮小引数: グローバルおよびノードごとのPicard反復が、局所的な条件 $\mu_h < 1$ のもとで収束することの証明。
• 単一グリッド残差保証: 任意の不正確なイテレートに対する、計算可能な成分ごとの誤差界。外挿ベースの推定器とは異なり、これは精緻化されたメッシュや漸近的誤差展開を必要としない。
• 誤差の分離: 非線形反復を終了することによって導入される誤差と、離散化誤差を明示的に分離する二段階の推定。
• メッシュに依存しないレゾルベント境界: レゾルベント境界がメッシュサイズに依存しないことを証明するブロック化Volterra推定。これにより収束証明を完結させている。

結果
数値実験により、4つのモデルにおける理論的知見が検証された。

1. 保証の信頼性: 証明された境界は、グローバルなリプシッツ定数 $S$ が大幅に1より大きい場合（例：$S \approx 7.7$）であっても、グローバルPicardイテレートに対して観測された代数誤差を一貫して上回っている。「有効性（effectivity）」、すなわち境界と誤差の比は、反復が進むにつれて向上する。
2. 収束率:
   • 滑らかなベンチマークにおいて、提案手法は2次収束を達成し、分数階Adams PECE法（次数 $\approx 1+\alpha$）や左端点則（次数 1）を凌駕する。
   • 初期特異性を持つ問題では、一様メッシュは減少した収束次数（例：$\approx 1.2$）を示すが、勾配メッシュを用いることで2次のレートを正常に回復できる。
3. 既知の解を持たない問題: 正解を知らない非線形問題に対しても、メッシュの精緻化に伴い反復回数が減少するなど、頑健に動作する。
4. 計算コスト: 暗黙的手法は、非線形求解と重みの構築が必要であるため、明示的なベースライン（Adams PECE、左端点則）よりも計算コストが高い。しかし、非線形ソルバーのための厳密な停止基準を提供し、左端点則よりも優れた精度を提供する。

意義と主張
著者らは、本論文の貢献を、新しい求積則の発見ではなく、既存の積分法に対する理論的およびアルゴリズム的な洗練として控えめに位置づけている。その意義は以下の通りである：

• 局所的なメッシュ特性を利用することで、グローバルな縮小が失敗する場合（$S > 1$）でも、暗黙的積分法の適用範囲を広げたこと。
• 漸近的な仮定や追加のメッシュ精緻化計算に依存しない、非線形ソルバーのための厳密で計算可能な停止基準を提供したこと。
• 代数誤差と離散化誤差を分離する統一されたフレームワークを提供し、非線形分数階Volterra方程式における信頼性の高い誤差推定を可能にしたこと。

著者らは、勾配メッシュ、平滑化変換、またはこれらの方程式に対する積分則を初めて使用したとは主張しておらず、むしろ、対角可解性テスト、重み付き縮小の証明、および単一グリッド残差保証の具体的な統合を主張している。