An Implicit Particle-in-Cell Method for Granular Materials

An implicit-in-time method for granular materials is described. The method combines the material point method, a first-order contact algorithm, and a Newton–Krylov equation solver to give improved energy conservation, stabilization of the finite-grid-instability, and the correct description of collisions between grains.