[论文解读] Maximum bound principles for a class of semilinear parabolic equations and exponential time differencing schemes

Qiang Du, Lili Ju+2 more|arXiv (Cornell University)|May 23, 2020

Advanced Mathematical Modeling in Engineering参考文献 93被引用 12

一句话总结

本文为半线性抛物方程中的最大界原理（MBP）建立了理论框架，证明当初始条件和边界条件满足此类有界性时，解在绝对值上保持一致有界。本文提出了无条件时间保持MBP的指数时间差分（ETD）格式，确保了稳定性和高阶精度，且对标量、向量和矩阵值系统均给出了严格的误差估计和能量稳定性证明。

ABSTRACT

The ubiquity of semilinear parabolic equations has been illustrated in their numerous applications ranging from physics, biology, to materials and social sciences. In this paper, we consider a practically desirable property for a class of semilinear parabolic equations of the abstract form $u_t=\mathcal{L}u+f[u]$ with $\mathcal{L}$ being a linear dissipative operator and $f$ being a nonlinear operator in space, namely a time-invariant maximum bound principle, in the sense that the time-dependent solution $u$ preserves for all time a uniform pointwise bound in absolute value imposed by its initial and boundary conditions. We first study an analytical framework for some sufficient conditions on $\mathcal{L}$ and $f$ that lead to such a maximum bound principle for the time-continuous dynamic system of infinite or finite dimensions. Then, we utilize a suitable exponential time differencing approach with a properly chosen generator of contraction semigroup to develop first- and second-order accurate temporal discretization schemes, that satisfy the maximum bound principle unconditionally in the time-discrete setting. Error estimates of the proposed schemes are derived along with their energy stability. Extensions to vector- and matrix-valued systems are also discussed. We demonstrate that the abstract framework and analysis techniques developed here offer an effective and unified approach to study the maximum bound principle of the abstract evolution equation that cover a wide variety of well-known models and their numerical discretization schemes. Some numerical experiments are also carried out to verify the theoretical results.

研究动机与目标

• 建立线性耗散算子 𝒮 和非线性算子 f 的充分解析条件，以保证抽象半线性抛物方程的时间不变最大界原理（MBP）成立。
• 构造一阶和二阶指数时间差分（ETD）格式，实现无条件时间保持MBP，且不依赖时间步长大小。
• 为所提出的ETD格式在标量和向量/矩阵值情形下，推导严格的误差估计并证明能量稳定性。
• 将抽象框架扩展至具有狄利克雷、诺伊曼和周期边界条件的系统，并分析不同矩阵范数（2-范数与Frobenius范数）下的MBP行为。
• 探讨在积分因子（IF）格式中无条件保持MBP的可行性，并与ETD方法进行比较。

提出的方法

• 基于收缩半群理论，推导线性算子 𝒮 和非线性算子 f 的解析条件，以确保解在所有时间保持有界。
• 利用收缩半群的生成元构造ETD格式，以在时间离散设置下无条件保持MBP。
• 应用指数积分方法，实现一阶和二阶时间精度，同时在无时间步长限制下保持MBP。
• 通过能量稳定性分析和误差估计，验证ETD格式在各种边界条件下的收敛性与鲁棒性。
• 将框架扩展至向量和矩阵值系统，分析 ℓ²、ℓ∞ 以及矩阵2-范数和Frobenius范数下的MBP行为。
• 将ETD与积分因子（IF）格式进行比较，表明尽管IF格式在时间步长约束下可有条件保持MBP，但ETD格式实现了无条件MBP保持。

实验结果

研究问题

• RQ1 在何种条件下，线性算子 𝒮 和非线性算子 f 使得半线性抛物方程在时间连续情形下满足最大界原理（MBP）？
• RQ2 能否构造出无条件时间保持MBP的指数时间差分（ETD）格式，且不依赖时间步长大小？
• RQ3 所提出的ETD格式在标量、向量和矩阵值系统中的误差估计与能量稳定性特性为何？
• RQ4 在矩阵值系统中，MBP在不同矩阵范数下（特别是2-范数与Frobenius范数）的行为如何？
• RQ5 积分因子（IF）格式能否无条件保持MBP？若不能，其在何种条件下可实现条件性保持？

主要发现

• 本文建立了 𝒮 和 f 的充分条件，证明只要初始和边界数据满足有界性，抽象半线性抛物方程的解在所有时间均保持绝对值有界。
• 构造了一阶和二阶ETD格式，实现无条件MBP保持，即对任意时间步长均保持有界性，无稳定性限制。
• 推导了理论误差估计，表明在解和数据具有适当正则性假设下，ETD格式达到最优收敛速率。
• 证明了ETD格式的能量稳定性，确保离散能量不随时间增长，从而支持长时间计算的数值可靠性。
• 对于矩阵值系统，ETD格式无条件保持了关于2-范数的MBP，但关于Frobenius范数的MBP保持问题仍为开放问题，原因在于难以控制非线性项的上界。
• 尽管IF格式可在时间步长约束下条件性保持MBP，但本文指出通过引入适当的稳定器，未来可能实现IF格式的无条件MBP保持，为后续研究指明方向。