[论文解读] Scalable and modular material point method for large-scale simulations

Krishna Kumar, Jeffrey Salmond+5 more|arXiv (Cornell University)|Sep 29, 2019

Fluid Dynamics Simulations and Interactions被引用 5

一句话总结

本文提出了一种可扩展且模块化的材料点法（MPM）框架，专为大规模连续介质力学模拟而设计，通过面向对象编程和动态负载均衡技术，高效处理滑坡和混凝土流动等大变形问题。该方法结合了混合欧拉-拉格朗日格式与优化的数据结构及并行化技术，在不同计算平台上实现了高性能与高可移植性。

ABSTRACT

In this paper, we describe a new scalable and modular material point method (MPM) code developed for solving large-scale problems in continuum mechanics. The MPM is a hybrid Eulerian-Lagrangian approach, which uses both moving material points and computational nodes on a background mesh. The MPM has been successfully applied to solve large-deformation problems such as landslides, failure of slopes, concrete flows, etc. Solving these large-deformation problems result in the material points actively moving through the mesh. Developing an efficient parallelisation scheme for the MPM code requires dynamic load-balancing techniques for both the material points and the background mesh. This paper describes the data structures and algorithms employed to improve the performance and portability of the MPM code. An object-oriented programming paradigm is adopted to modularise the MPM code. The Unified Modelling Language (UML) diagram of the MPM code structure is shown in Figure 1.

研究动机与目标

• 解决模拟滑坡和坡体失稳等大变形问题所面临的计算挑战。
• 开发一种在并行架构上可高效扩展的高性能 MPM 代码。
• 通过面向对象设计范式提升代码的可移植性与可维护性。
• 在并行 MPM 模拟中实现对材料点和背景网格的动态负载均衡策略。
• 通过模块化软件架构，实现对复杂地质力学和流体行为的可扩展模拟。

提出的方法

• MPM 采用混合欧拉-拉格朗日格式，利用材料点追踪材料属性，通过背景网格求解控制方程。
• 采用面向对象编程模型对代码进行模块化设计，UML 图表展示了类结构，以支持可扩展性与可维护性。
• 对材料点和计算节点均应用动态负载均衡技术，确保大变形模拟过程中工作负载均衡。
• 通过可移植的数据结构设计，支持并行执行中的高效内存访问与缓存一致性。
• 通过清晰的类接口支持本构模型、时间积分格式和求解器的模块化实现。
• 通过内存聚合与材料点与网格节点间高效数据移动实现性能优化。

实验结果

研究问题

• RQ1 如何使材料点法在地质力学与连续介质力学的大规模模拟中兼具可扩展性与模块化？
• RQ2 哪些面向对象设计模式能够在保持高性能的同时，实现可维护且可扩展的 MPM 代码？
• RQ3 在并行 MPM 模拟中，如何有效应用对材料点和背景网格的动态负载均衡？
• RQ4 在大规模 MPM 应用中，哪些数据结构与算法最能实现高性能与可移植性？
• RQ5 该模块化框架如何在统一的仿真环境中支持多样化本构模型与求解器的集成？

主要发现

• 所提出的 MPM 框架通过有效实现材料点与网格节点的动态负载均衡，在并行架构上实现了高可扩展性。
• 面向对象设计支持物理模型、时间积分器与求解器的模块化实现，显著提升了代码的可维护性与可扩展性。
• 使用优化的数据结构改善了内存访问模式，有助于提升缓存性能与计算效率。
• 该框架能够稳定且准确地模拟复杂地质力学问题，如滑坡与混凝土流动等大变形行为。
• 模块化与可移植的设计使代码可在不同 HPC 平台间轻松适配，仅需极少重构。
• 基于 UML 的软件架构为通过新增物理模型或数值方法扩展框架提供了清晰蓝图。