Three-dimensional Magnetohydrodynamic Model Research Articles

The Dynamic Response of Jovian Magnetotail Reconnection to Enhanced Solar Wind Ram Pressure

In this study, we employ the Grid Agnostic Magnetohydrodynamic (MHD) for Extended Research Applications (GAMERA), a high-resolving-power, three-dimensional global MHD model, to simulate magnetotail reconnection in Jupiter's magnetosphere. While previous satellite observations have provided initial statistics on magnetotail reconnection properties at Jupiter, they have been limited in spacetime coverage, leaving the dynamic process of Jovian magnetotail reconnection and its response to the solar wind (SW) poorly understood. Using MHD simulations, we quantitatively analyze the temporal evolution and spatial dependence of nightside reconnection in Jupiter's magnetotail under ideal quiet and enhanced SW conditions. Our results demonstrate that magnetotail reconnection tends to occur in the midnight and postmidnight sectors, with a low occurrence in the premidnight sector, consistent with both Galileo and Juno observations and predictions by Delamere & Bagenal. The magnetic local time (MLT)–radial distribution of magnetotail reconnection is broad, indicating that Jovian magnetotail reconnection is always dynamic rather than steady state. Enhanced SW ram pressure can decrease the MLT coverage of magnetotail reconnection by compressing Jupiter's magnetosphere. However, the occurrence of magnetotail reconnection near the midnight and postmidnight sectors is enhanced by SW compression beyond 60 R J, but is not significantly impacted by SW compression within 60 R J. Conversely, SW compression suppresses reconnection in the premidnight sector, leading to a stronger dawn–dusk asymmetry in the occurrence and location of magnetotail reconnection. This study validates the applicability of the GAMERA code for simulating Jupiter’s magnetosphere and provides complementary insights into the dynamic structure and the SW response of Jupiter’s magnetosphere.

The Dynamic Response of Jovian Magnetotail Reconnection to Enhanced Solar Wind Ram Pressure

In this study, we employ the Grid Agnostic Magnetohydrodynamic (MHD) for Extended Research Applications (GAMERA), a high-resolving-power, three-dimensional global MHD model, to simulate magnetotail reconnection in Jupiter's magnetosphere. While previous satellite observations have provided initial statistics on magnetotail reconnection properties at Jupiter, they have been limited in spacetime coverage, leaving the dynamic process of Jovian magnetotail reconnection and its response to the solar wind (SW) poorly understood. Using MHD simulations, we quantitatively analyze the temporal evolution and spatial dependence of nightside reconnection in Jupiter's magnetotail under ideal quiet and enhanced SW conditions. Our results demonstrate that magnetotail reconnection tends to occur in the midnight and postmidnight sectors, with a low occurrence in the premidnight sector, consistent with both Galileo and Juno observations and predictions by Delamere & Bagenal. The magnetic local time (MLT)–radial distribution of magnetotail reconnection is broad, indicating that Jovian magnetotail reconnection is always dynamic rather than steady state. Enhanced SW ram pressure can decrease the MLT coverage of magnetotail reconnection by compressing Jupiter's magnetosphere. However, the occurrence of magnetotail reconnection near the midnight and postmidnight sectors is enhanced by SW compression beyond 60 R J, but is not significantly impacted by SW compression within 60 R J. Conversely, SW compression suppresses reconnection in the premidnight sector, leading to a stronger dawn–dusk asymmetry in the occurrence and location of magnetotail reconnection. This study validates the applicability of the GAMERA code for simulating Jupiter’s magnetosphere and provides complementary insights into the dynamic structure and the SW response of Jupiter’s magnetosphere.

Latest Evolution of the X-Ray Remnant of SN 1987A: Beyond the Inner Ring

