Numerical relaxation of a 3D MHD Taylor - Woltjer state subject to abrupt expansion

Since the advent of Taylor-Woltjer theory [J B Taylor, PRL, 33, 1139 (1974), L Woltjer, PNAS, 44, 489 (1958)], it has been widely believed that situations with perfectly conducting boundaries and near ideal conditions, the final state of MHD system would be force-free Taylor-Woltjer states defined as $\vec{\nabla} \times \vec{B} = \alpha \vec{B}$ with $\alpha$ as a constant and $\vec{B}$ is the magnetic field defined over a volume $V$. These states are of fundamental importance in fusion plasmas [J B Taylor, RMP 58, 741 (1986)]. More recently, several new MHD models have been proposed - for example Reduced Multi-region relaxed MHD [S R Hudson {\it et al}, Phys. Plasmas, 19, 112502 (2012)] and arbitrary scale relaxation model to Taylor-Woltjer state [H Qin {\it et al}, PRL, 109, 235001 (2012)] to mention a few. In the present work, we use a 3D compressible MHD solver in cartesian geometry which can handle conducting or periodic as well has mixed boundary conditions to investigate numerically the arbitrary scale relaxation model proposed by Qin et al [H Qin {\it et al}, PRL, 109, 235001 (2012)]. For this purpose, we consider two volumes $V_{init}$ and $V_{final}$. We load the 3D MHD solver in the limit of zero compressibility with a Taylor-Woltjer state $B_{init}(x,y,z,t=0)$ and let it again a numerical evolve with conducting boundaries at $V_{init}$ to make sure that we have obtained a numerically steady Taylor - Woltjer state for volume $V_{init}$. Followed by this procedure, we "suddenly" relax the boundaries to a new volume $V_{final}$, such that $V_{init}$ $<$ $V_{final}$ and evaluate whether or not the system attains quasi-steady state. Details of the numerical method used, the protocol followed, the expansion technique and the novelty of this numerical experiment and details of our results have been presented in this paper.