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.
Before we begin, below is a gif of what running our simulation looks like: