Abstract
To study acoustic wave propagation and the corresponding energy deposition in partially ionized plasmas, we use a two-fluid computational model that treats neutrals and charged particles ...(electrons and ions) as two separate fluids. This two-fluid model takes into account the ion–neutral collisions, ionization, and recombination, allowing us to investigate both the collisional and reactive interactions between uncoupled ions and neutrals in the plasmas. In the present numerical simulations, the initial density is specified to reach hydrostatic equilibrium, and as a comparison, chemical equilibrium is also taken into account to provide a density profile that differs from typical hydrostatic equilibrium profiles. External velocity drivers are then imposed to generate monochromatic acoustic waves. As is well known, the upward propagating acoustic waves steepen in gravitationally stratified plasmas due to the exponentially decreasing density, and they heat the plasmas in the nonlinear regimes where kinetic energy is dissipated by shock waves and collisional interactions. In particular, the lower ionization fraction resulting from the present initial chemical equilibrium significantly enhances the heating efficiency. Moreover, the ionization process absorbs a significant amount of energy, and the decoupling between ions and neutrals is also enhanced while considering ionization and recombination. Therefore, simulations without considering ionization and recombination may overestimate the overall heating effects but also underestimate the energy dissipation. The results also suggest that a more accurate ionization and recombination model could be essential for improving the modeling of partially ionized plasmas.
Abstract
Many aspects of the 3D structure and evolution of interplanetary coronal mass ejections (ICMEs) remain unexplained. Here, we investigate two main topics: (1) the coherence scale of magnetic ...fields inside ICMEs, and (2) the dynamic nature of ICME magnetic complexity. We simulate ICMEs interacting with different solar winds using the linear force-free spheromak model incorporated into the EUHFORIA model. We place a swarm of ∼20,000 spacecraft in the 3D simulation domain and characterize ICME magnetic complexity and coherence at each spacecraft based on the simulated time series. Our simulations suggest that ICMEs retain a lower complexity and higher coherence along their magnetic axis, but that a characterization of their global complexity requires crossings along both the axial and perpendicular directions. For an ICME of initial half angular width of 45° that does not interact with other large-scale solar wind structures, global complexity can be characterized by as little as 7–12 spacecraft separated by 25°, but the minimum number of spacecraft rises to 50–65 (separated by 10°) if interactions occur. Without interactions, ICME coherence extends for 45°, 20°–30°, 15°–30°, and 0°–10° for
B
,
B
ϕ
,
B
θ
, and
B
r
, respectively. Coherence is also lower in the ICME west flank compared to the east flank due to Parker spiral effects. Moreover, coherence is reduced by a factor of 3–6 by interactions with solar wind structures. Our findings help constrain some of the critical scales that control the evolution of ICMEs and aid in the planning of future dedicated multispacecraft missions.
Abstract We investigate the role of perpendicular diffusion in shaping the energetic ion spectrum in corotating interaction regions (CIRs), focusing on its mass-to-charge ( A / Q ) ratio dependence. ...We simulate a synthetic CIR using the EUropean Heliospheric FORecasting Information Asset and model the subsequent ion acceleration and transport by solving the focused transport equation incorporating both parallel and perpendicular diffusion. Our results reveal distinct differences in ion spectra between scenarios with and without perpendicular diffusion. In the absence of perpendicular diffusion, ion spectra near CIRs show a strong ( A / Q ) ϵ dependence with ϵ depending on the turbulence spectral index, agreeing with theoretical predictions. In contrast, the incorporation of perpendicular diffusion, characterized by a weak A / Q dependence, leads to similar spectra for different ion species. This qualitatively agrees with observations of energetic particles in CIRs.
Abstract
A key challenge in space weather forecasting is accurately predicting the magnetic field topology of interplanetary coronal mass ejections (ICMEs), specifically the north–south magnetic ...field component (
B
z
) for Earth-directed CMEs. Heliospheric MHD models typically use spheromaks to represent the magnetic structure of CMEs. However, when inserted into the ambient interplanetary magnetic field, spheromaks can experience a phenomenon reminiscent of the condition known as the “spheromak tilting instability,” causing its magnetic axis to rotate. From the perspective of space weather forecasting, it is crucial to understand the effect of this rotation on predicting
B
z
at 1 au while implementing the spheromak model for realistic event studies. In this work, we study this by modeling a CME event on 2013 April 11 using the European Heliospheric Forecasting Information Asset. Our results show that a significant spheromak rotation up to 90° has occurred by the time it reaches 1 au, while the majority of this rotation occurs below 0.3 au. This total rotation resulted in poor predicted magnetic field topology of the ICME at 1 au. To address this issue, we further investigated the influence of spheromak density on mitigating rotation. The results show that the spheromak rotation is less for higher densities. Importantly, we observe a substantial reduction in the uncertainties associated with predicting
B
z
when there is minimal spheromak rotation. Therefore, we conclude that spheromak rotation adversely affects
B
z
prediction in the analyzed event, emphasizing the need for caution when employing spheromaks in global MHD models for space weather forecasting.
In order to address the growing need for more accurate space-weather predictions, a new model named EUHFORIA (EUropean Heliospheric FORecasting Information Asset) was recently developed. We present ...the first results of the performance assessment for the solar-wind modeling with EUHFORIA and identify possible limitations of its present setup. Using the basic EUHFORIA 1.0.4 model setup with the default input parameters, we modeled background solar wind (no coronal mass ejections) and compared the obtained results with
Advanced Composition Explorer
(ACE)
in-situ
measurements. For the purposes of statistical study we developed a technique of combining daily EUHFORIA runs into continuous time series. The combined time series were derived for the years 2008 (low solar activity) and 2012 (high solar activity), from which
in-situ
speed and density profiles were extracted. We find for the low-activity phase a better match between model results and observations compared to the high-activity time interval considered. The quality of the modeled solar-wind parameters is found to be rather variable. Therefore, to better understand the results obtained we also qualitatively inspected characteristics of coronal holes,
i.e.
the sources of the studied fast streams. We discuss how different characteristics of the coronal holes and input parameters to EUHFORIA influence the modeled fast solar wind, and suggest possibilities for the improvement of the model.
Abstract
The propagation of interplanetary coronal mass ejections (ICMEs) in the heliosphere is influenced by many physical phenomena, related to the internal structure of the ICME and its ...interaction with the ambient solar wind and magnetic field. As the solar magnetic field is modulated by the 11 yr dynamo cycle, our goal is to perform a theoretical exploratory study to assess the difference of propagation of an ICME in typical minimum and maximum activity backgrounds. We define a median representative CME at 0.1 au, using both observations and numerical simulations, and describe it using a spheromak model. We use the heliospheric propagator EUropean Heliospheric FORecasting Information Asset to inject the same ICME in two different background wind environments. We then study how the environment and the internal CME structure impact the propagation of the ICME toward Earth, by comparison with an unmagnetized CME. At minimum of activity, the structure of the heliosphere around the ecliptic causes the ICME to slow down, creating a delay with the polar parts of the ejecta. This delay is more important if the ICME is faster. At maximum of activity, a southern coronal hole causes a northward deflection. For these cases, we always find that the ICME at the maximum of activity arrives first, while the ICME at the minimum of activity is actually more geoeffective. The sign of the helicity of the ICME is also a crucial parameter, but at the minimum of activity only, since it affects the magnetic profile and the arrival time up to 8 hr.
Abstract
This paper is dedicated to the new implicit unstructured coronal code COCONUT, which aims at providing fast and accurate inputs for space-weather forecasting as an alternative to empirical ...models. We use all 20 available magnetic maps of the solar photosphere covering the date of 2019 July 2, which corresponds to a solar eclipse on Earth. We use the same standard preprocessing on all maps, then perform coronal MHD simulations with the same numerical and physical parameters. We conclude by quantifying the performance of each map using three indicators from remote-sensing observations: white-light total solar eclipse images for the streamers’ edges, EUV synoptic maps for coronal holes, and white-light coronagraph images for the heliospheric current sheet. We discuss the performance of space-weather forecasting and show that the choice of the input magnetic map has a strong impact. We find performances between 24% and 85% for the streamers’ edges, 24%–88% for the coronal hole boundaries, and a mean deviation between 4° and 12° for the heliospheric current sheet position. We find that the HMI runs perform better on all indicators, with GONG-ADAPT being the second-best choice. HMI runs perform better for the streamers’ edges, and GONG-ADAPT for polar coronal holes, HMI synchronic for equatorial coronal holes, and the streamer belt. We especially illustrate the importance of the filling of the poles. This demonstrates that the solar poles have to be taken into account even for ecliptic plane previsions.
Abstract Gradual solar energetic particle (SEP) events are associated with shocks driven by coronal mass ejections (CMEs). The merging of two CMEs (so-called cannibalistic CMEs) and the interaction ...of their associated shocks has been linked to some of the most powerful solar storms ever recorded. Multiple studies have focused on the observational aspects of these SEP events, yet only a handful have focused on modeling similar CME–CME interactions in the heliosphere using advanced magnetohydrodynamic (MHD) models. This work presents, to our knowledge, the first modeling results of a fully time-dependent 3D simulation that captures both the interaction of two CMEs and its effect on the acceleration and transport of SEPs. This is achieved by using an MHD model for the solar wind and CME propagation together with an integrated SEP model. We perform different simulations and compare the behavior of the energetic protons in three different solar wind environments, where a combination of two SEP-accelerating CMEs are modeled. We find that particle acceleration is significantly affected by the presence of both CMEs in the simulation. Initially, less efficient acceleration results in lower-energy particles. However, as the CMEs converge and their shocks eventually merge, particle acceleration is significantly enhanced through multiple acceleration processes between CME-driven shocks, resulting in higher particle intensities and energy levels.
Abstract
We developed a novel global coronal COCONUT (Coolfluid Corona Unstructured) model based on the COOLFluiD code. The steady-state model is predetermined by magnetograms set as boundary ...conditions, while inside the numerical domain the corona is described by MHD equations. This set of equations is solved with the use of an implicit solver on unstructured grids. Here we present numerically obtained results for two extremes of the solar activity cycle represented by CR 2161 and CR 2219 for solar maximum and minimum, respectively. We discuss the impact of reconstruction level on representative solar corona solutions and thus also the impact of small magnetic structures on the overall structure of the solar wind. Moreover, both cases correspond to particular solar eclipses, namely those in 2015 March and 2019 July, to allow us the direct comparison of simulations with observed coronal features. We use a validation scheme proposed by Wagner et al. (from less to more sophisticated methods, i.e., visual classification, feature matching, streamer direction and width, brute force matching, topology classification). The detailed comparison with observations reveals that our model recreates relevant features such as the position, direction, and shape of the streamers (by comparison with white-light images) and the coronal holes (by comparison with extreme ultraviolet images) for both cases of minimum and maximum solar activity. We conclude that an unprecedented combination of accuracy, computational speed and robustness even in the case of maximum activity is accomplished at this stage, with possible further improvements in a foreseeable perspective.