求解两类偏微分方程的弱Galerkin有限元法

发布时间：2018-06-26 22:17

  本文选题：弱Galerkin有限元法 + 有限元方法　； 参考：《吉林大学》2017年博士论文

【摘要】：弱 Galerkin 有限元法(Weak Galerkin Finite Element Method)是一种求解偏微分方程的高效数值方法,是对经典有限元方法的继承和发展.最早由王军平和叶秀在2011年提出,并用于求解二阶椭圆问题.它与经典有限元方法的不同之处在于定义了弱函数,进一步针对弱函数引入了弱微分(弱导数,弱梯度,弱散度,弱旋度等)的概念,并将其应用到变分形式中,然后对导出的变分形式进行数值求解.弱Galerkin有限元法一经提出,便受到了广泛关注.经过几年的发展,它的内容更加丰富,主要包括:弱Galerkin混合有限元的提出和稳定子的引入,使得弱Galerkin有限元法应用起来更加方便;杂交技术和自适应技术的应用,使得求解速度更快;将多尺度思想引入弱Galerkin有限元法,使得它能够处理来自物理和工程领域的多尺度问题等.此外,随着科研内容的不断完善,该方法能处理问题的类型也越来越多.研究者们已经将弱Galerkin有限元法用于处理包括二阶椭圆交界面问题,Bihar-monic方程求解问题,Helmholtz方程求解问题,Stokes方程求解问题,以及Maxwell方程求解问题等多个经典问题.本文主要研究一维光栅衍射问题和Navier-Stokes方程的求解问题.并针对两类问题各自的难点,给出一些处理技巧,设计相应的弱Galerkin有限元格式.光栅衍射问题存在两个难点:第一个难点是问题求解区域无界.如何对求解区域进行截断,并给出相应的边界条件将直接影响到计算结果.第二个难点是在交界面附近,电磁波的导数变化剧烈.经典的有限元方法不能够很好的模拟交界面附近的振荡或跳跃.因此,寻找更加稳定,精度更高的数值方法是十分必要的.对于第一个难点,在x1方向使用周期边界条件,在x2方向引入透射边界条件将区域截断成矩形区域[0,Λ]×[-b,b].其中Λ表示x1方向的周期,b是一个正常数.之后,在得到的有界区域上求如下Helmholtz方程(△α + k2)uα=，，0 在Ω 中,需要强调的是,uα在x1方向满足周期边界条件,在x2方向满足透射边界条件.周期边界条件较自然,这里主要给出透射边界条件,形式如下:其中,算子△α和Tj,j = 1,2分别在(3.3)和(3.8)中给出定义.针对光栅问题的第二个难点,我们采用带稳定子的弱Galerkin有限元法求解上述方程.稳定子的引入使得弱Galerkin有限元法的解是弱连续的,用它来捕捉电磁波在光栅交界面附近的振荡非常合适.基于上述讨论,我们分别在理论分析和数值实验上验证了算法的可靠性.在理论上,对带有透射边界条件和周期边界条件的弱Galerkin有限元格式进行了细致的分析,给出了解的存在性,唯一性和收敛性的详细证明,并导出了最优收敛阶,形式如下:定理令u∈Hm+1(Ω)和uh∈Vh分别代表方程(3.5)-(3.7)的真解和方程(3.15)的弱Galerκκin有限元解.设eh=Qhu-uh={ei,eb},并假设衍射问题的对偶问题有一个解ω∈H1+s其中s∈(0,1],那么存在一个常数h*0使得‖eh‖wh≤Chm‖u‖m+1,‖ei‖≤Chm+s‖u‖m+1,对任意的h∈(0,h*)均成立.现有关于弱Galerkin有限元法的收敛性分析多是针对Dirichlet边界条件展开的,而对于带透射边界条件的理论结果还很有限.本文的研究在理论上丰富了弱Galerkin有限元法的内容,对于其它类似的问题具有一定的参考价值.在数值方面,我们由简入繁地对三个光栅衍射问题进行了模拟,并从收敛性和光栅效率两方面验证算法的精度.关于收敛性,数值实验结果表明,误差关于‖·‖wh-范数和L2-范数分别达到了 1阶和2阶收敛速率,这与理论分析数据吻合.关于光栅效率,反射效率和透射效率之和稳定在1左右,这与事实吻合.另外,我们还给出了主要衍射级的光栅效率以及复杂光栅具有高振荡的电场强度函数图像.通过对三个光栅问题的模拟,我们发现弱Galerkin有限元法在求解光栅衍射问题时有以下优势:(1)无论对于简单光栅还是复杂光栅,该算法都非常稳定,高效且具有高精度;(2)弱Galerkin有限元法的解是弱连续的,它能较好地模拟电磁波在光栅交界面附近的振荡.Navier-Stokes方程是计算流体力学最主要的方程之一.在一个带有Lipschitz边值条件的有界区域Ω ∈R2内,不可压粘性流可以由与速度u和压力p相关的Navier-Stokes方程来描述,形式如下:ut-v%

本文编号：2071602