A simple solver
2D ideal MHD on a periodic grid: Rusanov fluxes, forward Euler, a projection for div B = 0. Only NumPy.
A simple solver
2D ideal MHD on a periodic grid: Rusanov fluxes, forward Euler, a projection for div B = 0. Only NumPy.
Test problems
The Orszag-Tang vortex, the MHD rotor, a magnetized blast wave and a multi-fluid run, as executed notebooks.
Room to grow
Higher-order reconstruction, HLL and HLLD fluxes, constrained transport: each a lesson with a check you can run.