论文部分内容阅读
电磁场问题是当今的一个热门的研究课题,其数值模拟有重要的应用价值.经典的电磁场问题一般由包含二阶curl算子的偏微分方程所描述,可采用有限元方法在H(curl)空间中进行求解,其中Nedelec元是最佳选择.求解此类向量型偏微分方程计算量无疑是很大的,因此其高效高精度的数值模拟是当前人们特别关注的问题.超收敛分析恰好是实现高效高精度的一种直接有效的方法.近年来,包含四阶curl算子的偏微分方程频频出现在磁流体动力学问题和透射特征值问题中,它们在很多领域都有广泛的应用.一般在H(cur12)空间中求解此类问题,考虑正方形上最低阶的curl-curl协调元(也称H(curl2)协调元),每个单元上有24个自由度,远远大于最低阶Nedelec元的4个自由度,因此对于这种curl-curl元而言,其超收敛意义是更大的.本论文第一部分内容即在于针对curl协调元和curl-curl协调元研究电磁场问题求解中的超收敛性质.我们选用与Nedelecc元等价的一组由Legendre多项式的组合而成的分层矢量基函数来求解电磁场问题,分析其在二维正方形、三维立方体剖分下的超收敛现象:首先由插值算子的定义以及Legendre多项式的定义可以得到电场E和磁场H的插值的超收敛性质;由Legendre多项式的正交性以及积分恒等式技巧得到E,H超逼近性质;进而得到它们的超收敛估计.值得注意的是,上述超收敛性质都是该单元的天然超收敛性质,也即点态超收敛.基于点态超收敛,我们应用PPR后处理技术进行多项式重构,最终得到E,H的基于L2度量的全局意义上的超收敛.据笔者所知,这是第一篇使用PPR技术获得全局超收敛的电磁场文章.进一步,考虑简化的quad-curl问题,为了应用curl-curl协调元求解且保证解的存在唯一性,需将原问题转化为鞍点问题,进而得到有限元的离散格式.由Babuska-Brezzi理论,本文证明了该有限元解在不同范数意义下的超逼近性质.由Bramble-Hilbert引理等,我们证明了该解在不同范数意义下的超收敛性质.本文第二部分内容则主要针对磁流体动力学问题和透射特征值问题中出现的quad-curl算子,构造高阶分层curl-curl协调的基函数,进而使用有限元方法求解该类问题.我们通过选用指标为-1,-2的广义Jacobi多项式,利用谱元思想,分别构造二维正方形和三维六面体上高阶分层矢量基函数—curl-curl协调元.这种基函数与传统Lagrange型基函数相比,更具有实用性,特别是对这种向量型基函数而言.首先,由多项式的正交性,得到的总体质量矩阵和刚度矩阵具有很好的稀疏性,便于线性方程组的求解;其次,由于其基函数的显示表达,我们可以较为容易的使用任意次多项式,进而可验证其hversion及p-version的收敛性.数值实验表明,我们构造的元可以达到谱精度,且存在很多超收敛现象.据笔者所知,六面体上的curl-curl协调元为本文原创性研究成果.