Institut für Physik und Astronomie
Refine
Document Type
- Article (14)
- Doctoral Thesis (1)
Language
- English (15)
Is part of the Bibliography
- yes (15)
Keywords
- magnetohydrodynamics (MHD) (15) (remove)
Institute
There is a recurring question in solar physics regarding whether or not electric currents are neutralized in active regions (ARs). This question was recently revisited using three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetic flux emergence into the solar atmosphere. Such simulations showed that flux emergence can generate a substantial net current in ARs. Other sources of AR currents are photospheric horizontal flows. Our aim is to determine the conditions for the occurrence of net versus neutralized currents with this second mechanism. Using 3D MHD simulations, we systematically impose line-tied, quasi-static, photospheric twisting and shearing motions to a bipolar potential magnetic field. We find that such flows: (1) produce both direct and return currents, (2) induce very weak compression currents-not observed in 2.5D-in the ambient field present in the close vicinity of the current-carrying field, and (3) can generate force-free magnetic fields with a net current. We demonstrate that neutralized currents are in general produced only in the absence of magnetic shear at the photospheric polarity inversion line-a special condition that is rarely observed. We conclude that. photospheric flows,. as magnetic flux emergence, can build up net currents in the solar atmosphere, in agreement with recent observations. These results thus provide support for eruption models based on pre-eruption magnetic fields that possess a net coronal current.
Solar-like stars maintain their magnetic fields thanks to a dynamo mechanism. The Babcock-Leighton dynamo is one possible dynamo that has the particularity to require magnetic flux tubes. Magnetic flux tubes are assumed to form at the bottom of the convective zone and rise buoyantly to the surface. A delayed dynamo model has been suggested, where the delay accounts for the rise time of the magnetic flux tubes; a time, that has been ignored by former studies.
The present thesis aims to study the applicability of the flux tube/Babcock-Leighton dynamo to other stars. To do so, we attempt to constrain the rise time of magnetic flux tubes thanks to the first fully compressible MHD simulations of rising magnetic flux tubes in stratified rotating spherical shells.
Such simulations are limited to an unrealistic parameter space, therefore, a scaling relation is required to scale the results to realistic physical regimes. We extended earlier works on 2D scaling relations and derived a general scaling law valid for both 2D and 3D. We then carried out two large series of numerical experiments and verified that the scaling law we have derived indeed applies to the fully non-linear case. It allowed us to extract a constraint for the rise time of magnetic flux tubes that is valid for any solar-like star. We finally introduced this constraint to a delayed dynamo model.
By carrying out simulations of a mean-field, delayed, flux tube/Babcock-Leighton dynamo, we were able to identify a new dynamo regime resulting from the delay. This regime requires delays about an entire cycle and exhibits subequipartition magnetic activity. Revealing this new regime shows that even for long delays the flux tube/Babcock-Leighton dynamo can still deliver non-decaying solutions and remains a good candidate for a wide range of solar-like stars.
Context. Extrapolations of solar photospheric vector magnetograms into three-dimensional magnetic fields in the chromosphere and corona are usually done under the assumption that the fields are force-free. This condition is violated in the photosphere itself and a thin layer in the lower atmosphere above. The field calculations can be improved by preprocessing the photospheric magnetograms. The intention here is to remove a non-force-free component from the data.
Aims. We compare two preprocessing methods presently in use, namely the methods of Wiegelmann et al. (2006, Sol. Phys., 233, 215) and Fuhrmann et al. (2007, A&A, 476, 349).
Methods. The two preprocessing methods were applied to a vector magnetogram of the recently observed active region NOAA AR 10 953. We examine the changes in the magnetogram effected by the two preprocessing algorithms. Furthermore, the original magnetogram and the two preprocessed magnetograms were each used as input data for nonlinear force-free field extrapolations by means of two different methods, and we analyze the resulting fields.
Results. Both preprocessing methods managed to significantly decrease the magnetic forces and magnetic torques that act through the magnetogram area and that can cause incompatibilities with the assumption of force-freeness in the solution domain. The force and torque decrease is stronger for the Fuhrmann et al. method. Both methods also reduced the amount of small-scale irregularities in the observed photospheric field, which can sharply worsen the quality of the solutions. For the chosen parameter set, the Wiegelmann et al. method led to greater changes in strong-field areas, leaving weak-field areas mostly unchanged, and thus providing an approximation of the magnetic field vector in the chromosphere, while the Fuhrmann et al. method weakly changed the whole magnetogram, thereby better preserving patterns present in the original magnetogram. Both preprocessing methods raised the magnetic energy content of the extrapolated fields to values above the minimum energy, corresponding to the potential field. Also, the fields calculated from the preprocessed magnetograms fulfill the solenoidal condition better than those calculated without preprocessing.
The onset of a solar eruption is formulated here as either a magnetic catastrophe or as an instability. Both start with the same equation of force balance governing the underlying equilibria. Using a toroidal flux rope in an external bipolar or quadrupolar field as a model for the current-carrying flux, we demonstrate the occurrence of a fold catastrophe by loss of equilibrium for several representative evolutionary sequences in the stable domain of parameter space. We verify that this catastrophe and the torus instability occur at the same point; they are thus equivalent descriptions for the onset condition of solar eruptions.
We quantitatively address the conjecture that magnetic helicity must be shed from the Sun by eruptions launching coronal mass ejections in order to limit its accumulation in each hemisphere. By varying the ratio of guide and strapping field and the flux rope twist in a parametric simulation study of flux rope ejection from approximately marginally stable force-free equilibria, different ratios of self- and mutual helicity are set and the onset of the torus or helical kink instability is obtained. The helicity shed is found to vary over a broad range from a minor to a major part of the initial helicity, with self helicity being largely or completely shed and mutual helicity, which makes up the larger part of the initial helicity, being shed only partly. Torus-unstable configurations with subcritical twist and without a guide field shed up to about two-thirds of the initial helicity, while a highly twisted, kink-unstable configuration sheds only about one-quarter. The parametric study also yields stable force-free flux rope equilibria up to a total flux-normalized helicity of 0.25, with a ratio of self- to total helicity of 0.32 and a ratio of flux rope to external poloidal flux of 0.94. These results numerically demonstrate the conjecture of helicity shedding by coronal mass ejections and provide a first account of its parametric dependence. Both self- and mutual helicity are shed significantly; this reduces the total initial helicity by a fraction of ∼0.4--0.65 for typical source region parameters.
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.
Force-free equilibria containing two vertically arranged magnetic flux ropes of like chirality and current direction are considered as a model for split filaments/prominences and filament-sigmoid systems. Such equilibria are constructed analytically through an extension of the methods developed in Titov & Demoulin and numerically through an evolutionary sequence including shear flows, flux emergence, and flux cancellation in the photospheric boundary. It is demonstrated that the analytical equilibria are stable if an external toroidal (shear) field component exceeding a threshold value is included. If this component decreases sufficiently, then both flux ropes turn unstable for conditions typical of solar active regions, with the lower rope typically becoming unstable first. Either both flux ropes erupt upward, or only the upper rope erupts while the lower rope reconnects with the ambient flux low in the corona and is destroyed. However, for shear field strengths staying somewhat above the threshold value, the configuration also admits evolutions which lead to partial eruptions with only the upper flux rope becoming unstable and the lower one remaining in place. This can be triggered by a transfer of flux and current from the lower to the upper rope, as suggested by the observations of a split filament in Paper I. It can also result from tether-cutting reconnection with the ambient flux at the X-type structure between the flux ropes, which similarly influences their stability properties in opposite ways. This is demonstrated for the numerically constructed equilibrium.
We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the timescale of magnetic-field growth depends on the shock strength.