Based on our Chandra imaging-spectroscopic observations, we present the latest evolution of the X-ray remnant of SN 1987A. Recent changes in the electron temperatures and volume emission measures suggest that the blast wave in SN 1987A is moving out of the dense inner ring structure, also called the equatorial ring (ER). The 0.5–2.0 keV X-ray light curve shows a linearly declining trend (by ∼4.5% yr−1) between 2016 and 2020 as the blast wave heats the hitherto unknown circumstellar medium (CSM) outside the ER. While the peak X-ray emission in the latest 0.3–8.0 keV image is still within the ER, the radial expansion rate in the 3.0–8.0 keV images suggests an increasing contribution of the X-ray emission from less dense CSM since 2012, at least partly from beyond the ER. It is remarkable that, since 2020, the declining soft X-ray flux has stabilized around ∼7 × 10−12 erg s−1 cm−2, which may signal a contribution from the reverse-shocked outer layers of ejecta as predicted by the three-dimensional magnetohydrodynamic models. In the latest ACIS spectrum of supernova remnant 1987A in 2022 we report a significant detection of the Fe K line at ∼6.7 keV, which may be due to changing thermal conditions of the X-ray emitting CSM and/or the onset of reverse shock interactions with the Fe ejecta.

Three-dimensional modeling of multi-component vacuum arc considering anode vapor in actual magnetic field

For vacuum circuit breakers, the anode will change to active mode when the current is high enough. The vapor produced by the anode will become another important source of plasma. In this paper, a three-dimensional (3D) magnetohydrodynamic model considering anode vapor and multiple components is developed to study the vacuum circuit breakers with cup-shaped CuCr contact. The characteristics of arc under the actual magnetic field with different anode temperature distribution and different cathode utilization are simulated. Moreover, the ionization and recombination processes between different kinds of components are considered. The simulation results show that the arc current concentrates toward the center. The current value is high between electrodes slots under the actual magnetic field, which significantly affects the arc characteristics. At the same time, higher anode temperatures will produce more neutral atom vapor area (NAVA). Current is concentrated outside NAVA due to the low conductivity of NAVA. As NAVA expands, the constriction of current on the outside is more severe. The simulation results in this paper are consistent with the experimental results.

Influence of an Upstream Dynamic Pressure Pulse on Venus–Solar Wind Interaction by Multispecies Magnetohydrodynamic Simulation

As an unmagnetized planet, Venus does not have a global dipole magnetic field like Earth. The solar wind interacts directly with the Venusian ionosphere. In this paper, the interaction between solar wind and the Venusian upper atmosphere is investigated by conducting a three-dimensional multispecies magnetohydrodynamic (MHD) model, in which the Venusian ionosphere is developed self-consistently by considering chemical reactions among four major ion species (H+, O2 +, O+, and CO2 +). Based on this developed MHD model, we mainly study the effects of the impulsive solar wind dynamic pressure on the interaction process. Simulation results reveal that the locations of the bow shock and magnetic pileup boundary show an abrupt compression, then expansion, and finally shrinkage to reach a new steady state with the influence of the impulsive dynamic pressure. The recovery time of the large-scale magnetic field is a few minutes. In addition, the peak escape fluxes are enhanced by about 1 order of magnitude due to the increased solar wind dynamic pressure pulse.

Influence of Crustal Magnetic Fields on Horizontal Plasma Transport and Ion Escape on Mars

Owing to its unevenly distributed crustal fields, Mars acts as a unique obstacle to the solar wind. In the presence of the crustal fields, the transport of the planetary ions on the dayside ionosphere exhibits north–south asymmetry. Additionally, the heavy-ion loss in the magnetotail is affected by the crustal fields. In this paper, a three-dimensional multispecies magnetohydrodynamic model is employed to simulate Mars–solar wind interactions. Numerical results indicate that the meridional transport is dominant in most areas on the dayside ionosphere. In the presence of the crustal fields, the meridional transport on the southern hemisphere (southward transport) is reduced by more than 70% above the strong crustal sources, and the zonal velocity shows local changes inside strong and weak crustal field regions. These effects result in an increase or decrease in the number density of the heavy ions reaching the terminator, thereby influencing the thickness of the ionosphere. Decreased southward velocity leads to a reduction in the heavy-ion loss on the southern magnetotail. The radial outward flux is reduced by more than 30% for O2 + and CO2 + and by 10% for O+. This study shows that in addition to the zonal transport, the meridional transport is important for the day-to-night transport on the dayside of Mars. Collectively, the horizontal plasma transport, controlled by crustal fields, is associated with the altered ionosphere structure and reduced heavy-ion loss in the magnetotail.

