  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.
