几种求解H(curl)与H(grad)型偏微分方程有限元离散系统的多水平快速算法

摘要

H(curl)型和H(grad)型偏微分方程组(PDEs)在电磁场计算和固体力学等领域具有广泛应用.本文针对频域电磁场PML逼近方程组、频域电磁场隐身超常材料模型和不定时谐Maxwell方程组等几类重要的H(curl)型偏微分方程组的棱有限元离散系统,以及Possion方程和线弹性力学方程组等几类重要的H(grad)型偏微分方程的间断Galerkin(DG)有限元离散系统和二次协调有限元离散系统,研究了相应的多水平快速算法.主要内容如下:
　　针对一类二维电磁场散射问题PML逼近方程的对称不定线性棱有限元离散系统,首次利用迭代两网格法的思想为其构造了一种迭代两网格法,该方法本质地将对原复杂离散系统的求解转化为三个更为简单的辅助系统的求解.同时为其中的两个辅助系统分别构造了高效预条件子,以及基于相应预条件子的PMINRES方法和PCG方法.数值实验表明,与经典的迭代法(如PMINRES法)相比,我们设计的迭代两网格法在迭代次数和CPU时间上均具有明显的优势.
　　通过将完美匹配层(PML)方法、电磁场理论以及隐身超常材料理论等相结合,给出了一种二维频域电磁场隐身超常材料模型,数值实验检验了模型的正确性.进一步,为一类电磁场散射问题PML逼近方程组和上述隐身超常材料模型方程组的线性棱有限元离散系统,结合HX预条件方法,构造了一种Schur补型的多水平预条件子,数值实验验证了算法的有效性.
　　针对不定时谐Maxwell方程组的棱有限元离散系统,构造了一种基于两水平方法和HX预条件子的多水平加性预条件子,本质性地将细网格空间上原不定问题的求解转化为粗网格空间和细网格上curl算子核空间上原不定问题的求解,以及在细网格空间上的对称正定问题的求解.当粗网格尺寸充分小时,证明了基于该预条件子的PGMRES法的一致收敛性.数值实验验证了理论结果的正确性.
　　针对一类二阶椭圆问题的一般对称内罚间断有限元离散系统,通过将DG有限元空间分解为“高频”部分和线性协调元空间,为其构造一种新的稳定的空间分解.利用该空间分解,首先为DG有限元离散系统设计了多水平预条件子.基于辅助空间预条件理论,我们证明了该预条件系统的条件数是一致有界的,同时为其中的线性元系统设计了BPX预条件子,从而得到一种多水平加性预条件子.接着还为DG有限元离散系统设计了一种两水平迭代法,并证明该迭代法是一致收敛的.数值实验验证了上述理论的正确性.
　　针对平面弹性力学问题二次有限元离散系统,给出了一种基于二分网格的局部多层网格(LMG)法.通过将二次有限元空间分解成‘高频’部分和线性元空间,结合二分网格和插值算子的特性,证明了该空间分解的稳定性和强Cauchy-Schwarz不等式成立,进而得到了LMG算法的一致收敛结果.数值实验验证了理论的正确性.

目录

封面

声明

中文摘要

英文摘要

目录

第一章 绪论

1.1 研究背景及现状

1.2 主要工作

1.3 预备知识和实验环境

第二章 一类电磁场PML方程有限元离散系统的迭代两网格法

2.1 模型问题与线性棱有限元方程

2.2 基于PMINRES方法的数值实验

2.3 一种迭代两网格法

2.4 算法实现与数值实验

第三章 求解两类Maxwell方程组的Schur补型多水平预条件子

3.1 模型问题和有限元方程

3.2 一种Schur补型多水平预条件子

3.3 数值实验

第四章 求解不定时谐Maxwell方程组的多水平加性预条件子

4.1 模型问题及有限元方程

4.2 NSPD问题的预条件子一般理论

4.3 一种多水平加性预条件子

4.4 数值实验

第五章 一类间断有限元方程的多水平快速算法

5.1 模型问题与间断有限元方程

5.2 空间分解的稳定性

5.3 一种多水平预条件子

5.4 一种两水平迭代法

第六章 一种求解线弹性力学问题自适应有限元离散系统的

6.1 模型问题和自适应有限元方程

6.2 LMG算法和收敛性分析

6.3 数值实验

总结与展望

参考文献

致谢

个人简历、在学期间发表的学术论文及研究成果