Modelling chromospheric evaporation in response to coronal heating

This thesis presents a new computationally efficient method for modelling the response of the solar corona to the release of energy. During impulsive heating events, the coronal temperature increases which leads to a downward heat flux into the transition region (TR). The plasma is unable to radiate this excess conductive heating and so the gas pressure increases locally. The resulting pressure gradient drives an upflow of dense material, creating an increase in the coronal density. This density increase is often called chromospheric evaporation. A process which is highly sensitive to the TR resolution in numerical simulations. If the resolution is not adequate, then the downward heat flux jumps over the TR and deposits the heat in the chromosphere, where it is radiated away. The outcome is that with an under-resolved TR, major errors occur in simulating the coronal density evolution. We address this problem by treating the lower transition region as a discontinuity that responds to changing coronal conditions through the imposition of a jump condition that is derived from an integrated form of energy conservation. In this thesis, it is shown that this method permits fast and accurate numerical solutions in both one-dimensional and multi-dimensional simulations. By modelling the TR with this appropriate jump condition, we remove the influence of poor numerical resolution and obtain the correct evaporative response to coronal heating, even when using resolutions that are compatible with multi-dimensional magnetohydrodynamic simulations