[论文解读] Analytical computation of moderate-degree fully-symmetric cubature rules on the triangle

Stefanos‐Aldo Papanicolopulos|arXiv (Cornell University)|Nov 16, 2011

Advanced Numerical Analysis Techniques参考文献 14被引用 4

一句话总结

本文提出了一种基于对称多项式的解析方法，用于计算三角形上最高至15次的完全对称求积规则，通过使用对称多项式简化矩方程组。该方法可实现高精度规则的精确、封闭形式求解，包括新发现的13次和15次规则，显著提升了三角形区域数值积分的精度与效率。

ABSTRACT

A method is developed to compute analytically fully symmetric cubature rules on the triangle by using symmetric polynomials to express the two kinds of invariance inherent in these rules. Rules of degree up to 15, some of them new and of good quality, are computed and presented.

研究动机与目标

• 开发一种用于计算三角形上完全对称求积规则的解析方法，避免依赖数值求解器。
• 将此前仅能数值求解的求积规则解析范围扩展至15次，尤其针对高次规则。
• 通过代数简化确保解的唯一性与精确性，从而提升求积规则的质量与可靠性。
• 提出新的、高质量的13次与15次求积规则，部分为此前未知的规则，适用于有限元与计算物理等实际应用。

提出的方法

• 该方法利用对称多项式表达完全对称求积规则的对称性性质，降低矩方程组的复杂度。
• 定义初等对称多项式时采用符号约定 (−1)^k，以简化代数运算与求解推导。
• 将矩方程组重新表述为面积坐标对称多项式的形式，形成更具结构性且可解的非线性方程组。
• 利用关联多项式方程 ∑_{j=0}^n ˜x_{n−j}x^j = 0 的根，从对称多项式不变量中恢复原始面积坐标。
• 该方法利用对称多项式基本定理，确保任意坐标对称函数均可表示为初等对称多项式的组合。
• 该方法可实现对给定规则的所有解的解析推导，确保解的唯一性，并避免数值求解器可能存在的收敛问题。

实验结果

研究问题

• RQ1 能否使用对称多项式在15次以内实现三角形上完全对称求积规则的解析计算？
• RQ2 与笛卡尔坐标或极坐标形式相比，使用对称多项式如何简化矩方程组？
• RQ3 能否通过解析方法推导出新的、高质量的13次与15次求积规则？其性能如何与现有规则比较？
• RQ4 解析计算能否保证求积规则解的唯一性与正确性，从而避免数值求解中的收敛陷阱？

主要发现

• 本文计算了最高至15次的完全对称求积规则，发现了多个新规则，包括高质量的13次与15次规则。
• 对于13次规则，提出了三种类型为[1,4,4]、共37个点的独立规则，其中一种质量为PI，三种质量为NI，表现出高精度。
• 该方法成功推导出此前仅能通过数值方法获得的规则的解析解，确保了解的唯一性与正确性。
• 使用对称多项式显著降低了矩方程组的复杂度，使系统化且简化的解析求解成为可能。
• 所计算的规则具有精确的权重与面积坐标，例如形心处权重为3.333333333333333e-01，对称点处权重为5.637138317907531e-02，其坐标如(0.1416807749137785, 0.4291596125431107, 0.4291596125431107)。
• 该方法可识别同一阶次与点配置下的多个非唯一解，例如13次时存在四种不同的[1,4,4]型规则，各自具有不同的权重与坐标集合。