[论文解读] Weakening the tight coupling between geometry and simulation in isogeometric analysis: from sub- and super- geometric analysis to Geometry Independent Field approximaTion (GIFT)

Elena Atroshchenko, Gang Xu+2 more|arXiv (Cornell University)|Jan 1, 2017

Advanced Numerical Analysis Techniques参考文献 12被引用 1

一句话总结

本文提出几何无关场近似（GIFT）方法，该方法在等几何分析中将几何参数化与解近似分离。通过使用不同的样条空间（如用NURBS表示几何，用PHT样条表示场），GIFT在不改变CAD几何的前提下实现了自适应局部加密，即使在非嵌套空间下仍能实现最优收敛率，尽管标准Patch测试失效。

ABSTRACT

This paper presents an approach to generalize the concept of isogeometric analysis (IGA) by allowing different spaces for parameterization of the computational domain and for approximation of the solution field. The method inherits the main advantage of isogeometric analysis, i.e. preserves the original, exact CAD geometry (for example, given by NURBS), but allows pairing it with an approximation space which is more suitable/flexible for analysis, for example, T-splines, LR-splines, (truncated) hierarchical B-splines, and PHT-splines. This generalization offers the advantage of adaptive local refinement without the need to re-parameterize the domain, and therefore without weakening the link with the CAD model. We demonstrate the use of the method with different choices of the geometry and field splines, and show that, despite the failure of the standard patch test, the optimum convergence rate is achieved for non-nested spaces.

研究动机与目标

• 本文旨在克服传统等几何分析（IGA）的刚性，即几何与场近似使用相同的样条空间。
• 解决基于NURBS的IGA因张量积结构而难以实现局部加密的局限性。
• 目标是构建一种框架，在保持精确CAD几何的同时，允许灵活且局部增强的近似空间，以提高分析效率。
• 该方法旨在保持与CAD的紧密集成，同时实现解场的自适应h-加密，且独立于几何参数化。
• 研究旨在验证：即使场空间与几何空间非嵌套，仍可实现最优收敛率。

提出的方法

• GIFT将计算域参数化（使用NURBS或其他几何样条）与解近似空间（使用T样条、LR样条、PHT样条或截断层级B样条）解耦。
• 该方法采用弱一致的Galerkin弱形式在相同参数域上求解PDE，但对场变量使用不同的近似空间。
• 仅对解近似空间（如PHT样条）进行局部加密，不修改几何参数化，也无需重新参数化。
• 基于每个参数单元的PDE残差计算残差误差指标，并通过标记策略根据误差阈值选择需加密的单元。
• 加密遵循PHT样条规则：在交叉点和T型顶点处创建新基函数，通过将新顶点处的Bézier系数重置为零来更新现有基函数。
• 解的新控制变量通过在细化后的近似空间中重新求解PDE获得，而非通过旧控制点更新。

实验结果

研究问题

• RQ1 在等几何分析中，几何与场近似之间的解耦方法是否能在非嵌套近似空间下仍保持最优收敛率？
• RQ2 GIFT中标准Patch测试的失效是否会影响方法的一致性与收敛性？
• RQ3 是否可以在不改变原始CAD几何或无需重参数化的情况下，对解场执行自适应局部加密？
• RQ4 使用不同样条类型（如用NURBS表示几何，PHT样条表示场）如何影响数值解的精度与效率？
• RQ5 何种误差指标与标记策略可确保在GIFT框架下实现有效且稳定的局部加密？

主要发现

• GIFT在2D与3D问题中均实现了最优收敛率，即使场与几何近似空间非嵌套亦然。
• 该方法在整个过程中保持精确CAD几何，保留了等几何分析的关键优势。
• 尽管标准Patch测试失效，但因采用弱形式与适当的误差指标，方法仍保持一致且收敛。
• PHT样条的使用使得解场可实现有效的局部h-加密，且不改变底层几何参数化。
• 基于残差的误差指标 eK = ||f + Δuh||_L2(K) × hK 可靠地用于标记需加密的单元。
• 混合标记策略——先按误差大小标记前ϵ%的单元，再应用均值标记——在误差分布高度不均匀时提升了鲁棒性。