MilsPD: a lightweight multiphysics-integrated peridynamics solver

This paper presents MilsPD, a highly efficient multiphysics-integrated peridynamics (PD) solver designed for fracture mechanics applications. It addresses a key computational challenge in the field: the expensive simulation of large-scale multiphysics problems involving evolving fractures and fluid–solid interactions. Central to its performance is a set of fused multiphysics algorithms that update both coupled thermo-mechanical and thermo-hydrodynamic governing equations in optimized passes over material points. Using massive parallelism, our GPU implementation achieves speedups exceeding 1200 × relative to the leading CPU-based codes. We demonstrate the practical utility and superior performance of MilsPD by simulating the large-scale quenching process and formation of columnar joints, a task that remains computationally infeasible for the existing CPU-based multiphysics PD approach. MilsPD thus provides a robust new framework that significantly enhances the capability for high-fidelity failure analysis on previously intractable problems.