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)
- corona (1)
- coronal mass ejections (1)
- coronal mass ejections (CMEs) (1)
- flares (1)
- flow (1)
- fluid shells (1)
- hopf-bifurcation (1)
- magnetoconvection (1)
- magnetohydrodynamics (1)
- nonsymmetric linear-systems (1)
- number (1)
- onset (1)
- sun (1)
- thermal-convection (1)
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.
We investigate the dynamo effect in a flow configuration introduced by G. O. Roberts in 1972. Based on a clear energetic hierarchy of Fourier components on the steady-state dynamo branch, an approximate model of interacting modes is constructed covering all essential features of the complete system but allowing simulations with a minimum amount of computation time. We use this model to study the excitation mechanism of the dynamo, the transition from stationary to time-dependent dynamo solutions and the characteristic properties of the latter ones.
We study a transition to hyperchaos in the two-dimensional incompressible Navier-Stokes equations with periodic boundary conditions and an external forcing term. Bifurcation diagrams are constructed by varying the Reynolds number, and a transition to hyperchaos (HC) is identified. Before the onset of HC, there is coexistence of two chaotic attractors and a hyperchaotic saddle. After the transition to HC, the two chaotic attractors merge with the hyperchaotic saddle, generating random switching between chaos and hyperchaos, which is responsible for intermittent bursts in the time series of energy and enstrophy. The chaotic mixing properties of the flow are characterized by detecting Lagrangian coherent structures. After the transition to HC, the flow displays complex Lagrangian patterns and an increase in the level of Lagrangian chaoticity during the bursty periods that can be predicted statistically by the hyperchaotic saddle prior to HC transition.
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
The multiplicity of stable convection patterns in a rotating spherical fluid shell heated from the inner boundary and driven by a central gravity field is presented. These solution branches that arise as rotating waves (RWs) are traced for varying Rayleigh number while their symmetry, stability, and bifurcations are studied. At increased Rayleigh numbers all the RWs undergo transitions to modulated rotating waves (MRWs) which are classified by their spatiotemporal symmetry. The generation of a third frequency for some of the MRWs is accompanied by a further loss of symmetry. Eventually a variety of MRWs, three-frequency solutions, and chaotic saddles and attractors control the dynamics for higher Rayleigh numbers.
It is shown that the ff effect of mean-field magnetohydrodynamics, which consists in the generation of a mean electromotive force along the mean magnetic field by turbulently fluctuating parts of velocity and magnetic field, is equivalent to the simultaneous generation of both turbulent and mean-field magnetic helicities, the generation rates being equal in magnitude and opposite in sign. In the particular case of statistically stationary and homogeneous fluctuations this implies that the ff effect can increase the energy in the mean magnetic field only under the condition that also magnetic helicity is accumulated there.