对流扩散问题多尺度方法-计算数学专业论文.docxVIP

中国科学技术大学学位论文原创性和授权使用声明 本人声明所呈交的学位论文，是本人在导师指导下进行研究工 作所取得的成果。除己特别加以标注和致谢的地方外，论文中不包含 任何他人己经发表或撰写过的研究成果。与我一同工作的同志对本 研究所做的贡献均己在论文中作了明确的说明。 本人授权中国科学技术大学拥有学位论文的部分使用权，即:学 校有权按有关规定向国家有关部门或机构送交论文的复印件和电子 版，允许论文被查阅和借阅，可以将学位论文编入有关数据库进行检 索，可以采用影印、缩印或扫描等复制手段保存、汇编学位论文。 保密的学位论文在解密后也遵守此规定。 作者签名: 年 月 日 第 1 第 1 章 绪论 第 1 第 1 章 绪论 PAGE PAGE 11 PAGE PAGE 10 第 1 章 绪论 在本文中，我们将研究两部分内容:第一部分我们构造了奇异摄动对流扩 散问题的指数拟合局部间断有限元方法;第二部分我们构造了多尺度对流扩散 问题的异质多尺度内部惩罚间断有限元方法。 在这一章中，我们首先简单的介绍间断有限元方法和多尺度方法。然后给 出本文所研究的问题和得到的主要结果。 1.1 间断有限元方法 间断有限元 (DG) 方法最早是 1973 年由 Reed 和 Hill [66] 提出，用于求解稳 态线性双曲中子输运方程: σu + ? · (au) = f, 2其中 σ 是一个实数，a(x) 是分片常数。在 1974 年，LeSaint 和 Raviart [61] 第一次 给出了这个方法的误差分析，并且证明了当使用 k 次分片多项式函数作为基函 数时，在一般的三角形网格上的 O(hk ) 收敛性和在 Cartesian 网格上的最佳收敛 阶 O(hk+1)。随后在 1986 年，Johnson 和 Pitkara¨nta [57] 将这个结论提高到一般 三角形网格上的 O(hk+ 1 ) 收敛，丽这个收敛阶在 1991 年由 Peterson [65] 通过数 值算例证实了是最佳的。我们先通过二维稳态线性对流方程来简单介绍这个方 法 [10]。考虑二维稳态线性对流问题: 2 aux + buy = f (x, y), 0 ≤ x, y ≤ 1, u(x, 0) = g1(x), 0 ≤ x ≤ 1, (1.1) u(0, y) = g2(y), 0 ≤ y ≤ 1, 其中对流系数 a 和 b 是正常数。求解区域 [0, 1]2 上的网格定义为 Ii,j = {(x, y) : xi? 1 ≤ x ≤ xi+ 1 , yj? 1 ≤ y ≤ yj+ 1 }, 1 ≤ i ≤ Nx, 1 ≤ j ≤ Ny , 2 2 2 2 22其中 2 2 0 = x 1 2  x 3 2  · · · xNx+ 1 = 1 和 0 = y 1 2  y 3 2  · · · yNy + 1 = 1 分别是 x 方向和 y 方向上的离散。在方程两端乘以连续的检验函数 v，并且在区 间 Ii,j 上积分，运用分部积分公式，可以得到下面的弱形式 r r ? (auvx + buvy )dxdy + a r yj+ 1  2u(xi+ 1 , y)v(xi+ 1 , y)dy 2 Ii,j 22r yj+ 1 2 2 yj? 1 2 2 2r xi+ 1 2 ?a u(xi? 1 , y)v(xi? 1 , y)dy + b u(x, yj+ 1 )v(x, yj+ 1 )dx 2 2 2yj? 1 2 2r xi+ 1 2 2 2 2xi? 1 2 r r ?b u(x, yj? 1 )v(x, yj? 1 )dx = f vdxdy. (1.2) 2xi? 1 2 2 2 Ii,j 间断有限元区别于一般有限元的主要之处就是在间断有限元空间 V kh = {v : v|Ii,j ∈ V k k (Ii,j ), 1 ≤ x ≤ Nx, 1 ≤ y ≤ Ny } (1.3) 上求解问题，其中 P k (Ii,j ) 表示的是单元 Ii,j 上次数不超过 k 次的多项式空间。 在这个空间里，并不要求函数在单元边界上连续。现在取解函数 uh 和检验函数 hvh 都位于间断有限元空间 V k 中。此时由于 uh 和 vh 在边界上间断，从丽边界项 h uh(xi+ 1 , y), uh(x, yj+ 1 ), uh(xi? 1 , y), uh(x, yj? 1 ) 2 2 2 2 和 vh(xi+ 1 , y), vh(x, yj+ 1 ), vh(xi? 1 , y), vh(x, yj? 1 ) 2 2 2 2 没有明确的定义。检验函数 vh 在单元 Ii,j 边界上的值通常选取为内部值: i+ 1