No Arabic abstract
Solar flares are an explosive phenomenon, where super-sonic flows and shocks are expected in and above the post-flare loops. To understand the dynamics of post-flare loops, a two-dimensional magnetohydrodynamic (2D MHD) simulation of a solar flare has been carried out. We found new shock structures in and above the post-flare loops, which were not resolved in the previous work by Yokoyama and Shibata 2001. To study the dynamics of flows along the reconnected magnetic field, kinematics and energetics of the plasma are investigated along selected field lines. It is found that shocks are crucial to determine the thermal and flow structures in the post-flare loops. On the basis of the 2D MHD simulation, we have developed a new post-flare loop model which we call the pseudo-2D MHD model. The model is based on the 1D MHD equations, where all the variables depend on one space dimension and all the three components of the magnetic and velocity fields are considered. Our pseudo-2D model includes many features of the multi-dimensional MHD processes related to magnetic reconnection (particularly MHD shocks), which the previous 1D hydrodynamic models are not able to include. We compare the shock formation and energetics of a specific field line in the 2D calculation with those in our pseudo-2D MHD model, and we found that they give similar results. This model will allow us to study the evolution of the post-flare loops in a wide parameter space without expensive computational cost and without neglecting important physics associated with magnetic reconnection.
We investigated the dynamic evolution of a 3-dimensional (3D) flux rope eruption and magnetic reconnection process in a solar flare, by simply extending 2-dimensional (2D) resistive magnetohydrodynamic simulation model of solar flares with low $beta$ plasma to 3D model. We succeeded in reproducing a current sheet and bi-directional reconnection outflows just below the flux rope during the eruption in our 3D simulations. We calculated four cases of a strongly twisted flux rope and a weakly twisted flux rope in 2D and 3D simulations. The time evolution of a weakly twisted flux rope in 3D simulation shows similar behaviors to 2D simulation, while a strongly twisted flux rope in 3D simulation shows clearly different time evolution from 2D simulation except for the initial phase evolution. The ejection speeds of both strongly and weakly twisted flux ropes in 3D simulations are larger than 2D simulations, and the reconnection rates in 3D cases are also larger than 2D cases. This indicates a positive feedback between the ejection speed of a flux rope and the reconnection rate even in the 3D simulation, and we conclude that the plasmoid-induced reconnection model can be applied to 3D. We also found that small scale plasmoids are formed inside a current sheet and make it turbulent. These small scale plasmoid ejections has role in locally increasing reconnection rate intermittently as observed in solar flares, coupled with a global eruption of a flux rope.
The suns chromosphere is a highly dynamic, partially-ionized region where spicules (hot jets of plasma) form. Here we present a two-fluid MHD model to study the chromosphere, which includes ion-neutral interaction and frictional heating. Our simulation recovers a magnetic canopy shape that forms quickly, but is also quickly disrupted by the formation of a jet. Our simulation produces a shock self-consistently, where the jet is driven by the frictional heating, which is much greater than the ohmic heating. Thus, our simulation demonstrates that the jet could be driven purely by thermal effects due to ion-neutral collisions and not by magnetic reconnection. We plan to improve the model to include photo-chemical effects and radiation.
We perform a zero-$beta$ magnetohydrodynamic simulation for the C7.7 class flare initiated at 01:18 UT on 2011 June 21 using the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). The initial condition for the simulation involves a flux rope which we realize through the regularized Biot-Savart laws, whose parameters are constrained by observations from the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO) and the Extreme Ultraviolet Imager (EUVI) on the twin Solar Terrestrial Relations Observatory (STEREO). This data-constrained initial state is then relaxed to a force-free state by the magneto-frictional module in MPI-AMRVAC. The further time-evolving simulation results reproduce the eruption characteristics obtained by SDO/AIA 94 A, 304 A, and STEREO/EUVI 304 A observations fairly well. The simulated flux rope possesses similar eruption direction, height range, and velocities to the observations. Especially, the two phases of slow evolution and fast eruption are reproduced by varying the density distribution in light of the filament material draining process. Our data-constrained simulations also show other advantages, such as a large field of view (about 0.76 solar radii). We study the twist of the magnetic flux rope and the decay index of the overlying field, and find that in this event, both the magnetic strapping force and the magnetic tension force are sufficiently weaker than the magnetic hoop force, thus allowing the successful eruption of the flux rope. We also find that the anomalous resistivity is necessary in keeping the correct morphology of the erupting flux rope.
The magnetohydrodynamics of active region NOAA 11283 is simulated using an initial non-force-free magnetic field extrapolated from its photospheric vector magnetogram. We focus on the magnetic reconnections at a magnetic null point that participated in the X2.1 flare on 2011 September 6 around 22:21 UT (SOL2011-09-06T22:21X2.1) followed by the appearance of circular flare ribbons and coronal dimmings. The initial magnetic field from extrapolation displays a three-dimensional (3D) null topology overlying a sheared arcade. Prior to the flare, magnetic loops rise due to the initial Lorentz force, and reconnect at the 3D null, leading to expansion and loss of confined plasma that produce the observed pre-flare coronal dimmings. Further, the simulated dynamics documents the transfer of twist from the arcade to the overlying loops through reconnections, developing a flux rope. The non-parallel field lines comprising the rope and lower-lying arcades form an X-type geometry. Importantly, the simultaneous reconnections at the 3D null and the X-type geometry can explain the observed circular and parallel flare ribbons. Reconnections at the 3D null transform closed inner spine field lines into open field lines of the outer spine. The footpoints of these open field lines correspond to a ring-shaped coronal dimming region, tracing the dome. Further, the flux rope bifurcates because of these reconnections which also results in the generation of open magnetic field lines. The plasma loss along the open field lines can potentially explain the observed coronal dimming.
We analyse the oscillatory properties of resonantly damped transverse kink oscillations in two-dimensional prominence threads. The fine structures are modelled as cylindrically symmetric magnetic flux tubes with a dense central part with prominence plasma properties and an evacuated part, both surrounded by coronal plasma. The equilibrium density is allowed to vary non-uniformly in both the transverse and the longitudinal directions.We examine the influence of longitudinal density structuring on periods, damping times, and damping rates for transverse kink modes computed by numerically solving the linear resistive magnetohydrodynamic (MHD) equations. The relevant parameters are the length of the thread and the density in the evacuated part of the tube, two quantities that are difficult to directly estimate from observations. We find that both of them strongly influence the oscillatory periods and damping times, and to a lesser extent the damping ratios. The analysis of the spatial distribution of perturbations and of the energy flux into the resonances allows us to explain the obtained damping times. Implications for prominence seismology, the physics of resonantly damped kink modes in two-dimensional magnetic flux tubes, and the heating of prominence plasmas are discussed.