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

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

Advanced Numerical Analysis Techniques被引用 14

一句话总结

PetIGA 是一个高性能、开源的等几何分析框架，利用 PETSc 实现基于 NURBS 离散化的偏微分方程的可扩展、并行求解。它在 4096 个核心下实现了强可扩展性，并通过显式计算雅可比矩阵，其效率比数值微分方法高出一个数量级。

ABSTRACT

We present PetIGA, a code framework to approximate the solution of partial differential equations using isogeometric analysis. PetIGA can be used to assemble matrices and vectors which come from a Galerkin weak form, discretized with Non-Uniform Rational B-spline basis functions. We base our framework on PETSc, a high-performance library for the scalable solution of partial differential equations, which simplifies the development of large-scale scientific codes, provides a rich environment for prototyping, and separates parallelism from algorithm choice. We describe the implementation of PetIGA, and exemplify its use by solving a model nonlinear problem. To illustrate the robustness and flexibility of PetIGA, we solve some challenging nonlinear partial differential equations that include problems in both solid and fluid mechanics. We show strong scaling results on up to 4096 cores, which confirm the suitability of PetIGA for large scale simulations.

研究动机与目标

• 开发一个基于 NURBS 基函数的可扩展、高性能等几何分析框架。
• 与 PETSc 集成，利用其并行线性代数和求解器支持大规模模拟。
• 通过高效的矩阵组装和雅可比矩阵计算，降低求解非线性 PDE 的计算成本。
• 通过模块化、可扩展的软件设计，促进科学计算中的代码重用和互操作性。
• 在固体力学和流体力学的复杂非线性问题上，展示框架的鲁棒性和性能。

提出的方法

• 使用 PETSc 作为基础库，提供并行线性代数、Krylov 类求解器和时间积分功能。
• 采用张量积 NURBS 基函数实现 Galerkin 加权弱形式，实现高阶、光滑的离散化。
• 支持基函数中任意连续性和多项式阶次，实现更高阶精度。
• 通过用户显式定义的残差和雅可比矩阵函数实现最佳性能，同时支持通过局部数值微分实现自动微分。
• 通过基于 MPI 的并行化实现强可扩展性，并通过 ParaView 支持可视化。
• 提供命令行接口，用于在运行时配置求解器、预条件子和问题参数。

实验结果

研究问题

• RQ1 能否构建一个高性能、开源的等几何分析框架，使其在大规模 HPC 系统上实现高效可扩展？
• RQ2 在等几何分析中，显式雅可比矩阵计算与数值微分技术相比，性能表现如何？
• RQ3 PetIGA 在现代超级计算机上对复杂非线性 PDE 的可扩展性达到何种程度？
• RQ4 该框架能否处理具有高阶连续性和复杂几何形状的固体力学和流体力学中的挑战性问题？
• RQ5 与 PETSc 的集成如何在科学模拟中同时实现高性能和灵活性？

主要发现

• PetIGA 在 4096 个核心下实现了强可扩展性，证实其适用于大规模模拟。
• 显式雅可比矩阵计算比基线局部数值微分慢约四倍，而着色有限差分法则慢逾十倍。
• 对于包含 128³ 个单元、p=3、k=2 的 3D Bratu 问题，使用局部数值微分的全局雅可比矩阵计算时间为 1.72 秒，而显式计算仅需 0.52 秒。
• 该框架成功求解了固体力学和流体力学中的复杂非线性问题，包括相场问题和大变形问题。
• 使用 PETSc 使得与高级求解器和预条件子的无缝集成成为可能，显著提升了性能和可扩展性。
• 该框架通过统一接口支持与外部库（如 Hypre、Trilinos/ML 和 MUMPS）的完全互操作性。