二维交错网格上的新型间断伽罗金有限元方法-计算数学专业论文.docxVIP

摘 要 近二十年来，间断伽罗金有限元方法飞速发展，广泛用于空气动力学、等离子体物 理、半导体等众多领域，各类新型间断伽罗金有限元方法应运而生。本文主要研究近期提 出的交错网格上的间断伽罗金有限元方法和交错网格上的中心间断伽罗金有限元方法， 并将其用于解决二维问题。数值结果表明:对于守恒律问题，交错网格上的间断伽罗金有 限元方法和交错网格上的中心间断伽罗金有限元方法都是k + 1阶的:对于扩散方程，交错 网格上的间断伽罗金有限元方法的运算量没有增加，且当k I= 1时，交错网格上的间断伽 罗金有限元方法是k+1阶的，而交错网格上的中心间断伽罗金有限元方法始终是k+1阶。 ii Abstract In the past twenty years, the discontinuous galerkin finite element method is widely used in aerodynamics, plasma physics, semiconductor and other fields. By the development of the method, some new discontinuous galerkin methods have been proposed. In this paper, we will discuss the discontinuous galerkin method and the central discontinuous galerkin method on overlapping cells and apply them to problems in 2D. The numerical results showed that: for conservation law, the convergence rate of the discontinuous galerkin method and the central discontinuous Galerkin method are k + 1; for diffusion equation, the convergence rate of the discontinuous galerkin method is k if k = 1, it’s k + 1 if k != 1, the convergence rate of the central discontinuous galerkin method on overlapping cells is k + 1 in both cases. iii 第一童 引 言 间断伽罗金有限元方法(Discontinuous Galerkin Finite element 方法，简称DG((方法〉 是利用完全间断的分片多项式空间作为近似解和试验函数空间的一种有限元方法，它最 初是在1973年由Reed和Hill提出是为了求解中子输运方程[1]: σu + ? · (au) = f, in ? 这里σ是一个实数，而a是一个常向量，当时Reed和Hill所才用的方法是方程两边同乘以试 验函数v ，并且在σ的任意子集K上积分，利利用分部积分原理，原方程写成: σ(u, v)K ? (u, a · ?v)K + (a · nK u, v)?K = (f, v)K , 这里nK 表示?K的单位外法向，且有: f f (u, v)K = uv dx, (w, v)?K = K wv ds K 接下来在?上建立三角网格剖分Th = {K}，用uh表示u的近似解，uh限定在每一个元K上 是至多k阶的多项式。近似解uh由下式唯一确定: σ(uh, v)K ? (uh, a · ?v)K + (h?, v)?K = (f, v)K ?K ∈ Th 其中v ∈ P k (K)，即v为K上至多k阶多项式，而h?是数值流通量由下式确定: h?(x) = a · nK (x) lim(x ? sa). s↓0 1此后LeSaint和Raviart[7]，Delfour，Hager和Trochu[9]，Schwb[10]等人先后对DG方法在 常微分方程中的应用和精度进行了分析。1974年，LeSaint和Raviart[7]首先给出了关于DG方 法的数学分析，并且证明了在L2(?)范数下对于一般的三角剖分的(?x)k 收敛速度，并且在 单变量的笛卡尔网格上对k阶多项式张量积具有(?x)k+1的收敛速度。在1986年，Johnson和 Pitkar¨anta[8]证明了对于一般的三角剖分，DG方法具有(?x)k+ 2 的收敛速度.1988年，Richter [11]在某些二维的非笛卡尔网格下证明