Abstract
This paper presents a set of computationally efficient, adaptive grids for magnetic reconnection phenomenon where the current density can develop large gradients in the reconnection region. Four-field extended MagnetoHydroDynamics (MHD) equations with hyperviscosity terms are transformed so that the curvilinear coordinates replace the Cartesian coordinates as the independent variables, and moving grids' velocities are also considered in this transformed system as a part of interpolating the physical solutions from the old grid to the new grid as time advances. The curvilinear coordinates derived from the current density through the Monge–Kantorovich (MK) optimization approach help to reduce the resolution requirements during the computation.