Hierarchical Optimization Time Integration for CFL-rate MPM Stepping

by   Xinlei Wang, et al.

We propose Hierarchical Optimization Time Integration (HOT) for efficient implicit time stepping of the Material Point method (MPM) irrespective of simulated materials and conditions. HOT is a MPM-specialized hierarchical optimization algorithm that solves nonlinear time step problems for large-scale MPM systems near the CFL-limit, e.g., with step sizes around 0.01s. HOT provides "out-of-the-box" convergent simulations across widely varying materials and computational resolutions without parameter tuning. As the first MPM solver enhanced by h-multigrid, HOT is highly parallelizable and, as we show in our analysis, robustly 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 10x performance speedup across a wide range of challenging benchmark test simulations.


page 12

page 13

page 19


A Gauss-Seidel projection method with the minimal number of updates for stray field in micromagnetic simulations

Magnetization dynamics in magnetic materials is often modeled by the Lan...

Do NHL goalies get hot in the playoffs? A multilevel logistic regression analysis

The hot-hand theory posits that an athlete who has performed well in the...

Rapid Detection of Hot-spots via Tensor Decomposition with applications to Crime Rate Data

We propose an efficient statistical method (denoted as SSR-Tensor) to ro...

Kinetic Simulation of Collisional Magnetized Plasmas with Semi-Implicit Time Integration

Plasmas with varying collisionalities occur in many applications, such a...

Combining p-multigrid and multigrid reduced in time methods to obtain a scalable solver for Isogeometric Analysis

Isogeometric Analysis (IgA) has become a viable alternative to the Finit...

A Transfer Learning Approach for Microstructure Reconstruction and Structure-property Predictions

Stochastic microstructure reconstruction has become an indispensable par...

Hot Roller Embossing for the Creation of Microfluidic Devices

We report on the hot roller embossing of polymer sheets for the creation...

Please sign up or login with your details

Forgot password? Click here to reset