J. Léorat

An interior penalty Galerkin method for the MHD equations in heterogeneous domains

The Maxwell equations in the magnetohydrodynamic (MHD) limit in heterogeneous domains composed of conducting and non-conducting regions are solved by using Lagrange finite elements and by enforcing continuities across interfaces using an Interior Penalty technique a la Baker [Finite element methods for elliptic equations using non-conforming elements, Math. Comp. 31 (137) (1977) 45-59]. The method is shown to be stable and convergent and is validated by convergence tests. It is used to compute Ohmic decay in various compact conducting domains and to simulate the kinematic dynamo action in two different geometries.

Nonlinear magnetohydrodynamics in axisymmetric heterogeneous domains using a Fourier/finite element technique and an interior penalty method

The Maxwell equations in the MHD limit in heterogeneous axisymmetric domains composed of conducting and non-conducting regions are solved by using a mixed Fourier/Lagrange finite element technique. Finite elements are used in the meridian plane and Fourier modes are used in the azimuthal direction. Parallelization is made with respect to the Fourier modes. Continuity conditions across interfaces are enforced using an interior penalty technique. The performance of the method is illustrated on kinematic and full dynamo configurations.

Nonlinear dynamo action in a precessing cylindrical container.

It is numerically demonstrated by means of a magnetohydrodynamics code that precession can trigger the dynamo effect in a cylindrical container. When the Reynolds number, based on the radius of the cylinder and its angular velocity, increases, the flow, which is initially centrosymmetric, loses its stability and bifurcates to a quasiperiodic motion. This unsteady and asymmetric flow is shown to be capable of sustaining dynamo action in the linear and nonlinear regimes. The magnetic field thus generated is unsteady and quadrupolar. These numerical evidences of dynamo action in a precessing cylindrical container may be useful for an experiment now planned at the Dresden sodium facility for dynamo and thermohydraulic studies in Germany.

Impact of Impellers on the Axisymmetric Magnetic Mode in the VKS2 Dynamo Experiment

In the von Kármán Sodium 2 (VKS2) successful dynamo experiment of September 2006, the observed magnetic field showed a strong axisymmetric component, implying that nonaxisymmetric components of the flow field were acting. By modeling the induction effect of the spiraling flow between the blades of the impellers in a kinematic dynamo code, we find that the axisymmetric magnetic mode is excited. The control parameters are the magnetic Reynolds number of the mean flow, the coefficient measuring the induction effect alpha, and the type of boundary conditions. We show that using realistic values of alpha, the observed critical magnetic Reynolds number, Rm;{c} approximately 32, can be reached easily with ferromagnetic boundary conditions. We conjecture that the dynamo action achieved in this experiment may not be related to the turbulence in the bulk of the flow, but rather to the alpha effect induced by the impellers.

A new Finite Element Method for magneto-dynamical problems: two-dimensional results

A mixed Lagrange finite element technique is used to solve the Maxwell equations in the magneto-hydrodynamic (MHD) limit in an hybrid domain composed of vacuum and conducting regions. The originality of the approach is that no artificial boundary condition is enforced at the interface between the conducting and the insulating regions and the non-conducting medium is not approximated by a weakly conducting medium as is frequently done in the literature. As a first evaluation of the performance of the method, we study two-dimensional (2D) configurations, where the flow streamlines of the conducting fluid are planar, i.e., invariant in one direction, and either the magnetic field (“magnetic scalar” case) or the electric field (“electric scalar” case) is parallel to the invariant direction. Induction heating, eddy current generation, and magnetic field stretching are investigated showing the usefulness of finite element methods to solve magneto-dynamical problems with complex insulating boundaries.

Electromagnetic induction in non-uniform domains

Kinematic simulations of the induction equation are carried out for different setups suitable for the von-Kármán-Sodium (VKS) dynamo experiment. The material properties of the flow driving impellers are modeled by means of high-conducting and high-permeability disks in a cylindrical volume filled with a conducting fluid. Two entirely different numerical codes are mutually validated by showing quantitative agreement on Ohmic decay and kinematic dynamo problems using various configurations and physical parameters. Field geometry and growth rates are strongly modified by the material properties of the disks even if the disks are thin. In contrast the influence of external boundary conditions remains small. Utilizing a VKS like mean fluid flow and high-permeability disks yield a reduction of the critical magnetic Reynolds number Rmc for the onset of dynamo action of the simplest non-axisymmetric field mode. However, this threshold reduction is not sufficient to fully explain the VKS experiment. We show that this reduction of Rmc is influenced by small variations in the flow configuration so that the observed reduction may be changed with respect to small modifications of setup and properties of turbulence.

Effects of discontinuous magnetic permeability on magnetodynamic problems

A novel approximation technique using Lagrange finite elements is proposed to solve magneto-dynamics problems involving discontinuous magnetic permeability and non-smooth interfaces. The algorithm is validated on benchmark problems and is used for kinematic studies of the Cadarache von Karman Sodium 2 (VKS2) experimental fluid dynamo.

Generation of axisymmetric modes in cylindrical kinematic mean-field dynamos of VKS type

In an attempt to understand why the dominating magnetic field observed in the von-Kármán-sodium (VKS) dynamo experiment is axisymmetric, we investigate in the present article the ability of mean field models to generate axisymmetric eigenmodes in cylindrical geometries. An α-effect is added to the induction equation and we identify reasonable and necessary properties of the α distribution so that axisymmetric eigenmodes are generated. The parametric study is done with two different simulation codes. We find that simple distributions of α-effect, either concentrated in the disk neighborhood or occupying the bulk of the flow, require unrealistically large values of the parameter α to explain the VKS observations.

Two spinning ways for precession dynamo.

It is numerically demonstrated by means of a magnetohydrodynamic code that precession can trigger dynamo action in a cylindrical container. Fixing the angle between the spin and the precession axis to be 1/2π, two limit configurations of the spinning axis are explored: either the symmetry axis of the cylinder is parallel to the spin axis (this configuration is henceforth referred to as the axial spin case), or it is perpendicular to the spin axis (this configuration is referred to as the equatorial spin case). In both cases, the centro-symmetry of the flow breaks when the kinetic Reynolds number increases. Equatorial spinning is found to be more efficient in breaking the centro-symmetry of the flow. In both cases, the average flow in the reference frame of the mantle converges to a counter-rotation with respect to the spin axis as the Reynolds number grows. We find a scaling law for the average kinetic energy in term of the Reynolds number in the axial spin case. In the equatorial spin case, the unsteady asymmetric flow is shown to be capable of sustaining dynamo action in the linear and nonlinear regimes. The magnetic field is mainly dipolar in the equatorial spin case, while it is is mainly quadrupolar in the axial spin case.

Nonlinear dynamo in a short Taylor–Couette setup

It is numerically demonstrated by means of a magnetohydrodynamics code that a short Taylor–Couette setup with a body force can sustain dynamo action. The magnetic threshold is comparable to what is usually obtained in spherical geometries. The linear dynamo is characterized by a rotating equatorial dipole. The nonlinear regime is characterized by fluctuating kinetic and magnetic energies and a tilted dipole whose axial component exhibits aperiodic reversals during the time evolution. These numerical evidences of dynamo action in a short Taylor–Couette setup may be useful for developing an experimental device.