求解对流扩散方程的全离散局部间断Galerkin方法

来源 :南京大学 | 被引量 : 2次 | 上传用户:ly_exe
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
局部间断Galerkin(LDG)方法自提出以来已经被广泛应用于求解对流扩散方程,并在高阶偏微分方程的数值求解方法研究中取得极大成功,然而相比于数值格式的快速发展和应用,理论研究则相对滞后;相关文献中关于全离散LDG方法的理论分析更是凤毛麟角。在实际计算中,求解发展型方程总是要借助于一定的时间离散方法,因此全离散方法的理论分析至关重要。本文我们将以对流扩散方程为主体,着重研究全离散LDG方法的稳定性和误差估计。本文的主要内容有三章:  第二章,我们考虑显式Runge-Kutta(RK)全离散LDG方法(简记为EXRK-LDG)求解带Dirichlet边界条件的对流扩散方程。此时有两个主要难点:一是区域边界处的数值流通量设置方式;二是RK方法每个中间时间层在边界处的设置。不适当的设置方法会影响格式的整体精度。本文将利用能量分析方法,建立三阶EXRK-LDG方法的误差估计,进而给出数值流通量的一种经济有效的选取方法和每个中间时间层的边界设置方法。在这样的设置方法下,可以证明,当时间步长(τ)满足CFL条件c(τ)/h<λc和d(τ)/h2≤λd时,三阶EXRK-LDG格式在时间和空间上都能达到L2模的最优阶收敛。这里c,d分别是对流项和扩散项系数,h是空间网格尺寸,λc和λd是给定的CFL数。  如果对流占优,则显格式的时空限制条件为(τ)=O(h),显式时间离散方法是一个很好的选择。但是对于扩散占优情形,显式时间离散对时间步长的限制为(τ)=O(h2),这个条件比较苛刻。为此,我们也将研究一类半隐半显式(Implicit-Explicit,简称IMEX)时间离散方法,对于对流项采用显式离散方式,而对于扩散项采用隐式离散方式。这类时间离散方法能克服显式时间离散小时间步长的限制,可以高效求解扩散占优的对流扩散问题,尤其是扩散部分是线性而对流部分是非线性的情形。  第三章,考虑RK型的IMEX全离散LDG格式(简记为IMEX-RK-LDG)。从一维线性对流扩散方程出发,通过建立数值解的梯度和跳量与梯度的数值解之间的重要关系,以及IMEX-RK-LDG格式满足的能量方程,我们将利用能量分析方法证明,在时间步长满足(τ≤<(τ)0的条件下,几个特殊的IMEX-RK-LDG格式是L2稳定的,这里(τ)0不依赖于空间尺寸h,只与对流项和扩散项系数有关。严格的分析表明,(τ)0与扩散项系数成正比,与对流项系数的平方成反比。在这个条件下,我们也将证明IMEX-RK-LDG格式具有最优的L2模收敛阶。  第四章,以一维和高维非线性对流扩散方程为模型,研究IMEX-RK-LDG方法和多步IMEX全离散LDG方法(简记为IMEX-MS-LDG)的稳定性和误差估计。其中,第一节通过建立与线性情形类似的LDG空间离散性质,并借助于先验误差假设,得到与第三章类似的结果。第二节将利用能量分析方法证明,在时间步长满足(τ)≤τ0的条件下,几个特殊的IMEX-MS-LDG格式满足能量范数稳定,且具有最优的L2模收敛阶。第三节通过建立高维空间上数值解的梯度和跳量与梯度的数值解之间的重要关系,得到与第三章类似的稳定性结果,同时,借助于间断有限元空间的椭圆投影,给出高维空间上IMEX-RK-LDG方法的最优L2模误差估计。
其他文献
该文采用时域有限差分方法(FDTD)分析了圆筒形进气道模型(包括终端开路、闭路和终端含有叶片)和S-形(或称后端倾斜方形)进气道模型散射特性,并分别给出了远区场和RCS结果.同
本文首先给出了张量范畴和拟Hopf代数的定义,然后证明了拟Hopf代数的表示范畴是一个张量范畴。接下来,我们把朱永昌等关于阶和指标的定义(见参考文献[32])推广到了张量范畴,并且证
该文通过对中国固定资产投资规模波动的实证分析,表明对其进行监测预警及调控的必要性与可行性.在此基础上,建立了固定资产投资规模监测预警的突变模型.随后,给出关于在建投
本硕士论文中集中了作者在攻读硕士学位期间的主要研究成果,主要研究的对象有:一元浅水波方程:修正的Novikov方程和Dullin-Gottwald-Holm方程;二元浅水波方程:二元的Camassa-H
图论是离散数学的一个重要的分支,它在生产管理,军事,交通运输,计算机科学与技术,通信工程等领域都有着广泛的应用.1736年,Euler[14]发表的关于K(o)gsberg七桥问题的论文是图论领域
国外某大型水泥公司生产线配备若干台高效胶带提升机,型号是BW-G800/400,产量为450~500 t/h,提升高度为91 m,提升速度为1.9 m/s。传动系统由主传动和辅助传动组成,主传动驱动
该文根据前不变凸函数的定义来定义了前不变凸集值映射,并且推广定义了各种广义不变凸集值映射,研究了它们的基本性质,及其之间的关系,在此基础上讨论了前不变凸集值映射在优