一阶双曲问题间断有限元的后验误差分析.pdfVIP

2012年5月 计算数学 第34卷第2期 SINICA May,2012 MATHEMATICANUMERICA V01．34，No．2 一阶舣曲问题问断有限元的后验误差分析q) 张 铁 李 铮 (东北大学数学系，沈阳110004) 摘 要 一阶双曲问题的有限元后验误差估计至今没有得到很好的解决．本文对d维区域上一阶双曲 问题的k次间断有限元逼近提出了一种新的后验误差分析方法，进而建立了间断有限元解在DG 范数下(强于L2范数)基于误差余量型的后验误差估计．数值计算验证了本文理论分析的有效性． 本文方法也适用于其他变分问题有限元逼近的后验误差分析， 关键词：一阶双曲问题；间断有限元方法；后验误差分析 MR(2000)主题分类：65N30；65N15 1．引 言 自适应计算是现代科学计算的重要组成部分，后验误差估计则是有限元自适应计算的理 论基础．目前有限元后验误差分析主要有两个途经：一是利用对偶论证技术获得基于误差余 量型的后验误差估计；二是利用后处理有限元解的超收敛性质得到渐近准确的后验误差估计． 对于椭圆问题连续或者间断有限元方法已经有丰富的后验误差估计研究结果11j引．然而对于 一阶双曲问题，由于它的对偶问题不具有所需要的正则性估计，目前也没有后处理有限元解的 超收敛结果，这使得上述两种方法不再有效．此外，由于一阶双曲问题的变分方程不具有日1 正定性，那些适用于椭圆问题H1范数后验误差估计的导数恢复技术也很难得到应用．近年来 虽然也有一些一阶双曲问题间断有限元后验误差估计的研究结果[3--6J，但这些后验误差估计 方法的有效性都存在问题．文[3，4]利用对偶论证技术给出了一种后验误差估计，然而所给出 的误差估计子中包含有对偶问题的未知解z，而对偶问题和原问题一样难于求解．文[3,4]是 采用数值方法对z近似求解，再将近似解代入估计子中．显然这样的后验误差估计子缺乏可 J：存 靠性和有效性．最近文【5]提出了另一种后验误差估计方法，它基于一个关键性的假设【7 元直径减半的剖分瓦／。上的间断有限元解．这个假设是非常牵强的，并且它的真实性无法得 到理论证明．文[6】曾在一种较弱范数给出了后验误差估计，但这种范数不具有实际意义并且 所得结果仅限于一维问题．本文提出了一种新的后验误差分析方法，从而对d维区域上一阶 双曲问题的七次间断有限元逼近建立了DG范数(强于L2范数)下基于误差余量型的，高效 的后验误差估计．本文方法的主要思想是引进一个高阶间断有限元解做为比较函数，并利用先 验误差估计结果把精确解的后验误差估计转化为两个间断有限元解之间的误差上界估计．数 值实验验证了本文后验误差估计的有效件．本文方法也适用于其他变分问题的有限元逼近的 后验误差分析． 本文中使用通常的Sobolev空间及其范数符号，并用C表示与剖分尺寸h无关的常数． +2011年12月7日收到 1)基金项目：国家自然科学基金资助项目(No 万方数据 216 计算数学 问题及其间断有限元逼近 考虑典型的一阶双曲问题 = Cu三卢·Vu+au (2．1) f，X∈Q u=g，37∈aQ一 其中QC E aQ=aQ+UaQ一，aQ士={X