Effect of numerical resolution on synthetic observables of simulated coronal loops

Abstract

Increasingly realistic simulations of the corona are used to predict synthetic observables for instruments onboard both existing and upcoming heliophysics space missions. Synthetic observables play an important role in constraining coronal heating theories. Choosing the spatial resolution of numerical simulations involves a trade-off between accuracy and computational cost. Since the numerical resolution affects not only the scale of structures that can be resolved, but also thermodynamic quantities such as the average coronal density, it is important to quantify the effect on synthesized observables. Using 3D radiative magnetohydrodynamic simulations of coronal loops at three different grid spacings, from 60 km down to 12 km, we find that changes in numerical resolution lead to differences in thermodynamic quantities and stratification as well as dynamic behaviour. Higher grid resolution results in a more complex and dynamic atmosphere. The resolution affects the emission intensity as well as the velocity distribution, thereby affecting synthetic spectra derived from the simulation. The distribution of synthetic coronal loop strand sizes changes as more fine-scale structure is resolved. A number of parameters, however, seem to start to saturate from our chosen medium grid resolution on. Our study shows that while choosing a sufficiently high resolution matters when comparing forward-modelled observables with data from current and future space missions, for most purposes not much is gained by further increasing the resolution beyond a grid spacing of 24 km, which seems to be adequate for reproducing bulk loop properties and forward-modelled emission, representing a good trade-off between accuracy and computational resource.