论文部分内容阅读
本文针对不同类型的非线性方程(包括磁流体动力学方程、带有阻尼项的Stokes方程、对流扩散方程、反应-扩散方程及非线性抛物方程等),分别从非协调有限元方法、最小二乘有限元方法、误差常数的精细估计、新的混合变分形式下的二重网格法、变网格法等不同角度出发,对混合有限元方法的构造,收敛性分析、超逼近和超收敛等方面进行深入系统的探讨.
首先考虑了一类低阶非协调单元(包括四面体元和六面体元),用于逼近一个三维的定常不可压、完全耦合的非线性磁流体动力学方程,在磁场分别属于H1(Ω)3和H(curl;Ω)时,采用非协调混合有限元法来分析文献[14]和[27]中的方程,证明了离散问题解的存在惟一性并给出了相应未知量的最优误差估计.而且采用一种新的方法证明了离散的Poincaré-Friedrichs不等式.
其次用带约束的非协调旋转Q1元和分片常数元来逼近定常的、不可压带有阻尼项的Stokes方程的速度和压力.证明了逼近解的存在惟一性.再利用精确解和逼近解的先验估计,并恰当选择方程中出现的参数α,ν和Υ,得到了最优误差估计及超逼近结果.最后,通过插值后处理技术,导出了速度的H1-模和压力L2-模的O(h2)阶的整体超收敛.
再次研究了对流扩散方程的最小二乘非协调有限元格式及其两种修正格式,用矩形(EQrot1)元和零阶R-T元分别来逼近位移和应力,借助于单元本身的特殊性质,给出了逼近问题解的存在惟一性,得到了位移H1-模和应力H(div)-模的O(h)阶的误差估计.同时,在直角三角形网格下,采用最小二乘有限元法,用零阶的R-T元和P1元去离散该方程,给出了应力的H(div)-模及位移的H1-模的误差常数的精细估计,并给出了数值算例验证了理论分析的正确性.
而后构造了一种新的混合变分形式,用最低阶矩形协调混合元去逼近半线性反应扩散方程,利用椭圆投影的特殊性质,得到了新的混合变分全离散形式下的未知量L2-模的误差收敛阶为O(Δt+h2)比[136]传统的混合变分形式下的误差提高.了一阶.随后,采用二重网格算法,迭代两步得到了网格比H=O(h1/3)及收敛阶O(△t+h+H3)的结果,这一结果文献[136]需迭代三步才能达到.而本章迭代三步后,得到了网格比是H=O(h2/9),收敛阶为O(△t+h+H9/2).和[136]相比,达到了迭代同样步数可以增加网格比而减少计算量的目的.
最后利用非协调(EQrat1)元和零阶R-T元,对一类非线性抛物方程构造了一种新的混合变网格格式.根据该单元相容误差在能量模意义下比插值误差高一阶的特殊性质,给出了收敛性分析并得到了最优阶误差估计.本章的结果可推广到协调混合有限元逼近的任意收敛阶情形.