论文部分内容阅读
本文研究了扭曲网格上扩散方程的数值求解问题,主要内容包括各种线性精确格式的构造、分析和数值试验。线性精确方法是一种在扭曲网格上构造扩散方程高精度有限体积格式的启示性途径,它要求当扩散方程的解析解是关于自变量的线性函数并且扩散系数是常数时每一个离散步骤都是精确成立的。在此框架下我们首先利用离散泛函分析工具从理论上证明了结构四边形网格上基于网格节点双线性插值九点格式的数值稳定性和收敛性。接下来我们利用线性精确方法改进了结构四边形网格上基于变分原理构造的有限体积格式,新的格式突破了边界条件的限制,可以在统一的框架下求解各类边值问题;在此基础上我们证明了改进后的格式与支撑算子格式的某种等价性。最后我们讨论了非结构网格上扩散方程的数值求解问题,将前述基于变分原理的线性精确格式推广到非结构三角网格情形。在上述各部分内容中,我们都通过一些算例考察了线性精确格式在扭曲网格上的数值表现,获得了满意的数值结果。我们的分析及数值试验表明,在构造扭曲网格上扩散方程的高精度有限体积格式方面,线性精确方法不失为一条有效的途径。