杂交间断有限元方法多边形网格上解二阶椭圆问题_张玲丽.pdf

摘要：本文提出一种杂交间断有限元方法，用于解二阶椭圆问题。首先通过网格单元内部的自由度以及单元边界上的自由度，重构一个新的多项式基函数，然后采用间断有限元加罚方法进行加罚以保证数值格式的稳定性。该方法不仅能够降低间断伽辽金方法的计算量，而且能够适应多边形网格。同时，本文进行了理论分析，首先证明了能量范数下的误差估计，然后通过Aubin-Nitsche论证给出了L2范数下的误差估计。结果表明，两种范数下的收敛阶均是最优的。关键词：弱伽辽金；弱梯度；杂交有限元方法；椭圆问题；误差估计；多边形网格中图分类号：O242.21文献标志码：A文章编号：2096-854X（2022）060106-06A Hybridizable Discontinuous Finite Element Method forSolving Second Order Elliptic Problems on Polygonal MeshesZhang Lingli，Wang Huijuan，Xu Shipeng*（School of Mathematics and Computer Science,Jiangxi Science and Technology Normal University,Nanchang 330038,Jiangxi,P.R.China）Abstract:In this paper,a hybridizable discontinuous finite element method（Hybridizable FEM）is proposed tosolve second-order elliptic problems.Firstly,a new polynomial basis function is reconstructed through the degrees offreedom inside the grid elements and the degrees of freedom on the element boundaries,and then the discontinuousfinite element penalty method is used to ensure the stability of the numerical scheme.This method can not only reducethe computational complexity of the discontinuous Galerkin method,but also adapt to the polygonal meshes.At thesame time,this paper makes a theoretical analysis,first proves the error estimation under the energy norm,and thengives the error estimation under the L2norm through Aubin-Nitsche arguement.The results show that both the orders ofconvergence under the two norms are optimal.Key words:Weak galerkin;weak gradient;hybridizable FEM;elliptic problem;error estimate;polygonal mesh杂交间断有限元方法多边形网格上解二阶椭圆问题张玲丽，王慧娟，徐世鹏*（江西科技师范大学数学与计算机科学学院，江西 南昌330038）【数学计算】收稿日期：2022-06-16最终修回日期：2022-08-14接受日期：2022-08-15基金项目：江西省教育厅项目（GJJ201124）、江西科技师范大学博士启动基金项目（2018BSQD009）作者简介：张玲丽，女，在读硕士研究生，研究方向：偏微分方程数值解；王慧娟，女，在读硕士研究生，研究方向：偏微分方程数值解；*徐世鹏（通讯作者），男，讲师，博士，研究方向：偏微分方程数值解，E-mail:。江西科技师范大学学报Journal of Jiangxi Science&Technology Normal University第6期Issue 62022年12月Dec.20221前言众所周知，间断伽辽金（DG）有限元方法已经被用来模拟各种各样的偏微分方法，得到了广泛的应用。然而，与协调有限元方法相比，间断伽辽金方法应用于扩散问题时有两个缺点：一是未知的自由度数量明显增加，二是所得线性系统的稀疏性要小得多。为了克服未知的自由度数量多的缺点，许多研究员进行了大量地研究，例如，自适应间断有限元方法、重构间断有限元方法、杂交间断有限元方法、2022年高阶杂交有限元方法、虚拟有限元方法、弱间断有限元方法（WG-FEM）等等。在本文中，为了降低间断伽辽金方法的计算量，同时保留其好的特征，从而引入杂交间断伽辽金方法。类似于杂交间断有限元方法与弱有限元方法，方法的自由度存在于网格单元内部以及网格边界上，但是文中提出的方法的基本思想是，通过网格单元内部的自由度以及单元边界上的自由度，重构一个新的多项式基函数，然后间断有限元加罚方法进行加罚保证数值格式的稳定性，这点与已有的有限元方法不同。该方法不仅能够降低间断伽辽金方法的计算量，而且能够适应多边形网格。本文考虑求解一般系数和混合边界条件的椭圆型模型问题：（Au）=fin，（1.1a）Aun=gNon N，（1.1b）ugDon D，（1.1c）其中奂Rd（d2，3）是边界为的有界多边形或多面体开区域。定义：ND，且ND。假设fL2（），gNL2（N）且gDH1/2（D）。对称矩阵系数A：（aij（x）dd奂L（）dd满足一致椭圆性条件，即存在两个常数，0使得：TTA（x）T，Rd，x（1.2）现在，（1.1）式的弱形式是找到uH1（），使得在D上满足ugD，并且：Auvdx=fvdx+NgNvds，vH1D（）（1.3）其中H1D（）：vH1（）：v|D=0。为了完整地表述论文中的结果，我们引入类似于2，定义1.9中Th的匹配单纯形子网格。单纯形子网格Mh=（Th，Fh）满足以下条件：1、形状正则性；2、接触正则性。间断伽辽金方法文献3，对椭圆型问题和抛物型问题进行全面的分析，在对离散化参数作较温和的假设，可以保证不含内部惩罚项的完全稳定。对于右侧的每个单元是分段常数的特殊情况，证明了冒泡稳定的对称格式、冒泡稳定的非对称格式、使用拉格朗日乘子施加连续性的间断伽辽金方法和使用最低阶Raviart-Thomas单元的标准混合有限元方法和用于初始变量的分段常数单元之间的等价性。所有这些不同的公式都由相同的有限元解来满足。最后，将间断有限元方法与混合方法联系起来。最近，由王军平教授与叶秀教授引入的WG-FEM的方法6，7是通过使用离散弱梯度算子设计的，该算子应用于具有一定形状规则的任意多面体的有限元剖分上的不连续分段多项式。弱有限元方法的中心思想是引入弱函数、弱空间、弱算子（弱导数、弱散度、弱旋度等）到变分形式中，同时引入相应的离散弱函数、离散弱空间、离散弱算子（离散弱导数、离散弱散度、离散弱旋度等）到离散数值格式中，然后在接下来的程序中按照有限元框架使用它们。在实际应用中，弱有限元方法具有一些非常好的性质，例如，保留局部/整体质量守恒性质、不需数值格式加罚参数、能够采用多边形网格加密、逼近函数空间容易构造、凝聚技术减少自由度数量等，因此，被用来模拟各种各样的偏微分方程，例如，椭圆界面问题8、双调和问题9、Stokes问题10、Maxwell问题11等等。本文的组织结构如下：第二节介绍一些基本的知识，提出我们的数值方法，为理论分析作准备。第三节将进行理论分析，将证明能量范数与L2范数下的收敛阶是否达到最优。最后，第四节以总结本文的工作作为结束。2准备知识在这节中，将给出一些记号，定义变分形式与有限元空间，引入数值算法，为后面的误差分析作准备。Eh表示单元的所有边的集合，子集EN，ED和EI分别表示在边界N、D上和内部边的集合。对于每个单元TTh，hT表示T的直径。h表示Th上最大单元的直径，即h=maxTThhT。设T是Th中边界为T的单元。T上的弱函数为张玲丽，王慧娟，徐世鹏：杂交间断有限元方法多边形网格上解二阶椭圆问题107江西科技师范大学学报第6期v=v0，vb，其中v0L2（T）和vbL2（T）。对于每个TTh，其弱空间为：V（T）：=v=v0，vb：v0L2（T），vbL2（T），其中v0不一定是vb的迹。图1定义形状正则的多边形。左图：三角剖分网格。右图：对应的多边形网格考虑包含映射iV：H1（T）V（T）通过iVq：q|T，q|T，qH1（T）。通过这个映射，可以把Sobolev空间H1（T）嵌入到V（T）中，因此通过用iVq来作用每个qH1（T），可以认为H1（T）是V（T）的一个子空间。用TTh上局部空间V（T）定义全局网格Th上的弱空间：V：v=v0，vb：v0|TL2（T），vb|EL2（E），TTh，EEh。值得注意的是，弱空间V中每个v=v0，vb的分量vb在Eh上是单值的，即v01，vb1：=v0，vb|T1V（T1），v02，vb2：=v0，vb|T2V（T2）且E=T1T2，则vb1|E=vb2|E。设Pk（D）和Pk（D?）分别是区域DRd和D?Rd-1上的次数不超过k0的多项式的集合。定义每个TTh上的局部弱离散空间，其中k0，Vh（T）：=vh=v0，vbV（T）：v0|TPk（T），vb|TPk+1（T）定义Th上的全局弱离散空间：Vh（T）：=vh=v0，vbV：v0|TPk（T），vb|EPk+1（E），TTh，EEh和V0h：vh=v0，vbVh：vb|TD=0，TTh，EEh。对于每个TTh，给定vv0，vbV（T），引入一个投影算子R1，（2.8a）-（2.8b），在KHh（T）上定义|K：v|KPk+1（E），使得：K（v0R|Kv）dx=0Pk（K）（2.1a）E（v0R|Kv）ds=0Pk1（K）（2.1b）其中k0且E：KTEh，即E是Eh上三角形单元K的边或面。由于投影算子R在本文的分析中起着重要的作用，将利用4，5中的技巧在下面的引理中给出一种证明方法，证明它的存在性和近似性。引理1：投影算子R是适定的，并且具有逼近性质：如果uHk2（K），则有：u-Ru0，KhK1/2u-Ru0，KhK（u-Ru）0，KChk+2Ku Hk+2（K），（2.2）其中C是与单元KHh（T）无关的常数。证明：（2.1a）-（2.1b）中的方程数满足d+k()deqn（2.1a）+d+kd-()1eqn（2.1b）=d+k+1()ddimPk+1（K），这意味着我们只需要证明投影的唯一性，即如果u=0，则Ru=0。设i（1i3）是K的重心坐标，在E上1=0。根据（2.1b），在E上有u=0，Ru=0，因此对于一些pPk（K）有Ru=1p。在（2.1a）中取=p，我们得到：K1p2dx=0。考虑到当xK时，10，得到p=0，因此Ru=0。显然，0，K?和下面的形式supv?Pk（K?）（，v?）K?v?0，K?+supw?Pk+1（E?）（，w?）E?w?0，E?1082022年都是Pk+1（K?）上的范数，其中K?表示参考单元，因此由于是有限维的，所以等价于下列式子：u0，K?CsupvPk（K?）（u，v）K?v0，K?0，K?+supwPk+1（E?）（u，w）E?w0，E?()，uPk+1（K?）。（2.3）设P：L2（K）Pk+1（K）是L2投影：K（u-Pu）vdx=0，vPk+1（K）有如下性质：若uHk+2（K），则u-Ru0，KhK1/2u-Ru0，KhK（u-Pu）0，KChk+2KuHk+2（K）。（2.4）利用（2.1