[论文解读] PetIGA: High-Performance Isogeometric Analysis

Nathan Collier, Lisandro Dalcín+1 more|arXiv (Cornell University)|May 20, 2013

Advanced Numerical Analysis Techniques参考文献 10被引用 38

一句话总结

PetIGA 是一个基于 PETSc 构建的高性能等几何分析框架，可使用基于 NURBS 的 Galerkin 空间离散化高效求解偏微分方程。该框架支持复杂问题，如大变形超弹性材料与湍流 Navier-Stokes 流动，具备可扩展的求解器和用于等参数映射的高阶导数精确计算能力。

ABSTRACT

In this paper we present a code framework for isogeometric analysis which we call PetIGA. We base this framework heavily on components from PETSc, a scientific library geared toward implementing scalable solvers needed for approximating solutions to partial differential equations. Our open source library can be used to assemble matrices and vectors which, for this release, come from a Galerkin weak form, discretized with the Non-Uniform Rational B-spline (NURBS) basis. We also detail the computation of these basis functions as well as how higher-order derivatives may be pushed forward for use in isoparametric mappings. Finally, we show performance results on a variety of nonlinear and time dependent problems illustrating the flexibility and efficiency of our framework. We solve large deformation hyper-elasticity, the Cahn-Hilliard equation, the Navier-Stokes-Korteweg equation, as well as Navier-Stokes in context of turbulent flow modeled with the variational multiscale method.

研究动机与目标

• 开发一个可扩展的开源等几何分析框架，以高效求解偏微分方程。
• 将 PETSc 的可扩展求解器集成到等几何分析中，实现高性能计算。
• 实现基于 NURBS 基函数及其高阶导数的精确计算，用于等参数映射。
• 在具有挑战性的非线性和时间依赖性问题上展示该框架的能力。
• 支持先进应用，如 Cahn-Hilliard 方程与 Navier-Stokes-Korteweg 方程，结合变分多尺度建模。

提出的方法

• 该框架利用 PETSc 的可扩展线性代数与求解器基础设施，从 Galerkin 加权弱形式中组装刚度矩阵与载荷向量。
• NURBS 基函数以高精度计算，支持任意多项式次数与节点区间。
• 推导并使用 NURBS 基函数的高阶导数，以计算等参数映射中的几何雅可比行列式。
• 该框架支持基于变分原理推导的弱形式的等几何空间离散化，采用 Galerkin 方法。
• 通过隐式时间积分与牛顿-Krylov 方法，实现时间依赖性与非线性 PDE 的求解。
• 实现模块化且可扩展，支持与先进变分多尺度方法集成，用于湍流流动建模。

实验结果

研究问题

• RQ1 如何利用现有的高性能科学计算库，将等几何分析高效扩展至大规模问题？
• RQ2 基于 NURBS 的等几何离散化在具有挑战性的非线性和时间依赖性 PDE 上的性能如何？
• RQ3 NURBS 基函数的高阶导数在等参数映射中能够以多高的精度计算并有效利用？
• RQ4 该框架能否处理复杂的物理模型，如超弹性材料与相场方程？
• RQ5 使用 PETSc 在多大程度上实现了等几何分析中的可扩展性与性能提升？

主要发现

• PetIGA 通过高阶 NURBS 离散化与可扩展求解器，成功求解了大变形超弹性材料问题。
• 该框架在 Cahn-Hilliard 方程上表现出稳健性能，能够精确捕捉相分离动力学。
• Navier-Stokes-Korteweg 模拟结果表明，该框架能稳定且准确地解析多相流中的尖锐界面。
• 采用变分多尺度方法的湍流流动模拟结果，证实了该框架在复杂流体动力学问题中的求解能力。
• 性能结果验证了该框架在大规模问题上的可扩展性与高效性，证实了 PETSc 与等几何分析集成的有效性。