Modeling and simulation study on the ablation characteristics of DC contactor contact

The DC contactor can frequently break the high current control circuit, and the arc motion characteristic in the current breaking process is closely related to its electrical life and reliability. In this paper, a three-dimensional magnetohydrodynamic arc model of the DC contactor considering contact motion and external circuit characteristics is established. On this basis, the influence of the vapor produced by the contact ablation on the arc behavior is studied, and the ablation rate of the cathode in a single breaking process is calculated. The results show that the arc temperature is reduced, the motion speed is slowed down, and the arc burning time is shortened due to the presence of copper vapor produced by ablation. The ablation rate of the cathode during a single arc is 8.9 µg/C.

Influence of the interstellar magnetic field and 11-year cycle of solar activity on the heliopause nose location

Context. The heliosphere is formed by the interaction between the solar wind (SW) plasma emanating from the Sun and a magnetised component of local interstellar medium (LISM) inflowing on the Sun. A separation surface called the heliopause (HP) forms between the SW and the LISM. Aims. In this article, we define the nose of the HP and investigate the variations in its location. These result from a dependence on the intensity and direction of the interstellar magnetic field (ISMF), which is still not well known but has a significant impact on the movement of the HP nose, as we try to demonstrate in this paper. Methods. We used a parametric study method based on numerical simulations of various forms of the heliosphere using a time-dependent three-dimensional magnetohydrodynamic (3D MHD) model of the heliosphere. Results. The results confirm that the nose of the HP is always in a direction that is perpendicular to the maximum ISMF intensity directly behind the HP. The displacement of the HP nose depends on the direction and intensity of the ISMF, with the structure of the heliosphere and the shape of the HP depending on the 11-year cycle of solar activity. Conclusions. In the context of the planned space mission to send the Interstellar Probe (IP) to a distance of 1000 AU from the Sun, our study may shed light on the question as to which direction the IP should be sent. Further research is needed that introduces elements such as current sheet, reconnection, cosmic rays, instability, or turbulence into the models.

Heliospheric 3-D MHD ENLIL simulations of multi-CME and multi-spacecraft events

Interpreting multi-spacecraft heliospheric observations of the evolving solar wind (SW) streams with propagating and interacting coronal mass ejections (CMEs) is challenging. Numerical simulations can provide global context and suggest what may and may not be observed. The heliospheric three-dimensional (3D) magnetohydrodynamic (MHD) ENLIL model can provide a near-real-time prediction of heliospheric space weather, and it is used at NASA Community Coordinated Modeling Center (CCMC), NOAA Space Weather Prediction Center (SWPC), and UK Meteorological Office (MetOffice). However, this version does not show its full potential, especially in the case of multi-CME events observed by various spacecraft. We describe tools developed to interpret remote observations and in-situ measurements better and apply them to multi-CME events observed by ACE, STEREO-A, Parker Solar Probe (PSP), BepiColombo, and Solar Orbiter. We present some results on 1) global structures of the SW speed and density at the ecliptic, 2) the evolution of SW parameters at the spacecraft, 3) magnetic field connectivity at the spacecraft, 4) automatic detection of shock parameters and alert plots, and 5) synthetic white-light (WL) imaging. This paper is not on model initialization or analyzing specific CME events, but it describes features not used at space weather prediction centers and provided by NASA/CCMC Run-On-Request service. This paper advertises new tools and shows their benefits when applied to selected heliospheric space weather events observed at near-Earth, PSP, Solar Orbiter, and STEREO-A spacecraft.

3D physico-chemical model of a pre-stellar core

