Create Your Own Constrained Transport Magnetohydrodynamics Simulation (With Python)

In today’s recreational coding exercise, we will simulate the Orszag-Tang vortex — a demonstration of supersonic 2D magnetohydrodynamics (MHD) turbulence — using a constrained transport finite volume algorithm. This tutorial builds on our Finite Volume tutorial for the compressible Euler equations. Here, we will add a magnetic field coupled to the fluid flow.

You may find the accompanying Python code on github.