Linking computational models to follow the evolution of heated coronal plasma

Abstract

A ‘proof of principle’ is presented, whereby the Ohmic and viscous heating determined by a three-dimensional (3D) MHD model of a coronal avalanche is used as the coronal heating input for a series of field-aligned, one-dimensional (1D) hydrodynamic models. Three-dimensional coronal MHD models require large computational resources. For current numerical parameters, it is difficult to model both the magnetic field evolution and the energy transport along field lines for coronal temperatures much hotter than 1 MK, because of severe constraints on the time step from parallel thermal conduction. Using the 3D MHD heating derived from a simulation and evaluated on a single field line, the 1D models give coronal temperatures of 1 MK and densities 1014–1015 m−3 for a coronal loop length of 80 Mm. While the temperatures and densities vary smoothly along the field lines, the heating function leads to strong asymmetries in the plasma flows. The magnitudes of the velocities in the 1D model are comparable with those seen in 3D reconnection jets in our earlier work. Advantages and drawbacks of this approach for coronal modelling are discussed.