A Positivity-Preserving Expectation Scheme for Hamilton--Jacobi--Bellman Equations with Oblique Robin Boundary Conditions

技術要約：斜めロビン境界条件を持つHJB方程式のための正値性を保持する期待値スキーム

問題設定
本論文は、斜めロビン境界条件（oblique Robin boundary conditions）に従うハミルトン・ヤコビ・ベルマン（HJB）方程式の数値解法を扱う。問題は、領域 $O_T = [0, T) \times O$ 上の後退形式で定式化される：
$$\begin{cases} -\partial_t u + H(t, x, Du, D^2u, u) = 0, & (t, x) \in O_T, \\ L(t, x, Du, u) = 0, & (t, x) \in [0, T) \times \partial O, \\ u(T, x) = \Psi(x), & x \in O, \end{cases}$$
ここで、ハミルトニアン $H$ は異方性または退化的な拡散（$\sigma$ によって制御される）およびドリフト（$\mu$）を含む。境界作用素 $L$ は、斜微分項とゼロ次項 $k(t,x,b)u$（ロビン項）およびソース $g$ を含む。標準的な有限差分法や有限要素法は、厳密な対角優位性や特定のメッシュ条件が満たされない限り、異方性拡散において非負性（単調性）を保持できないことが多く、これは退化的なケースにおいて問題となる。

手法
著者らは、代数的なステンシル分解ではなく、ファインマン・カッツ表現から導出された条件付き期待値フレームワークに基づく完全離散スキームを提案している。

1. 内部ノード（反射ファインマン・カッツ恒等式）：

   • スキームは、$m$ 次元のラデマッハー・ベクトルによって駆動される 1 ステップ弱オイラー近似を用いて拡散過程を離散化し、$P=2^m$ 個の等確率な枝（branches）を生成する。
   • 鏡面反射（Mirror Reflection）： 試行点 $X^*$ が領域の外側に落ちた場合、境界上のその斜投影 $p_{\gamma_b}(X^*)$ に関して反射させる。反射点 $\tilde{X}$ は $\tilde{X} = X^* - 2d_{\gamma_b}(X^*)\gamma_b$ となる。
   • 境界局所時間（Boundary Local Time）： 往復距離 $e_d = 2d_{\gamma_b}(X^*)$ が、境界局所時間の離散増分として機能する。
   • 減衰（Attenuation）： ロビン係数 $k$ は、減衰因子 $e^{-k e_d}$（継続値に対して）および $e_d e^{-k e_d/2} g$（境界ソースに対して）を通じてのみスキームに寄与する。決定的なことに、枝の確率は正確に $1/P$ のままであり、対角優位性の条件を必要とせずに重みの正値性が保持される。
   • 補間： 反射点における値は、$P1$（ piecewise linear）補間を用いて空間メッシュ上で回収される。

2. 境界ノード（Same-Level Closure）：

   • 内部ノードとは異なり、境界ノードは空間オフセット $\ell_h \asymp h$ にわたる別個の same-level 方程式を用いて更新される。
   • 境界条件は以下のように離散化される：
     $$\frac{U_i - I_h[U](x_i - \ell_h \gamma(x_i))}{\ell_h} + k(t_n, x_i)U_i - g(t_n, x_i) = 0.$$
   • この系は一般に暗黙的である。固定された線形データに対しては、疎で非特異なM行列を与える。制御問題の場合、境界写像は縮小写像となり、解の存在を保証する。

3. 制御の取り扱い：

   • スキームは、枝の平均を形成した後に、制御集合 $A \times B$ に対する下限（または上限、定式化による）を取ることによって、制御を扱う。制御は各ノードおよび各タイムステップで新たに選択される。

主な貢献

1. 対角優位性なしでの正値性： 本スキームは、確率的な構成（等確率および非負の減衰因子）から直接、非負の係数を生成する。これにより、コンパクトなステンシルによる対角優位性の条件や、$\Delta t$ と $h$ の間のCFL型の関係を課すことなく、異方性および退化的な拡散を扱うことができる。
2. 幾何学的ロビン処理： 本手法は、内側へのシフト係数を必要とせずに（従来のセミ・ラグラン的手法とは異なり）、鏡面反射と離散局所時間を用いて斜めロビン条件を処理する。ロビンデータは減衰を通じてのみ結合し、枝の確率は変化しない。
3. 単調性と安定性： スキームは順序保存型（単調）であり、$L^\infty$ 安定である。境界閉鎖は縮小写像を形成し、内部更新は明示的である。
4. 収束理論： 精密化条件 $\Delta t, h \to 0$ かつ $h^2/\Delta t \to 0$ の下で、本スキームは半緩和極限引数（Barles–Souganidis フレームワーク）を介して一意な粘性解に収束する。
5. 誤差評価： 均一な $C^{1,1}$ 境界残差を持つ古典的な解 $u \in C^3_b$ に対して、誤差界は以下のように導かれる：
   $$\max |U^n_i - u(t_n, x_i)| \le C \left( \Delta t^{1/2} + \frac{h^2}{\Delta t} \right).$$
   項のバランスをとるために $\Delta t \asymp h^{4/3}$ とすると $O(h^{2/3})$ のオーダーが得られ、標準的な選択 $\Delta t \asymp h$ では $O(h^{1/2})$ となる。

数値結果
論文では、単位円盤および混合領域の矩形における4つの2次元実験を報告している：

• 比較 ($k=0$)： $k=0$ の問題について、較正された参照実装（Calzolaら）と比較している。提案手法は、調整可能な内側オフセットパラメータを必要とせずに、競争力のある精度（次数 $\approx 0.96$）を達成し、細かいメッシュ上でもよりバランスの取れた挙動を示す。
• 可変ロビン ($k \neq 0$)： 可変ロビン係数を用いたテストにより、スキームの安定性が確認された。「バランスの取れた」タイムステップ・スケジュール（$\Delta t \propto h^{4/3}$）は、$\Delta t = h$ の場合に比べて誤差を約55%減少させ、理論的な $O(h^{2/3})$ の境界よりも鋭い、観測された次数 $1.242\text{--}1.261$（前漸近的）を得た。
• 角度感度： 角度制御数 $K$ に対する誤差の減衰は、$O(K^{-2})$ と一致している。
• 非滑らかな/混合ケース： 円錐型の粘性解（退化拡散）および混合ロビン–ディリクレ問題の実験では、単調な誤差減衰を示したが、収束率は低く、これらは滑らかな領域の定理の対象外である。

意義と主張
著者らは、本研究が斜めロビン境界条件を含めることで、彼らの以前の期待値ベースのフレームワークを拡張したものであると主張しており、これは界面現象（熱対流、電気化学など）を伴う問題において重要なステップであるとしている。主な意義は、正値性の要件をメッシュ制約から切り離したことにある。確率的な反射メカニズムを用いることで、本スキームは、標準的な有限差分法で通常必要とされる制限的な幾何学的条件を課すことなく、退化的な異方性拡散に対して単調性と安定性を維持できる。

論文では、滑らかなテストにおいて観測された前漸近的な次数が理論的境界を超えることもあると控えめに述べているが、理論的な保証は粘性解のフレームワークと特定の正則性仮定に限定されている。本手法は、特に斜め境界条件や退化的な拡散が存在する場合において、標準的なステンシルが失敗するHJB方程式に対する、堅牢で正値性を保持する代替案として提示されている。