杂交间断有限元方法多边形网格上解二阶椭圆问题_张玲丽.pdf
《杂交间断有限元方法多边形网格上解二阶椭圆问题_张玲丽.pdf》由会员分享,可在线阅读,更多相关《杂交间断有限元方法多边形网格上解二阶椭圆问题_张玲丽.pdf(7页珍藏版)》请在咨信网上搜索。
1、摘要:本文提出一种杂交间断有限元方法,用于解二阶椭圆问题。首先通过网格单元内部的自由度以及单元边界上的自由度,重构一个新的多项式基函数,然后采用间断有限元加罚方法进行加罚以保证数值格式的稳定性。该方法不仅能够降低间断伽辽金方法的计算量,而且能够适应多边形网格。同时,本文进行了理论分析,首先证明了能量范数下的误差估计,然后通过Aubin-Nitsche论证给出了L2范数下的误差估计。结果表明,两种范数下的收敛阶均是最优的。关键词:弱伽辽金;弱梯度;杂交有限元方法;椭圆问题;误差估计;多边形网格中图分类号:O242.21文献标志码:A文章编号:2096-854X(2022)060106-06A H
2、ybridizable 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 hy
3、bridizable 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 discont
4、inuousfinite 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
5、 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;polyg
6、onal mesh杂交间断有限元方法多边形网格上解二阶椭圆问题张玲丽,王慧娟,徐世鹏*(江西科技师范大学数学与计算机科学学院,江西 南昌330038)【数学计算】收稿日期:2022-06-16最终修回日期:2022-08-14接受日期:2022-08-15基金项目:江西省教育厅项目(GJJ201124)、江西科技师范大学博士启动基金项目(2018BSQD009)作者简介:张玲丽,女,在读硕士研究生,研究方向:偏微分方程数值解;王慧娟,女,在读硕士研究生,研究方向:偏微分方程数值解;*徐世鹏(通讯作者),男,讲师,博士,研究方向:偏微分方程数值解,E-mail:。江西科技师范大学学报Journa
7、l of Jiangxi Science&Technology Normal University第6期Issue 62022年12月Dec.20221前言众所周知,间断伽辽金(DG)有限元方法已经被用来模拟各种各样的偏微分方法,得到了广泛的应用。然而,与协调有限元方法相比,间断伽辽金方法应用于扩散问题时有两个缺点:一是未知的自由度数量明显增加,二是所得线性系统的稀疏性要小得多。为了克服未知的自由度数量多的缺点,许多研究员进行了大量地研究,例如,自适应间断有限元方法、重构间断有限元方法、杂交间断有限元方法、2022年高阶杂交有限元方法、虚拟有限元方法、弱间断有限元方法(WG-FEM)等等。在本
8、文中,为了降低间断伽辽金方法的计算量,同时保留其好的特征,从而引入杂交间断伽辽金方法。类似于杂交间断有限元方法与弱有限元方法,方法的自由度存在于网格单元内部以及网格边界上,但是文中提出的方法的基本思想是,通过网格单元内部的自由度以及单元边界上的自由度,重构一个新的多项式基函数,然后间断有限元加罚方法进行加罚保证数值格式的稳定性,这点与已有的有限元方法不同。该方法不仅能够降低间断伽辽金方法的计算量,而且能够适应多边形网格。本文考虑求解一般系数和混合边界条件的椭圆型模型问题:(Au)=fin,(1.1a)Aun=gNon N,(1.1b)ugDon D,(1.1c)其中奂Rd(d2,3)是边界为的
9、有界多边形或多面体开区域。定义: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,对椭圆型问题和抛物型问题进
10、行全面的分析,在对离散化参数作较温和的假设,可以保证不含内部惩罚项的完全稳定。对于右侧的每个单元是分段常数的特殊情况,证明了冒泡稳定的对称格式、冒泡稳定的非对称格式、使用拉格朗日乘子施加连续性的间断伽辽金方法和使用最低阶Raviart-Thomas单元的标准混合有限元方法和用于初始变量的分段常数单元之间的等价性。所有这些不同的公式都由相同的有限元解来满足。最后,将间断有限元方法与混合方法联系起来。最近,由王军平教授与叶秀教授引入的WG-FEM的方法6,7是通过使用离散弱梯度算子设计的,该算子应用于具有一定形状规则的任意多面体的有限元剖分上的不连续分段多项式。弱有限元方法的中心思想是引入弱函数、
11、弱空间、弱算子(弱导数、弱散度、弱旋度等)到变分形式中,同时引入相应的离散弱函数、离散弱空间、离散弱算子(离散弱导数、离散弱散度、离散弱旋度等)到离散数值格式中,然后在接下来的程序中按照有限元框架使用它们。在实际应用中,弱有限元方法具有一些非常好的性质,例如,保留局部/整体质量守恒性质、不需数值格式加罚参数、能够采用多边形网格加密、逼近函数空间容易构造、凝聚技术减少自由度数量等,因此,被用来模拟各种各样的偏微分方程,例如,椭圆界面问题8、双调和问题9、Stokes问题10、Maxwell问题11等等。本文的组织结构如下:第二节介绍一些基本的知识,提出我们的数值方法,为理论分析作准备。第三节将进
12、行理论分析,将证明能量范数与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(
13、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
14、,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)
15、,在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()
16、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?表示参考单元,因此由于是有限维的,
17、所以等价于下列式子: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.1a),(2.3)和柯西-施瓦茨不等式,得到:Ru-Pu0,K?CsupwPk+1(E?)(Ru-Pu,w)E?w0,E?CRu-Pu0,E?。利用一个简单的标度引理和逆不等式,得到:Ru-Pu0,KC
18、h1/2KRu-Pu0,ECh1/2KR(u-Pu)0,ECh1/2Ku-Pu0,E和(Ru-Pu)0,KCh1KRu-Pu0,ECh1/2Ku-Pu0,E。最后,利用性质(2.4)、三角不等式和迹不等式可以得出证明。如果EKK,则当E奂,v|E=v|E,q|E=q|E,v|E=v|En,q|E=q|En时,标量函数v和向量函数q的跳跃和平均值为:v|E=12(v+v-),v|E=v+n+v-n-,q|E=12(q+q-),q|E=q+n+q-n-。EThKqnvdx=EE0hEq+n+v+q-n-v-ds+EE0hEqnvds=EEhEqv+qvds=EEhEqvds。于是,给出杂交间断有限
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 杂交 间断 有限元 方法 多边形 网格 上解二阶 椭圆 问题 张玲丽
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。