Numerical simulation of reconnection in an emerging magnetic flux region

Numerical solutions in two dimensions are presented for the resistive MHD equations of an initial boundary value problem, simulating reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. This numerical model displays four main phases, which are interpreted in terms of (1) a slowly evolving quasi-steady phase during which most of the magnetic flux emerges, with approximate equilibrium between magnetic and pressure forces; (2) an impulsive phase in which either the reconnection or continuing emergence of the first phase comes to disrupt the equilibrium, and extensive acceleration occurs as the high pressure region within the emerging region drives the fluid upwards and outwards; (3) a second quasi-steady phase; and (4) a potential-static phase in which continuing reconnection, ohmic dissipation, and fluid transport through boundaries depletes the system of all currents and flows.