Context. Pre-stellar cores represent the earliest stage of the formation process of stars and planets. By characterizing the physical and chemical structure of these cores, we can establish the initial conditions for star and planet formation and determine to what degree the chemical composition of pre-stellar cores is inherited by the later stages. Aims. We aim to determine the underlying causes of spatial chemical segregation observed in pre-stellar cores and study the effects of the core structure and external environment on the chemical structure of pre-stellar cores. Methods. A three-dimensional (3D) magnetohydrodynamic model of a pre-stellar core embedded in a dynamic star-forming cloud was post-processed with a sequentially continuum radiative transfer, a gas-grain chemical model, and a line-radiative transfer model. The results were analyzed and compared to observations of CH3OH and c-C3H2 in L1544. We compared nine different chemical models to the observations to determine which initial conditions are compatible with the observed chemical segregation in the prototypical pre-stellar core L1544. Results. The model is able to reproduce several aspects of the observed chemical differentiation in L1544. Extended methanol emission is shifted towards colder and more shielded regions of the core envelope, while c-C3H2 emission overlaps with the dust continuum, which is consistent with the observed chemical structure. Furthermore, these results are consistent across a broad spectrum of chemical models. Increasing the strength of the interstellar radiation field or the cosmic-ray ionization rate with respect to the typical values assumed in nearby star-forming regions leads to synthetic maps that are inconsistent with the observed chemical structure. Conclusions. Our model shows that the observed chemical dichotomy in L1544 can arise as a result of uneven illumination due to the asymmetrical structure of the 3D core and the environment within which the core has formed. This highlights the importance of the 3D structure at the core-cloud transition on the chemistry of pre-stellar cores. The reported effect is likely to affect later stages of the formation process of stars and planets through chemical inheritance.

Influence of the Martian crustal magnetic fields on the Mars-solar wind interaction and plasma transport

The plasma transport process is important for the ionosphere of Mars, which controls the structure of the ionosphere above an altitude of 200 km. Plasma transport from the dayside ionosphere is crucial for producing the nightside ionosphere on Mars. The alteration in dayside plasma transport in the presence of crustal fields may influence the distribution of Martian ionospheric plasma and plasma escape in the magnetotail. This study employed a three-dimensional multispecies magnetohydrodynamic (MHD) model to simulate Mars-solar wind interactions. We show the magnetic field distribution and plasma velocity variation on the Martian day-side. The results indicate that the ion transport from low- to high-solar-zenith-angle areas in the south is inhibited by crustal fields, leading to a reduction in the ion number density and a thinner ionosphere near the southern terminator. Many heavy ions remain in the southern dayside ionosphere rather than moving to the nightside. In addition, the maximum reduction in the tailward flux of the planetary ions calculated by the MHD simulation is more than 50% at the southern terminator, indicating an inhibitory effect of the crustal fields on day-to-night transport. These effects may lead to a reduction in ion number density in the southern nightside ionosphere. Finally, we demonstrate a decrease in the global heavy-ion loss rate.

Three-dimensional core-collapse supernovae with complex magnetic structures – II. Rotational instabilities and multimessenger signatures

ABSTRACT The gravitational collapse of rapidly rotating massive stars can lead to the onset of the low T/|W| instability within the central proto-neutron star (PNS), which leaves strong signatures in both the gravitational wave (GW) and neutrino emission. Strong large-scale magnetic fields are usually invoked to explain outstanding stellar explosions of rapidly rotating progenitors, but their impact on the growth of such instability has not yet been cleared. We analyse a series of three-dimensional magnetohydrodynamic models to characterize the effects of different magnetic configurations on the development of the low T/|W| and the related multimessenger features. In the absence of magnetic fields, we observe the growth on dynamical time-scales of the low T/|W|, associated with a strong burst of GW and a correlated modulation of the neutrino emission. However, models with a strong magnetic field show a quenching of the low T/|W|, due to a flattening of the rotation profile in the first ∼100 ms after shock formation caused by the magnetic transport of angular momentum. The associated GW emission is weakened by an order of magnitude, exhibits a broader spectral shape, and has no dominant feature associated with the PNS large-scale oscillation modes. Neutrino luminosities are damped along the equatorial plane due to a more oblate PNS, and the only clear modulation in the signal is due to Standing Accretion Shock Instability activity. Finally, magnetized models produce lower luminosities for νe than for $\bar{\nu }_e$, which is connected to a higher concentration of neutron-rich material in the PNS surroundings.

