The simulation of the dust and gas environment around cometary nuclei is a challenging computational problem. We develop scalable algorithms to advance physical models of comets to take into account the detailed shape, rotation, and physical properties of comets. The tool suit runs on the HLRN supercomputer and is used to predict the activity of comet 67P/Churyumov-Gerasimenko, which is observed by the ESA Rosetta mission.