首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Abstract

We investigate the influence of differential rotation on magnetic instabilities for an electrically conducting fluid in the presence of a toroidal basic state of magnetic field B 0 = BMB0(r, θ)1 φ and flow U0 = UMU0 (r, θ)1φ, [(r, θ, φ) are spherical polar coordinates]. The fluid is confined in a rapidly rotating, electrically insulating, rigid spherical container. In the first instance the influence of differential rotation on established magnetic instabilities is studied. These can belong to either the ideal or the resistive class, both of which have been the subject of extensive research in parts I and II of this series. It was found there, that in the absence of differential rotation, ideal modes (driven by gradients of B 0) become unstable for Ac ? 200 whereas resistive instabilities (generated by magnetic reconnection processes near critical levels, i.e. zeros of B0) require Ac ? 50. Here, Λ is the Elsasser number, a measure of the magnetic field strength and Λc is its critical value at marginal stability. Both types of instability can be stabilised by adding differential rotation into the system. For the resistive modes the exact form of the differential rotation is not important whereas for the ideal modes only a rotation rate which increases outward from the rotation axis has a stabilising effect. We found that in all cases which we investigated Λc increased rapidly and the modes disappeared when Rm ≈ O(ΛC), where the magnetic Reynolds number Rm is a measure of the strength of differential rotation. The main emphasis, however, is on instabilities which are driven by unstable gradients of the differential rotation itself, i.e. an otherwise stable fluid system is destabilised by a suitable differential rotation once the magnetic Reynolds number exceeds a certain critical value (Rm )c. Earlier work in the cylindrical geometry has shown that the differential rotation can generate an instability if Rm ) ?O(Λ). Those results, obtained for a fixed value of Λ = 100 are extended in two ways: to a spherical geometry and to an analysis of the effect of the magnetic field strength Λ on these modes of instability. Calculations confirm that modes driven by unstable gradients of the differential rotation can exist in a sphere and they are in good agreement with the local analysis and the predictions inferred from the cylindrical geometry. For Λ = O(100), the critical value of the magnetic Reynolds number (Rm )c Λ 100, depending on the choice of flow U0 . Modes corresponding to azimuthal wavenumber m = 1 are the most unstable ones. Although the magnetic field B 0 is itself a stable one, the field strength plays an important role for this instability. For all modes investigated, both for cylindrical and spherical geometries, (Rm )c reaches a minimum value for 50 ≈ Λ ≈ 100. If Λ is increased, (Rm )c ∝ Λ, whereas a decrease of Λ leads to a rapid increase of (Rm )c, i.e. a stabilisation of the system. No instability was found for Λ ≈ 10 — 30. Optimum conditions for instability driven by unstable gradients of the differential rotation are therefore achieved for ≈ Λ 50 — 100, Rm ? 100. These values lead to the conclusion that the instabilities can play an important role in the dynamics of the Earth's core.  相似文献   

2.
Linear magnetoconvection in a model of a non-uniformly stratified horizontal rotating fluid layer with a toroidal magnetic field is investigated for no-slip and finitely electrically conductive boundaries and with very thin stably stratified upper sublayer. The basic parabolic temperature profile is determined by the temperature difference between the boundaries and by the homogeneous heat source distribution in the layer. This results in a density pattern, in which a stably stratified upper sublayer is present. The developed diffusive perturbations (modes) are strongly affected by the complicated coupling of viscous, thermal and magnetic diffusive processes. The calculations were performed for various values of Roberts number (q ≪ 1 and q = O(1)). The mean electromotive force produced by the developed hydromagnetic instabilities is investigated to find the modes, which can be appropriate for creating the α-effect. It was found that the azimuthal part of the EMF is dominant for westward modes when the Elsasser number Λ ≲ O(1).  相似文献   

3.
Abstract

