椭圆型偏微分方程的弱有限元方法.pdfVIP

中国科学: 数学 2015 年 第45 卷 第7 期: 1061 1092 综 述 椭圆型偏微分方程的弱有限元方法 献给林群教授 华诞 ∗ 王军平 王春梅 Division of Mathematical Sciences, National Science Foundation, Arlington, VA 22230, USA; School of Mathematics, Georgia Institute of Technology, Atlanta, GA 30332-0160, USA; 南京师范大学泰州学院数学系, 泰州 225300 E-mail: jwang@nsf.gov, cwang462@math.gatech.edu 收稿日期: 2014-11-03; 接受日期: 2015-04-15; * 通信作者 The research of Wang Junping (王军平) was supported by the National Science Foundation IR/D Program 摘要 弱有限元方法(weak Galerkin finite element methods), 简称WG 方法, 是求解偏微分方程的一 种全新且高效的数值方法. WG 有限元方法的主要思想是, 对间断函数引入广义弱微分算子, 并将其 应用于通常的变分形式中以对相应的偏微分方程进行数值求解, 数值解的连续性则通过稳定子以弱形 式来实现. 本文以二阶椭圆问题为例, 详细介绍弱有限元方法的原理和基础, 并给出相关的理论分析. 此外, 本文简单介绍其他椭圆问题的弱有限元方法的离散格式. 弱有限元方法的最大特点是, (1) 有限 元剖分允许任意多边形或多面体; (2) 总体刚度矩阵可通过单元刚度矩阵叠加而得; (3) 逼近函数构造 简单, 且极易满足相应的稳定性条件; (4) 格式可做杂交处理以并行消去单元内部自由度. 关键词 弱有限元方法 弱梯度 弱散度 弱旋度 稳定子 杂交弱有限元方法 椭圆方程 主题分类 65N15, 65N30, 35J50, 76D07, 35B45 引言 有限元方法是数值求解椭圆边值问题的一种方法. 此方法首先于 20 世纪 50 年代初由工程师们 提出, 并用于求解简单的结构问题. 有限元方法作为一种系统的数值方法, 并奠定其数学基础, 则是20 世纪 60 年代中期, 由以冯康先生为代表的中国学者与西方学者独立完成. 从历史上看, 还应该提到 Courant, 他在1943 年已经提出过在三角形网格上用分片线性函数去逼近二阶椭圆问题的Dirichlet 问 题, 这是有限元方法最原始的思想. 至今, 作为一类实用有效的计算工具, 有限元方法已广泛应用于科 学研究领域众多问题的数值求解; 同时在数学上已建立了一套完整的理论, 成为计算数学、应用数学 研究者和工程师们的学习和研究科目. 有限元方法以及与它相关的各种新的计算方法的构造和分析仍 在层出不穷地涌现, 目前仍处在百花齐放的阶段. 用一个简单的模型问题说明有限元方法的主要思想. 考虑Poisson 方程的齐次Dirichlet 边值问题: − ∆u = f, 在Ω, 英文引用格式 Wang J P, Wang C M. Weak Galerkin finite element methods for elliptic PDEs (in Chinese). Sci Sin Math, 2015, 45: 1061–1092, doi: 10.1360/N012014-00233 王军平等: 椭圆型偏微分方程的弱有限元方法 u = 0, 在∂Ω, 其中Ω 是平面2 中的一个有界区域. 其相应的变分问题是, 求u ∈ H 1 (Ω), 使得 0 (∇u, ∇v) = (f, v), ∀ v ∈ H 1 (Ω). (1.1)