不连续伽辽金中的数值求积

计算科学 正交 双曲-pde 不连续-galerkin 一体化
2021-11-29 04:48:15

我想知道哪种方法是数值积分勒让德多项式的最佳方法。我正在构建一个不连续的 Galerkin CFD 代码,其中勒让德多项式用作光谱投影中的基函数。

在以前的代码中,我习惯于使用 Gauss-Legendre 数值求积,其中只有k允许获得度多项式的精确积分的点2k1. 现在,由于这样的交点是勒让德多项式的根,我需要使用k+1积分点k-th Legendre 多项式(否则多项式在其根处求值,所有高阶积分为零)。还有其他有效的方法吗?

非常感谢您的帮助。

1个回答

一般而言,您为有限元空间使用哪个基对于求积无关紧要。例如,如果您使用次数多项式k(无论是拉格朗日基、勒让德基、单项式还是其他任何东西)那么质量矩阵的被积函数具有多项式次数2k,你需要一个高斯正交k+1点正确整合。这是因为这样的高斯求积可以积分任何次数的多项式2k+1确切地说,无论您实际上如何根据基数表示该多项式。