[论文解读] 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$的双二次样条，实现四阶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.

研究动机与目标

• 在无需几何限制的非结构化四边形网格上实现适用于分析的等几何离散化。
• 解决从CAD数据重建适用于高阶PDE的分析友好模型的挑战。
• 开发一种低阶样条构造方法，支持四阶问题的$H^2$-非协调方法。
• 确保样条基函数非奇异、局部支持，并具备类似B样条的性质，如单位分解和非负性。
• 通过Bézier提取提供一种实用且显式的实现框架。

提出的方法

• 在完全非结构化的四边形网格（包括具有特殊顶点的网格）上构建双二次样条。
• 在所有顶点和大多数边上强制实现$C^1$连续性，其余边实现近乎-$C^1$连续性。
• 采用基于Bézier提取的框架，显式表示样条基函数，以支持高效实现。
• 确保基函数满足关键性质：非负性、局部支持和单位分解。
• 设计样条空间为$H^2$-非协调，适用于四阶PDE。
• 通过数值实验验证最优收敛率和良好条件的系统。

实验结果

研究问题

• RQ1 在非结构化四边形网格上的双二次样条能否在保持四阶问题良好条件基函数的同时实现近乎-$C^1$连续性？
• RQ2 所提出的样条构造在高阶PDE的逼近精度和条件性方面表现如何？
• RQ3 该方法在无几何限制的情况下，对包含任意数量特殊顶点的任意四边形网格的处理能力如何？
• RQ4 Bézier提取框架是否能实现样条基函数的高效且稳定的实现？
• RQ5 该方法在各种网格类型下对四阶问题的收敛行为如何？

主要发现

• 所提出的近乎-$C^1$双二次样条在非结构化四边形网格上对四阶PDE实现了最优逼近阶。
• 样条基函数非奇异，并在数值测试中表现出良好条件性。
• 该方法支持任意非结构化网格，包括具有任意数量特殊顶点的网格。
• 基函数保持了类似B样条的性质，如非负性、局部支持和单位分解。
• 数值实验确认了在各种测试问题中均实现最优收敛率和鲁棒的条件性。
• Bézier提取框架实现了样条基函数的直接且高效的实现。