间断有限元方法.doc

《间断有限元方法.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、2016年夏季学期研究生课程考核（读书报告、研究报告）考核科目:间断有限元方法及其应用学生所在院（系）:理学院数学系学生所在学科:学生姓名:学号学生类别考核结果阅卷人1．引言间断Galerkin(DG)方法兼有有限元与有限体积方法的特征。如同一般有限元方法那样，DG方法利用单元多项式空间作为近似解和检验函数空间，但是与传统的有限元方法不同,有限元函数空间基函数都是完全间断的分片多项式，各个单元之间的通信也需要像有限体积方法那样通过在单元边界上构造合适的数值流通量来实现。因此DG方法既保持了一般有限元方法

2、和有限体积方法的优点，又克服了各自的不足。该方法可采用局部高阶插值的方法构造基函数，具有灵活处理边界条件以及可显式求解间断问题的能力，克服了一般有限元方法不适于间断问题的缺点,以及有限体积方法必须通过扩大模板进行重构来提高精度的不足。因此间断Galerkin(DG)方法的出现拓展了传统有限元方法的应用范围,改善了人们对传统有限元方法的认识。2．DG的基本概念间断Galerkin方法最早由Reed和Hill在1973年为解决中子输运方程问题而提出。随后众多学者对间断有限元方法提出了改进和发展特别是90年代

3、以来,以Cockbum和舒其望为代表提出了Runge-Kutta间断Galerkin(RKDG)方法，该方法结合TVD(TVD:TotalVariationDiminishing)Runge-Kutta时间离散方法和间断有限元求解一维双曲守恒律方程（组）以至于高维双曲守恒律方程（组），能够适合复杂计算区域和边界条件，可以精确的捕捉激波和接触间断。它不但在光滑区域可以保证高精度，而且在间断区域可以保持数值无振荡，分辨率高，可以证明收敛到熵解。这些优点使得RKDG成为计算流体力学流行的方法之一，并被广泛应用

4、到气象学、海洋学、湍流、电磁学、石油勘探、水动力学等离子物理和图像处理等领域。同样是在20世纪70年代,内惩罚(IP:InteriorPenalty)类方法被独立地提出来求解摘圆和抛物方程。内惩罚方法后来也被归为间断Galerkin方法一种，本文记为内惩罚间断Galerkin(IPDG)方法。内惩罚间断有限元的发展与同时代求解双曲守恒律的间断有限元方法保持相对对立，该方法的侧重点在于选择合适的惩罚项保持格式的稳定性，而不在于如何构造数值流通量。基于DG方法求解双曲守恒律的巨大成功，许多学者考虑运用DG方

5、法的思想求解扩散方程，但如果只是简单地将DG方法推广到扩散方程得到的数值格式并不准确。例如考虑一维热传导方程。（1）将求解区域剖分为网格其中心点的坐标为，网格步长，记为和为在处的左右极限，即，定义单元端点处的跳跃和均值分别为：，。注意,在不引起混淆的情况下,我们仍然用表示数值解。选取有限元空间为（2）为在单元上的次多项式。定义方程（1）的DG弱解形式为：求，对，使得（3）为数值流量，舒其望己经证明数值流量，简单的取在端点两侧导数均值会导致解不稳定,数值解与精确解是不相容的，与真解有O(1)的误差，称之为

6、“subtleinconsistency”。3．求解扩散方程的各种DG方法的构造3.1LDG方法LDG方法是Cockbum和舒其望在1998年提出的,其思想来源于F.Bassi和S.Rebay求解可压缩Navier-Stokes方程的文章。LDG通过引入辅助变量将含有高阶导数的微分方程写成只含有一阶导数的偏微分方程组,然后用DG方法进行空间离散。引入辅助变量将方程(1)重新写为，然后应用DG得到下列格式：求，对，使得（4）正确的设计数值流量，是得到稳定和高精度方法的关键，可以证明交替地选取左右极限值，或

7、者均可以保证格式稳定而且达到最优收敛阶。当LDG格式(4)中基函数取为每个单元上的局部基函数时，辅助变量是局部可解的，这正是该方法被称为“局部”间断Galerkin方法的由来。3.2Baumann-OdenDG方法Baumann-OdenDG方法并不引入辅助变量,而是通过在弱形式(3)的单元边界上添加惩罚项以保证稳定性。其格式为求，对，使得（5）这里的数值流量直接取为端点两侧导数均值就可以保证该格式的稳定性。实际上Baumann-Oden方法也可以认为是一种内惩罚类方法。3.3dGRPDG方法Gassn

8、er等将双曲守恒律中的黎曼问题推广到扩散方程。对方程(1)两端同时乘以检验函数，在时空区域求积分,通过分部积分可得dGRPDG格式，求对，使得（6）其中，是数值流量,通过求解扩散方程的广义黎曼问题而得到。3.4Cheng-ShuDG方法Cheng-ShuDG方法在热传导方程(1)两边同时乘以检验函数，在单元上积分，通过多次的分部积分将近似解上的导数转移到检验函数上，在界面上设计合适的数值流量得到Cheng-ShuDG格式：求，对，使得（7）