Advancements in Peridynamics.jl for Dynamic Fracture Simulations | Partmann | JuliaCon 2024

Advancements in Peridynamics.jl for Dynamic Fracture Simulations by Kai Partmann
PreTalx: pretalx.com/ju...
Slides: pretalx.com/me...
Peridynamics is a nonlocal continuum formulation that uses an integro-differential equation to describe the relative displacement and forces between material points. This equation is also fulfilled for material discontinuities, allowing for modeling of crack propagation, fragmentation, and many other effects. Peridynamics has been shown to correctly map complex fracture problems with different materials, making it suitable for various engineering applications.
Our ongoing development aims to address the package's critical limitations and needed functionality enhancements. The newest version of the package acknowledges these challenges, particularly the limited support for high-performance computing. To improve this, we have rewritten the package from scratch, incorporating a dual strategy: multithreading via Threads.jl for an interactive user experience and MPI.jl integration for scaling simulations on HPC clusters. This effort involved rethinking the underlying data structure and implementing a customized halo-exchange approach to communicate between threads and MPI ranks. Additionally, a new API focuses on early error detection, addressing past issues where wrong user input could only be caught at the start of a simulation. Inspired by user feedback, the new API offers a more streamlined experience for single and multi-body simulations.
The bond-based model is the first peridynamics formulation with significant limitations and a high sensitivity to discretization errors. New material models that overcome these limitations exist. In the new version of the package, we introduce some of them, like the correspondence formulation, allowing for advanced simulations with higher accuracy and fewer discretization errors. The talk will feature simple tutorials on example simulations, highlighting the functionality of the new API and the expanded capabilities due to new material models and MPI integration.