Issue |
A&A
Volume 676, August 2023
|
|
---|---|---|
Article Number | A83 | |
Number of page(s) | 13 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/202346788 | |
Published online | 11 August 2023 |
The role of plasma β in global coronal models: Bringing balance back to the force
1
Centre for Mathematical Plasma Astrophysics, KU Leuven,
Celestijnenlaan 200B,
3001
Leuven, Belgium
e-mail: michaela.brchnelova@kuleuven.be
2
Shenzhen Key Laboratory of Numerical Prediction for Space Storm, Institute of Space Science and Applied Technology, Harbin Institute of Technology,
Shenzhen
518055, PR China
3
Key Laboratory of Solar Activity and Space Weather, National Space Science Center, Chinese Academy of Sciences,
Beijing
100190, PR China
4
Institute of Physics, University of Maria Curie-Skłodowska,
ul. Radziszewskiego 10,
20-031
Lublin, Poland
Received:
2
May
2023
Accepted:
8
June
2023
Context. COolfluid COrona uNstrUcTured (COCONUT) is a global coronal magnetohydrodynamic (MHD) model that was recently developed and will soon be integrated into the ESA Virtual Space Weather Modelling Centre (VSWMC). In order to achieve robustness and fast convergence to a steady state for numerical simulations with COCONUT, several assumptions and simplifications were made during its development, such as prescribing filtered photospheric magnetic maps to represent the magnetic field conditions in the lower corona. This filtering leads to smoothing and lower magnetic field values at the inner boundary (i.e. the solar surface), resulting in an unrealistically high plasma β (greater than 1 in a large portion of the domain).
Aims. We aim to examine the effects of prescribing such filtered magnetograms in global coronal simulations and formulate a method for achieving more realistic plasma β values and improving the resolution of electromagnetic features without losing computational performance.
Methods. We made use of the newly developed COCONUT solver to demonstrate the effects of the highly pre-processed magnetic maps set at the inner boundary and the resulting high plasma β on the features in the computational domain. Then, in our new approach, we shifted the inner boundary to 2 R⊙ from the original 1.01 R⊙ and preserved the prescribed highly filtered magnetic map. With the shifted boundary, the boundary density and pressure were also naturally adjusted to better represent the considered physical location. This effectively reduces the prescribed plasma β and leads to a more realistic setup. The method was applied on a magnetic dipole, a minimum (2008) and a maximum (2012) solar activity case, to demonstrate its effects.
Results. 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 β at the boundary and inside the computational domain. This is also demonstrated via synthetic white light imaging (WLI) and with the validation against tomography data. 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, as its capacity to resolve features or structures is superior, the operational performance can be additionally increased if needed.
Conclusions. The newly developed method is thus deemed as a good potential replacement of the original setup for operational purposes, providing higher physical detail of the resolved profiles while preserving a good convergence and robustness of the solver.
Key words: magnetohydrodynamics (MHD) / methods: numerical / Sun: corona
© The Authors 2023
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.