On the coupling of fluid dynamics and electromagnetism at the top of the earth's core

Abstract

A kinematic approach to short-term geomagnetism has recently been based upon pre-Maxwell frozen-flux electromagnetism. A complete dynamic theory requires coupling fluid dynamics to electromagnetism.

A geophysically plausible simplifying assumption for the vertical vorticity balance, namely that the vertical Lorentz torque is negligible, is introduced and its consequences are developed. The simplified coupled magnetohydrodynamic system is shown to conserve a variety of magnetic and vorticity flux integrals. These provide costraints on eligible models for the geomagnetic main field, its secular variation, and the horizontal fluid motions at the top of the core, and so permit a number of tests of the underlying assumptions.     

The effect of time-dependent γ-pumping on buoyant magnetic structures

ABSTRACT

In this paper, we explore for the first time the interactions of the net downward, time-dependent, γ-pumping overlying an imposed layer of magnetic fluid, in a polytropic atmosphere. Our calculations show that an equipartition of energy, between the magnetic and kinetic components, must be reached for buoyancy-driven magnetic structures to rise into the pumping region. However, structures do not rise unhindered, as in a previous investigation. We show that the evolution and other features of the emerging magnetic flux structures are significantly affected by the temporal variation of the γ-pumping. The rate of emerging structures, the strength of magnetic concentrations and the extent to how far magnetic field can travel were all found to depend on the timescale of the γ-pumping.     

The rapid dissipation of magnetic fields in highly conducting fluids

Abstract

This paper treats the dynamical conditions that obtain when long straight parallel twisted flux tubes in a highly conducting fluid are packed together in a broad array. It is shown that there is generally no hydrostatic equilibrium. In place of equilibrium there is a dynamical nonequilibrium, leading to neutral point reconnection and progressive coalescence of neighboring tubes (with the same sense of twisting), forming tubes of larger diameter and reduced twist. The magnetic energy in the twisting of each tube declines toward zero, dissipated into small-scale motions of the fluid and thence into heat.

The physical implications are numerous. For instance, it has been suggested that the subsurface magnetic field of the sun is composed of close-packed twisted flux tubes. Any such structures are short lived, at best.

The footpoints of the filamentary magnetic fields above bipolar magnetic regions on the sun are continually shuffled and rotated by the convection, so that the fields are composed of twisted rubes. The twisting and mutual wrapping is converted directly into fluid motion and heat by the dynamical nonequilibrium, so that the work done by the convection of the footpoints goes directly into heating the corona above. This theoretical result is the final step, then, in understanding the assertion by Rosner, Tucker, and Valana, and others, that the observed structure of the visible corona implies that it is heated principally by direct dissipation of the supporting magnetic field. It is the dynamical nonequilibrium that causes the dissipation, in spite of the high electrical conductivity. It would appear that any bipolar magnetic field extending upward from a dense convective layer into a tenuous atmosphere automatically produces heating, and a corona of some sort, in the sun or any other convective star.     

Tangential discontinuities and the optical analogy for stationary fields. VII. Effects of resistivity

Abstract

The relatively large resistivity in the solar photosphere and chromosphere softens the ideal tangential discontinuities of magnetostatic equilibrium into continuous transitions in field direction over scales of 0.1–10 km. This softening is communicated upward at the Alfvén speed into the active solar corona. The degree of softening is a vital part of the theory of magnetic heat input to the active X-ray corona, because the very low resistivity of the coronal gas provides effective dissipation only if the current sheets are reduced to a thickness of 10^?2km.

A close examination of the problem shows that the Alfvén transit time up into the corona is large compared to the characteristic time of 1 sec in which the coronal tangential discontinutities are formed. It also shows that the principal effect of the resistivity is to create a thin surface layer of fluid on adjacent flux bundles, which causes a general drift of the flux but does not directly broaden the current sheets higher up in the field. In fact the motions of the surface layers do not extend upwards beyond the first winding pattern at each end of a coronal loop.

It appears that the photospheric and chromospheric resistivity is without striking consequences for magnetic heating in the corona.     

Convection in a rotating cylinder with its sidewall having a finite thermal conductance

Abstract

