[論文レビュー] High order gradient, curl and divergence conforming spaces, with an application to NURBS-based IsoGeometric Analysis

René R. Hiemstra, R.H.M. Huijsmans+1 more|arXiv (Cornell University)|Sep 9, 2012

Advanced Numerical Analysis Techniques参考文献 22被引用数 6

ひとこと要約

本稿では、NURBSに基づく高次勾配、回転、発散適合有限要素空間を導入し、離散的にガウスの定理やストークスの定理といった保存則を正確に満たす。部分単位の分割から基底関数を構築することで、曲がったメッシュ上でもストークス流れにおける正確な発散なしの速度場を保証し、多項式次数にかかわらず最適収束性と頑健性を達成する。

ABSTRACT

Conservation laws, in for example, electromagnetism, solid and fluid mechanics, allow an exact discrete representation in terms of line, surface and volume integrals. We develop high order interpolants, from any basis that is a partition of unity, that satisfy these integral relations exactly, at cell level. The resulting gradient, curl and divergence conforming spaces have the property that the conservation laws become completely independent of the basis functions. This means that the conservation laws are exactly satisfied even on curved meshes. As an example, we develop high order gradient, curl and divergence conforming spaces from NURBS - non uniform rational B-splines - and thereby generalize the compatible spaces of B-splines developed by Buffa et al.[1]. We give several examples of 2D Stokes flow calculations which result, amongst others, in a point wise divergence free velocity field.

研究の動機と目的

• 曲がったメッシュ上でも、勾配、回転、発散に関する基本的積分保存則（ガウスの定理、ストークスの定理、微分積分学の基本定理）を正確に満たす高次有限要素空間の構築。
• 等長補間スプライン空間をNURBSに一般化し、等長解析（IGA）における応用を可能にする。
• 基底関数の選択に依存しない離散的保存性を保証することで、安定性と物理的一致性を向上。
• 非圧縮性ストークスおよびナビエ＝ストークス問題の解法における有効性を示し、正確な発散なしの速度場を実現。
• ベンチマーク問題（テイラー＝コウエット流れ、リフティングキャビティ流れなど）を通じた検証により、最適収束性と頑健性を示す。

提案手法

• 任意の部分単位基底から高次補間関数を構築し、ガウスの定理、ストークスの定理、微分積分学の基本定理といった基本的積分定理を正確に満たす。
• 幾何的実体（節点、辺、面、体積）に関連する自由度を定義し、離散的設定において保存則の物理的構造を保持。
• NURBSにこのフレームワークを適用し、非一様有理Bスプラインを用いて正確な幾何表現を可能にする。
• 混合ガレルキン法を用い、速度、圧力、フラックス変数の混合有限要素空間を採用し、インフラ・スタビリティを構造的に保証。
• 非一様重みを有するNURBSパッチを用いて、曲がった領域（例えば、中空領域）を正確にモデル化。
• h-型細分化の研究を実施し、さまざまな流れ配置における速度および圧力の最適収束率を検証。

実験結果

リサーチクエスチョン

• RQ1 任意の部分単位基底から構築した高次有限要素空間は、セル単位で積分保存則を正確に満たすことができるか？
• RQ2 NURBSに基づく等長解析は、非圧縮流れにおける正確な発散なしの速度場をサポートするようにどのように拡張できるか？
• RQ3 本手法は、曲がった非一様メッシュ上でも最適収束率および離散的保存性を維持するか？
• RQ4 リフティングキャビティのような特異性を有する流れ問題において、特別な処理なしに物理的一致性を保ちながら扱えるか？
• RQ5 本手法の保存性特性が、複雑な流体動力学シミュレーションにおける数値的安定性および精度をどの程度向上させるか？

主な発見

• 2次元ストークス流れにおいて、領域全体にわたり点ごとの発散が正確にゼロとなる、正確な発散なしの速度場を達成。
• テイラー＝コウエット流れのテストケースにおいて、多項式次数p = 1からp = 4の範囲で、L2ノルムにおける速度の最適h収束率p+1が観察された。
• 圧力はすべての細分化レベルおよび多項式次数において定数であり、機械精度まで正確に保たれた。
• リフティングキャビティ流れのシミュレーションでは、粗い9×9メッシュおよび低次近似でも、組み込みの保存性特性のおかげでベンチマーク結果と優れた一致を示した。
• NURBSパッチにおける非一様重み（例：コーナーで0.5√2）が存在しても、最適収束性および正確な幾何表現が維持された。
• 物理的構造を離散的空間に組み込むことで、基底関数の選択に依存せず、インフラ・スタビリティが自動的に保証された。