  Magnetic reconnection is an important process in astrophysical environments,
as it re-configures magnetic field topology and converts magnetic energy into
thermal and kinetic energy. In extreme astrophysical systems, such as black
hole coronae and pulsar magnetospheres, radiative cooling modifies the energy
partition by radiating away internal energy, which can lead to the radiative
collapse of the layer. In this paper, we perform 2D & 3D simulations to model
the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access
cooling rates in the laboratory necessary to investigate reconnection in a
previously unexplored radiatively-cooled regime. These simulations are
performed in GORGON, an Eulerian resistive magnetohydrodynamic code, which
models the experimental geometry comprising two exploding wire arrays driven by
20 MA of current on the Z machine (Sandia National Laboratories). Radiative
losses are implemented using non-local thermodynamic equilibrium tables
computed using the atomic code Spk, and we probe the effects of radiation
transport by implementing both a local radiation loss model and P$_{1/3}$
multi-group radiation transport. The load produces highly collisional,
super-Alfv\'enic $(M_{A} \approx 1.5)$, supersonic $(M_S \approx 4-5)$ plasma
flows which generate a reconnection layer ($L/{\delta} \approx 100, S_L \approx
400$). The reconnection layer undergoes radiative collapse when the radiative
losses exceed Ohmic and compressional heating ${\tau}_{cool}^{-1}/{\tau}_A^{-1}
\approx 100$; this generates a cold strongly compressed current sheet, leading
to an accelerated reconnection rate, consistent with theoretical predictions.
Finally, the current sheet is unstable to the plasmoid instability, but the
magnetic islands are extinguished by strong radiative cooling before ejection
from the layer.
