Hierarchical Optimization Time Integration for CFL-Rate MPM Stepping

Abstract

We propose Hierarchical Optimization Time Integration (HOT) for efficient implicit timestepping of the material point method (MPM) irrespective of simulated materials and conditions. HOT is an MPM-specialized hierarchical optimization algorithm that solves nonlinear timestep problems for large-scale MPM systems near the CFL limit. HOT provides convergent simulations out of the box across widely varying materials and computational resolutions without parameter tuning. As an implicit MPM timestepper accelerated by a custom-designed Galerkin multigrid wrapped in a quasi-Newton solver, HOT is both highly parallelizable and robustly convergent. As we show in our analysis, HOT maintains consistent and efficient performance even as we grow stiffness, increase deformation, and vary materials over a wide range of finite strain, elastodynamic, and plastic examples. Through careful benchmark ablation studies, we compare the effectiveness of HOT against seemingly plausible alternative combinations of MPM with standard multigrid and other Newton-Krylov models. We show how these alternative designs result in severe issues and poor performance. In contrast, HOT outperforms existing state-of-the-art, heavily optimized implicit MPM codes with an up to 10× performance speedup across a wide range of challenging benchmark test simulations.

Lagrangian–Eulerian multidensity topology optimization with the material point method

“…However, when the initial configuration is far from a local optimum, which is often true in static solves, Newton's method may fail to provide a proper search direction as the Hessian can be indefinite. 65,66,67 Teran et al 68 proposed a positive definite fix to project the Hessian to a symmetric positive definite form to guarantee that a descent direction can be found. This method is referred to as projected Newton (PN) throughout the paper and is applied in the static solve.…”

A momentum-conserving implicit material point method for surface tension with contact angles and spatial gradients

“…Hyde et al [2020] provide a thorough discussion of the state of the art. Our approach utilizes the particle-based MPM [ de Vaucorbeil et al 2020;Sulsky et al 1994] PIC technique, largely due to its natural ability to handle self collision [Fei et al 2018[Fei et al , 2017Guo et al 2018;, topology change Wolper et al 2020Wolper et al , 2019, diverse materials [Daviet and Bertails-Descoubes 2016;Klár et al 2016;Ram et al 2015;Schreck and Wojtan 2020;Stomakhin et al 2013;Wang et al 2020c;Yue et al 2015] as well as implicit time stepping with elasticity Fei et al 2018;Stomakhin et al 2013;Wang et al 2020b]. We additionally use the APIC method [Fu et al 2017;Jiang et al 2015 for its conservation properties and beneficial suppression of noise.…”

An implicit updated lagrangian formulation for liquids with large surface energy

“…Müller et al combine SPH with surface tracking [Müller 2009]. Lastly, although our approach is the first fully implicit MPM discretization of surface tension, there are many existing fully implicit MPM discretizations of elastoplastic materials [Fei et al 2018;Stomakhin et al 2013;Wang et al 2020].…”