[論文レビュー] 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

ひとこと要約

本稿は、楕円型PDEのNURBSに基づくアイソジオメトリック離散化において生じる剛性行列および質量行列の条件数に対する厳密な理論的上限を提供する。$h$-refinementでは、剛性行列の条件数は$h^{-2}$に比例するが、質量行列は一様に有界である。$p$-refinementでは、両行列とも指数関数的に増大し、それぞれ剛性行列の上界が$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$-refinement）および多項式次数の増加（$p$-refinement）におけるアイソジオメトリック行列の条件数の解析を行うこと。
• 2階楕円型PDEのNURBSに基づくアイソジオメトリック離散化における条件数推定の厳密な証明と数値的検証を提供すること。
• 特に多項式次数$p$の増加に伴う条件数の挙動について、先行研究における$p$-refinementの理論的上限の欠如を補完すること。
• アイソジオメトリック法における行列の条件数を定量的に評価することで、効率的な反復解法の設計と解析を支援すること。

提案手法

• Bスプライン基底関数の性質とその安定性定数を用いて、剛性行列および質量行列の最小固有値および最大固有値の上界を導出する。
• Lipschitz領域$\Omega \subset \mathbb{R}^d$（$d=2,3$）における楕円型PDEのガラーキン弱形式に基づく理論的解析を実施し、$C^{p-m}$-連続なNURBS基底関数を用いる。
• Bスプラインの安定性定数を$p2^p$として推定し、最小固有値の上限を導出する。de Boorの予想（$2^p$）が成り立つ場合、より鋭い上限が得られることに注意する。
• $p$-refinementにおける極値固有値および条件数の理論的推定を、粗いメッシュおよび細かいメッシュの両方で数値実験により検証する。
• 積分不等式、Bスプラインの性質、行列の条件数理論を組み合わせた厳密な数学的証明を用いて、上限を導出する。
• 古典的FEMの挙動と比較することで、アイソジオメトリック法の性能、特に行列の条件数の観点から文脈を整理する。

実験結果

リサーチクエスチョン

• RQ1 アイソジオメトリック離散化における$h$-refinementの際、剛性行列の条件数はメッシュサイズ$h$に対してどのようにスケーリングされるか？
• RQ2 質量行列の条件数は$h$-refinement下でどのように振る舞い、$h$に依存せずに一様に有界であるか？
• RQ3 $p$-refinementにおける剛性行列の条件数の増大率は何か？また、$p$の関数として上界を導出できるか？
• RQ4 最小固有値および最大固有値の理論的上限を導出可能か？また、Bスプラインの安定性定数にどのように依存するか？
• RQ5 数値的結果は、$p$-refinementにおける極値固有値および条件数の理論的推定とどのように一致するか？

主な発見

• $h$-refinementでは、剛性行列の条件数は上界$C h^{-2}$で抑えられ、古典的FEMと一致し、数値的にも確認されている。
• $h$-refinementでは、質量行列の条件数は$h$に依存せず一様に有界であり、これは反復解法にとって好都合である。
• $p$-refinementでは、剛性行列の条件数は指数関数的に増大し、上界として$p^{2d+2}4^{pd}$（$d=2,3$）で抑えられる。
• $p$-refinementでは、質量行列の条件数も指数関数的に増大し、上界として$p^{2d}4^{pd}$（$d=2,3$）で抑えられる。
• 数値的結果は理論的推定を裏付けている：$p$が増加するに従い、剛性行列の最小固有値は減少するが、細かいメッシュでは最大固有値は$p$に比例して増大する。
• $p > 10$では、丸め誤差が数値的精度に影響を与え、実験では$p=10$を超える計算は信頼性を欠く。