解一维双曲守恒律方程和抛物方程的间断有限元法的逐点和-Math.PDFVIP

中国科学 数学 年 第 卷 第 期 自然科学基金项目进展专栏 解一维双曲守恒律方程和抛物方程的间断 有限元法的逐点和区间平均值误差估计 献给林群教授 华诞 ∗ 曹外香 张智民 北京计算科学研究中心, 北京 100094; Department of Mathematics, Wayne State University, Detroit, MI 48202, USA E-mail: wxcao@csrc.ac.cn, zzhang@math.wayne.edu 收稿日期: 2014-10-12; 接受日期: 2015-03-18; * 通信作者 国家自然科学基金(批准号:和 资助项目 摘要 本文研究求解一维双曲守恒律方程和抛物方程的间断有限元法的超收敛性质. 具体来说, 对于 双曲守恒律方程和抛物方程, 当分别选择迎风和交替的数值流量时, 本文证明在合适的初始化条件下, 间断有限元解在迎风点上(双曲方程) 或数值迹在节点上(抛物方程) 的逐点误差和区间平均值误差均 以2k + 1 阶的速度收敛, 其中k 是间断有限元空间多项式的次数. 这个结果是对Cao 等人(2014) 以 及Cao 和Zhang (2014) 的超收敛结果的进一步改进和完善, 即把逐点误差的收敛阶由原来的2k + 1 2 提高到2k + 1, 而区间平均值误差的收敛阶则从2k 提高到2k + 1, 从而与 Cao 等人 (2014)、Cao 和 Zhang (2014) 以及Yang 和Shu (2012) 的数值试验结果完全吻合. 关键词 超收敛 间断有限元法 双曲方程 抛物方程 区间平均值 逐点误差估计 主题分类 引言 间断有限元或间断Galerkin (discontinuous Galerkin, DG) 方法是一类使用完全不连续的分片多 项式空间作为解空间和检验函数空间的有限元方法. 它最早是 1973 年 Reed 和 Hill [1] 在研究中子 输运方程时提出来的. 随后, Cockburn 和 Shu 采用间断有限元方法进行空间离散去求解与时间有关 的非线性方程(参见文献 [2–6]), 进一步推动了这个方法的发展. 局部间断有限元 (LDG) 方法是间断 有限元法的一个延伸 它主要用于求解各种含高阶导数的方程 由于数值逼近在区间边界上的不连续 , ( 性, 间断有限元法不便于直接求解高阶方程). 局部间断有限元方法的基本思想就是把含高阶导数的 方程写成只含一阶导数的方程组形式, 然后再用间断有限元方法一一求解. 为了确保数值格式的稳定 性, 选取合适的数值流量是至关重要的. Cockburn 和 Shu [7] 首次用局部间断有限元方法处理了与时 间有关的对流扩散方程. 此后, 局部间断有限元方法被应用到各种各样的高阶方程, 如含三阶导数的 Korteweg-de Vries (KdV) 方程 (参见文献 [8])、含时间项的双调和方程、完全非线性的k(n, n, n) 方程 和Cahn-Hilliard 类型的方程等 (参见文献 [9–12]). 英文引用格式 Cao W X, Zhang Z M. Point-wise and cell average error estimates of the DG and LDG methods for 1D hyperbolic and parabolic equations (in Chinese). Sci Sin Math, 2015, 45: 1115–1132, doi: 10.1360/N012014-00214 曹外香等 解一维双曲守恒律方程和抛物方程的间断有限元法的逐点和区间平均值误差估计 间断有限元法解偏微分方程的超收敛性质也是近年来学者们比较感兴趣的研究领域, 这方面的内 容可参见文献 [13–20]. 最近, Yang 和Shu [21] 在研究一维双