This paper presents the first attempt to examine the stability of a poloidal magnetic field in a rapidly rotating spherical shell of electrically conducting fluid. We find that a steady axisymmetric poloidal magnetic field loses its stability to a non-axisymmetric perturbation when the Elsasser number A based on the maximum strength of the field exceeds a value about 20. Comparing this with observed fields, we find that, for any reasonable estimates of the appropriate parameters in planetary interiors, our theory predicts that all planetary poloidal fields are stable, with the possible exception of Jupiter. The present study therefore provides strong support for the physical relevance of magnetic stability analysis to planetary dynamos. We find that the fluid motions driven by magnetic instabilities are characterized by a nearly two-dimensional columnar structure attempting to satisfy the Proudman-Taylor theorm. This suggests that the most rapidly growing perturbation arranges itself in such a way that the geostrophic condition is satisfied to leading order. A particularly interesting feature is that, for the most unstable mode, contours of the non-axisymmetric azimuthal flow are closely aligned with the basic axisymmetric poloidal magnetic field lines. As a result, the amplitude of the azimuthal component of the instability is smaller than or comparable with that of the poloidal component, in contrast with the instabilities generated by toroidal decay modes (Zhang and Fearn, 1994). It is shown, by examining the same system with and without fluid inertia, that fluid inertia plays a secondary role when the magnetic Taylor number Tm ? 105. We find that the direction of propagation of hydromagnetic waves driven by the instability is influenced strongly by the size of the inner core.  相似文献   

4.
The equations for an electrically conducting fluid in cylindrical coordinates are linearized assuming that the inertial terms in the momentum equation can be ignored (small Rossby number), and that the ratio of the Elsasser number and magnetic Reynolds number is one. After these assumptions, the governing equations are linearized about an ambient solution which vanishes at the the equator. Upon assuming large Elsasser and magnetic Reynolds number, the solutions to the linearized equations are approximated by wave trains having very short wave length (relative to the core radius) but which vary slowly (on a scale of the core radius). The period of the waves is much longer than a day but much shorter than the period of the slow hydromagnetic oscillations. These waves are found to be trapped in a region about the equator and away from the axis of rotation. The waves break at a latitudinal wave region boundary, in the sense that the waves become exponentially large in a boundary layer, having as an exponent some positive power of the large azimuthal wave number. This behavior is amplified as the Elsasser number becomes smaller while still remaining relatively large. Waves in more Earth-like parameter regimes are discussed briefly.  相似文献   

5.
Abstract

We propose a method of derivation of global asymptotic solutions of the hydromagnetic dynamo problem at large magnetic Reynolds number. The procedure reduces to matching the local asymptotic forms for the magnetic field generated near individual extrema of generation strength. The basis of the proposed method, named here the Maximally-Efficient-Generation Approach (MEGA), is the assertion that properties of global asymptotic solutions of the kinematic dynamo are determined by the distribution of the generation strength near its leading extrema and by the number and distribution of the extrema.

The general method is illustrated by the global asymptotic solution of the α2-dynamo problem in a slab. The nature of oscillatory solutions revealed earlier in numerical simulations and the reasons for the dominance of even magnetic modes in slab geometry are clarified.

Applicability of the asymptotic solutions at moderate values of the asymptotic parameter is also discussed. We confirm this applicability using comparisons with complementary asymptotic expansions and numerical simulations. In particular, this justifies application of the MEGA solutions to estimation of the generation threshold.  相似文献   

6.
Abstract

Magnetic instabilities play an important role in the understanding of the dynamics of the Earth's fluid core. In this paper we continue our study of the linear stability of an electrically conducting fluid in a rapidly rotating, rigid, electrically insulating spherical geometry in the presence of a toroidal basic state, comprising magnetic field BMB O(r, θ)1ø and flow UMU O(r, θ)1ø The magnetostrophic approximation is employed to numerically analyse the two classes of instability which are likely to be relevant to the Earth. These are the field gradient (or ideal) instability, which requires strong field gradients and strong fields, and the resistive instability, dependent on finite resistivity and the presence of a zero in the basic state B O(r,θ). Based on a local analysis and numerical results in a cylindrical geometry we have established the existence of the field gradient instability in a spherical geometry for very simple basic states in the first paper of this series. Here, we extend the calculations to more realistic basic states in order to obtain a comprehensive understanding of the field gradient mode. Having achieved this we turn our attention to the resistive instability. Its presence in a spherical model is confirmed by the numerical calculations for a variety of basic states. The purpose of these investigations is not just to find out which basic states can become unstable but also to provide a quantitative measure as to how strong the field must become before instability occurs. The strength of the magnetic field is measured by the Elsasser number; its critical value c describing the state of marginal stability. For the basic states which we have studied we find c 200–1000 for the field gradient mode, whereas for the resistive modes c 50–160. For the field gradient instability, c increases rapidly with the azimuthal wavenumber m whereas in the resistive case there is no such pronounced difference for modes corresponding to different values of m. The above values of c indicate that both types of instability, ideal and resistive, are of relevance to the parameter regime found inside the Earth. For the resistive mode, as is increased from c, we find a shortening lengthscale in the direction along the contour BO = 0. Such an effect was not observable in simpler (for example, cylindrical) models.  相似文献   