Rapid blue- and redshifted excursions in Hα line profiles synthesized from realistic 3D magnetohydrodynamic simulations

Context. Rapid blue- and redshifted excursions (RBEs and RREs) may play an important role in mass-loading and heating the solar corona, but their nature and origin are still debatable. Aims. We aim to model these features to learn more about their properties, formation, and origin. Methods. We created a realistic three-dimensional (3D) magnetohydrodynamic model of a solar plage region. Synthetic Hα spectra were generated and the spectral signatures of these features identified. The magnetic field lines associated with these events were traced, and the underlying dynamic was studied. Results. The model reproduces many properties of RBEs and RREs well, such as spatial distribution, lateral movement, length, and lifetimes. Synthetic Hα line profiles, similarly to observed ones, show a strong blue- or redshift as well as asymmetries. These line profiles are caused by the vertical component of velocities higher than 30 − 40 km s−1, which mostly appear in the height range 2 − 4 Mm. By tracing magnetic field lines, we show that the vertical velocity that causes the appearance of RBEs or RREs is always associated with the component of velocity perpendicular to the magnetic field lines. Conclusions. The study confirms the hypothesis that RBEs and RREs are signs of Alfvénic waves with, in some cases, a significant contribution from slow magneto-acoustic modes.

Nonlinear Fast Magnetosonic Waves in Solar Prominence Pillars

We investigate the properties of nonlinear fast magnetosonic (NFM) waves in a solar prominence, motivated by recent high-resolution and high-cadence Hinode/Solar Optical Telescope (SOT) observations of small-scale oscillations in a prominence pillar. As an example, we analyze the details of the 2012 February 14 Hinode/SOT observations of quasi-periodic propagating features consistent with NFM waves, imaged in emission in Ca ii and in the far blue wing of Hα. We perform wavelet analysis and find oscillations in the 1–3 minutes period range. Guided by these observations, we model the NFM waves with a three-dimensional magnetohydrodynamics (3D MHD) model, extending previous 2.5D MHD studies. The new model includes the structure of the high-density, low-temperature material of the prominence pillar embedded in the hot corona, in both potential and non-force-free sheared magnetic field configurations. The nonlinear model demonstrates the effects of mode coupling and the propagating density compressions associated with linear and NFM waves. The guided fast magnetosonic waves, together with density compressions and currents, are reproduced in the 3D pillar structure. We demonstrate for the first time the dynamic effects of the Lorentz force due to the magnetic shear in the non-force-free field on the pillar structure and on the propagation of the waves. The insights gained from the 3D MHD modeling are useful for improving the coronal seismology of prominence structures that exhibit fast MHD wave activity.

X-Ray Morphology Due to Charge-exchange Emissions Used to Study the Global Structure around Mars

Soft X-ray emissions induced by solar wind ions that collide with neutral material in the solar system have been detected around planets, and were proposed as a remote probe for solar wind interaction with the Martian exosphere. A multi-fluid three-dimensional magnetohydrodynamic model is adopted to derive the global distributions of solar wind particles. Spherically symmetric exospheric H, H2, He, O, and CO2 density profiles and a sophisticated hybrid model that includes charge-exchange and proton–neutral excitation processes are used to study the low triplet line ratio (0.77 ± 0.58) of O vii and the total X-ray luminosity around Mars. We further calculate the emission factor α-value with different neutrals over wide ion-abundance and velocity ranges. Our results are in good agreement with those of previous reports. The evolution of the charge stage of solar wind ions shows that sequential recombination due to charge-exchange can be negligible in the interaction region. This only appears below an altitude of 400 km. The anonymous low disk G-ratio can be easily explained by the collisional quenching effect at neutral densities higher than 1011 cm−3. However, the quenching contribution is small in Mars’ exosphere and only appears below 400 km. Charge-exchange with H2 and N2 is still the most likely reason for this low G-ratio. X-ray emissivity maps in collisions with different neutrals differ from each other. A clear bow shock arising from the collision with all the neutrals is in accordance with previous reports. The resulting total X-ray luminosity of 6.55 MW shows better agreement with the XMM-Newton observation of 12.8 ± 1.4 MW than that of previous predictions.

