[论文解读] Discrete maximum principle and a Delaunay-type mesh condition for linear finite element approximations of two-dimensional anisotropic diffusion problems

Weizhang Huang|arXiv (Cornell University)|Aug 3, 2010

Advanced Numerical Methods in Computational Mathematics被引用 4

一句话总结

本文提出了一种针对二维各向异性扩散问题的线性有限元方法的Delaunay型网格条件，确保满足离散最大值原理（DMP）。该条件通过从反扩散矩阵导出的度量，将经典的Delaunay准则推广至各向异性扩散情形，且被证明弱于现有的各向异性非钝角条件，从而可在更粗、更灵活的网格上实现单调解。

ABSTRACT

The finite element solution of two-dimensional anisotropic diffusion problems is considered. A Delaunay-type mesh condition is developed for linear finite element approximations to satisfy a discrete maximum principle. The condition is shown to be weaker than the existing anisotropic non-obtuse angle condition. It reduces to the well known Delaunay condition for the special case with the identity diffusion matrix. Numerical results are presented to verify the theoretical findings.

研究动机与目标

• 开发一种网格条件，确保线性有限元近似在二维各向异性扩散问题中满足离散最大值原理（DMP）。
• 通过将扩散矩阵纳入几何准则，将经典Delaunay条件推广至各向异性扩散情形。
• 建立一种弱于各向异性非钝角条件的网格条件，从而在保持单调性的同时允许更灵活、更粗的网格。
• 基于全局刚度矩阵而非局部单元分析，为该条件提供理论基础。
• 通过在具有不同几何和各向异性特性的网格上进行数值实验，验证理论结果。

提出的方法

• 通过分析线性有限元公式化的全局刚度矩阵，而非依赖于局部单元刚度矩阵，推导出新的网格条件。
• 引入基于反扩散矩阵 $\mathbb{D}^{-1}$ 的度量相关角度条件，定义在 $\mathbb{D}^{-1}$-范数下边之间的夹角。
• 提出一种Delaunay型条件，要求在 $\mathbb{D}^{-1}$-范数下，对于所有单元对，共用边对角的夹角和不超过 $\pi$。
• 建立该条件为线性有限元近似在各向异性扩散问题中满足DMP的充分条件。
• 采用具有一致权重和节点的数值积分，近似弱形式，确保与理论公式的兼容性。
• 利用已知的在 $\mathbb{D}^{-1}$-度量下的Delaunay网格生成技术，通过边交换过程生成满足新条件的网格。

实验结果

研究问题

• RQ1 能否为线性有限元方法制定一种Delaunay型网格条件，以确保在各向异性扩散问题中满足离散最大值原理？
• RQ2 与现有的各向异性非钝角条件相比，该新网格条件在网格灵活性和单调性保持方面有何强度差异？
• RQ3 当扩散矩阵为单位矩阵时，该条件是否退化为经典Delaunay条件？
• RQ4 该条件能否用于指导边交换算法的设计，以在各向异性扩散中实现单调有限元解？
• RQ5 当网格不满足新条件时，数值解中欠冲和过冲的收敛行为如何？

主要发现

• 所提出的Delaunay型网格条件严格弱于各向异性非钝角条件，允许更广泛的网格类别满足DMP。
• 对于常数扩散矩阵，该条件简化为 $\alpha_{ij,\mathbb{D}^{-1}}^{K} + \alpha_{ij,\mathbb{D}^{-1}}^{K'} \leq \pi$，这推广了经典Delaunay条件。
• 数值结果表明，满足新条件的网格（如图4b、4c所示）产生的有限元解始终在0和1之间，无欠冲或过冲。
• 违反该条件的网格（如图4a、4d所示）即使在网格加密后，仍会产生显著的欠冲和过冲。
• 当单元数 $N$ 增加时，欠冲和过冲的幅度以 $O(N^{-0.5})$ 的速率缓慢衰减，证实了该条件被违反时单调性的缺失。
• 该条件具有几何可解释性，可作为边交换算法的准则，用于生成单调且具备各向异性感知能力的网格。