[论文解读] Stochastic Solution of Elliptic and Parabolic Boundary Value Problems for the Spectral Fractional Laplacian

Mamikon Gulian, Guofei Pang|arXiv (Cornell University)|Dec 4, 2018

Advanced Mathematical Modeling in Engineering参考文献 23被引用 5

一句话总结

本文提出了针对带有非零狄利克雷边界条件的谱分数阶拉普拉斯算子 $(-\Delta_{\Omega,g})^{\alpha/2}$ 的椭圆与抛物型边值问题的新型随机解公式（费曼-卡茨型）。该方法利用次级停止的布朗运动，通过蒙特卡罗模拟提供了一套概率框架，实现了在有界区域内的高效、局部并行计算，并具备严格的收敛性与正则性分析。

ABSTRACT

We prove and implement stochastic solution (or Feynman-Kac) formulas for boundary value problems involving the spectral fractional Laplacian with nonzero Dirichlet boundary condition. The main tools used in the proofs are the abstract Cauchy problem for Feller semigroups together with Balakrishnan's theory of fractional powers. We show the the spectral fractional Laplacian with nonzero Dirichlet boundary conditions is the generator of an appropriate Feller semigroup for subordinate stopped Brownian motion, and obtain a stochastic solution formula for the fractional heat equation in a bounded domain. We then obtain precise regularity and steady-state convergence properties of the parabolic problem using the eigenfunction expansion of the classical solution, which leads to estimates for the survival probability of subordinate stopped Brownian motion. These results allow us to take the long-time limit in the parabolic formula to establish a stochastic solution formula for the Dirichlet boundary value problem. These stochastic solution formulas for the spectral fractional Laplacian with nonzero boundary conditions are novel, and allow for efficient, embarrassingly parallel local solution of the above boundary value problems. We discuss the discretization of these formulas, and verify them in dimensions two and three with benchmark examples. We study the effect of the number of path samples and the time step size for path discretization on the accuracy of the solution.

研究动机与目标

• 开发并验证带有非零狄利克雷边界条件的谱分数阶拉普拉斯算子的随机解公式。
• 通过费勒半群与下放理论，建立谱分数阶拉普拉斯算子的概率解释。
• 通过基于路径的蒙特卡罗采样，实现分数阶PDE的高效、天然并行数值求解。
• 利用特征函数展开分析分数阶热方程的正则性与稳态收敛性。
• 通过二维与三维区域中的数值基准测试验证所提公式的准确性。

提出的方法

• 基于巴拉克里斯南的分数阶幂理论，将谱分数阶拉普拉斯算子表示为费勒半群的生成元。
• 构造过程 $X^{\Omega,\alpha}_{t} = \sqrt{2}B_{T_{\alpha/2}(t) \wedge \tau_{\Omega}}$，即布朗运动经 $\alpha/2$-稳定下放过程时间变换，并在离开区域时停止。
• 建立非齐次分数阶柯西问题的费曼-卡茨公式：$u(t,x) = \mathbb{E}_x\left[ f(X_t) \chi_{\tau_\Omega > T_{\alpha/2}(t)} + g(X_t) \chi_{\tau_\Omega \leq T_{\alpha/2}(t)} + \int_0^t r(t-s,X_s) \chi_{\tau_\Omega > T_{\alpha/2}(s)} ds \right]$。
• 通过取 $t \to \infty$ 将解推广至稳态狄利克雷问题，得到 $u(x) = \mathbb{E}_x\left[ g(X_{T^{-1}_{\alpha/2}(\tau_\Omega)}) \right] + \mathbb{E}_x\left[ \int_0^{T^{-1}_{\alpha/2}(\tau_\Omega)} r(X_s) ds \right]$。
• 利用特征函数展开推导次级停止过程的正则性与存活概率估计。
• 通过在二维与三维区域中基于路径采样的数值实现与验证，分析路径数量与时间步长对收敛性的影响。

实验结果

研究问题

• RQ1 如何通过随机过程表示带有非零狄利克雷边界条件的谱分数阶拉普拉斯算子？
• RQ2 次级停止的布朗运动与分数阶抛物型及椭圆型PDE解之间有何关联？
• RQ3 次级停止过程的存活概率如何影响解的正则性与收敛性？
• RQ4 所提出的随机解公式能否实现保证收敛性与并行可扩展性的数值实现？
• RQ5 蒙特卡罗离散化在路径数量与时间步长方面的精度与收敛行为如何？

主要发现

• 证明了谱分数阶拉普拉斯算子 $(-\Delta_{\Omega,g})^{\alpha/2}$ 是与次级停止布朗运动相关联的费勒半群的生成元。
• 基于抽象柯西问题与下放理论，推导并证明了分数阶热方程的随机解公式。
• 通过在抛物型公式中取 $t \to \infty$，获得狄利克雷问题的稳态解，其形式为涉及过程退出时间的路径期望。
• 二维与三维区域的数值实验表明，随着路径数量增加与时间步长减小，解呈现收敛性，验证了方法的准确性。
• 该方法实现了局部并行、无需网格的分数阶PDE求解，在高维下表现出良好的可扩展性。
• 通过特征函数展开推导出正则性与存活概率估计，证实了该解框架的适定性。