[论文解读] Condition number estimates for matrices arising in NURBS based isogeometric discretizations of elliptic partial differential equations

K.P.S. Gahalaut, Satyendra Tomar+1 more|arXiv (Cornell University)|Jun 26, 2014

Advanced Numerical Analysis Techniques被引用 8

一句话总结

本文为基于NURBS的等几何离散化椭圆PDE问题中刚度矩阵和质量矩阵的条件数提供了严格的理论界。研究证明，在$h$-型加密下，刚度矩阵的条件数按$h^{-2}$增长，而质量矩阵则保持一致有界；在$p$-型加密下，两类矩阵均表现出指数增长，其上界分别为$p^{2d+2}4^{pd}$（刚度）和$p^{2d}4^{pd}$（质量），其中$d=2,3$。这些结果对分析等几何分析中迭代求解器的性能至关重要。

ABSTRACT

We derive bounds for the minimum and maximum eigenvalues and the spectral condition number of matrices for isogeometric discretizations of elliptic partial differential equations in an open, bounded, simply connected Lipschitz domain $Ω\subset \mathbb{R}^d$, $d\in\{2,3\}$. We consider refinements based on mesh size $h$ and polynomial degree $p$ with maximum regularity of spline basis functions. For the $h$-refinement, the condition number of the stiffness matrix is bounded above by a constant times $ h^{-2}$ and the condition number of the mass matrix is uniformly bounded. For the $p$-refinement, the condition number grows exponentially and is bounded above by $p^{2d+2}4^{pd}$ and $p^{2d}4^{pd}$ for the stiffness and mass matrices, respectively. Rigorous theoretical proofs of these estimates and supporting numerical results are provided.

研究动机与目标

• 推导等几何分析中刚度与质量矩阵的最小与最大特征值及谱条件数的理论界。
• 分析在$h$-型加密（网格加密）与$p$-型加密（多项式阶数加密）下等几何矩阵的条件性。
• 为基于NURBS的二阶椭圆PDE等几何离散化提供条件数估计的严格证明与数值验证。
• 解决先前文献中缺乏$p$-型加密下条件数界的理论空白，特别是针对多项式阶数$p$增加时条件数行为的研究。
• 通过量化等几何方法中矩阵的条件性，支持高效迭代求解器的设计与分析。

提出的方法

• 利用B样条基函数的性质及其稳定性常数，推导刚度与质量矩阵最小与最大特征值的上界。
• 基于Lipschitz区域$\Omega \subset \mathbb{R}^d$（$d=2,3$）中椭圆PDE的Galerkin弱形式，采用$C^{p-m}$-连续的NURBS基函数进行理论分析。
• 利用B样条的稳定性常数（估计为$p2^p$）来界定最小特征值；指出若de Boor猜想（$2^p$）成立，则可获得更紧的界。
• 在最粗与最细的网格上进行数值实验，以验证$p$-型加密下极值特征值与条件数的理论估计。
• 采用涉及积分不等式、B样条性质与矩阵条件数理论的严格数学证明，推导出相关界。
• 将结果与经典有限元法（FEM）的行为进行比较，以定位等几何方法在矩阵条件性方面的性能表现。

实验结果

研究问题

• RQ1 在等几何离散化的$h$-型加密中，刚度矩阵的条件数如何随网格尺寸$h$变化？
• RQ2 在$h$-型加密下，质量矩阵的条件数行为如何？是否一致有界？
• RQ3 在$p$-型加密下，刚度矩阵的条件数增长速率如何？能否以$p$表示其上界？
• RQ4 能否为等几何矩阵的最小与最大特征值推导出理论界？其与B样条稳定性常数的关系如何？
• RQ5 在$p$-型加密下，数值结果与极值特征值及条件数的理论估计相比如何？

主要发现

• 在$h$-型加密下，刚度矩阵的条件数上界为$C h^{-2}$，与经典FEM一致，且经数值验证。
• 在$h$-型加密下，质量矩阵的条件数一致有界，与$h$无关，这对迭代求解器有利。
• 在$p$-型加密下，刚度矩阵的条件数呈指数增长，其上界为$p^{2d+2}4^{pd}$，其中$d=2,3$。
• 在$p$-型加密下，质量矩阵的条件数也呈指数增长，其上界为$p^{2d}4^{pd}$，同样适用于$d=2,3$。
• 数值结果验证了理论估计：随着$p$增加，刚度矩阵的最小特征值减小，而最大特征值在精细网格上随$p$线性增长。
• 当$p > 10$时，舍入误差开始影响数值精度，导致实验中可靠计算的$p$值上限约为10。