[論文レビュー] Using exact geometry information in finite element computations

Luca Heltai, Wolfgang Bangerth+2 more|arXiv (Cornell University)|Jan 1, 2019

Computational Geometry and Mesh Generation被引用数 3

ひとこと要約

本論文は、高度なシミュレーション技術に必要な2つの基本的な幾何クエリを特定することで、有限要素計算に正確な幾何情報を取り込むフレームワークを提案する。この研究では、点が幾何内にあるか否かの判定（point-in-geometry）と、最近接点への投影（closest point projection）という2つのプリミティブが、標準的な産業用幾何表現においても効率的にサポート可能であることを示し、メッシュベースの近似に依存せず、より優れた自己適合メッシュの細分化と高次法を実現可能であることを明らかにする。

ABSTRACT

The traditional workflow in continuum mechanics simulations is that a geometry description -- obtained using Constructive Solid Geometry or Computer Aided Design tools -- forms the input for a mesh generator. The mesh is then used as the sole input for the finite element, finite volume, and finite difference solver, which at this point no longer has access to the original geometry. However, many more modern techniques -- for example, adaptive mesh refinement and the use of higher order geometry approximation methods -- really do need information about the underlying geometry to realize their full potential. We have undertaken an exhaustive study of where typical finite element codes use geometry information, with the goal of determining what information geometry tools would have to provide. Our study shows that all geometry needs inside the simulators can be satisfied by just two ``primitives'': elementary queries posed by the simulation software to the geometry description. We then show that it is possible to provide these primitives in all of the commonly used ways in which geometries are described in common industrial workflows. We illustrate our solutions using examples from adaptive mesh refinement for complex geometries.

研究の動機と目的

• 現代の有限要素ソルバが正確な幾何情報を利用するために必要な最小限の幾何クエリのセットを特定すること。
• メッシュ生成後に幾何データが破棄される従来のFEMワークフローの制限を是正することにより、高度なシミュレーション技術の実現を妨げる要因を解消すること。
• 元の幾何記述へのアクセスを維持することで、高次幾何近似と自己適合メッシュ細分化を可能にすること。
• 提案された幾何プリミティブが、CADやCSGワークフローで一般的に使われる産業用幾何表現と互換性を持つことを示すこと。

提案手法

• シミュレーション中に幾何情報がアクセスされるすべてのケースを抽出するために、有限要素コードベースを包括的に分析する。
• FEMにおけるすべての幾何依存演算が、2つの基本的クエリに還元可能であることを同定する：点が幾何内にあるか否かの判定（point-in-geometry）と、最近接点への投影（closest point projection）。
• NURBS、CAD、CSGなど、下位の幾何表現に依存しない形で、これらの2つのプリミティブを公開する幾何抽象化レイヤーを設計する。
• 自己適合メッシュ細分化を用いた複雑な幾何形状におけるフレームワークの実装と検証を行い、一貫した性能向上と精度向上を示す。
• 既存の有限要素ソルバに変更を加えずに、幾何プリミティブを統合することで、後方互換性と効率性を確保する。

実験結果

リサーチクエスチョン

• RQ1 自己適合メッシュ細分化や高次法といった高度な有限要素技術を支援するために、最小限で十分な幾何クエリのセットは何か？
• RQ2 従来のメッシュベースのパイプラインを損なわず、正確な幾何情報が既存の有限要素ワークフロー内でどのように効率的にアクセス・利用できるか？
• RQ3 提案された幾何プリミティブは、CAD、CSG、NURBSといった多様な産業用幾何表現において実装可能か？
• RQ4 正確な幾何情報を統合することで、複雑な幾何形状における自己適合メッシュ細分化の精度と収束性はどの程度向上するか？

主な発見

• 有限要素シミュレーションにおけるすべての幾何依存演算は、単に2つの基本的クエリに還元可能である：点が幾何内にあるか否かの判定（point-in-geometry）と、最近接点への投影（closest point projection）。
• 提案された幾何プリミティブは、CAD、CSG、NURBSを含むすべての主要な産業用幾何表現と互換性があり、広範な適用可能性を有する。
• 正確な幾何情報の統合により、特に高歪みや複雑なトポロジーを持つ領域において、より正確で効率的な自己適合メッシュ細分化が可能になる。
• フレームワークは計算効率を維持しており、コア有限要素ソルバの変更を必要とせず、既存のシミュレーションパイプラインへの実用的導入を可能にする。