[論文レビュー] Almost-$C^1$ splines: Biquadratic splines on unstructured quadrilateral meshes and their application to fourth order problems

Thomas Takacs, Deepesh Toshniwal|arXiv (Cornell University)|Jan 27, 2022

Advanced Numerical Analysis Techniques被引用数 4

ひとこと要約

本稿では、非構造的四角形メッシュ上にほぼ-$C^1$な双二次スプラインを導入し、4次偏微分方程式（PDE）に対する$H^2$-非適合なアイソジオメトリック離散化を可能にする。この手法は、すべての辺および頂点でほぼ$C^1$連続性を達成し、任意のメッシュ配置をサポートしており、最適な近似特性を持つ、良好に条件付けられた、Bスプラインに類似た基底を提供する。

ABSTRACT

Isogeometric Analysis generalizes classical finite element analysis and intends to integrate it with the field of Computer-Aided Design. A central problem in achieving this objective is the reconstruction of analysis-suitable models from Computer-Aided Design models, which is in general a non-trivial and time-consuming task. In this article, we present a novel spline construction, that enables model reconstruction as well as simulation of high-order PDEs on the reconstructed models. The proposed almost-$C^1$ are biquadratic splines on fully unstructured quadrilateral meshes (without restrictions on placements or number of extraordinary vertices). They are $C^1$ smooth almost everywhere, that is, at all vertices and across most edges, and in addition almost (i.e. approximately) $C^1$ smooth across all other edges. Thus, the splines form $H^2$-nonconforming analysis-suitable discretization spaces. This is the lowest-degree unstructured spline construction that can be used to solve fourth-order problems. The associated spline basis is non-singular and has several B-spline-like properties (e.g., partition of unity, non-negativity, local support), the almost-$C^1$ splines are described in an explicit Bézier-extraction-based framework that can be easily implemented. Numerical tests suggest that the basis is well-conditioned and exhibits optimal approximation behavior.

研究の動機と目的

• 幾何的制限なしに、非構造的四角形メッシュ上で解析に適したアイソジオメトリック離散化を可能にすること。
• 高次偏微分方程式（PDE）のための、CADデータから解析に適したモデルを再構築する課題に対処すること。
• 4次問題に適した$H^2$-非適合法をサポートする低次のスプライン構成を開発すること。
• スプライン基底が非特異的で、局所的にサポートされ、分割単位および非負性といったBスプラインに類似た性質を有することを保証すること。
• Bézier抽出を用いた実装フレームワークを明示的かつ実用的に提供すること。

提案手法

• 完全に非構造的な四角形メッシュ上に双二次スプラインを構築し、特異頂点を含むものも含む。
• すべての頂点および大多数の辺で$C^1$連続性を強制し、残りの辺ではほぼ$C^1$連続性を達成する。
• Bézier抽出に基づくフレームワークを用いて、スプライン基底を明示的に表現し、効率的な実装を可能にする。
• 基底が非負性、局所的サポート、分割単位といった重要な性質を満たすことを保証する。
• スプライン空間を$H^2$-非適合に設計し、4次PDEに適したものとする。
• 最適な収束率および良好に条件付けられたシステムであるかを数値的に検証する。

実験結果

リサーチクエスチョン

• RQ1 非構造的四角形メッシュ上に定義された双二次スプラインは、4次問題に対して良好に条件付けられた基底を維持しながら、ほぼ$C^1$連続性を達成できるか？
• RQ2 提案されたスプライン構成は、高次PDEの近似精度および条件数の観点でどの程度の性能を示すか？
• RQ3 幾何的制限なしに、任意の数の特異頂点を含む任意の四角形メッシュをどの程度適応できるか？
• RQ4 Bézier抽出フレームワークは、スプライン基底の効率的かつ安定した実装を可能にするか？
• RQ5 さまざまなメッシュタイプにおいて、4次問題に対する収束挙動はどのようになるか？

主な発見

• 提案されたほぼ-$C^1$な双二次スプラインは、非構造的四角形メッシュ上での4次PDEに対して最適な近似次数を達成する。
• スプライン基底は非特異的であり、数値実験で良好に条件付けられた挙動を示す。
• 本手法は、任意の非構造的メッシュ、任意の数の特異頂点を含むものも含めてサポートする。
• 基底は非負性、局所的サポート、分割単位といったBスプラインに類似た性質を維持する。
• 数値実験により、最適な収束率およびさまざまなテスト問題におけるロバストな条件付けが確認された。
• Bézier抽出フレームワークにより、スプライン基底の直感的で効率的な実装が可能となった。