An investigation is made of steady thermal convection of a Boussinesq fluid confined in a vertically-mounted rotating cylinder. The top and bottom endwall disks are thermal conductors at temperatures T_t and T_b with δT = T_t ? T_b >0. The vertical sidewall has a finite thermal conductance. A Newtonian heat flux condition is adopted at the sidewall. The Rayleigh number of the fluid system is large to render a boundary layer-type flow. Finite-difference numerical solutions to the full Navier-Stokes equations are obtained. The vertical motions within the buoyancy layer along the sidewall induce weak meridional flows in the interior. Because of the Coriolis acceleration, the meridional flows give rise to azimuthal flows relative to the rotating container. Strong vertical gradients of azimuthal flows exist in the regions near the endwalls. As the stratification effect increases, concentration of flow gradients in thin endwall boundary layers becomes more pronounced. The azimuthal flow field exhibits considerable horizontal gradients. The temperature field develops horizontal variations superposed on the dominant vertical distribution. As either the sidewall thermal conductance or the stratification effect decreases, the temperature distribution tends to the profile varying linearly with height. Comparisons of the sizes of the dynamic effects demonstrate that, in the bulk of flow field, the vertical shear of azimuthal velocity is supported by the horizontal temperature gradient, resulting in a thermal-wind relation.     

Topological semi-dynamos

Abstract

Magnetic field generation in a continuous medium in processes without self-excitation—the so-called semi-dynamo, involving as essential elements both magnetohydrodynamic processes and the presence of an impressed e.m.f.—has been studied for the case of the topological pumping effect on the magnetic field generation by an impressed e.m.f. operating in a three-dimensional Bénard convection layer.

Under conditions of interest for astrophysical applications the magnetic flux produced can exceed substantially that excited by the e.m.f. in the absence of motion.

The results obtained permitted an evaluation of the general quasi-steady magnetic field of the Sun generated by an azimuthal Coriolis e.m.f. which is active in the outermost layers of the convective envelope, taking into account small-scale convective and turbulent motions. In the polar regions of the Sun this field can reach ～10^?1 G.     

Magnetic buoyancy instabilities incorporating rotation

Abstract

Recent calculations suggest that the bulk of the solar toroidal field may be stored in a thin, convectively stable region situated between the convection zone proper and the radiative zone. Determining the stability properties of such a field is therefore important with implications for both the generation and escape of magnetic flux. The plane layer, linear stability analysis of Hughes (1985) is extended to incorporate the effects of uniform rotation. Detailed studies are made of interchange, or “axisymmetric” modes and of undular, or wavelike, motions, considering modes of both low and high frequency. The force due to rotation acts to constrain the fluid motions, a feature which is strongly stabilizing for direct modes, but can, in certain circumstances, be destabilizing for oscillatory motions.

For the interchange modes we show that the instability discussed at length by Hughes (1985), driven by fields increasing with height, is still present and indeed may be enhanced by rotational effects. We also study the more conventional instabilities, discussing the transformation between direct and oscillatory modes and considering in detail some peculiar properties of the oscillatory instabilities.

The more relevant instabilities in an astrophysical context are likely to be undular modes. Previous studies of low frequency modes driven by top heavy field gradients are extended to consider modes of various frequencies for a wide range of parameter values. Of particular interest is the occurrence of two distinct modes of instability for bottom heavy field gradients. We also exhibit some of the peculiar stability boundaries which can result when none of the competing influences in the problem is dominant.     

Magnetoconvection

Abstract

Cowling investigated the effect of an imposed magnetic field on convection in order to explain the origin of sunspots. After summarizing the classical linear theory of Boussinesq magnetoconvection, this review proceeds to more recent nonlinear results. Weakly nonlinear theory is used to establish the relevant bifurcation structure, which involves steady, oscillatory and chaotic solutions. Behaviour found in numerical experiments can then be related to these analytical results. Thereafter, attention is focused on the astrophysically relevant problem of fully compressible magnetoconvection. Steady two-dimensional nonlinear solutions show two important effects: stratification introduces an asymmetry between rising and falling fluid, while compressibility leads to evacuated magnetic flux sheets. Time-dependent behaviour includes transitions between standing waves and travelling waves, as well as changes in horizontal scale, leading to the development of more complicated spatial structures. Work on three-dimensional models, which is now in progress, will lead to a better understanding of the structure of a sunspot.     

