几类有限体积元及有限体积格式的数值分析

论文题目: 几类有限体积元及有限体积格式的数值分析

论文类型: 博士论文

论文专业: 计算数学

作者: 杨旻

导师: 袁益让

关键词: 有限体积元方法,有限体积方法,高精度,对称格式,半正定问题,误差估计

文献来源: 山东大学

发表年度: 2005

论文摘要: 有限体积元（FVE）和有限体积（FV）方法作为守恒型问题的离散方法具有较长的历史,这类方法被广泛的应用于科学工程计算领域中,如流体力学,热传导,物质运移以及石油工业等,在早期的文献中[2-6],这类方法亦被称为积分有限差分方法,盒式方法和广义差分方法,并且大部分相关的数值结果是关于线性问题的。有限体积元和有限体积格式的建立是基于所谓的”平衡”方法:首先,将问题在每一个离散单元（称为控制体积或对偶单元）上表示为局部平衡形式;其次,利用散度公式,将该平衡形式表示为积分守恒形式;最后,根据计算要求采用不同的技巧将守恒形式离散。一般而言,我们可以将有限体积元和有限体积方法看作是介于有限元和有限差分方法之间的第三类离散方法。它们一方面具有有限元方法的灵活性,适于处理复杂区域及边界问题,另外一方面,类似有限差分方法离散格式简单易于计算。有限体积元和有限体积方法的主要区别在于有限体积元方法涉及两个空间,其中解空间为初始剖分上的分片多项式函数空间,检验函数空间为对偶剖分上的分片常数函数空间,它类似有限元方法采用Galerkin技巧先形成变分形式再求解。而有限体积方法直接利用差分技巧在每个控制体积上对平衡方程离散后再求解。有限体积元和有限体积方法被广泛应用源于其格式构造的简便以及保持数值流量的局部守恒性这两个优点,尤其是第二个优点很重要,它使得的该方法在估计诸如流体力学,半导体模拟,热传导等流量为主的问题时非常有效。尽管如此,这类方法也有一些比较严重的缺点:1.基本上,大部分FVE和FV格式关于空间均只有一阶收敛精度;2.对于FVE方法而言,由于解空间和检验函数空间不相同,导致了高维问题即便是定常系数,其离散格式的系数矩阵往往也不对称,这无论给问题的理论分析还是实际计算都带来了较大的困难,因为在实际工程计算中,很多常用的算法,如共轭梯度法,都依赖于系数矩阵的对称性。为了获得高精度的格式可以采用两种方法。1.选用规则或对称的网格剖分。Cai,Ewing等人证明了对于椭圆和积分微分方程,当采用对称网格时,相应的FVE格式的H1模误差估计具有二阶空间精度。Vassilevski,Petrova和Lazarov基于三角形网格提出了一类FV格式,当网格规则时,获得了超收敛的结果,即O（h2）。Weiser和Wheeler研究了非一致张量积网格上的FV方法以及超收敛性。Lazarov,Mishev和Vassilevski基于正方形网格剖分提出了两类迎风FV格式,包括修正迎风格式和Il’in格式,这些格式都具有O(hm-1),3/2≤m≤3的收敛精度。但是通过第一种方法构造高精度格式,一方面对网格剖分有较大的限制,另外要求处理的问题比较简单,因此在实际应用中有一定的局限性。2.对于FVE格式,还可以通过选取高次有限元空间作为解空间以获得高精度收敛结果。田和陈以Poisson方程为模型,基于三角形网格上的二次Lagrange元和三次Hermite元,首先提出了两类高精度FVE格式并证得了最优的H1模误差估计,而后陈又研究了常系数

论文目录:

Abstract

Abstract (In Chinese)

Chapter 1 High Accuracy Finite Volume Element Methods On Triangular Meshes

§1.1 Quadratic Finite Volume Element Methods for Nonlinear Parabolic Systems

§1.1.1 Meshes and notation

§1.1.2 Some auxiliary results

§1.1.3 FVE scheme and error estimates

§1.1.4 Numerical experiment

§1.2 A Multistep Finite Volume Element Scheme Along Characteristics for Nonlinear Convection Diffusion Problems

§1.2.1 Meshes and notation

§1.2.2 Multistep FVE scheme and error estimates

§1.2.3 Numerical experiment

§1.3 Cubic Finite Volume Element Methods for Elliptic Equations With Variable Coefficients

§1.3.1 Notation and method

§1.3.2 Some auxiliary results

§1.3.3 Error estimates

§1.3.4 Numerical experiment

Caapter 2 Symmetric Finite Volume Element Methods Based On Linear Tensor Product Elements for 2-D And 3-D Convection Diffusion Problems

§2.1 Symmetric Finite Volume Element Methods Along Characteristics for 3-D Convection Diffusion Problems

§2.1.1 Meshes and notation

§2.1.2 Some auxiliary results

§2.1.3 Symmetric FVE scheme and error estimates

§2.1.4 Numerical experiment

§2.2 Symmetric Multistep Finite Volume Element Methods Along Characteristics for 2-D Nonlinear Convection Diffusion Problems

§2.2.1 Meshes and notation

§2.2.2 Some auxiliary results

§2.2.3 Symmetric multistep FVE scheme and error estimates

§2.2.4 Numerical experiment

Chapter 3 A Symmetric Finite Volume Element Scheme With Second Order Accuracy for Nonlinear Convection Diffusion Problems

§3.1 Introduction

§3.2 Meshes and notation

§3.3 Some auxiliary results

§3.4 Symmetric high accuracy FVE scheme and error estimates

§3.5 Numerical experiment

Chapter 4 Multistep Finite Volume Methods for Positive Semidefinite Two-Phase Incompressible Flow In Porous Media

§4.1 Introduction

§4.2 Multistep Finite Volume Scheme

§4.3 Error Estimates

§4.4 Numerical Experiment

References

Acknowledgments (In Chinese)

Curriculum Vitae (In Chinese)

学位论文评阅及答辩情况表

发布时间: 2005-10-17

参考文献

[1].波方程中一些新的能量守恒有限体积元方法[D]. 闫金亮.南京师范大学2016
[2].两类界面问题的有限体积元方法[D]. 朱玲.南京师范大学2015
[3].界面问题的有限体积元法研究[D]. 高艳妮.吉林大学2016
[4].发展型方程的混合有限体积元方法及数值模拟[D]. 方志朝.内蒙古大学2013
[5].流体力学中几类波方程的有限体积元方法[D]. 王全祥.南京师范大学2013
[6].混合有限体积元法[D]. 田万福.吉林大学2010
[7].热传导型半导体瞬态问题的数值解法和分析[D]. 陈传军.山东大学2006
[8].电阻抗成像的数值模拟和分析[D]. 李久平.山东大学2008
[9].解二阶椭圆型方程的高次有限体积元法的若干研究[D]. 丁玉琼.吉林大学2010
[10].几类随机波方程的有限体积元算法研究[D]. 朱泉涌.南京师范大学2014

几类有限体积元及有限体积格式的数值分析

猜你喜欢