Refine
Has Fulltext
- no (61)
Year of publication
Document Type
- Article (61) (remove)
Is part of the Bibliography
- yes (61) (remove)
Keywords
- magnetohydrodynamics (MHD) (3)
- instabilities (2)
- Sun (1)
- Sun: atmosphere (1)
- Sun: corona (1)
- Sun: coronal mass ejections (CMEs) (1)
- Sun: flares (1)
- Sun: magnetic topology (1)
- corona (1)
- coronal mass ejections (CMEs) (1)
Institute
- Institut für Physik und Astronomie (61) (remove)
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.
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.
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.
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.
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 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.
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 numerically investigate nonlinear asymmetric square patterns in a horizontal convection layer with up-down reflection symmetry. As a novel feature we find the patterns to appear via the skewed varicose instability of rolls. The time-independent nonlinear state is generated by two unstable checkerboard (symmetric square) patterns and their nonlinear interaction. As the bouyancy forces increase the interacting modes give rise to bifurcations leading to a periodic alternation between a nonequilateral hexagonal pattern and the square pattern or to different kinds of standing oscillations.
Motivated by the successful Karlsruhe dynamo experiment, a relatively low-dimensional dynamo model is proposed. It is based on a strong truncation of the magnetohydrodynamic (MHD) equations with an external forcing of the Roberts type and the requirement that the model system satisfies the symmetries of the full MHD system, so that the first symmetry-breaking bifurcations can be captured. The backbone of the Roberts dynamo is formed by the Roberts flow, a helical mean magnetic field and another part of the magnetic field coupled to these two by triadic mode interactions. A minimum truncation model (MTM) containing only these energetically dominating primary mode triads is fully equivalent to the widely used first-order smoothing approximation. However, it is shown that this approach works only in the limit of small wave numbers of the excited magnetic field or small magnetic Reynolds numbers ($Rm ll 1$). To obtain dynamo action under more general conditions, secondary mode
We present a straightforward comparison of model calculations for the alpha-effect, helicities, and magnetic field line twist in the solar convection zone with magnetic field observations at atmospheric levels. The model calculations are carried out in a mixing-length approximation for the turbulence with a profile of the solar internal rotation rate obtained from helioseismic inversions. The magnetic field data consist of photospheric vector magnetograms of 422 active regions for which spatially-averaged values of the force-free twist parameter and of the current helicity density are calculated, which are then used to determine latitudinal profiles of these quantities. The comparison of the model calculations with the observations suggests that the observed twist and helicity are generated in the bulk of the convection zone, rather than in a layer close to the bottom. This supports two-layer dynamo models where the large-scale toroidal field is generated by differential rotation in a thin layer at the bottom while the alpha-effect is operating in the bulk of the convection zone. Our previous observational finding was that the moduli of the twist factor and of the current helicity density increase rather steeply from zero at the equator towards higher latitudes and attain a certain saturation at about 12 - 15 degrees. In our dynamo model with algebraic nonlinearity, the increase continues, however, to higher latitudes and is more gradual. This could be due to the neglect of the coupling between small-scale and large-scale current and magnetic helicities and of the latitudinal drift of the activity belts in the model
We study possible interrelations between the 300-year record of the yearly sunspot numbers and the solar inertial motion (SIM) using the recently developed technique of synchronization analysis. Phase synchronization of the sunspot cycle and the SIM is found and statistically confirmed in three epochs (1734-1790, 1855-1875 and 1907-1960) of the whole period 1700-2000. These results give quantitative support to the hypothesis that there is a weak interaction between the solar activity and the SIM.
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.