The role of plasma beta in global coronal models: Bringing balance back to the force

COCONUT is a global coronal magnetohydrodynamic model recently developed. In order to achieve robustness and fast convergence to steady-state, several assumptions have been made during its development, such as prescribing filtered photospheric magnetic maps for representing the magnetic field in the lower corona. This filtering leads to smoothing and lower magnetic field values at the inner boundary, resulting in an unrealistically high plasma beta.In this paper, we examine the effects of prescribing such filtered magnetograms and formulate a method for achieving more realistic plasma beta values without losing computational performance. We demonstrate the effects of the highly pre-processed magnetic maps and the resulting high plasma beta on the features in the domain. Then, in our new approach, we shift the inner boundary to 2 Rs and preserve the prescribed highly filtered magnetic map. This effectively reduces the prescribed plasma beta and leads to a more realistic setup. The method is applied on a magnetic dipole, a minimum (2008) and a maximum (2012) solar activity case, to demonstrate its effects. The results obtained with the proposed approach show significant improvements in the resolved density and radial velocity profiles, and far more realistic values of the plasma \{beta} at the boundary and inside the computational domain. This is also demonstrated via synthetic white light imaging and with the validation against tomography. The computational performance comparison shows similar convergence to a limit residual on the same grid when compared to the original setup. Considering that the grid can be further coarsened with this new setup, the operational performance can be additionally increased if needed. The newly developed method is thus deemed as a good potential replacement of the original setup for operational purposes.
Comment: accepted for publication, 13 pages, 17 figures, 1 table