Flows within the Ekman layer during spin-up of a thermally stratified fluid

Abstract

Finite-difference numerical solutions were obtained to present the flow and temperature field details within the transient Ekman layer during spin-up of a thermally stratified fluid in a cylinder. This complements the earlier studies on stratified spin-up which examined the flows in the interior core region. As the stratification increases, the following changes in the flow field are noticeable. The radial velocity in the Ekman layer decreases in magnitude. The azimuthal flows adjust smoothly from the interior region to the endwall boundary, and the Ekman layer in the azimuthal flow field fades. Vertical motions are inhibited, resulting in a weakened Ekman pumping. The axial vorticity field behaves similarly to the azimuthal flows. The temperature deviation from the equilibrium profile decreases, and the heat transfer flux from the endwall to the fluid decreases. The thickness of the thermal layer is larger than the velocity layer thickness. Illustrative comparisons of the relative sizes of the terms in the governing equations are conducted in order to assess the stratification effect in the adjustment process of the fluid.     

Resistive instabilities in rapidly rotating fluids: Linear theory of the convective modes

Abstract

This paper analyzes the linear stability of a rapidly-rotating, stratified sheet pinch in a gravitational field, g, perpendicular to the sheet. The sheet pinch is a layer (O ? z ? d) of inviscid, Boussinesq fluid of electrical conductivity σ, magnetic permeability μ, and almost uniform density ρ _o; z is height. The prevailing magnetic field. B _o(z), is horizontal at each z level, but varies in direction with z. The angular velocity, Ω, is vertical and large (Ω ? V_A/d, where V_A = B₀√(μρ₀) is the Alfvén velocity). The Elsasser number, Λ = σB² ₀/2Ωρ₀, measures σ. A (modified) Rayleigh number, R = gβd²/ρ₀V² _A, measures the buoyancy force, where β is the imposed density gradient, antiparallel to g. A Prandtl number, P_K = μσK, measures the diffusivity, k, of density differences.     

Flow of a stratified fluid bounded by walls of finite thermal conductance

Abstract

A study is made of the behavior of a thermally stratified fluid in a container when the non-horizontal boundaries have finite thermal conductance. The theory of Rahm and Walin is briefly recounted. Numerical solutions to the Navier-Stokes equations for a Boussinesq fluid in a cylinder, adopting a Newtonian heat flux condition at the vertical sidewall, are presented. Results on the details of flow and temperature fields are given over ranges of the Rayleigh number Ra, the container aspect ratio H, and the sidewall conductance S. As S increases, the isotherms in the meridional plane are horizontal at small radii but they diverge at large radii. This creates temperature nonuniformilies in the horizontal direction, and convective motions result. The salient features of the interior temperature profiles are captured by the theoretical model. The velocity field is characterized by two oppositely-directed circulations. As Ra or S varies, the qualitative circulation patterns remain substantially unchanged, but the magnitudes of the convective flows differ by large amounts. The effects of the externally-imposed parameters on the flow and temperature structures are examined.     

Plume capture by divergent plate motions: implications for the distribution of hotspots,geochemistry of mid-ocean ridge basalts,and estimates of the heat flux at the core–mantle boundary

The coexistence of stationary mantle plumes with plate-scale flow is problematic in geodynamics. We present results from laboratory experiments aimed at understanding the effects of an imposed large-scale circulation on thermal convection at high Rayleigh number (10⁶≤Ra≤10⁹) in a fluid with a temperature-dependent viscosity. In a large tank, a layer of corn syrup is heated from below while being stirred by large-scale flow due to the opposing motions of a pair of conveyor belts immersed in the syrup at the top of the tank. Three regimes are observed, depending on the ratio V of the imposed horizontal flow velocity to the rise velocity of plumes ascending from the hot boundary, and on the ratio λ of the viscosity of the interior fluid to the viscosity of the hottest fluid in contact with the bottom boundary. When V≪1 and λ≥1, large-scale circulation has a negligible effect on convection and the heat flux is due to the formation and rise of randomly spaced plumes. When V>10 and λ>100, plume formation is suppressed entirely, and the heat flux is carried by a sheet-like upwelling located in the center of the tank. At intermediate V, and depending on λ, established plume conduits are advected along the bottom boundary and ascending plumes are focused towards the central upwelling. Heat transfer across the layer occurs through a combination of ascending plumes and large-scale flow. Scaling analyses show that the bottom boundary layer thickness and, in turn, the basal heat flux q depend on the Peclet number, Pe, and λ. When λ>10, q∝Pe^1/2 and when λ→1, q∝(Peλ)^1/3, consistent with classical scalings. When applied to the Earth, our results suggest that plate-driven mantle flow focuses ascending plumes towards upwellings in the central Pacific and Africa as well as into mid-ocean ridges. Furthermore, plumes may be captured by strong upwelling flow beneath fast-spreading ridges. This behavior may explain why hotspots are more abundant near slow-spreading ridges than fast-spreading ridges and may also explain some observed variations of mid-ocean ridge basalt (MORB) geochemistry with spreading rate. Moreover, our results suggest that a potentially significant fraction of the core heat flux is due to plumes that are drawn into upwelling flows beneath ridges and not observed as hotspots.     

