Abstract. Numerical experiments of the emergence of magnetic flux from the uppermost layers of the solar interior to the photosphere and its further eruption into the low atmosphere and corona are carried out. We use idealized models for the initial stratification and magnetic field distribution below the photosphere similar to those used for multidimensional flux emergence experiments in the literature. The energy equation is adiabatic except for the inclusion of ohmic and viscous dissipation terms, which, however, become important only at interfaces and reconnection sites. Three-dimensional experiments for the eruption of magnetic flux both into an unmagnetized corona and into a corona with a preexisting ambient horizontal field are presented. The shocks preceding the rising plasma present the classical structure of nonlinear Lamb waves. The expansion of the matter when rising into the atmosphere takes place preferentially in the horizontal directions: a flattened (or oval) low plasma-β ball ensues, in which the field lines describe loops in the corona with increasing inclination away from the vertical as one goes toward the sides of the structure. Magnetograms and velocity field distributions on horizontal planes are presented simultaneously for the solar interior and various levels in the atmosphere. Since the background pressure and density drop over many orders of magnitude with increasing height, the adiabatic expansion of the rising plasma yields very low temperatures. To avoid this, the entropy of the rising fluid elements should be increased to the high values of the original atmosphere via heating mechanisms not included in the present numerical experiments. The eruption of magnetic flux into a corona with a preexisting magnetic field pointing in the horizontal direction yields a clear case of essentially three-dimensional reconnection when the upcoming and ambient field systems come into contact. The coronal ambient field is chosen at time t = 0 perpendicular to the direction of the tube axis and thus, given the twist of the magnetic tube, almost anti-parallel to the field lines at the upper boundary of the rising plasma ball. A thin, dome-shaped current layer is formed at the interface between the two flux systems, in which ohmic dissipation and heating are taking place. The reconnection proceeds by merging successive layers on both sides of the reconnection site; however, this occurs not only at the cusp of the interface, but, also, gradually along its sides in the direction transverse to the ambient magnetic field. The topology of the magnetic field in the atmosphere is thereby modified: the reconnected field lines typically are part of the flanks of the tube below the photosphere but then join the ambient field system in the corona and reach the boundaries of the domain as horizontal field lines.

Solar coronal loops are observed to be remarkably stable structures. A magnetohydrodynamic stability analysis of a model loop by the energy method suggests that the main reason for stability is the fact that the ends of the loop are anchored in the dense photosphere. In addition to such line-tying, the effect of a radial pressure gradient is incorporated in the analysis.Two-ribbon flares follow the eruption of an active region filament, which may lie along a magnetic flux tube. It is suggested that the eruption is caused by the kink instability, which sets in when the amount of magnetic twist in the flux tube exceeds a critical value. This value depends on the aspect ratio of the loop, the ratio of the plasma to magnetic pressure and the detailed transverse magnetic structure. For a force-free field of uniform twist the critical twist is 3.3 ~r, and for other fields it is typically between 2~" and 6~r. Occasionally active region loops may become unstable and give rise to small loop flares, which may also be a result of the kink instability.

Aims. We investigate the effects of varying the magnetic field strength and the twist of a flux tube as it rises through the solar interior and emerges into the atmosphere. Methods. Using a 3D numerical MHD code, we consider a simple stratified model, comprising of one solar interior layer and three overlying atmospheric layers. We set a horizontal, twisted flux tube in the lowest layer. The specific balance of forces chosen results in the tube being fully buoyant and the temperature is decreased in the ends of the tube to encourage the formation of an Ω-shape along the tube's length. We vary the magnetic field strength and twist independently of each other so as to give clear results of the individual effects of each parameter. Results. We find a self-similar evolution in the rise and emergence of the flux tube when the magnetic field strength of the tube is modified. During the rise through the solar interior, the height of the crest and axis, the velocity of the crest and axis, and the decrease in the magnetic field strength of the axis of the tube are directly dependent upon the initial magnetic field strength given to the tube. No such self-similarity is evident when the twist of the flux tube is changed, due to the complex interaction of the tension force on the rise of the tube. For low magnetic field strength and twist values, we find that the tube cannot fully emerge into the atmosphere once it reaches the top of the interior since the buoyancy instability criterion cannot be fulfilled. For those tubes that do advance into the atmosphere, when the magnetic field strength has been modified, we find further self-similar behaviour in the amount of tube flux transported into the atmosphere. For the tubes that do emerge, the variation in the twist results in the buoyancy instability, and subsequent emergence, occurring at different locations along the tube's length.

Context. This paper extends the models of Craig & McClymont (1991, ApJ, 371, L41) and McLaughlin & Hood (2004, A&A, 420, 1129 to include finite β and nonlinear effects. Aims. We investigate the nature of nonlinear fast magnetoacoustic waves about a 2D magnetic X-point. Methods. We solve the compressible and resistive MHD equations using a Lagrangian remap, shock capturing code (Arber et al. 2001, J. Comp. Phys., 171, 151) and consider an initial condition in u × B ·ẑ (a natural variable of the system). Results. We observe the formation of both fast and slow oblique magnetic shocks. The nonlinear wave deforms the X-point into a "cusp-like" point which in turn collapses to a current sheet. The system then evolves through a series of horizontal and vertical current sheets, with associated changes in connectivity, i.e. the system exhibits oscillatory reconnection. Our final state is non-potential (but in force balance) due to asymmetric heating from the shocks. Larger amplitudes in our initial condition correspond to larger values of the final current density left in the system. Conclusions. The inclusion of nonlinear terms introduces several new features to the system that were absent from the linear regime.

No abstract

No abstract

The remarkable magnetohydrodynamic stability of solar coronal loops has been attributed lo the anchoring of the ends of loops in the dense photosphere. However, all the previous analyses of such line-tying have been approximate, in the sense that they give only upper or lower bounds on the critical amount of twist (or the critical looplength) required for the breakdown of stability. The object of the present paper is to remove these approximations and determine the exact value for the critical twist. When it is exceeded the magnetic field becomes kink unstable and a flare may be initiated.A simple analytic stability calculation is described for an idealised loop. This is followed by the development of a general numerical technique for any loop profile, which involves solving the partial differential equations of motion. It is found, for example, that a force-free field of uniform twist possesses a critical twist of 2.4971. by comparison with the previous bounds of 2n, for stability, and 3.3 n, for instability.

Abstract. The nature of fast magnetoacoustic and Alfvén waves is investigated in a zero β plasma. This gives an indication of wave propagation in the low β solar corona. It is found that for a two-dimensional null point, the fast wave is attracted to that point and the front of the wave slows down as it approaches the null point, causing the current density to accumulate there and rise rapidly. Ohmic dissipation will extract the energy in the wave at this point. This illustrates that null points play an important role in the rapid dissipation of fast magnetoacoustic waves and suggests the location where wave heating will occur in the corona. The Alfvén wave behaves in a different manner in that the wave energy is dissipated along the separatrices. For Alfvén waves that are decoupled from fast waves, the value of the plasma β is unimportant. However, the phenomenon of dissipating the majority of the wave energy at a specific place is a feature of both wave types.