[論文レビュー] Improved approximation of eigenvalues in isogeometric methods for multi-patch geometries and Neumann boundaries

Thomas Horger, Alessandro Reali+2 more|arXiv (Cornell University)|Jan 23, 2017

Advanced Numerical Analysis Techniques被引用数 2

ひとこと要約

本稿では、マルチパッチ幾何構造における界面での弱い連続性と法線方向微分のペナルティを課すことで、高次のペナルティ技法をアイソジオメトリックモアター法に導入する。モアター結合と法線方向微分のペナルティを組み合わせることで、第四階の問題および固有値問題に対する安定的かつ高次の空間離散化が可能となり、第二階の固有値問題におけるスペクトル汚染（「外れ値」）が効果的に排除される。複雑な幾何構造でも精度を維持する。

ABSTRACT

We present a systematic study on higher-order penalty techniques for isogeometric mortar methods. In addition to the weak-continuity enforced by a mortar method, normal derivatives across the interface are penalized. The considered applications are fourth order problems as well as eigenvalue problems for second and fourth order equations. The hybrid coupling enables the discretization of fourth order problems in a multi-patch setting as well as a convenient implementation of natural boundary conditions. For second order eigenvalue problems, the pollution of the discrete spectrum - typically referred to as 'outliers' - can be avoided. Numerical results illustrate the good behaviour of the proposed method in simple systematic studies as well as more complex multi-patch mapped geometries for linear elasticity and Kirchhoff plates.

研究の動機と目的

• アイソジオメトリック解析における第二階問題の離散固有値近似におけるスペクトル汚染（「外れ値」）を解消すること。
• 複雑な界面を持つマルチパッチ幾何構造上での第四階問題の安定的かつ高精度な離散化を可能とすること。
• 連続性および法線方向微分制約を弱く強制する高次のペナルティ技法の体系的フレームワークを提供すること。
• アイソジオメトリック法における自然境界条件（ノイマン境界条件）の実装を容易にすること。

提案手法

• Nitsche法を用いたペナルティ項を備えたモアター法により、パッチ界面を越えて弱い連続性を強制する。
• 界面を越える法線方向微分に対して高次のペナルティ項を導入し、法線方向フラックスの連続性を向上させる。
• モアター連続性と法線方向微分ペナルティをハイブリッド的に組み合わせることで、第四階問題の一貫性のある離散化が可能となる。
• 本手法は第二階および第四階の固有値問題に適用され、特にスペクトル汚染の回避に注力されている。
• 本手法はマルチパッチ幾何構造をサポートしており、線形弾性およびキルヒホッフ板問題に対して検証されている。
• 境界フラックスの弱い強制を通じて、境界条件が自然に組み込まれる。

実験結果

リサーチクエスチョン

• RQ1 高次のペナルティ項は、第二階問題のアイソジオメトリック固有値近似におけるスペクトル汚染を効果的に排除できるか？
• RQ2 複雑な界面を持つマルチパッチアイソジオメトリック設定において、第四階問題をどのように高精度に離散化できるか？
• RQ3 法線方向微分をペナルティ化することで、アイソジオメトリックモアター法の精度および安定性がどの程度向上するか？
• RQ4 本手法は、複雑なマッピングを伴うマルチパッチ幾何構造上でも高次の収束率を維持できるか？
• RQ5 ハイブリッドモアター-ペナルティ法は、条件数およびスペクトル特性の観点から、標準的なアイソジオメトリック法と比較してどのように異なるか？

主な発見

• 提案手法は、第二階問題の離散固有値近似におけるスペクトル汚染（「外れ値」）を効果的に排除した。
• 本手法は、複雑なマッピングを伴うマルチパッチ幾何構造上でも、第四階問題に対して安定的かつ高精度な解を得た。
• 数値実験により、線形弾性問題およびキルヒホッフ板問題の両方で高次の収束率が確認された。
• ハイブリッドモアター-ペナルティ式は、ノイマン境界条件の自然な実装を可能とした。
• 体系的な研究および複雑な幾何構造のテストを通じて、離散系のロバストな性能と良好な条件数が示された。