Best approach for implementing particle breakage in DEM

[tianzhen12345]

I am interested in extending PhasicFlow to account for particle breakage (fragmentation) during high-energy impacts or compression. Currently, particles are treated as rigid bodies with constant mass.
Does the current solver support dynamic addition/deletion of particles during runtime?
Are there any existing hooks in the contact laws where breakage criteria (e.g., Tavares model) could be integrated?

[PhasicFlow]

Regarding the first question: yes. it is possible to do it through particles object. Everything is done automatically. See, how particles insertion mechanism works, you can so the same.
Regrading the second question: No, you need to develop a new contact law and mark those particles which are about to break due to high energy impacts. You need to implement a new contact force model, a new interaction class model for this. For a better help, see how we implemented coarse graining in the phasicFlow framework.