二维超材料中麦克斯韦方程组棱有限元的超收敛分析-计算数学专业论文.docxVIP

摘 要 本文运用混合有限元方法求解了二维超材料中时域麦克斯韦方程组的控制模 型。本文所使用的混合有限元空间是第一类线性 ? ?′??′??? 元。尽管内边中点不是自 然的超收敛点，但是经过平均，我们在三角形网格的内边中点得到了超收敛结果。 经计算，我们发现了将有限元解在两个相邻的三角形单元作平均以后，真解在 内边中点的值与有限元解在该点的平均值之间的误差是 ?(?2)，这比有限元解和真 解的 ?2 误差 ?(?) 高出了一阶，于是我们在三角形网格的内边中点得到了超收敛 结果。基于这个发现，我们给出了相应的理论证明。 本文还得到结论: 等边三角形网格的中心不是超收敛点，并且给出了相应的理 论证明和数值算例。 本文进一步地将所得到的结果推广到了线性的 ? ?′??′??? 矩形元，矩形网格的中 心点是超收敛点是已经存在的结果。经过平均技术，我们在矩形网格内边中点也得 到了超收敛结果，并且给出了相应的理论证明和数值算例。 关键词：麦克斯韦方程组；超材料；? ?′??′??? 元；超收敛 Abstract In this paper, mixed finite element method is used to solve the governing model of time dependent Maxwell’s equations involving metamaterial in 2-D. The mixed finite element spaces here are the first family of linear ? ?′??′??? elements. Though there exits no natural superconvergence points at the midpoints of the interior edges for the triangular mesh, superconvergence happens here after a simple average tech- nique. In our computational experiment, we find that after a simple average of the finite element solution over two over the two neighboring elements, an ?(?2) error can be obtained between the analytic solution and the average value of the finite element solutions at the midpoints of interior edges over the two neighboring ele- ments and this error is one order higher than the standard ?2 error ?(?) between the analytic solution and the finite element solution. Thus, We obtain the super- convergence results at the midpoints of the interior edges and we also provide the corresponding comprehensive analysis. In this paper, we also present a result that the centers of the equilateral trian- gular meshes are points without the phenomenon of superconvergence. The corre- sponding analysis and computational experiments are also carried out. Furthermore, we extend our result to the lowest order ? ?′??′??? rectangular mesh. We’ve already known that superconvergence happens at the centers of the rectangular mesh. In this paper, after a simple average technique, we find that superconvergence also happens at