Refine
Has Fulltext
- no (31)
Year of publication
Document Type
- Article (31) (remove)
Language
- English (31) (remove)
Is part of the Bibliography
- yes (31)
Keywords
- flow (1)
- fluid shells (1)
- hopf-bifurcation (1)
- magnetoconvection (1)
- nonsymmetric linear-systems (1)
- number (1)
- onset (1)
- thermal-convection (1)
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.
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.
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 study Hamiltonian chaos generated by the dynamics of passive tracers moving in a two-dimensional fluid flow and describe the complex structure formed in a chaotic layer that separates a vortex region from the shear flow. The stable and unstable manifolds of unstable periodic orbits are computed. It is shown that their intersections in the Poincare map as an invariant set of homoclinic points constitute the backbone of the chaotic layer. Special attention is paid to the finite time properties of the chaotic layer. In particular, finite time Lyapunov exponents are computed and a scaling law of the variance of their distribution is derived. Additionally, the box counting dimension as an effective dimension to characterize the fractal properties of the layer is estimated for different duration times of simulation. Its behavior in the asymptotic time limit is discussed. By computing the Lyapunov exponents and by applying methods of symbolic dynamics, the formation of the layer as a function of the external forcing strength, which in turn represents the perturbation of the originally integrable system, is characterized. In particular, it is shown that the capture of KAM tori by the layer has a remarkable influence on the averaged Lyapunov exponents. (C) 2004 Elsevier Ltd. All rights reserved
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 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.
Symmetry-breaking bifurcations have been studied for convection in a nonrotating spherical shell whose outer radius is twice the inner radius, under the influence of an externally applied central force field with a radial dependence proportional to 1/r(5). This work is motivated by the GeoFlow experiment, which is performed under microgravity condition at the International Space Station where this particular central force can be generated. In order to predict the observable patterns, simulations together with path-following techniques and stability computations have been applied. Branches of axisymmetric, octahedral, and seven-cell solutions have been traced. The bifurcations producing them have been identified and their stability ranges determined. At higher Rayleigh numbers, time-periodic states with a complex spatiotemporal symmetry are found, which we call breathing patterns.
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.
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.
Bifurcations of dynamos in rotating and buoyancy-driven spherical Rayleigh-Benard convection in an electrically conducting fluid are investigated numerically. Both nonmagnetic and magnetic solution branches comprised of rotating waves are traced by path-following techniques, and their bifurcations and interconnections for different Ekman numbers are determined. In particular, the question of whether the dynamo branches bifurcate super- or sub-critically and whether a direct link to the primary pure convective states exists is answered.
Bifurcations in rotating spherical shell convection under the influence of differential rotation
(2021)
The bifurcations of thermal convection in a rotating spherical shell heated from the inner sphere and driven by the buoyancy of a central gravity field are studied numerically. This model of spherical Rayleigh-Benard convection describes large-scale convection in planets and in the outer zones of celestial bodies. In this work, the influence of an additionally imposed differential rotation of the inner sphere with respect to the outer one on the heat transfer and, more generally, on the whole bifurcation structure is investigated. In addition to numerical simulations, path-following techniques are applied in order to compute both stable and unstable solution branches. The dynamics and the heat transfer are essentially determined by a global bifurcation, which we have identified as a homoclinic bifurcation that consists of a collision of a stable modulated rotating with an unstable rotating wave.