The onset of magnetoconvection at large Prandtl number in a rotating layer I. Finite magnetic diffusion

Abstract

This paper develops further a convection model that has been studied several times previously as a very crude idealization of planetary core dynamics. A plane layer of electrically-conducting fluid rotates about the vertical in the presence of a magnetic field. Such a field can be created spontaneously, as in the Childress—Soward dynamo, but here it is uniform, horizontal and externally-applied. The Prandtl number of the fluid is large, but the Ekman, Elsasser and Rayleigh numbers are of order unity, as is the ratio of thermal to magnetic diffusivity. Attention is focused on the onset of convection as the temperature difference applied across the layer is increased, and on the preferred mode, i.e., the planform and time-dependence of small amplitude convection. The case of main interest is the layer confined between electrically-insulating no-slip walls, but the analysis is guided by a parallel study based on illustrative boundary conditions that are mathematically simpler.     

Global aspects of solar wind–ionosphere interactions

Recent observations have quantified the auroral wind O⁺ outflow in response to magnetospheric inputs to the ionosphere, notably Poynting energy flux and precipitating electron density. For moderate to high activity periods, ionospheric O⁺ is observed to become a significant or dominant component of plasma pressure in the inner plasma sheet and ring current regions. Using a global circulation model of magnetospheric fields and its imposed ionospheric boundary conditions, we evaluate the global ionospheric plasma response to local magnetospheric conditions imposed by the simulation and evaluate magnetospheric circulation of solar wind H⁺, polar wind H⁺, and auroral wind O⁺. We launch and track the motions of millions of test particles in the global fields, launched at randomly distributed positions and times. Each particle is launched with a flux weighting and perpendicular and parallel energies randomly selected from defined thermal ranges appropriate to the launch point. One sequence is driven by a two-hour period of southward interplanetary magnetic field for average solar wind intensity. A second is driven by a 2-h period of enhanced solar wind dynamic pressure for average interplanetary field. We find that the simulated ionospheric O⁺ becomes a significant plasma pressure component in the inner plasma sheet and outer ring current region, particularly when the solar wind is intense or its magnetic field is southward directed. We infer that the reported empirical scalings of auroral wind O⁺ outflows are consistent with a substantial pressure contribution to the inner plasma sheet and plasma source surrounding the ring current. This result violates the common assumption that the ionospheric load is entirely confined to the F layer, and shows that the ionosphere is often an important dynamic element throughout the magnetosphere during moderate to large solar wind disturbances.     

Nonlinear convection in an imposed horizontal magnetic field

Abstract

Nonlinear two-dimensional magnetoconvection, with a Boussinesq fluid driven across the field-lines, is taken as a model for giant-cell convection in the sun and late-type stars. A series of numerical experiments shows the sensitivity of the horizontal scale of convection to the applied field and to the Rayleigh number R. Overstable oscillations occur in cells as broad as they are deep, but increasing R leads to steady motions of much greater wavelength. Purely geometrical effects can cause oscillation: this work implies that strong horizontal field will in general lead to time-dependent convection.     

The onset of magnetoconvection at large Prandtl number in a rotating layer II. Small magnetic diffusion

Abstract