A magnetohydrodynamic simulation of the dayside magnetic reconnection between the solar wind and the Martian crustal field

Using a three-dimensional multispecies magnetohydrodynamic model, we study the effects of the orientation of the interplanetary magnetic field (IMF), solar wind dynamic pressure (Pd), and the location of the intense crustal field, on the dayside magnetic reconnection between the solar wind and the Martian crustal field. Our main results are as follows: (1) Different IMF orientations result in different magnetic field configurations and reconnection conditions on the Martian dayside. When the intense crustal field is located on the dayside, the dayside magnetic reconnection tends to occur in the region with solar zenith angles (SZA) ≈45° in the southern hemisphere for the IMF with a southward component. When the IMF has a northward component, the magnetic field lines are piled up in the same place and the Martian magnetic pileup boundary (MPB) appears as a local bulged “mini-magnetopause”. Under the pure radial IMF, the magnetic reconnection is absent, which might be due to the presence of additional outward magnetic tension and kinetic effects. (2) Dayside reconnection can change the shape of the Martian MPB, while the bow shock is weakly affected. When the IMF has a southward component, the dayside magnetic reconnection happens and the MPB is located closer to Mars with a “cusp” shape. When the IMF has a northward component, the Martian MPB expands with a local bulged “mini-magnetopause”. For the pure radial IMF condition, the subsolar region of the MPB is located closer to Mars than that under other IMF directions. The influence of the IMF cone angles on the Martian bow shock is much less than that on the MPB, and the bow shock locations are very close to the model results of another author found in the literature. (3) With increasing Pd, the size of the crustal field region decreases and the draped fields correspondingly move to lower altitudes where the IMF and crustal field have the same direction. When the IMF has a southward component and the magnetic reconnection occurs at SZA ≈ 45°, the reconnection site, the region of the closed topology of the crustal field, and the draped IMF, do not change much with increasing Pd. We suggest that the multipolar crustal magnetic fields can protect the solar wind IMF from further reconnecting with the crustal field to a lower altitude when Pd is enhanced.

Diverse Polarimetric Features of AGN Jets from Various Viewing Angles: Towards a Unified View

Here, we demonstrate that polarization properties show a wide diversity depending on viewing angles. To simulate images of a supermassive black hole and surrounding plasma, we performed a full-polarimetric general relativistic radiative transfer based on three-dimensional general relativistic magnetohydrodynamics models with moderate magnetic strengths. Under an assumption of a hot-jet and cold-disk in the electron temperature prescription, we confirmed a typical scenario where polarized synchrotron emissions from the funnel jet experience Faraday rotation and conversion in the equatorial disk. Further, we found that linear polarization vectors are inevitably depolarized for edge-on-like observers, whereas a portion of vectors survive and reach the observers in face-on-like cases. We also found that circular polarization components have persistent signs in the face-on cases, and changing signs in the edge-on cases. It is confirmed that these features are smoothly connected via intermediate viewing-angle cases. These results are due to Faraday rotation/conversion for different viewing angles, and suggest that a combination of linear and circular polarimetry can give a constraint on the inclination between the observer and black hole’s (and/or disk’s) rotating-axis and plasma properties in the jet–disk structure. These can also lead to a more statistical and unified interpretation for a diversity of emissions from active galactic nuclei.

Predicting solar wind streams from the inner-heliosphere to Earth via shifted operator inference

