论文部分内容阅读
本文主要针对带间断系数的扩散问题和奇异摄动反应-扩散问题,在各向异性网格条件下推导健壮的局部重构型后验误差估计,进而实现各向异性自适应计算。针对带间断系数的扩散问题,推导了适用于各向异性网格的局部重构型后验误差估计。基于平衡通量的局部守恒假设,一个抽象的误差上界首先被建立。剩下的任务就是构造满足假设的平衡通量。这里,我们给出了平衡通量的两种构造方式:直接法和混合有限元法。其中直接法是最容易想到的,对网格各向异性的影响不需要特别的处理,而混合有限元法需要将网格的各向异性尺寸引入到平衡通量的构造中。这类后验估计关于扩散系数是健壮的,且可应用于有限体积法,有限差分法和有限元法。我们分别在能量范数和对偶范数意义下表明了所构造的误差估计子是可靠的、有效的。事实上,这些估计在各向同性网格条件下与Vohralík(J Sci Comput46:397–438,2011)得到的各向同性后验估计等价,因此从某种意义上来说,我们推广了Vohral′k的结果。 本文针对奇异摄动反应-扩散问题的有限体积法,推导了适用于各向异性网格的局部重构型后验误差估计。类似地,基于平衡通量的局部守恒假设,我们首先在能量范数意义下建立了一个抽象的误差上界。这个平衡通量的构造采用直接法。我们在能量范数意义下表明了所构造的误差估计子是可靠的、有效的且关于反应系数是健壮的。事实上,这些估计在各向同性网格条件下与Cheddadi等人(ESAIM-Math Model Numer Anal43:867–888,2009)得到的各向同性后验估计等价,因此从某种意义上来说,我们推广了Cheddadi等人的结果。基于所得的各向异性后验误差估计,我们推导出相应的各向异性误差指示子,从而用于引导网格的各向异性自适应加密。事实上,原始形式的各向异性后验估计并不能直接用于引导网格的各向异性自适应加密,因为它不能表示出误差在各个方向的贡献,且对齐测度中含有未知的真解。为了得到各向异性误差指示子,我们改写了对齐测度的格式,并通过后处理的方式得到它的一个逼近。最后,我们给出了一个各向异性网格自适应加密算法,能够利用所得到的各向异性误差指示子生成一个合适的度量场,从而生成一个新的各向异性网格。