COCONUT-MF: Two-fluid ion-neutral global coronal modelling

The global coronal model COCONUT was originally developed to replace models such as the WSA model in space weather forecasting to improve the physical accuracy of the predictions. This model has, however, several simplifications implemented in its formulation to allow for rapid convergence, one of which includes a single-fluid treatment. In this paper, we have two goals. Firstly, we aim to introduce a novel multi-fluid global coronal model and validate it with simple cases as well as with real data-driven applications. Secondly, we aim to investigate to what extent considering a single-fluid plasma in the global coronal model might affect the resulting plasma dynamics, and thus whether the assumptions on which the single-fluid coronal model is based are justified. We developed a multi-fluid global coronal model, COCONUT-MF, which resolves the ion and neutral fluid equations separately. While this model is still steady-state and thus does not resolve unsteady processes, it can account for charge exchange, chemical and collisional contributions. We present the results of the ion-neutral modelling for a dipole, a minimum of solar activity, and a solar maximum. We demonstrate the higher accuracy of the applied AUSM+ scheme compared to HLL. Subsequently, we also evaluate the effects of the considered ion-neutral coupling terms on the resulting plasma dynamics. Despite the very low concentration of neutrals, these terms still affect the flow field to a limited but non-negligible extent (up to 5 to 10% locally). Even though the coronal plasma is generally assumed to be collisionless, our results show that there is sufficient collisionality in it to couple the two fluids. Follow-up work will include extension of the model to lower atmospheric layers of the Sun and inclusion of more advanced physical terms such as heating and radiation.