[论文解读] Error analysis of linearized semi-implicit Galerkin finite element methods for nonlinear parabolic equations

Buyang Li, Weiwei Sun|arXiv (Cornell University)|Aug 23, 2012

Advanced Numerical Methods in Computational Mathematics参考文献 19被引用 108

一句话总结

本文提出了一种针对非线性抛物方程（特别是3D焦耳加热系统）的线性化半隐式伽辽金有限元方法的新型误差分析框架。通过将误差分解为时间与空间分量，并利用时间离散系统，作者在无任何时间步长限制的条件下，建立了无条件的最优 $L^2$ 和 $H^1$ 误差估计，克服了先前工作中为保证强范数有界性而必须施加严格时间步长条件的关键局限。

ABSTRACT

This paper is concerned with the time-step condition of commonly-used linearized semi-implicit schemes for nonlinear parabolic PDEs with Galerkin finite element approximations. In particular, we study the time-dependent nonlinear Joule heating equations. We present optimal error estimates of the semi-implicit Euler scheme in both the $L^2$ norm and the $H^1$ norm without any time-step restriction. Theoretical analysis is based on a new splitting of the error and precise analysis of a corresponding time-discrete system. The method used in this paper can be applied to more general nonlinear parabolic systems and many other linearized (semi)-implicit time discretizations for which previous works often require certain restriction on the time-step size $τ$.

研究动机与目标

• 解决在非线性抛物PDE的线性化半隐式格式误差分析中通常施加的时间步长限制问题。
• 克服依赖于逆不等式和 $L^∞$-范数有界性假设的归纳法，这些假设需要时间步长约束。
• 在 $L^2$ 和 $H^1$ 范数下建立无任何时间步长大小 $\tau$ 限制的最优误差估计。
• 开发一个适用于更广泛非线性抛物系统（超越焦耳加热模型）的一般性框架。

提出的方法

• 通过时间离散抛物系统，引入一种新的数值误差分解方法，将其划分为时间分量与空间分量。
• 利用仅依赖于空间网格尺寸 $h$ 的界，分析时间离散系统的伽辽金有限元误差，而不依赖于 $\tau$。
• 通过归纳法与逆不等式控制误差的 $L^\infty$-范数，无需先验有界性假设，从而避免时间步长限制。
• 应用格朗沃尔不等式，推导出 $L^2$ 与 $H^1$ 范数下的最优误差估计。
• 采用 $W^{1,p}$ 估计与插值误差界，控制涉及电导率 $\sigma(u)$ 的非线性项。
• 通过归纳法与与网格相关的逆不等式，建立离散解在强范数下的有界性，从而实现无条件误差控制。

实验结果

研究问题

• RQ1 能否在不施加时间步长限制的条件下，为非线性抛物方程的线性化半隐式伽辽金有限元格式导出最优误差估计？
• RQ2 是否可以通过一种新的误差分解策略，消除误差分析中对数值解 $L^\infty$-范数有界的依赖？
• RQ3 如何为3D非线性焦耳加热系统无条件地获得 $H^1$-范数误差估计？
• RQ4 所提出的框架能否推广至其他非线性抛物系统及更高阶时间格式？
• RQ5 在新分析框架下，误差对空间网格尺寸 $h$ 与时间步长 $\tau$ 的精确依赖关系是什么？

主要发现

• 无条件地实现了最优 $L^2$ 误差估计，阶数为 $O(h^{3/2})$，且对时间步长大小 $\tau$ 无任何限制。
• 无条件地建立了最优 $H^1$ 误差估计，阶数为 $O(h)$，且无时间步长限制。
• 通过归纳法与格朗沃尔不等式，证明了 $L^2$ 误差界 $\|e_h^n\|_{L^2} \leq Ch^{3/2}$，其仅依赖于 $h$-相关界。
• 通过使用时间离散系统与逆不等式，避免了对误差 $L^\infty$-范数的控制，从而消除了对时间步长条件的需求。
• 该分析框架具有通用性，可推广至其他非线性抛物系统及更高阶时间格式。
• 其核心创新在于误差分解策略以及利用时间离散系统以解耦时间与空间误差贡献，从而实现无条件稳定性与最优性。