PION: Simulating bow shocks and circumstellar nebulae around massive stars

Expanding nebulae are produced by mass loss from stars, especially during late stages of evolution. We describe the algorithms and methods implemented in the radiation-magnetohydrodynamics (MHD) code PION for highly scalable simulations using static mesh-refinement. We present results from 3D MHD simulations of bow shocks around runaway massive stars, and of the expansion of a fast wind from a Wolf-Rayet star into the slow wind from a previous red supergiant phase of evolution. PION is free software that can be downloaded from https://www.pion.ie/