This paper develops further a convection model that has been studied several times previously as a very crude idealization of planetary core dynamics. A plane layer of electrically-conducting fluid rotates about the vertical in the presence of a magnetic field. Such a field can be created spontaneously, as in the Childress-Soward dynamo, but here it is uniform, horizontal and externally-applied. The Prandtl number of the fluid is large, but the Ekman, Elsasser and Rayleigh numbers are of unit order. In Part I of this series, it was also supposed that the ratio thermal diffusivity diffusivity/magnetic diffusivity is O(1), but here we suppose that this ratio is large. The character of the solution is changed in this limit. In the case of main interest, when the layer is confined between electrically-insulating no-slip walls, the solution is significantly different from the solution when the mathematically simpler, illustrative boundary conditions also considered in Part I are employed. As in Part I, attention is focussed on the onset of convection as the temperature difference applied across the layer is increased, and on the preferred mode, i.e., the planform and time-dependence of small amplitude convection.     

Intermittent convection

Abstract

A theoretical analysis of pseudo two-dimensional, finite-amplitude, thermal convection is made for an infinite Prandtl number fluid which is subjected to a constant heat flux out of the top boundary and insulated at the bottom. For large Rayleigh numbers the convective flow becomes intermittent and the system is characterized by the following cyclic process: the formation of a thermal boundary layer by diffusion, the instability of this layer when it becomes sufficiently thick, the destruction of the layer by the convective flow, the dying down of the convection, and the reforming of the thermal boundary layer by diffusion. The periodicity and the horizontal wave number of the intermittent convective flow are found to be independent of the depth of the fluid layer but depend on the rate of cooling and the properties of the fluid.     

Thermal and magnetically driven instabilities in a non-constantly stratified fluid layer

Abstract

The linear hydromagnetic stability of a non-constantly stratified horizontal fluid layer permeated by an azimuthal non-homogeneous magnetic field is investigated for various widths of the stably stratified part of the layer in the geophysical limit q→0 (q is the ratio of thermal and magnetic diffusivities). The choice of the strength of the magnetic field B_o is as in Soward (1979) (see also Soward and Skinner, 1988) and the equations for the disturbances are treated as in Fearn and Proctor (1983). It was found that convection is developed in the whole layer regardless of the width of its stably stratified part. The thermal instability penetrates essentially from the unstably stratified part of the layer into the stably stratified part for A ～ 1 (A characterises the ratio of the Lorentz and Coriolis forces). When the magnetic field is strong (A>1) the thermal convection is suppressed in the stably stratified part of the layer. However, in this case, it is replaced by the magnetically driven instability; which is fully developed in the whole layer. The thermal instabilities always propagate westward and exist for all the modes m. The magnetically driven instabilities propagate either westward or eastward according to the width of the stably and unstably stratified parts and exist only for the mode m=1.     

The destabilising nature of differential rotation

Abstract

In a rapidly rotating, electrically conducting fluid we investigate the thermal stability of the fluid in the presence of an imposed toroidal magnetic field and an imposed toroidal differential rotation. We choose a magnetic field profile that is stable. The familiar role of differential rotation is a stabilising one. We wish to examine the less well known destabilising effect that it can have. In a plane layer model (for which we are restricted to Roberts number q = 0) with differential rotation, U = sΩ(z)1 _?, no choice of Ω(z) led to a destabilising effect. However, in a cylindrical geometry (for which our model permits all values of q) we found that differential rotations U = sΩ(s)1 _? which include a substantial proportion of negative gradient (dΩ/ds ≤ 0) give a destabilising effect which is largest when the magnetic Reynolds number R _m = O(10); the critical Rayleigh number, Ra _c, is about 7% smaller at minimum than at R_m = 0 for q = 10⁶. We also find that as q is reduced, the destabilising effect is diminished and at q = 10^?6, which may be more appropriate to the Earth's core, the effect causes a dip in the critical Rayleigh number of only about 0.001%. This suggests that we see no dip in the plane layer results because of the q = 0 condition. In the above results, the Elsasser number A = 1 but the effect of differential rotation is also dependent on A. Earlier work has shown a smooth transition from thermal to differential rotation driven instability at high A [A = O(100)]. We find, at intermediate A [A = O(10)], a dip in the Ra_c vs. R_m curve similar to the A = 1 case. However, it has Ra_c ≤ 0 at its minimum and unlike the results for high A, larger values of R_m result in a restabilisation.