We present a new simulation setup using the MURaM radiative MHD code that allows the study of the formation of collisional polarity inversion lines (cPILs) in the photosphere and the coronal response including flares. In this scheme, we start with a bipolar sunspot configuration and set the spots on collision course by imposing the appropriate velocity field at the footpoints in the subphotospheric boundary. We produce different setups with the same initial spot separation by varying physical parameters such as the collision speed and minimum collision distance. While all setups lead to the formation of an EUV and X-ray sigmoid structure, only the cases with a close passing of the spots cause flares and mass eruptions. The energy release is in the 1–2 × 1031erg range, putting the simulated flares into the upper C-class to lower M-class range of GOES X-ray 1–8 Å flux. While the setup with the more distant passing of the spots does not lead to a flare, the corona is nonetheless substantially heated, suggesting noneruptive energy-release mechanisms. We focus our discussion on two particular setups that differ in spot coherence and resulting cPIL length persistence. We find different timings in the transition from a sheared magnetic arcade to magnetic flux rope (MFR); the setup with a large length but shorter duration cPIL produces a MFR during the eruption, while the MFR is preexisting in the setup with a large length and longer duration cPIL. While both result in flares of comparable strength and the eruption of a coronal mass ejection, the setup with preexisting MFR (and embedded filament) leads to an MFR eruption with a larger mass content.
Observations have shown a clear association of filament/prominence eruptions with the emergence of magnetic flux in or near filament channels. Magnetohydrodynamic (MHD) simulations have been employed to systematically study the conditions under which such eruptions occur. These simulations to date have modeled filament channels as 2D flux ropes or 3D uniformly sheared arcades. Here we present MHD simulations of flux emergence into a more realistic configuration consisting of a bipolar active region containing a line-tied 3D flux rope. We use the coronal flux-rope model of Titov et al. as the initial condition and drive our simulations by imposing boundary conditions extracted from a flux emergence simulation by Leake et al. We identify three mechanisms that determine the evolution of the system: (i) reconnection displacing footpoints of field lines overlying the coronal flux rope, (ii) changes of the ambient field due to the intrusion of new flux at the boundary, and (iii) interaction of the (axial) electric currents in the preexisting and newly emerging flux systems. The relative contributions and effects of these mechanisms depend on the properties of the preexisting and emerging flux systems. Here we focus on the location and orientation of the emerging flux relative to the coronal flux rope. Varying these parameters, we investigate under which conditions an eruption of the latter is triggered.
more » « less- NSF-PAR ID:
- 10490998
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 962
- Issue:
- 2
- ISSN:
- 0004-637X
- Format(s):
- Medium: X Size: Article No. 149
- Size(s):
- ["Article No. 149"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Radiative magnetohydrodynamic simulation includes sufficiently realistic physics to allow for the synthesis of remote sensing observables that can be quantitatively compared with observations. We analyze the largest flare in a simulation of the emergence of large flare-productive active regions described by Chen et al. The flare releases 4.5 × 1031erg of magnetic energy and is accompanied by a spectacular coronal mass ejection. Synthetic soft X-ray flux of this flare reaches M2 class. The eruption reproduces many key features of observed solar eruptions. A preexisting magnetic flux rope is formed along the highly sheared polarity inversion line between a sunspot pair and is covered by an overlying multipole magnetic field. During the eruption, the progenitor flux rope actively reconnects with the canopy field and evolves to the large-scale multithermal flux rope that is observed in the corona. Meanwhile, the magnetic energy released via reconnection is channeled down to the lower atmosphere and gives rise to bright soft X-ray post-flare loops and flare ribbons that reproduce the morphology and dynamic evolution of observed flares. The model helps to shed light on questions of where and when the a flux rope may form and how the magnetic structures in an eruption are related to observable emission properties.
-
Abstract We present both the observation and the magnetohydrodynamics (MHD) simulation of the M2.4 flare (SOL2017-07-14T02:09) of NOAA active region (AR) 12665 with a goal to identify its initiation mechanism. The observation by the Atmospheric Image Assembly (AIA) on board the Solar Dynamics Observatory (SDO) shows that the major topology of the AR is a sigmoidal configuration associated with a filament/flux rope. A persistent emerging magnetic flux and the rotation of the sunspot in the core region were observed with Magnetic Imager (HMI) on board the SDO on the timescale of hours before and during the flare, which may provide free magnetic energy needed for the flare/coronal mass ejection (CME). A high-lying coronal loop is seen moving outward in AIA EUV passbands, which is immediately followed by the impulsive phase of the flare. We perform an MHD simulation using the potential magnetic field extrapolated from the measured pre-flare photospheric magnetic field as initial conditions and adopting the observed sunspot rotation and flux emergence as the driving boundary conditions. In our simulation, a sigmoidal magnetic structure and an overlying magnetic flux rope (MFR) form as a response to the imposed sunspot rotation, and the MFR rises to erupt like a CME. These simulation results in good agreement with the observation suggest that the formation of the MFR due to the sunspot rotation and the resulting torus and kink instabilities were essential to the initiation of this flare and the associated coronal mass ejection.more » « less
-
Abstract Two new schemes for identifying field lines involved in eruptions, the r -scheme and q -scheme, are proposed to analyze the eruptive and confined nature of solar flares, as extensions to the original r m scheme proposed in Lin et al. Motivated by three solar flares originating from NOAA Active Region 12192 that are misclassified by r m , we introduce refinements to the r -scheme employing the “magnetic twist flux” to approximate the force balance acting on a magnetic flux rope (MFR); in the q -scheme, the reconnected field is represented by those field lines that anchor in the flare ribbons. Based on data obtained by the Solar Dynamics Observatory/Helioseismic and Magnetic Imager, the coronal magnetic field for 51 flares larger than M5.0 class, from 29 distinct active regions, is constructed using a nonlinear force-free field extrapolation model. Statistical analysis based on linear discriminant function analysis is then performed, revealing that despite both schemes providing moderately successful classifications for the 51 flares, the coronal mass ejection-eruptivity classification for the three target events can only be improved with the q -scheme. We find that the highly twisted field lines and the flare-ribbon field lines have equal average force-free constant α , but all of the flare-ribbon-related field lines are shorter than 150 Mm in length. The findings lead us to conclude that it is challenging to distinguish the MFR from the ambient magnetic field using any quantity based on common magnetic nonpotentiality measures.more » « less
-
more » « less
Abstract We report on the solar and interplanetary (IP) causes of the third largest geomagnetic storm (26 August 2018) in solar cycle 24. The underlying coronal mass ejection (CME) originating from a quiescent filament region becomes a 440 km/s magnetic cloud (MC) at 1 au after ∼5 days. The prolonged CME acceleration (for ∼24 hr) coincides with the time profiles of the post‐eruption arcade intensity and reconnected flux. Chen et al. (2019,
https://doi.org/10.3847/1538-4357/ab3f36 ) obtain a lower speed since they assumed that the CME does not accelerate after ∼12 hr. The presence of multiple coronal holes near the filament channel and the high‐speed wind from them seem to have the combined effect of producing complex rotation in the corona and IP medium resulting in a high‐inclination MC. The Dst time profile in the main phase steepens significantly (rapid increase in storm intensity) coincident with the density increase (prominence material) in the second half of the MC. Simulations using the Comprehensive Inner Magnetosphere‐Ionosphere model show that a higher ring current energy results from larger dynamic pressure (density) in MCs. Furthermore, the Dst index is highly correlated with the main‐phase time integral of the ring current injection that includes density, consistent with the simulations. A complex temporal structure develops in the storm main phase if the underlying MC has a complex density structure during intervals of southward IP magnetic field. We conclude that the high intensity of the storm results from the prolonged CME acceleration, complex rotation of the CME flux rope, and the high density in the 1‐au MC.
