CMHD Orszag-Tang Turbulence Problem: Magnetic Pressure (Rusanov WENO)

Computational Magnetohydrodynamics.
Rusanov Riemann's solver.
Resolution: 400x400.
Euler Time stepping method.
Time integration from t=0 to t=1.0 with timestep = 0.0001.
WENO reconstruction scheme to determine fluxes at cell boundaries.
Periodic boundary conditions.
Divergence of the magnetic field kept zero with constrained transport algorithm. (Poisson Equation is iteratively solved.)

Initial conditions:

rho = 25/(36*pi)
u = -sin(2*pi*y)
v = sin(2*pi*x)
p = 5/(12*pi)
BX = -sin(2*pi*Y)/sqrt(4*pi)
BY = sin(4*pi*X)/sqrt(4*pi)

As the evolution time increases, the vortex flow pattern becomes increasingly complicated due to the nonlinear interactions of waves.
As the simulation progresses, a number of strong waves have formed and passed through each other, creating turbulent flow features at all spatial scales.

Code was written in Julia from scratch.