[论文解读] 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.

研究动机与目标

• 开发高阶有限元空间，使其在曲面网格上精确满足梯度、旋度和散度的积分守恒定律。
• 将相容B样条空间推广至NURBS，以适用于等几何分析（IGA）。
• 确保离散守恒性质独立于基函数选择，从而提升稳定性与物理一致性。
• 通过数值实验验证该方法在求解不可压缩斯托克斯与纳维-斯托克斯问题中的有效性，实现严格无散的速度场。
• 通过基准问题验证该方法，包括泰勒-库埃特流和顶盖驱动腔流，结果表明其具有最优收敛率和鲁棒性。

提出的方法

• 从任意统一划分基函数构造高阶插值函数，确保精确满足基本积分定理（高斯定理、斯托克斯定理、微积分基本定理）。
• 将自由度与几何实体（节点、边、面、体）关联，保持离散设置中守恒定律的物理结构。
• 将该框架应用于NURBS，将相容B样条空间扩展至非均匀有理B样条，以实现精确几何表示。
• 采用混合伽辽金弱形式，结合速度、压力与通量变量的混合有限元空间，通过构造保证inf-sup稳定性。
• 使用具有非均匀权重的NURBS网格片精确建模曲面域，包括环形区域等复杂几何形状。
• 实施h型超细化研究，验证在各种流动构型中速度与压力的最优收敛率。

实验结果

研究问题

• RQ1 能否从任意统一划分基函数构造高阶有限元空间，使其在单元层面精确满足积分守恒定律？
• RQ2 如何将基于NURBS的等几何分析扩展至支持不可压缩流中严格无散的速度场？
• RQ3 所提方法在曲面、非均匀网格上是否保持最优收敛率与离散守恒性质？
• RQ4 该方法是否无需特殊处理即可处理流动问题中的奇点（如顶盖驱动腔流），同时保持物理一致性？
• RQ5 该方法的守恒性质在复杂流体动力学模拟中在多大程度上提升了数值稳定性与精度？

主要发现

• 该方法在2D斯托克斯流中实现了严格无散的速度场，全场点的散度精确为零。
• 在泰勒-库埃特流测试案例中，速度在L2范数下对多项式阶次p = 1至p = 4均观测到p+1阶的最优h-收敛率。
• 压力在所有细化级别与多项式阶次下均为常数，且精度达到机器精度。
• 尽管使用粗化的9×9网格与低阶近似，顶盖驱动腔流模拟仍与基准结果高度一致，归因于其内在的守恒性质。
• 即使在NURBS网格片中使用非均匀权重（如角点处权重为0.5√2），方法仍保持最优收敛率与精确几何表示。
• 该方法通过将物理结构编码至离散空间，自动保证inf-sup稳定性，且不依赖于基函数的选择。