The following is provided to button up this thread.
The setup ~1.8M particles and all ~50 species ran to completion, simulating 1 second of operation, using a single NVIDIA H100 GPU. The animation shows the volume fraction field and particles, colored by log-scale velocity magnitude.