Two-way coupled MHD-PIC simulations of magnetic reconnection in magnetic island coalescence

Abstract

We present simulations of magnetic reconnection with a newly developed coupled MHD-PIC code. In this work a global magnetohydrodynamic (MHD) simulation receives kinetic feedback within an embedded region that is modeled by a kinetic particle-in-cell (PIC) code. The PIC code receives initial and boundary conditions from the MHD simulation, while the MHD solution is updated with the PIC state. We briefly describe this coupling mechanism. This method is suitable for simulating magnetic reconnection problems, as we show with the example of reconnection in the coalescence of magnetic islands. We compare the MHD, Hall-MHD, fully PIC and coupled MHD-PIC simulations of the magnetic island coalescence. We find that the kinetic simulations are very different from the MHD and Hall-MHD results, while the coupled MHD-PIC simulations can remedy this discrepancy while saving computing time. The diffusion region is well resolved in the kinetic simulations, which is also captured by the coupled MHD-PIC model. The coupled simulation also reproduces the kinetic Hall magnetic fields correctly. We calculate the reconnection rates and find differences between the MHD and kinetic results. We find that the coupled MHD-PIC code can reasonably reproduce the kinetic reconnection rate when a larger PIC feedback region is used, while still saving significant computing time.