We present 3D zero-beta ideal MHD simulations of the solar flare/CME event that occurred in Active Region 11060 on 2010 April 8. The initial magnetic configurations of the two simulations are stable nonlinear force-free field and unstable magnetic field models constructed by Su et al. (2011) using the flux rope insertion method. The MHD simulations confirm that the stable model relaxes to a stable equilibrium, while the unstable model erupts as a CME. Comparisons between observations and MHD simulations of the CME are also presented.

Understanding the magnetic configuration of the source regions of coronal mass ejections (CMEs) is vital in order to determine the trigger and driver of these events. Observations of four CME productive active regions are presented here, which indicate that the pre-eruption magnetic configuration is that of a magnetic flux rope. The flux ropes are formed in the solar atmosphere by the process known as flux cancellation and are stable for several hours before the eruption. The observations also indicate that the magnetic structure that erupts is not the entire flux rope as initially formed, raising the question of whether the flux rope is able to undergo a partial eruption or whether it undergoes a transition in specific flux rope configuration shortly before the CME.

In this paper, we address the formation of a magnetic flux rope (MFR) that erupted on 2012 July 12 and caused a strong geomagnetic storm event on July 15. Through analyzing the long-term evolution of the associated active region observed by the Atmospheric Imaging Assembly and the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory, it is found that the twisted field of an MFR, indicated by a continuous S-shaped sigmoid, is built up from two groups of sheared arcades near the main polarity inversion line a half day before the eruption. The temperature within the twisted field and sheared arcades is higher than that of the ambient volume, suggesting that magnetic reconnection most likely works there. The driver behind the reconnection is attributed to shearing and converging motions at magnetic footpoints with velocities in the range of 0.1-0.6 km s(-1). The rotation of the preceding sunspot also contributes to the MFR buildup. Extrapolated three-dimensional non-linear force-free field structures further reveal the locations of the reconnection to be in a bald-patch region and in a hyperbolic flux tube. About 2 hr before the eruption, indications of a second MFR in the form of an S-shaped hot channel are seen. It lies above the original MFR that continuously exists and includes a filament. The whole structure thus makes up a stable double-decker MFR system for hours prior to the eruption. Eventually, after entering the domain of instability, the high-lying MFR impulsively erupts to generate a fast coronal mass ejection and X-class flare; while the low-lying MFR remains behind and continuously maintains the sigmoidicity of the active region.

The structure of the coronal magnetic field prior to eruptive processes and the conditions for the onset of eruption are important issues that can be addressed through studying the magnetohydrodynamic (MHD) stability and evolution of nonlinear force-free field (NLFFF) models. This paper uses data-constrained NLFFF models of a solar active region (AR) that erupted on 2010 April 8 as initial conditions in MHD simulations. These models, constructed with the techniques of flux rope insertion and magnetofrictional relaxation (MFR), include a stable, an approximately marginally stable, and an unstable configuration. The simulations confirm previous related results of MFR runs, particularly that stable flux rope equilibria represent key features of the observed pre-eruption coronal structure very well, and that there is a limiting value of the axial flux in the rope for the existence of stable NLFFF equilibria. The specific limiting value is located within a tighter range, due to the sharper discrimination between stability and instability by the MHD description. The MHD treatment of the eruptive configuration yields a very good agreement with a number of observed features, like the strongly inclined initial rise path and the close temporal association between the coronal mass ejection and the onset of flare reconnection. Minor differences occur in the velocity of flare ribbon expansion and in the further evolution of the inclination; these can be eliminated through refined simulations. We suggest that the slingshot effect of horizontally bent flux in the source region of eruptions can contribute significantly to the inclination of the rise direction. Finally, we demonstrate that the onset criterion, formulated in terms of a threshold value for the axial flux in the rope, corresponds very well to the threshold of the torus instability in the considered AR.

The rotation of erupting filaments in the solar corona is addressed through a parametric simulation study of unstable, rotating flux ropes in bipolar force-free initial equilibrium. The Lorentz force due to the external shear-field component and the relaxation of tension in the twisted field are the major contributors to the rotation in this model, while reconnection with the ambient field is of minor importance, due to the field's simple structure. In the low-beta corona, the rotation is not guided by the changing orientation of the vertical field component's polarity inversion line with height. The model yields strong initial rotations which saturate in the corona and differ qualitatively from the profile of rotation vs. height obtained in a recent simulation of an eruption without preexisting flux rope. Both major mechanisms writhe the flux rope axis, converting part of the initial twist helicity, and produce rotation profiles which, to a large part, are very similar within a range of shear-twist combinations. A difference lies in the tendency of twist-driven rotation to saturate at lower heights than shear-driven rotation. For parameters characteristic of the source regions of erupting filaments and coronal mass ejections, the shear field is found to be the dominant origin of rotations in the corona and to be required if the rotation reaches angles of order 90 degrees and higher; it dominates even if the twist exceeds the threshold of the helical kink instability. The contributions by shear and twist to the total rotation can be disentangled in the analysis of observations if the rotation and rise profiles are simultaneously compared with model calculations. The resulting twist estimate allows one to judge whether the helical kink instability occurred. This is demonstrated for the erupting prominence in the "Cartwheel CME" on 9 April 2008, which has shown a rotation of a parts per thousand aEuro parts per thousand 115(a similar to) up to a height of 1.5 R (aS (TM)) above the photosphere. Out of a range of initial equilibria which include strongly kink-unstable (twist I broken vertical bar=5 pi), weakly kink-unstable (I broken vertical bar=3.5 pi), and kink-stable (I broken vertical bar=2.5 pi) configurations, only the evolution of the weakly kink-unstable flux rope matches the observations in their entirety.