The Material Point Method (MPM) has become a cornerstone of physics-based simulation, widely used in geomechanics and computer graphics for modeling phenomena such as granular flows, viscoelasticity, fracture mechanics, etc. Despite its versatility, the original MPM suffers from cell-crossing instabilities caused by discontinuities in particle-grid transfer kernels. Existing solutions mitigate these issues by adopting smoother shape functions, but at the cost of increased computational overhead due to larger kernel support. In this paper, we propose a novel $C^2$-continuous compact kernel for MPM that achieves a unique balance between stability and computational efficiency. Our method integrates seamlessly with Affine Particle-In-Cell (APIC) and Moving Least Squares (MLS) MPM, while only doubling the number of grid nodes associated with each particle compared to linear kernels. At its core is an innovative dual-grid framework, which associates particles with grid nodes exclusively within the cells they occupy on two staggered grids, ensuring consistent and stable force computations. To further accelerate performance, we present a GPU-optimized implementation inspired by state-of-the-art massively parallel MPM techniques, achieving an additional $2\times$ speedup in G2P2G transfers over quadratic B-spline MPM. Comprehensive validation through unit tests, comparative studies, and stress tests demonstrates the efficacy of our approach in conserving both linear and angular momentum, handling stiff materials, and scaling efficiently for large-scale simulations. Our results highlight the transformative potential of compact, high-order kernels in advancing MPM's capabilities for stable, high-performance simulations, paving the way for more computationally efficient applications in computer graphics and beyond.
翻译:暂无翻译