[論文レビュー] PetIGA: High-Performance Isogeometric Analysis

Nathan Collier, Lisandro Dalcín+1 more|arXiv (Cornell University)|May 20, 2013

Advanced Numerical Analysis Techniques参考文献 10被引用数 38

ひとこと要約

PetIGA は、PETSc に基づく高性能な等幾何解析フレームワークであり、NURBS を用いたガレルキン離散化により偏微分方程式を効率的に解ける。大変形超弾性や乱流ナビエ＝ストークス流れといった複雑な問題を、スケーラブルなソルバーと等パラメトリック写像の高階微分の高精度計算で扱える。

ABSTRACT

In this paper we present a code framework for isogeometric analysis which we call PetIGA. We base this framework heavily on components from PETSc, a scientific library geared toward implementing scalable solvers needed for approximating solutions to partial differential equations. Our open source library can be used to assemble matrices and vectors which, for this release, come from a Galerkin weak form, discretized with the Non-Uniform Rational B-spline (NURBS) basis. We also detail the computation of these basis functions as well as how higher-order derivatives may be pushed forward for use in isoparametric mappings. Finally, we show performance results on a variety of nonlinear and time dependent problems illustrating the flexibility and efficiency of our framework. We solve large deformation hyper-elasticity, the Cahn-Hilliard equation, the Navier-Stokes-Korteweg equation, as well as Navier-Stokes in context of turbulent flow modeled with the variational multiscale method.

研究の動機と目的

• 偏微分方程式を効率的に解くためのスケーラブルでオープンソースの等幾何解析フレームワークの開発。
• 高性能コンピューティングを実現するため、PETSc のスケーラブルなソルバーを等幾何解析に統合すること。
• 等パラメトリック写像に用いる NURBS 基底関数およびその高階微分の正確な計算を可能とすること。
• 非線形および時間に依存する問題に対するフレームワークの能力を実証すること。
• 変分的マルチスケールモデリングを用いた Cahn-Hilliard 方程式や Navier-Stokes-Korteweg 方程式といった高度な応用を支援すること。

提案手法

• フレームワークは、PETSc のスケーラブルな線形代数およびソルバーインfraストラクチャを活用し、ガレルキン弱形式から剛性行列と荷重ベクトルを構築する。
• NURBS 基底関数は、任意の多項式次数および knot span をサポートする高精度で計算される。
• NURBS 基底関数の高階微分が導出され、等パラメトリック写像の幾何的ヤコビ行列の計算に用いられる。
• 変分原理から導かれる弱形式を用いたガレルキン法による等幾何離散化がサポートされる。
• 陰的時間積分法およびニュートン＝クリコフ法を用いて、時間に依存する非線形 PDE の解法が可能となる。
• 実装はモジュラで拡張可能であり、乱流流れモデリングに向けた高度な変分的マルチスケール手法との統合を可能としている。

実験結果

リサーチクエスチョン

• RQ1 既存の高性能科学計算ライブラリを用いて、等幾何解析を大規模問題にスケーラブルに拡張する方法は何か？
• RQ2 NURBS を用いた等幾何離散化は、非線形および時間に依存する PDE に対してどの程度の性能を示すか？
• RQ3 NURBS 基底関数の高階微分は、等パラメトリック写像にどの程度正確に計算され、活用されるか？
• RQ4 フレームワークは、超弾性性や相場面方程式のような複雑な物理モデルを処理できるか？
• RQ5 PETSc の使用は、等幾何解析におけるスケーラビリティとパフォーマンスをどの程度向上させるか？

主な発見

• PetIGA は、高次 NURBS 離散化とスケーラブルなソルバーを用いて、大変形超弾性問題を効果的に解いている。
• フレームワークは Cahn-Hilliard 方程式に対して堅牢な性能を示し、相分離ダイナミクスを正確に捉えている。
• Navier-Stokes-Korteweg のシミュレーションでは、多相流における鋭い界面の安定的かつ高精度な解像が確認された。
• 変分的マルチスケール法を用いた乱流シミュレーションは、フレームワークの複雑な流体力学への対応能力を裏付けた。
• パフォーマンス結果は、大規模問題においてスケーラビリティと効率性が確認され、PETSc と等幾何解析の統合が正当化された。