This repository contains python code that are used to simulate 3 magnetohydrodynamic (MHD) experiments. The python code utilizes an open-source environment called dedalus. All python code is run on the PSC-Bridges supercomputing cluster.
The directory "Hartmann" contains python code used to simulate the MHD Hartmann flow problem. It is a problem with well known analytical soluations. A typical simulation costs around 5 CPU-hours.
The directory ssxluo contains python code that were used to simulate
(1) the relaxation of a spheromak into a B-force-free equilibrium named the Taylor state (see original paper)
(2) the merging of two Taylor state plasmas due to pressure gradient (magnetic reconnection takes place in this experiment).
A taypical simulation run costs around 5,000-10,000 CPU-hours with a resolution of 28x24x360. The x and y axis of the solution domain can be distributively solved.