Phlegethon: a fully compressible magnetohydrodynamic code for simulations in stellar astrophysics

We present PHLEGETHON, a fully compressible, Eulerian magnetohydrodynamic (MHD) code designed for multidimensional simulations in stellar astrophysics. The code uses a time-explicit, second-order, finite-volume method optimized to model a wide range of dynamical processes in stars, from very low-Mach-number turbulent convection in the cores of massive stars to supersonic flows in subsurface convection zones. PHLEGETHON employs low-dissipation Riemann solvers and a well-balanced method to accurately capture slow flows arising from strongly stratified media. The induction equation is solved using a staggered constrained-transport method to ensure divergence-free evolution of the magnetic field. The MHD equations are coupled to arbitrary nuclear reaction networks solved in a time-implicit approach, together with super-time-stepping for efficient treatment of thermal diffusion. Equations of state appropriate for stellar plasmas are available, accounting for partial ionization, electron degeneracy, and electron-positron pair production. The code is implemented in a compact and user-friendly manner, and it scales to tens of thousands of CPU cores using MPI-based domain decomposition. We perform several verification tests to demonstrate the accuracy and versatility of the code, and present simulations of magnetoconvection in a core-collapse supernova progenitor star. The rich variety of physical effects and numerical methods implemented in PHLEGETHON enables the code to model diverse multidimensional processes that play a crucial role in stellar-interior dynamics, such as reactive convection, convective boundary mixing, internal-wave excitation, and magnetic-field amplification mechanisms. Within a single framework, these phenomena can be investigated across a wide range of stellar evolutionary stages, from main-sequence stars to supernova progenitors. PHLEGETHON is publicly accessible online.