[论文解读] Benchmark of Polygon Quality Metrics for Polytopal Element Methods

Marco Attene, Silvia Biasotti+7 more|arXiv (Cornell University)|May 27, 2019

Advanced Numerical Methods in Computational Mathematics参考文献 20被引用 4

一句话总结

本文提出了一套系统性基准，用于评估多边形网格质量度量与多胞体有限元方法（PEM）性能之间的相关性，重点针对虚拟元方法（VEM）。通过测试各种多边形形状及度量指标——包括面积、角度、边长和内核特性——识别出在这些较弱的形状正则性准则下，PEM求解器仍能保持稳定与准确，揭示出当前理论约束过于严格，且VEM在高度不规则的多边形上仍表现可靠。

ABSTRACT

Polytopal Element Methods (PEM) allow to solve differential equations on general polygonal and polyhedral grids, potentially offering great flexibility to mesh generation algorithms. Differently from classical finite element methods, where the relation between the geometric properties of the mesh and the performances of the solver are well known, the characterization of a good polytopal element is still subject to ongoing research. Current shape regularity criteria are quite restrictive, and greatly limit the set of valid meshes. Nevertheless, numerical experiments revealed that PEM solvers can perform well on meshes that are far outside the strict boundaries imposed by the current theory, suggesting that the real capabilities of these methods are much higher. In this work, we propose a benchmark to study the correlation between general 2D polygonal meshes and PEM solvers. The benchmark aims to explore the space of 2D polygonal meshes and polygonal quality metrics, in order to identify weaker shape-regularity criteria under which the considered methods can reliably work. The proposed tool is quite general, and can be potentially used to study any PEM solver. Besides discussing the basics of the benchmark, in the second part of the paper we demonstrate its application on a representative member of the PEM family, namely the Virtual Element Method, also discussing our findings.

研究动机与目标

• 弥合多胞体有限元方法（PEM）中理论形状正则性约束与实际求解器性能之间的差距。
• 系统评估二维多边形单元各种几何特性对PEM求解器精度与条件数的影响。
• 识别出仍能确保可靠性能的宽松质量度量标准，从而实现更灵活的网格生成。
• 提供一个模块化、可重用的基准框架，适用于任何PEM求解器，包括虚拟元方法（VEM）。

提出的方法

• 该基准对每个单元评估12种多边形质量度量，并汇总为最小值、最大值和平均值，形成36种网格级度量。
• 通过参数化多边形族生成方式，逐步施加对各项度量的应力，同时结合随机多边形以避免偏差。
• 采用虚拟元方法（VEM）作为代表性PEM求解器，以解的精度和线性系统条件数作为性能指标。
• 通过顶点自由度计算L∞和L2范数下的相对误差，利用∇-投影算子近似真实误差。
• 设计了一个模块化框架，可扩展至其他PEM方法和额外的质量度量。
• 研究采用在多边形域上的模型泊松问题，测试在不同网格质量下的收敛性与稳定性。

实验结果

研究问题

• RQ1 哪些多边形质量度量与PEM求解器（尤其是虚拟元方法）的精度和条件数相关性最强？
• RQ2 在多胞体方法中，形状正则性约束可在多大程度上放宽而不降低求解器性能？
• RQ3 极端多边形形状（如非星形或高度畸变的单元）如何影响PEM求解器的收敛性与稳定性？
• RQ4 诸如内切圆/外接圆、内核大小和边长比等几何特征能否用于预测求解器的可靠性？
• RQ5 在何种质量极限之下，PEM性能会彻底失效？

主要发现

• 虚拟元方法在高度不规则的多边形上仍能保持最优收敛率和良好的条件数，即使这些多边形违反了标准的星形性与形状正则性假设。
• 诸如内切圆半径与直径之比、最小角等度量与求解器性能具有强相关性，提示其可用于网格质量评估。
• 具有复杂几何结构（如雪花状）的非星形单元仍能产生精确解，表明当前理论约束过于保守。
• 基准揭示，内切圆半径小而直径大是导致条件数差的关键预测因子，而边长均匀性与顶点分布则有助于系统稳定。
• 研究发现，内核大小与外接圆半径与内切圆半径之比是最具预测力的度量指标之一，用于判断求解器可靠性。
• 该框架能提前检测到性能退化，从而实现对问题网格单元的早期识别。