Stability Analysis of Explicit MPM

Bai, Song; Schroeder, Craig

Stability Analysis of Explicit MPM

Files

v41i8pp019-030.pdf (37.96 MB)

movie.mp4 (119.32 MB)

tech-doc.pdf (580.7 KB)

Date

2022

Authors

Bai, Song
Schroeder, Craig

Publisher

The Eurographics Association and John Wiley & Sons Ltd.

Abstract

In this paper we analyze the stability of the explicit material point method (MPM). We focus on PIC, APIC, and CPIC transfers using quadratic and cubic splines in two and three dimensions. We perform a fully three-dimensional Von Neumann stability analysis to study the behavior within the bulk of a material. This reveals the relationship between the sound speed, CFL number, and actual time step restriction and its dependence on discretization options. We note that boundaries are generally less stable than the interior, with stable time steps generally decreasing until the limit when particles become isolated. We then analyze the stability of a single particle to derive a novel time step restriction that stabilizes simulations at their boundaries. Finally, we show that for explicit MPM with APIC or CPIC transfers, there are pathological cases where growth is observed at arbitrarily small time steps sizes. While these cases do not necessarily pose a problem for practical usage, they do suggest that a guarantee of stability may be theoretically impossible and that necessary but not sufficient time step restrictions may be a necessary and practical compromise.

CCS Concepts: Computing methodologies --> Physical simulation

        @article{10.1111:cgf.14620
,
journal = {Computer Graphics Forum},
title = {{Stability Analysis of Explicit MPM
}},
author = {Bai, Song and 
Schroeder, Craig
},
year = {2022
},
publisher = {The Eurographics Association and John Wiley & Sons Ltd.
},
ISSN = {1467-8659
},
DOI = {10.1111/cgf.14620
}
}