Thermal convection in Earth's inner core with phase change at its boundary

Inner core translation, with solidification on one hemisphere and melting on the other, provides a promising basis for understanding the hemispherical dichotomy of the inner core, as well as the anomalous stable layer observed at the base of the outer core - the F-layer - which might be sustained by continuous melting of inner core material. In this paper, we study in details the dynamics of inner core thermal convection when dynamically induced melting and freezing of the inner core boundary (ICB) are taken into account. If the inner core is unstably stratified, linear stability analysis and numerical simulations consistently show that the translation mode dominates only if the viscosity $\eta$ is large enough, with a critical viscosity value, of order $3 10^{18}$ Pas, depending on the ability of outer core convection to supply or remove the latent heat of melting or solidification. If $\eta$ is smaller, the dynamical effect of melting and freezing is small. Convection takes a more classical form, with a one-cell axisymmetric mode at the onset and chaotic plume convection at large Rayleigh number. [...] Thermal convection requires that a superadiabatic temperature profile is maintained in the inner core, which depends on a competition between extraction of the inner core internal heat by conduction and cooling at the ICB. Inner core thermal convection appears very likely with the low thermal conductivity value proposed by Stacey & Davis (2007), but nearly impossible with the much higher thermal conductivity recently put forward. We argue however that the formation of an iron-rich layer above the ICB may have a positive feedback on inner core convection: it implies that the inner core crystallized from an increasingly iron-rich liquid, resulting in an unstable compositional stratification which could drive inner core convection, perhaps even if the inner core is subadiabatic.Comment: 25 pages, 12 figure

Bound of dissipation on a plane Couette dynamo

Variational turbulence is among the few approaches providing rigorous results in turbulence. In addition, it addresses a question of direct practical interest, namely the rate of energy dissipation. Unfortunately, only an upper bound is obtained as a larger functional space than the space of solutions to the Navier-Stokes equations is searched. Yet, in general, this upper bound is in good agreement with experimental results in terms of order of magnitude and power law of the imposed Reynolds number. In this paper, the variational approach to turbulence is extended to the case of dynamo action and an upper bound is obtained for the global dissipation rate (viscous and Ohmic). A simple plane Couette flow is investigated. For low magnetic Prandtl number $P_m$ fluids, the upper bound of energy dissipation is that of classical turbulence (i.e. proportional to the cubic power of the shear velocity) for magnetic Reynolds numbers below $P_m^{-1}$ and follows a steeper evolution for magnetic Reynolds numbers above $P_m^{-1}$ (i.e. proportional to the shear velocity to the power four) in the case of electrically insulating walls. However, the effect of wall conductance is crucial : for a given value of wall conductance, there is a value for the magnetic Reynolds number above which energy dissipation cannot be bounded. This limiting magnetic Reynolds number is inversely proportional to the square root of the conductance of the wall. Implications in terms of energy dissipation in experimental and natural dynamos are discussed.Comment: In this new version, amistake (in equation 23 of the first version) is correcte

Remarks on compressible convection in Super-Earths

The radial density of planets increases with depth due to compressibility, leading to impacts on their convective dynamics. To account for these effects, including the presence of a quasi-adiabatic temperature profile and entropy sources due to dissipation, the compressibility is expressed through a dissipation number, $\mathcal{D}$, proportional to the planet's radius and gravity. In Earth's mantle, compressibility effects are moderate, but in large rocky or liquid exoplanets (Super-Earths), the dissipation number can become very large. This paper explores the properties of compressible convection when the dissipation number is significant. We start by selecting a simple Murnaghan equation of state that embodies the fundamental properties of condensed matter at planetary conditions. Next, we analyze the characteristics of adiabatic profiles and demonstrate that the ratio between the bottom and top adiabatic temperatures is relatively small and probably less than 2. We examine the marginal stability of compressible mantles and reveal that they can undergo convection with either positive or negative superadiabatic Rayleigh numbers. Lastly, we delve into simulations of convection performed using the exact equations of mechanics, neglecting inertia (infinite Prandtl number case), and examine their consequences for Super-Earths dynamics.Comment: 31 pages, 15 figure