7.
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.  相似文献   

8.
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.  相似文献   

9.
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 (Ω ? VA/d, where VA = B0√(μρ0) is the Alfvén velocity). The Elsasser number, Λ = σB2 0/2Ωρ0, measures σ. A (modified) Rayleigh number, R = gβd2/ρ0V2 A, measures the buoyancy force, where β is the imposed density gradient, antiparallel to g. A Prandtl number, PK = μσK, measures the diffusivity, k, of density differences.  相似文献   

10.
Abstract

Dynamic interaction between magnetic field and fluid motion is studied through a numerical experiment of nonlinear three-dimensional magnetoconvection in a rapidly rotating spherical fluid shell to which a uniform magnetic field parallel to its spin axis is applied. The fluid shell is heated by internal heat sources to maintain thermal convection. The mean value of the magnetic Reynolds number in the fluid shell is 22.4 and 10 pairs of axially aligned vortex rolls are stably developed. We found that confinement of magnetic flux into anti-cyclonic vortex rolls was crucial on an abrupt change of the mode of magnetoconvection which occurred at Δ = 1 ～ 2, where A is the Elsasser number. After the mode change, the fluid shell can store a large amount of magnetic flux in itself by changing its convection style, and the magnetostrophic balance among the Coriolis, Lorentz and pressure forces is established. Furthermore, the toroidal/poloidal ratio of the induced magnetic energy becomes less than unity, and the magnetized anti-cyclones are enlarged due to the effect of the magnetic force. Using these key ideas, we investigated the causes of the mode change of magnetoconvection. Considering relatively large magnetic Reynolds number and a rapid rotation rate of this model, we believe that these basic ideas used to interpret the present numerical experiment can be applied to the dynamics in the Earth's and other planetary cores.  相似文献   

11.
Abstract

An inviscid, electrically conducting fluid is contained between two rigid horizontal planes and bounded laterally by two vertical walls. The fluid is permeated by a strong uniform horizontal magnetic field aligned with the side wall boundaries and the entire system rotates rapidly about a vertical axis. The ratio of the magnitudes of the Lorentz and Coriolis forces is characterized by the Elsasser number, A, and the ratio of the thermal and magnetic diffusivities, q. By heating the fluid from below and cooling from above the system becomes unstable to small perturbations when the adverse density gradient as measured by the Rayleigh number, R, is sufficiently large.

With the viscosity ignored the geostrophic velocity, U, which is aligned with the applied magnetic field, is independent of the coordinate parallel to the rotation axis but is an arbitrary function of the horizontal cross-stream coordinate. At the onset of instability the value of U taken ensures that Taylor's condition is met. Specifically the Lorentz force, which results from marginal convection must not cause any acceleration of the geostrophic flow. It is found that the critical Rayleigh number characterising the onset of instability is generally close to the corresponding value for the usual linear problem, in which Taylor's condition is ignored and U is chosen to vanish. Significant differences can occur when q is small owing to a complicated flow structure. There is a central interior region in which the local magnetic Reynolds number, Rm , based on U is small of order q and on exterior region in which Rm is of order unity.  相似文献   

12.
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 Rm = 0 for q = 106. 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 Rac vs. Rm curve similar to the A = 1 case. However, it has Rac ≤ 0 at its minimum and unlike the results for high A, larger values of Rm result in a restabilisation.  相似文献   

