[论文解读] Speed of convergence of Chernoff approximations for two model examples: heat equation and transport equation

Pavel S. Prudnikov|arXiv (Cornell University)|Dec 17, 2020

advanced mathematical theories参考文献 22被引用 5

一句话总结

本文通过解析估计和数值实验，分析了切尔诺夫逼近在两个模型偏微分方程——热方程和输运方程——上的收敛速度。研究表明，收敛速率取决于初始条件的正则性，当初始数据不在生成元定义域之外时收敛更慢，并为不同切尔诺夫函数和初始条件提供了定量误差估计和图形证据。

ABSTRACT

Paul Chernoff in 1968 proposed his approach to approximations of one-parameter operator semigroups while trying to give a rigorous mathematical meaning to Feynman's path integral formulation of quantum mechanics. In early 2000's Oleg Smolyanov noticed that Chernoff's theorem may be used to obtain approximations to solutions of initial-value problems for linear partial differential equations (LPDEs) of evolution type with variable coefficients, including parabolic equations, Schrödinger equation, and some other. Chernoff expressions are explicit formulas containing variable coefficients of LPDE and the initial condition, hence they can be used as a numerical method for solving LPDEs. However, the speed of convergence of such approximations at the present time is understudied which makes it risky to employ this class of numerical methods. In the present paper we take two equations with known solutions (heat equation and transport equation) and study both analytically and numerically the speed of decay of the norm of the difference between Chernoff approximations and exact solutions. We also provide graphical illustrations of convergence and its rate. These model examples, being relatively simple, allow to demonstrate general properties of Chernoff approximations. The observations obtained build a base for the future employment of the approach based on Chernoff's theorem to the problem of construction of new numerical methods for solving initial-value problem for parabolic LPDEs with variable coefficients.

研究动机与目标

• 研究具有变系数的线性抛物型PDE的切尔诺夫逼近收敛速度。
• 评估初始条件正则性对收敛速率的影响。
• 为模型方程提供收敛速率的解析估计与数值验证。
• 为基于切尔诺夫定理发展适用于变系数抛物型PDE的新数值方法奠定基础。

提出的方法

• 推导切尔诺夫逼近应用于输运方程和热方程时的精确解与误差界。
• 应用切尔诺夫定理构造生成柯西问题解逼近的算子值函数。
• 利用$C_0$-半群理论定义演化方程的解并分析收敛性。
• 使用MATLAB对不同初始条件和切尔诺夫函数进行数值实验。
• 采用对数-对数图与非线性回归估计收敛阶数。
• 比较不同初始条件（包括$u_0(x) = \sin(x)$和$u_0(x) = \exp(-|x|)$）下的收敛速率。

实验结果

研究问题

• RQ1 对于具有光滑与非光滑初始数据的热方程，切尔诺夫逼近的收敛速率是多少？
• RQ2 切尔诺夫函数的选择如何影响输运方程的收敛速度？
• RQ3 当初始条件不在生成元的定义域中时，收敛速率是否会退化？
• RQ4 数值实验能否验证不同初始条件下收敛速度的理论估计？
• RQ5 基于观测到的收敛速率，给定初始条件属于哪个逼近子空间的阶？

主要发现

• 对于输运方程，切尔诺夫逼近的收敛速率为$O(1/n)$，且该速率与初始条件的正则性无关。
• 对于初始条件$u_0(x) = \sin(x)$的热方程，收敛速率为$O(1/n)$，与理论估计一致。
• 当$u_0(x) = \exp(-|x|)$（该函数不在生成元的定义域中）时，收敛速率仍为$O(1/n)$，表明即使对于正则性较低的初始数据，逼近子空间的阶仍被保持。
• 数值实验表明，不同切尔诺夫函数（包括$\gamma = 1/2$、$1/3$和$1/6$）下的收敛速率具有鲁棒性，均呈现一致的$O(1/n)$行为。
• 利用对数-对数图的图形分析确认了$O(1/n)$的收敛速率，并可通过回归方法估计收敛阶数。
• 研究表明，即使对于非光滑初始条件，收敛行为仍与理论预期一致，支持切尔诺夫逼近在实际数值方案中的应用。