中文

时间分数阶 Fokker-Planck 方程的二阶精度数值格式

数值分析 2021-06-29 v1 数值分析

摘要

研究了一种求解阶为 α(0,1)\alpha \in (0, 1) 且含一般驱动力的分数阶 Fokker-Planck 方程的时间步进 L1L1 格式。对 α(1/2,1)\alpha\in(1/2,1) 通过一种新颖简洁的方法获得了半离散解的稳定性界。我们的稳定性估计是 α\alpha 鲁棒的,在于其对 α\alpha 趋近于 11 的极限情形(此时模型退化为经典 Fokker-Planck 方程)依然有效,该极限具有实际重要性。在误差分析方面,我们对 α(1/2,1)\alpha\in(1/2,1) 获得了最优的二阶精度估计。采用时间分级网格以补偿连续解在原点附近的奇异行为。L1L1 格式与标准空间 Galerkin 有限元离散相结合,从数值上支撑我们的理论贡献。我们采用所得全离散可计算数值格式进行了一些数值测试。这些测试表明所施加的时间分级网格假设可进一步放宽,且我们观察到即便对 α(0,1/2]\alpha\in(0,1/2] 的情形(即超出理论覆盖的范围)也具有二阶精度。

关键词

引用

@article{arxiv.2106.14146,
  title  = {A second-order accurate numerical scheme for a time-fractional Fokker-Planck equation},
  author = {Kassem Mustapha and Omar M. Knio and Olivier P. Le Maître},
  journal= {arXiv preprint arXiv:2106.14146},
  year   = {2021}
}