13.
Abstract

We study the nonlinear stability of MHD waves propagating in a two-dimensional, compressible, highly magnetized, viscous plasma. These waves are driven by a weak, shear body force which could be imposed by large scale internal fluctuations present in the solar atmosphere.

The effects of anisotropic viscosity (leading to a cubic damping) and of the nonlinear coupling of the Alfven and the magnetoacoustic waves are analysed using Galerkin and multiple-scale analysis: the MHD equations are reduced to a set of nonlinear ordinary differential equations which is then suitably truncated to give a model dynamical system, representing the interaction of two complex Galerkin modes.

For propagation oblique to the background magnetic field, analytical integration shows that the low-wavenumber mode is physically unstable. For propagation parallel to the background magnetic field the high-wavenumber wave can undergo saddlenode bifurcations, in way that is similar to the van der Pol oscillator; these bifurcations lead to the appearance of a hysteresis cycle.

A numerical integration of the dynamical system shows that a sequence of Hopf bifurcations takes place as the Reynolds number is increased, up to the onset of nonperiodic behaviour. It also shows that energy can be transferred from the low- wavenumber to the high-wavenumber mode.  相似文献   

14.
Abstract

A spherical αω-dynamo is studied for small values of the viscous coupling parameter ε ～ v1/2, paying attention particularly to large dynamo numbers. The present study is a follow-up of the work by Hollerbach et al. (1992) with their choice of α-effect and Archimedean wind including also the constraint of magnetic field symmetry (or antisymmetry) due to equatorial plane. The magnetic field scaled by ε1/2 is independent of ε in the solutions for dynamo numbers smaller than a certain value of D b (the Ekman state) which are represented by dynamo waves running from pole to equator or vice-versa. However, for dynamo numbers larger than D b the solution bifurcates and subsequently becomes dependent on ε. The bifurcation is a consequence of a crucial role of the meridional convection in the mechanism of magnetic field generation. Calculations suggest that the bifurcation appears near dynamo number about 33500 and the solutions for larger dynamo numbers and ε = 0 become unstable and fail, while the solutions for small but non-zero ε are characterized by cylindrical layers of local maximum of magnetic field and sharp changes of geostrophic velocity. Our theoretical analysis allows us to conclude that our solution does not take the form of the usual Taylor state, where the Taylor constraint should be satisfied due to the special structure of magnetic field. We rather obtained the solution in the form of a “weak” Taylor state, where the Taylor constraint is satisfied partly due to the amplitude of the magnetic field and partly due to its structure. Calculations suggest that the roles of amplitude and structure are roughly fifty-fifty in our “weak” Taylor state solution and thus they can be called a Semi-Taylor state. Simple estimates show that also Ekman state solutions can be applicable in the geodynamo context.  相似文献   

15.
The velocity, pressure, perturbation magnetic field, helicity and electromotive force driven by an isolated buoyant parcel in an unbounded, rapidly rotating, electrically conducting fluid in the limit of small Elsasser number and very small Ekman number are calculated, visualized and analyzed. On the scale of the parcel, the solution is identical to that obtained in the limit of small Ekman number and zero Elsasser number. On the scale of the Taylor-column, it is elongated in the direction of the applied magnetic field and compressed in the direction perpendicular to it. The α-effect calculated by averaging the electromotive force on planes normal to rotation is strongly anisotropic: near the parcel and in the inner part of the Taylor-column it is strongest when the applied magnetic field is perpendicular to rotation and gravity; in the outer part of the Taylor-column it is strongest when the applied magnetic field is in the same plane as rotation and gravity.  相似文献   

