[论文解读] The LifeV library: engineering mathematics beyond the proof of concept

Luca Bertagna, Simone Deparis+3 more|arXiv (Cornell University)|Oct 18, 2017

Advanced Numerical Methods in Computational Mathematics参考文献 111被引用 24

一句话总结

本文介绍了 LifeV，一个基于 C++ 的开源有限元库，专为在 1D、2D 和 3D 中求解偏微分方程而设计，重点聚焦于实际工程应用。该库利用 MPI、Trilinos 和域分解方法，实现了在计算血流动力学、流固耦合、多孔介质以及冰盖动力学等领域的可扩展并行模拟，证明其已超越概念验证研究，具备实际部署能力。

ABSTRACT

LifeV is a library for the finite element (FE) solution of partial differential equations in one, two, and three dimensions. It is written in C++ and designed to run on diverse parallel architectures, including cloud and high performance computing facilities. In spite of its academic research nature, meaning a library for the development and testing of new methods, one distinguishing feature of LifeV is its use on real world problems and it is intended to provide a tool for many engineering applications. It has been actually used in computational hemodynamics, including cardiac mechanics and fluid-structure interaction problems, in porous media, ice sheets dynamics for both forward and inverse problems. In this paper we give a short overview of the features of LifeV and its coding paradigms on simple problems. The main focus is on the parallel environment which is mainly driven by domain decomposition methods and based on external libraries such as MPI, the Trilinos project, HDF5 and ParMetis. Dedicated to the memory of Fausto Saleri.

研究动机与目标

• 弥合计算力学中新型数值方法与其实现在工程问题中的实际部署之间的差距。
• 提供一个稳健、可扩展且支持多物理场问题先进数值方法的并行有限元框架。
• 使尖端数值技术（如单体 FSI、降基方法和 XFEM）能够应用于生产级模拟。
• 通过将该库整合到本科课程和继续教育项目中，支持教育与培训。
• 通过将代码库演进至现代 C++ 标准（C++11/14）和云兼容执行环境，确保长期可维护性和高性能。

提出的方法

• LifeV 采用基于表达式模板和泛型编程的模块化 C++ 架构，以实现高层级、可读性强且高效的有限元公式化。
• 并行化通过 MPI 与 Trilinos 生态系统的集成实现，包括 Epetra、AztecOO、ML、Belos 和 Zoltan，用于线性求解器、预处理和域分解。
• 该库以域分解方法为核心并行化策略，可在 HPC 和云平台中实现高效的负载均衡与可扩展性。
• 支持高级数值方法，包括稳定化有限元、降基方法、用于断裂介质的 XFEM 以及间断伽辽金格式。
• 代码库设计为可扩展和可重用，设有专门用于流体动力学、结构力学以及流固耦合等多物理场耦合的模块。
• 该库持续维护中，正在进行向 C++11/14 的迁移，以提升类型安全性、性能和开发人员生产力。

实验结果

研究问题

• RQ1 如何设计有限元库，使其既能支持数值方法的前沿研究，又能满足实际工程模拟的需求？
• RQ2 哪些设计模式和软件工程实践能够实现 C++ 中高效、可扩展且可维护的并行有限元计算？
• RQ3 像 LifeV 这类面向研究的库在多大程度上可成功应用于复杂的真实世界问题，如心脏力学和流固耦合？
• RQ4 Trilinos 与基于 MPI 的域分解集成如何在 HPC 和云架构上实现高性能、可扩展的解决方案？
• RQ5 现代 C++（C++11/14）在提升大规模科学软件的可维护性、可读性和性能方面发挥何种作用？

主要发现

• LifeV 已成功应用于计算血流动力学中的实际问题，包括流固耦合、心脏电机械学以及血管系统的多尺度建模。
• 该库在非定常纳维-斯托克斯问题上表现出良好的可扩展性和收敛性，其性能已在高性能计算集群上得到验证。
• Trilinos 与 MPI 的集成实现了大规模 PDE 问题的高效并行求解，域分解方法提供了有效的负载均衡与可扩展性。
• LifeV 已用于多孔介质流、冰盖动力学（正问题与反问题）以及三维问题的降基方法的生产级模拟。
• 该库已被应用于教学，包括本科课程和继续教育项目，证实其在研究原型之外的可用性。
• 向 C++11/14 的持续迁移正在提升代码清晰度、可维护性和性能，未来计划迁移到 BSD 许可证以促进采纳与集成。