Search

Context. Solar wind modelling has become a crucial area of study due to the increased dependence of modern society on technology, navigation, and power systems. Accurate space weather forecasts can predict upcoming threats to Earth's geospace and allow for harmful socioeconomic impacts to be mitigated. Coronal and heliospheric models must be as realistic as possible to achieve successful predictions. In this study, we examine a novel full magnetohydrodynamic (MHD) chain from the Sun to Earth. Aims. The goal of this study is to demonstrate the capabilities of the full MHD modelling chain from the Sun to Earth by finalising the implementation of the full MHD coronal model into the COolfluid COroNa UnsTructured (COCONUT) model and coupling it to the MHD heliospheric model Icarus. The resulting coronal model has significant advantages compared to the pre-existing polytropic alternative, as it includes more physics and allows for a more realistic modelling of bi-modal wind, which is crucial for heliospheric studies. In particular, we examine different empirical formulations for the heating terms in the MHD equations to determine an optimal one that would be able to mimic a realistic solar wind configuration most accurately. Methods. New heating source terms were implemented into the MHD equations of the pre-existing polytropic COCONUT model. A realistic specific heat ratio was applied. In this study, only thermal conduction, radiative losses, and approximated coronal heating function were considered in the energy equation. Multiple approximated heating profiles were examined to see the effect on the solar wind. The output of the coronal model was used to onset the 3D MHD heliospheric model Icarus. A minimum solar activity case was chosen as the first test case for the full MHD model. The numerically simulated data in the corona and the heliosphere were compared to observational products. First, we compared the density data to the available tomography data near the Sun and then the modelled solar wind time series in Icarus was compared to OMNI 1-min data at 1 AU. Results. A range of approximated heating profiles were used in the full MHD coronal model to obtain a realistic solar wind configuration. The bi-modal solar wind was obtained for the corona when introducing heating that is dependent upon the magnetic field. The modelled density profiles are in agreement with the tomography data. The modelled wind in the heliosphere is in reasonable agreement with observations. Overall, the density is overestimated, whereas the speed at 1 AU is more similar to OMNI 1-min data. The general profile of the magnetic field components is modelled well, but its magnitude is underestimated. Conclusions. We present a first attempt to obtain the full MHD chain from the Sun to Earth with COCONUT and Icarus. The coronal model has been upgraded to a full MHD model for a realistic bi-modal solar wind configuration. The approximated heating functions have modelled the wind reasonably well, but simple approximations are not enough to obtain a realistic density-speed balance or realistic features in the low corona and farther, near the outer boundary. The full MHD model was computed in 1.06 h on 180 cores of the Genius cluster of the Vlaams Supercomputing Center, which is only 1.8 times longer than the polytropic simulation. The extended model gives the opportunity to experiment with different heating formulations and improves the approximated function to model the real solar wind more accurately. [ABSTRACT FROM AUTHOR]

Context. Coronal mass ejections (CMEs) stand as intense eruptions of magnetized plasma from the Sun, and they play a pivotal role in driving significant changes of the heliospheric environment. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space weather forecasting. Aims. The primary objective of this paper is to establish a connection between CMEs and their progenitors in solar source regions, enabling us to infer the magnetic structures of CMEs before their full development. Methods. We created a dataset comprising a magnetic flux rope series with varying projection shapes (S-, Z-, and toroid-shaped), sizes, and toroidal fluxes using the Regularized Biot-Savart Laws (RBSL). These flux ropes were inserted into solar quiet regions with the aim of imitating the eruptions of quiescent filaments. Thereafter, we simulated the propagation of these flux ropes from the solar surface to a distance of 25 R⊙ with our global coronal magnetohydrodynamic (MHD) model COCONUT. Results. Our parametric survey revealed significant impacts of source flux ropes on the consequent CMEs. Regarding the flux-rope morphology, we find that the projection shape (e.g., sigmoid or torus) can influence the magnetic structures of CMEs at 20 R⊙, albeit with minimal impacts on the propagation speed. However, these impacts diminish as source flux ropes become fat. In terms of toroidal flux, our simulation results demonstrate a pronounced correlation with the propagation speed of CMEs as well as the successfulness in erupting. Conclusions. This work builds the bridge between the CMEs in the outer corona and their progenitors in solar source regions. Our parametric survey suggests that the projection shape, cross-section radius, and toroidal flux of source flux ropes are crucial parameters in predicting magnetic structures and the propagation speed of CMEs, providing valuable insights for space weather prediction. On the one hand, the conclusion drawn here could be instructive in identifying the high-risk eruptions with the potential to induce stronger geomagnetic effects (Bz and propagation speed). On the other hand, our findings hold practical significance for refining the parameter settings of launched CMEs at 21.5 R⊙ in heliospheric simulations, such as with EUHFORIA, based on observations for their progenitors in solar source regions. [ABSTRACT FROM AUTHOR]

It has been suggested that under solar coronal conditions, drift waves may contribute to coronal heating. Specific properties of the drift waves to be expected in the solar corona have, however, not yet been determined using more advanced numerical models. We investigate the linear properties of density-gradient-driven drift waves in the solar coronal plasma using gyrokinetic ion–electron simulations with the gyrokinetic code Gene, solving the Vlasov–Maxwell equations in five dimensions assuming a simple slab geometry. We determine the frequencies and growth rates of the coronal density gradient-driven drift waves with changing plasma parameters, such as the electron β , the density gradient, the magnetic shear, and additional temperature gradients. To investigate the influence of the finite Larmor radius effect on the growth and structure of the modes, we also compare the gyrokinetic simulation results to those obtained from drift-kinetics. In most of the investigated conditions, the drift wave has positive growth rates that increase with increasing density gradient and decreasing β. In the case of increasing magnetic shear, we find that from a certain point, the growth rate reaches a plateau. Depending on the considered reference environment, the frequencies and growth rates of these waves lie on the order of 0.1 mHz–1 Hz. These values correspond to the observed solar wind density fluctuations near the Sun detected by WISPR, currently of unexplained origin. As a next step, nonlinear simulations are required to determine the expected fluctuation amplitudes and the plasma heating resulting from this mechanism. [ABSTRACT FROM AUTHOR]

We have developed COCONUT -- a novel MHD solar corona model as a basis for EUHFORIA 2.0 space weather forecasting. The steady-state solution in a numerical domain extending from the solar surface up to 0.1AU is computed. Here the pre-processed magnetic maps prescribed at the inner boundary determine the magnetic topology. These magnetic maps consist of radial components of Helioseismic and Magnetic Imager magnetograms projected on spherical harmonics with a selected maximum frequency. Within the numerical domain a set of ideal MHD equations with gravity is solved implicitly on an unstructured mesh. This choice of grid allows us to omit the issue of singularities presented at the poles in typically used structured grids. Through heavy parallelization, an optimized CFL profile and other numerical acceleration techniques, we achieved superior convergence times with high robustness and reliability even for numerically challenging cases of maximum solar activity. The MHD solution is generally obtained within a few hours of computation, which is crucial for space weather forecasting systems. Our model passed a series of tests and was further validated using detailed comparison with observations. It successfully recreated magnetic structures (streamers, pseudo streamers, coronal holes) observed during the solar eclipses in March 2015 and in July 2019, thus corresponding to maximum and minimum of solar activity. The COCONUT code will be next extended to include the heating and radiation terms for a more realistic physical description. In the future, a multispecies formulation will be implemented allowing us to consider also the impact of electrons as well as of neutrals on the solution. [ABSTRACT FROM AUTHOR]