16.
The magnetoconvection problem under the magnetostrophic approximation is investigated as the nonlinear regime is entered. The model consists of a fluid filled sphere, internally heated, and rapidly rotating in the presence of a prescribed, axisymmetric, toroidal magnetic field. For simplicity only a dipole parity and a single azimuthal wavenumber (m = 2) is considered here. The leading order nonlinearity at small amplitude is the geostrophic flow U g which is introduced to the previously linear model (Walker and Barenghi, 1997a, b). Walker and Barenghi (1997c) considered parameter space above critical and found that U g acts as an equilibration mechanism for moderately supercritical solutions. However, for solutions well above critical a Taylor state is approached and the system can no longer equilibrate. More importantly though, in the context of this paper, is that subcritical solutions were found. Here subcritical solutions are considered in more detail. It was found that, at is strongly dependent on . ( is the critical value of the modified Rayleigh number is a measure of the maximum amplitude of the generated geostrophic flow while , the Elsasser number, defines the strength of the prescribed toroidal field.) Rm at proves to be the key measure in determining how far into the subcritical regime the system can advance.  相似文献   

17.
  相似文献   

18.
Abstract

The stratification profile of the Earth's magnetofluid outer core is unknown, but there have been suggestions that its upper part may be stably stratified. Braginsky (1984) suggested that the magnetic analog of Rossby (planetary) waves in this stable layer (the ‘H’ layer) may be responsible for a portion of the short-period secular variation. In this study, we adopt a thin shell model to examine the dynamics of the H layer. The stable stratification justifies the thin-layer approximations, which greatly simplify the analysis. The governing equations are then the Laplace's tidal equations modified by the Lorentz force terms, and the magnetic induction equation. We linearize the Lorentz force in the Laplace's tidal equations and the advection term in the magnetic induction equation, assuming a zeroth order dipole field as representative of the magnetic field near the insulating core-mantle boundary. An analytical β-plane solution shows that a magnetic field can release the equatorial trapping that non-magnetic Rossby waves exhibit. A numerical solution to the full spherical equations confirms that a sufficiently strong magnetic field can break the equatorial waveguide. Both solutions are highly dissipative, which is a consequence of our necessary neglect of the induction term in comparison with the advection and diffusion terms in the magnetic induction equation in the thin-layer limit. However, were one to relax the thin-layer approximations and allow a radial dependence of the solutions, one would find magnetic Rossby waves less damped (through the inclusion of the induction term). For the magnetic field strength appropriate for the H layer, the real parts of the eigenfrequencies do not change appreciably from their non-magnetic values. We estimate a phase velocity of the lowest modes that is rather rapid compared with the core fluid speed typically presumed from the secular variation.  相似文献   

19.
Abstract

Finite amplitude solutions for convection in a rotating spherical fluid shell with a radius ratio of η=0.4 are obtained numerically by the Galerkin method. The case of the azimuthal wavenumber m=2 is emphasized, but solutions with m=4 are also considered. The pronounced distinction between different modes at low Prandtl numbers found in a preceding linear analysis (Zhang and Busse, 1987) is also found with respect to nonlinear properties. Only the positive-ω-mode exhibits subcritical finite amplitude convection. The stability of the stationary drifting solutions with respect to hydrodynamic disturbances is analyzed and regions of stability are presented. A major part of the paper is concerned with the growth of magnetic disturbances. The critical magnetic Prandtl number for the onset of dynamo action has been determined as function of the Rayleigh and Taylor numbers for the Prandtl numbers P=0.1 and P=1.0. Stationary and oscillatory dynamos with both, dipolar and quadrupolar, symmetries are close competitors in the parameter space of the problem.  相似文献   

20.
Abstract

The weak-field Benard-type dynamo treated by Soward is considered here at higher levels of the induced magnetic field. Two sources of instability are found to occur in the intermediate field regime M ～ T 1/12, where M and T are the Hartmann and Taylor numbers. On the time scale of magnetic diffusion, solutions may blow up in finite time owing to destabilization of the convection by the magnetic field. On a faster time scale a dynamic instability related to MAC-wave instability can also occur. It is therefore concluded that the asymptotic structure of this dynamo is unstable to virtual increases in the magnetic field energy.

In an attempt to model stabilization of the dynamo in a strong-field regime we consider two approximations. In the first, a truncated expansion in three-dimensional plane waves is studied numerically. A second approach utilizes an ad hoc set of ordinary differential equations which contains many of the features of convection dynamos at all field energies. Both of these models exhibit temporal intermittency of the dynamo effect.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号