Solar wind conditions are predominantly predicted via three-dimensional numerical magnetohydrodynamic (MHD) models. Despite their ability to produce highly accurate predictions, MHD models require computationally intensive high-dimensional simulations. This renders them inadequate for making time-sensitive predictions and for large-ensemble analysis required in uncertainty quantification. This paper presents a new data-driven reduced-order model (ROM) capability for forecasting heliospheric solar wind speeds. Traditional model reduction methods based on Galerkin projection have difficulties with advection-dominated systems—such as solar winds—since they require a large number of basis functions and can become unstable. A core contribution of this work addresses this challenge by extending the non-intrusive operator inference ROM framework to exploit the translational symmetries present in the solar wind caused by the Sun's rotation. The numerical results show that our method can adequately emulate the MHD simulations and is more accurate than a reduced-physics surrogate model, the Heliospheric Upwind Extrapolation model.

ICARUS, a new inner heliospheric model with a flexible grid

Context. Simulating the propagation and predicting the arrival time of coronal mass ejections (CMEs) in the inner heliosphere with a full three-dimensional (3D) magnetohydrodynamic (MHD) propagation model requires a significant amount of computational time. For CME forecasting purposes, multiple runs may be required for different reasons such as ensemble modeling (uncertainty on input parameters) and error propagation. Moreover, higher resolution runs may be necessary, which also requires more CPU time, for example for the prediction of solar energetic particle acceleration and transport or in the framework of more in-depth studies about CME erosion and/or deformation during its evolution. Aims. In this paper we present ICARUS, a new inner heliospheric model for the simulation of a steady background solar wind and the propagation and evolution of superposed CMEs. This novel model has been implemented within the MPI-AMRVAC framework which enables the use of stretched grids and solution adaptive mesh refinement (AMR). The usefulness and efficiency (speed-up) of these advanced features are explored. In particular, we model a typical solar wind with ICARUS and then launch a simple cone CME and follow its evolution. We focus on the effect of radial grid stretching and two specific methods or criteria to trigger solution AMR on this typical simulation run. Methods. For the solar background wind simulation run, we limited the mesh refinement to the area(s) of interest, in this case a co-rotating interaction region (CIR). For the CME evolution run, on the other hand, we apply AMR where the CME is located by the use of a tracing function. As such, the grid is coarsened again after the CME has passed. Results. The implemented AMR is flexible and only refines the mesh in a particular sector of the computational domain, for example around the Earth or a single CIR, and/or for a particular feature such as CIR or CME shocks. Radial grid stretching alone yields speed-ups of up to 4 and more, depending on the resolution. Combined with solution adaptive mesh refinement, the speed-ups can be much larger depending on the complexity of the simulation (e.g., number of CIRs in the background wind, number of CMEs) and on the chosen AMR criteria, thresholds and the number of refinement levels. Conclusions. The ICARUS model implemented in the MPI-AMRVAC framework is a new inner heliospheric 3D MHD model that uses grid stretching as well as AMR techniques. The flexibility in the grid and its resolution allows an optimization of the computational time required for CME propagation simulations for both scientific and forecasting purposes.

The origin of underdense plasma downflows associated with magnetic reconnection in solar flares

Magnetic reconnection is a universal process that powers explosive energy-release events such as solar flares, geomagnetic substorms and some astrophysical jets. A characteristic feature of magnetic reconnection is the production of fast reconnection outflow jets near the plasma Alfvén speeds1,2. In eruptive solar flares, dark finger-shaped plasma downflows moving toward the flare arcade have been commonly regarded as the principal observational evidence for such reconnection-driven outflows3,4. However, they often show a speed much slower than that expected in reconnection theories5,6, challenging the reconnection-driven energy-release scenario in standard flare models. Here we present a three-dimensional magnetohydrodynamics model of solar flares. By comparing the model predictions with the observed plasma downflow features, we conclude that these dark downflows are self-organized structures formed in a turbulent interface region below the flare termination shock where the outflows meet the flare arcade, a phenomenon analogous to the formation of similar structures in supernova remnants. This interface region hosts a myriad of turbulent flows, electron currents and shocks, crucial for flare energy release and particle acceleration. Supra-arcade downflows are dark finger-like structures associated with solar flares. Simulations now suggest that these plasma downflows form within a turbulent region that may control turbulent flows, electron currents and shockwaves.