[论文解读] lifex: a flexible, high performance library for the numerical solution of complex finite element problems

Pasquale Claudio Africa|arXiv (Cornell University)|Jul 29, 2022

Advanced Numerical Methods in Computational Mathematics被引用 10

一句话总结

lifex 是一个基于 deal.II 的开源、高性能 C++ 库，专为求解复杂的多物理场、多尺度和多域有限元问题而设计，重点在于易用性、可扩展性和性能。它可在从笔记本电脑到百亿亿次 HPC 系统的各类架构上实现高效、可移植的仿真，达到接近理想的并行加速比，并支持任意有限元、复杂耦合以及运行时参数调优，且代码量极少。

ABSTRACT

Numerical simulations are ubiquitous in mathematics and computational science. Several industrial and clinical applications entail modeling complex multiphysics systems that evolve over a variety of spatial and temporal scales. This study introduces the design and capabilities of lifex, an open source C++ library for high performance finite element simulations of multiphysics, multiscale, and multidomain problems. lifex meets the emerging need for versatile, efficient computational tools that are easily accessed by users and developers. We showcase its flexibility and effectiveness on a number of illustrative examples and advanced applications of use and demonstrate its parallel performance up to thousands of cores.

研究动机与目标

• 应对在不同科学领域中对多功能、高效且用户友好的计算框架日益增长的需求，以模拟复杂的多物理场系统。
• 通过在统一环境中支持虚拟实验、敏感性分析、优化和不确定性量化，将软件开发与耗时的仿真工作流解耦。
• 提供一个高性能、可移植且可扩展的库，简化复杂有限元求解器的实现，同时保持数值精度和效率。
• 支持与标准文件格式（如 VTK）、网格类型（六面体和四面体）以及第三方库（PETSc、Trilinos、Boost）的无缝互操作性。
• 通过良好文档化、持续维护的代码库，推动可持续的开放科学生态系统，实现可重现、可重用和透明的虚拟实验。

提出的方法

• 利用现代 C++17 特性及简洁抽象的编程接口，简化复杂 PDE 及其耦合方案的定义。
• 以 deal.II 有限元库作为核心后端，复用其在自适应、高阶和 hp 精细化离散化方面的稳健基础设施。
• 采用基于 MPI 的分布式内存并行计算，通过域分解和负载均衡在数千个核心上实现高可扩展性。
• 支持运行时选择线性代数后端（PETSc、Trilinos 或 deal.II 原生）和求解器，实现无需修改代码的性能调优。
• 与 VTK 集成以支持可视化和网格 I/O，并支持以标准格式导入/导出网格和求解数据。
• 使用高级部署工具（Spack、Docker、mk、life x-env）确保在 Linux、Windows 和 macOS 系统上的可移植性和可重现性。

实验结果

研究问题

• RQ1 如何设计一个有限元库，以在复杂多物理场问题中实现高性能、易用性和可扩展性的平衡？
• RQ2 现代 C++ 抽象和模块化设计在多大程度上可以降低实现复杂性，同时不牺牲计算效率？
• RQ3 在异构 HPC 架构上，灵活且用户友好的有限元框架能实现多高的并行可扩展性？
• RQ4 该库在真实科学应用中，对任意有限元、复杂耦合方案和运行时参数配置的支持效率如何？
• RQ5 一个面向研究的库能否在保持对新手用户和专家开发者均友好同时，实现生产级性能和可移植性？

主要发现

• lifex 在最多数千个核心上实现了接近理想的并行加速比，展现出极高的计算效率和强可扩展性。
• 在性能关键组件中，库的计算开销可忽略不计，证实其抽象层不会影响执行速度。
• 通过 deal.II 后端支持任意有限元，可实现对各类 PDE 类型和离散化策略的精确模拟。
• 与 VTK 及标准网格格式的互操作性确保了与可视化和前后处理工具之间的无缝数据交换。
• 使用现代 C++17 和可移植部署栈（Docker、Spack、mk）可在 Linux、Windows 和 macOS 平台上实现可靠、可重现的构建。
• 库的模块化设计和运行时可配置性使得复杂求解器和耦合方案可用极少代码实现，显著提升了开发效率和代码可维护性。