Axisymmetric dynamo action produced by differential rotation, with anisotropic electrical conductivity and anisotropic magnetic permeability

The effect on dynamo action of an anisotropic electrical conductivity conjugated to an anisotropic magnetic permeability is considered. Not only is the dynamo fully axisymmetric, but it requires only a simple differential rotation, which twice challenges the well-established dynamo theory. Stability analysis is conducted entirely analytically, leading to an explicit expression of the dynamo threshold. The results show a competition between the anisotropy of electrical conductivity and that of magnetic permeability, the dynamo effect becoming impossible if the two anisotropies are identical. For isotropic electrical conductivity, Cowling's neutral point argument does imply the absence of an azimuthal component of current density, but does not prevent the dynamo effect as long as the magnetic permeability is anisotropic.Comment: 19 pages, 6 figure

A case for variational geomagnetic data assimilation: insights from a one-dimensional, nonlinear, and sparsely observed MHD system

Secular variations of the geomagnetic field have been measured with a continuously improving accuracy during the last few hundred years, culminating nowadays with satellite data. It is however well known that the dynamics of the magnetic field is linked to that of the velocity field in the core and any attempt to model secular variations will involve a coupled dynamical system for magnetic field and core velocity. Unfortunately, there is no direct observation of the velocity. Independently of the exact nature of the above-mentioned coupled system -- some version being currently under construction -- the question is debated in this paper whether good knowledge of the magnetic field can be translated into good knowledge of core dynamics. Furthermore, what will be the impact of the most recent and precise geomagnetic data on our knowledge of the geomagnetic field of the past and future? These questions are cast into the language of variational data assimilation, while the dynamical system considered in this paper consists in a set of two oversimplified one-dimensional equations for magnetic and velocity fields. This toy model retains important features inherited from the induction and Navier-Stokes equations: non-linear magnetic and momentum terms are present and its linear response to small disturbances contains Alfvén waves. It is concluded that variational data assimilation is indeed appropriate in principle, even though the velocity field remains hidden at all times; it allows us to recover the entire evolution of both fields from partial and irregularly distributed information on the magnetic field. This work constitutes a first step on the way toward the reassimilation of historical geomagnetic data and geomagnetic forecast

Geomagnetic Dipole Changes and Upwelling/Downwelling at the Top of the Earth's Core

The convective state of the top of Earth's outer core is still under debate. Conflicting evidence from seismology and geomagnetism provides arguments for and against a thick stably stratified layer below the core-mantle boundary. Mineral physics and cooling scenarios of the core favor a stratified layer. However, a non-zero secular variation of the total geomagnetic energy on the core-mantle boundary is evidence for the presence of radial motions extending to the top of the core. We compare the secular variation of the total geomagnetic energy with the secular variation of the geomagnetic dipole intensity and tilt. We demonstrate that both the level of cancellations of the sources and sinks of the dipole intensity secular variation, as well as the level of cancellations of the sources and sinks of the dipole tilt secular variation, are either larger than or comparable to the level of cancellations of the sources and sinks of the total geomagnetic energy secular variation on the core-mantle boundary, indicating that the latter is numerically significant hence upwelling/downwelling reach the top of the core. Radial motions below the core-mantle boundary are either evidence for no stratified layer or to its penetration by various dynamical mechanisms, most notably lateral heterogeneity of core-mantle boundary heat flux

Transition to Turbulence in the Hartmann Boundary Layer

The Hartmann boundary layer is a paradigm of magnetohydrodynamic (MHD) flows. Hartmann boundary layers develop when a liquid metal flows under the influence of a steady magnetic field. The present paper is an overview of recent successful attempts to understand the mechanisms by which the Hartmann layer undergoes a transition from laminar to turbulent flow. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/56018/1/125_ftp.pd