Refine
Year of publication
Document Type
- Article (62)
- Preprint (17)
- Monograph/Edited Volume (9)
- Doctoral Thesis (1)
- Postprint (1)
Keywords
- magnetohydrodynamics (MHD) (3)
- instabilities (2)
- magnetic fields (2)
- MHD (1)
- Sun (1)
- Sun: atmosphere (1)
- Sun: corona (1)
- Sun: coronal mass ejections (CMEs) (1)
- Sun: flares (1)
- Sun: magnetic topology (1)
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.
In this study, we model a sequence of a confined and a full eruption, employing the relaxed end state of the confined eruption of a kink-unstable flux rope as the initial condition for the ejective one. The full eruption, a model of a coronal mass ejection, develops as a result of converging motions imposed at the photospheric boundary, which drive flux cancellation. In this process, parts of the positive and negative external flux converge toward the polarity inversion line, reconnect, and cancel each other. Flux of the same amount as the canceled flux transfers to a flux rope, increasing the free magnetic energy of the coronal field. With sustained flux cancellation and the associated progressive weakening of the magnetic tension of the overlying flux, we find that a flux reduction of approximate to 11% initiates the torus instability of the flux rope, which leads to a full eruption. These results demonstrate that a homologous full eruption, following a confined one, can be driven by flux cancellation.
We have numerically studied the bifurcation properties of a sheet pinch with impenetrable stress-free boundaries. An incompressible, electrically conducting fluid with spatially and temporally uniform kinematic viscosity and magnetic diffusivity is confined between planes at x1=0 and 1. Periodic boundary conditions are assumed in the x2 and x3 directions and the magnetofluid is driven by an electric field in the x3 direction, prescribed on the boundary planes. There is a stationary basic state with the fluid at rest and a uniform current J=(0,0,J3). Surprisingly, this basic state proves to be stable and apparently to be the only time-asymptotic state, no matter how strong the applied electric field and irrespective of the other control parameters of the system, namely, the magnetic Prandtl number, the spatial periods L2 and L3 in the x2 and x3 directions, and the mean values B¯2 and B¯3 of the magnetic-field components in these directions.
Context: Most solar and stellar dynamo models use the alpha-Omega scenario where the magnetic field is generated by the interplay between differential rotation (the Omega effect) and a mean electromotive force due to helical turbulent convection flows (the alpha effect). There are, however, turbulent dynamo mechnisms that may complement the alpha effect or may be an alternative to it. Aims: We investigate models of solar-type dynamos where the alpha effect is completely replaced by two other turbulent dynamo mechanisms, namely the Omega x J effect and the shear- current effect, which both result from an inhomogeneity of the mean magnetic field. Methods: We studied axisymmetric mean-field dynamo models containing differential rotation, the Omega x J and shear-current effects, and a meridional circulation. The model calculations were carried out using the rotation profile of the Sun as obtained from helioseismic measurements and radial profiles of other quantities according to a standard model of the solar interior. Results: Without meridional flow, no satisfactory agreement of the models with the solar observations can be obtained. With a sufficiently strong meridional circulation included, however, the main properties of the large-scale solar magnetic field, namely, its oscillatory behavior, its latitudinal drift towards the equator within each half cycle, and its dipolar parity with respect to the equatorial plane, are correctly reproduced. Conclusions: We have thereby constructed the first mean-field models of solar-type dynamos that do not use the alpha effect.
The equilibrium states of electrically conducting fluids or plasmas have been a subject of intense study for a long time, motivated in particular by the interest in controlled thermonuclear fusion, as well as that in space and astrophysical phenomena such as plasma loops in the solar corona. If high temperatures prohibit solid walls, a conducting fluid can be held together by the action of an electric current passing through it with the pressure gradients being balanced by the Lorentz force. The resultant configuration is known as a pinch. In this paper we report on studies of the pinch in the geometry of a plane sheet.
A numerical bifurcation analysis of the electrically driven plane sheet pinch is presented. The electrical conductivity varies across the sheet such as to allow instability of the quiescent basic state at some critical Hartmann number. The most unstable perturbation is the two-dimensional tearing mode. Restricting the whole problem to two spatial dimensions, this mode is followed up to a time-asymptotic steady state, which proves to be sensitive to three- dimensional perturbations even close to the point where the primary instability sets in. A comprehensive three- dimensional stability analysis of the two-dimensional steady tearing-mode state is performed by varying parameters of the sheet pinch. The instability with respect to three-dimensional perturbations is suppressed by a sufficiently strong magnetic field in the invariant direction of the equilibrium. For a special choice of the system parameters, the unstably perturbed state is followed up in its nonlinear evolution and is found to approach a three-dimensional steady state.
A numerical bifurcation analysis of the electrically driven plane sheet pinch is presented. The electrical conductivity varies across the sheet such as to allow instability of the quiescent basic state at some critical Hartmann number. The most unstable perturbation is the two-dimensional tearing mode. Restricting the whole problem to two spatial dimensions, this mode is followed up to a time-asymptotic steady state, which proves to be sensitive to three-dimensional perturbations even close to the point where the primary instability sets in. A comprehensive three-dimensional stability analysis of the two-dimensional steady tearing-mode state is performed by varying parameters of the sheet pinch. The instability with respect to three-dimensional perturbations is suppressed by a sufficiently strong magnetic field in the invariant direction of the equilibrium. For a special choice of the system parameters, the unstably perturbed state is followed up in its nonlinear evolution and is found to approach a three-dimensional steady state.
Electrically conducting fluids in motion can act as self-excited dynamos. The magnetic fields of celestial bodies like the Earth and the Sun are generated by such dynamos. Their theory aims at modeling and understanding both the kinematic and dynamic aspects of the underlying processes. Kinematic dynamo models, in which for a prescribed flow the linear induction equation is solved and growth rates of the magnetic field are calculated, have been studied for many decades. But in order to get consistent models and to take into account the back-reaction of the magnetic field on the fluid motion, the full nonlinear system of the magnetohydrodynamic (MHD) equations has to be studied. It is generally accepted that these equations, i.e. the Navier-Stokes equation (NSE) and the induction equation, provide a theoretical basis for the explanation of the dynamo effect. The general idea is that mechanical energy pumped into the fluid by heating or other mechanisms is transferred to the magnetic field by nonlinear interactions. For two special helical flows which are known to be effective kinematic dynamos and which can be produced by appropriate external mechanical forcing, we review the nonlinear dynamo properties found in the framework of the full MHD equations. Specifically, we deal with the ABC flow (named after Arnold, Beltrami and Childress) and the Roberts flow (after G.~O. Roberts). The appearance of generic dynamo effects is demonstrated. Applying special numerical bifurcation-analysis techniques to high-dimensional approximations in Fourier space and varying the Reynolds number (or the strength of the forcing) as the relevant control parameter, qualitative changes in the dynamics are investigated. We follow the bifurcation sequences until chaotic states are reached. The transitions from the primary flows with vanishing magnetic field to dynamo-active states are described in particular detail. In these processes the stagnation points of the flows and their heteroclinic connections play a promoting role for the magnetic field generation. By the example of the Roberts flow we demonstrate how the break up of the heteroclinic lines after the primary bifurcation leads to a complicated intersection of stable and unstable manifolds forming a chaotic web which is in turn correlated with the spatial appearance of the dynamo.
We report on bifurcation studies for the incompressible magnetohydrodynamic equations in three space dimensions with periodic boundary conditions and a temporally constant external forcing. Fourier reprsentations of velocity, pressure and magnetic field have been used to transform the original partial differential equations into systems of ordinary differential equations (ODE), to which then special numerical methods for the qualitative analysis of systems of ODE have been applied, supplemented by the simulative calculation of solutions for selected initial conditions. In a part of the calculations, in order to reduce the number of modes to be retained, the concept of approximate inertial manifolds has been applied. For varying (incereasing from zero) strength of the imposed forcing, or varying Reynolds number, respectively, time-asymptotic states, notably stable stationary solutions, have been traced. A primary non-magnetic steady state loses, in a Hopf bifurcation, stability to a periodic state with a non-vanishing magnetic field, showing the appearance of a generic dynamo effect. From now on the magnetic field is present for all values of the forcing. The Hopf bifurcation is followed by furhter, symmetry-breaking, bifurcations, leading finally to chaos. We pay particular attention to kinetic and magnetic helicities. The dynamo effect is observed only if the forcing is chosen such that a mean kinetic helicity is generated; otherwise the magnetic field diffuses away, and the time-asymptotic states are non-magnetic, in accordance with traditional kinematic dynamo theory.
We report on bifurcation studies for the incompressible magnetohydrodynamic equations in three space dimensions with periodic boundary conditions and a temporally constant external forcing. Fourier representations of velocity, pressure and magnetic field have been used to transform the original partial differential equations into systems of ordinary differential equations (ODE), to which then special numerical methods for the qualitative analysis of systems of ODE have been applied, supplemented by the simulative calculation of solutions for selected initial conditions. In a part of the calculations, in order to reduce the number of modes to be retained, the concept of approximate inertial manifolds has been applied. For varying (increasing from zero) strength of the imposed forcing, or varying Reynolds number, respectively, time-asymptotic states, notably stable stationary solutions, have been traced. A primary non- magnetic steady state loses, in a Hopf bifurcation, stability to a periodic state with a non-vanishing magnetic field, showing the appearance of a generic dynamo effect. From now on the magnetic field is present for all values of the forcing. The Hopf bifurcation is followed by further, symmetry-breaking, bifurcations, leading finally to chaos. We pay particular attention to kinetic and magnetic helicities. The dynamo effect is observed only if the forcing is chosen such that a mean kinetic helicity is generated; otherwise the magnetic field diffuses away, and the time-asymptotic states are non-magnetic, in accordance with traditional kinematic dynamo theory.
Three-dimensional bouyancy-driven convection in a horizontal fluid layer with stress-free boundary conditions at the top and bottom and periodic boundary conditions in the horizontal directions is investigated by means of numerical simulation and bifurcation-analysis techniques. The aspect ratio is fixed to a value of 2√2 and the Prandtl number to a value of 6.8. Two-dimensional convection rolls are found to be stable up to a Rayleigh number of 17 950, where a Hopf bifurcation leads to traveling waves. These are stable up to a Rayleigh number of 30 000, where a secondary Hopf bifurcation generates modulated traveling waves. We pay particular attention to the symmetries of the solutions and symmetry breaking by the bifurcations.
We have studied the bifurcation structure of the incompressible two-dimensional Navier-Stokes equations with a special external forcing driving an array of 8×8 counterrotating vortices. The study has been motivated by recent experiments with thin layers of electrolytes showing, among other things, the formation of large-scale spatial patterns. As the strength of the forcing or the Reynolds number is raised the original stationary vortex array becomes unstable and a complex sequence of bifurcations is observed. The bifurcations lead to several periodic branches, torus and chaotic solutions, and other stationary solutions. Most remarkable is the appearance of solutions characterized by structures on spatial scales large compared to the scale of the forcing. We also characterize the different dynamic regimes by means of tracers injected into the fluid. Stretching rates and Hausdorff dimensions of convected line elements are calculated to quantify the mixing process. It turns out that for time-periodic velocity fields the mixing can be very effective.
We report on bifurcation studies for the incompressible Navier-Stokes equations in two space dimensions with periodic boundary conditions and an external forcing of the Kolmogorov type. Fourier representations of velocity and pressure have been used to approximate the original partial differential equations by a finite-dimensional system of ordinary differential equations, which then has been studied by means of bifurcation-analysis techniques. A special route into chaos observed for increasing Reynolds number or strength of the imposed forcing is described. It includes several steady states, traveling waves, modulated traveling waves, periodic and torus solutions, as well as a period-doubling cascade for a torus solution. Lyapunov exponents and Kaplan-Yorke dimensions have been calculated to characterize the chaotic branch. While studying the dynamics of the system in Fourier space, we also have transformed solutions to real space and examined the relation between the different bifurcations in Fourier space and toplogical changes of the streamline portrait. In particular, the time-dependent solutions, such as, e.g., traveling waves, torus, and chaotic solutions, have been characterized by the associated fluid-particle motion (Lagrangian dynamics).
Two-dimensional bouyancy-driven convection in a horizontal fluid layer with stress-free boundary conditions at top and bottom and periodic boundary conditions in the horizontal direction is investigated by means of numerical simulation and bifurcation-analysis techniques. As the bouyancy forces increase, the primary stationary and symmetric convection rolls undergo successive Hopf bifurcations, bifurcations to traveling waves, and phase lockings. We pay attention to symmetry breaking and its connection with the generation of large-scale horizontal flows. Calculations of Lyapunov exponents indicate that at a Rayleigh number of 2.3×105 no temporal chaos is reached yet, but the system moves nonchaotically on a 4-torus in phase space.
The dynamics and bifurcations of convective waves in rotating and buoyancy-driven spherical Rayleigh-Benard convection are investigated numerically. The solution branches that arise as rotating waves (RWs) are traced by means of path-following methods, by varying the Rayleigh number as a control parameter for different rotation rates. The dependence of the azimuthal drift frequency of the RWs on the Ekman and Rayleigh numbers is determined and discussed. The influence of the rotation rate on the generation and stability of secondary branches is demonstrated. Multistability is typical in the parameter range considered.
The dynamics and bifurcations of convective waves in rotating and buoyancy-driven spherical Rayleigh-Benard convection are investigated numerically. The solution branches that arise as rotating waves (RWs) are traced by means of path-following methods, by varying the Rayleigh number as a control parameter for different rotation rates. The dependence of the azimuthal drift frequency of the RWs on the Ekman and Rayleigh numbers is determined and discussed. The influence of the rotation rate on the generation and stability of secondary branches is demonstrated. Multistability is typical in the parameter range considered.
The bifurcations in a three-dimensional incompressible, electrically conducting fluid with an external forcing of the Roberts type have been studied numerically. The corresponding flow can serve as a model for the convection in the outer core of the Earth and is realized in an ongoing laboratory experiment aimed at demonstrating a dynamo effect. The symmetry group of the problem has been determined and special attention has been paid to symmetry breaking by the bifurcations. The nonmagnetic, steady Roberts flow loses stability to a steady magnetic state, which in turn is subject to secondary bifurcations. The secondary solution branches have been traced until they end up in chaotic states.
We study the dynamo properties of asymmetric square patterns in Boussinesq Rayleigh-B'enard convection in a plane horizontal layer. Cases without rotation and with weak rotation about a vertical axis are considered. There exist different types of solutions distinguished by their symmetry, among them such with flows possessing a net helicity and being capable of kinematic dynamo action in the presence as well as in the absence of rotation. In the nonrotating case these flows are, however, always only kinematic, not nonlinear dynamos. Nonlinearly the back-reaction of the magnetic field then forces the solution into the basin of attraction of a roll pattern incapable of dynamo action. But with rotation added parameter regions are found where the Coriolis force counteracts the Lorentz force in such a way that the asymmetric squares are also nonlinear dynamos.
The bifurcation behaviour of the 3D magnetohydrodynamic equations has been studied for external forcings of varying degree of helicity. With increasing strength of the forcing a primary non-magnetic steady state loses stability to a magnetic periodic state if the helicity exceeds a threshold value and to different non-magnetic states otherwise.
Using the incompressible magnetohydrodynamic equations, we have numerically studied the dynamo effect in electrically conducting fluids. The necessary energy input into the system was modeled either by an explicit forcing term in the Navier-Stokes equation or fully selfconsistently by thermal convection in a fluid layer heated from below. If the fluid motion is capable of dynamo action, the dynamo effect appears in the form of a phase transition or bifurcation at some critical strength of the forcing. Both the dynamo bifurcation and subsequent bifurcations that occur when the strength of the forcing is further raised were studied, including the transition to chaotic states. Special attention was paid to the helicity of the flow as well as to the symmetries of the system and symmetry breaking in the bifurcations. The magnetic field tends to be accumulated in special regions of the flow, notably in the vicinity of stagnation points or near the boundaries of convection cells.
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.
It has been suggested that coronal mass ejections (CMEs) remove the magnetic he-licity of their coronal source region from the Sun. Such removal is often regarded to be necessary due to the hemispheric sign preference of the helicity, which inhibits a simple annihilation by reconnection between volumes of opposite chirality. Here we monitor the relative magnetic he-licity contained in the coronal volume of a simulated flux rope CME, as well as the upward flux of relative helicity through horizontal planes in the simulation box. The unstable and erupting flux rope carries away only a minor part of the initial relative helicity; the major part remains in the volume. This is a consequence of the requirement that the current through an expanding loop must decrease if the magnetic energy of the configuration is to decrease as the loop rises, to provide the kinetic energy of the CME.
We investigate numerically the appearance of heteroclinic behavior in a three-dimensional, buoyancy-driven fluid layer with stress-free top and bottom boundaries, a square horizontal periodicity with a small aspect ratio, and rotation at low to moderate rates about a vertical axis. The Prandtl number is 6.8. If the rotation is not too slow, the skewed-varicose instability leads from stationary rolls to a stationary mixed-mode solution, which in turn loses stability to a heteroclinic cycle formed by unstable roll states and connections between them. The unstable eigenvectors of these roll states are also of the skewed-varicose or mixed-mode type and in some parameter regions skewed-varicose like shearing oscillations as well as square patterns are involved in the cycle. Always present weak noise leads to irregular horizontal translations of the convection pattern and makes the dynamics chaotic, which is verified by calculating Lyapunov exponents. In the nonrotating case the primary rolls lose, depending on the aspect ratio, stability to traveling waves or a stationary square pattern. We also study the symmetries of the solutions at the intermittent fixed points in the heteroclinic cycle.
We investigate numerically the appearance of heteroclinic behavior in a three-dimensional, buoyancy-driven fluid layer with stress-free top and bottom boundaries, a square horizontal periodicity with a small aspect ratio, and rotation at low to moderate rates about a vertical axis. The Prandtl number is 6.8. If the rotation is not too slow, the skewed-varicose instability leads from stationary rolls to a stationary mixed-mode solution, which in turn loses stability to a heteroclinic cycle formed by unstable roll states and connections between them. The unstable eigenvectors of these roll states are also of the skewed-varicose or mixed-mode type and in some parameter regions skewed-varicose like shearing oscillations as well as square patterns are involved in the cycle. Always present weak noise leads to irregular horizontal translations of the convection pattern and makes the dynamics chaotic, which is verified by calculating Lyapunov exponents. In the nonrotating case, the primary rolls lose, depending on the aspect ratio, stability to traveling waves or a stationary square pattern. We also study the symmetries of the solutions at the intermittent fixed points in the heteroclinic cycle.