首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
An essentially arbitrary function V(θ, λ) defined on the surface of a sphere can be expressed in terms of spherical harmonics V(θ, Λ) = a n=1m=0npmn (cos θ) (gmn cos mΛ + hmn sin mΛ) where the Pnm are the seminormalized associated Legendre polynomials used in geomagnetism, normalized so that 〈[Pmn(cos θ) cos mΛ]〉2 =1/(2n+1) The angular brackets denote an average over the sphere. The class of functions V(θ, λ) under consideration is that normally of interest in physics and engineering. If we consider an ensemble of all possible orientations of our coordinate system relative to the sphere, then the coefficients gnm and hnm will be functions of the particular coordinate system orientation, but 〈:(gmn)2〉) = 〈(hmn)2 = Sn/(2n=1) where Sn =m=0n [(gmn)2 + (hmn)2] for any orientation of the coordinate system (Sn is invariant under rotation of the coordinate system). The averages are over all orientations of the system relative to the sphere. It is also shown that 〈gmnglp〉 = 〈hmnhlp〉 = 0 for lm or pn and 〈gmnhlp〉 = 0 fro all n, m, p, l.  相似文献   

2.
The correlation discovered by Hide and Malin between the variable parts of the Earth's gravitational field and magnetic field (suitably displaced in longitude) was tentatively and qualitatively explained by them in terms of the influence on both fields of irregularities (or “surface bumps”) at the core-mantle interface. In this paper, a quantitative analysis of this phenomenon is developed, through study of an idealised problem in which conducting fluid occupying the region z < η(x) flows over the surface z = η(x) in the presence of a magnetic field (B0,0,0), the whole system rotating with angular velocity (0,0,Ω). It is assumed that |η′(x)| « 1 so that perturbation methods are applicable. Determination of the magnetic potential in the “mantle” region z < η(x) requires solution of the full hydromagnetic problem in the fluid. It is shown that three wave modes are excited, two of which (for values of the parameters of the problem of geophysical interest) have a boundary layer character. Phase interactions between these modes lead to a shift and a distortion of the magnetic pattern relative to the gravitational pattern. The correlation between the gravitational potential and the magnetic potential (shifted by a distance x0) is determined on the plane z = d (d a? |η|) as a function of x0/d and the curves obtained are qualitatively similar to that based on the observed data; the maximum correlation obtained varies between 0.67 and 1, depending on values of the parameters of the problem, and is about 0.72 for reasonable estimates of these parameters in the geophysical context.  相似文献   

3.
A rigorous singular perturbation theory is developed to estimate the electric field E produced in the mantle M by the core dynamo when the electrical conductivity σ in M depends only on radius r, and when |r?rln σ| ? 1 in most of M. It is assumed that σ has only one local minimum in M, either (a) at the Earth's surface ?V, or (b) at a radius b inside the mantle, or (c) at the core-mantle boundary ?K. In all three cases, the region where σ is no more than e times its minimum value constitutes a thin critical layer; in case (a), the radial electric field Er ≈ 0 there, while in cases (b) and (c), Er is very large there. Outside the critical layer, Er ≈ 0 in all three cases. In no case is the tangential electric field ES small, nearly toroidal, or nearly calculable from the magnetic vector potential A as ??tAS. The defects in Muth's (1979) argument which led him to contrary conclusions are identified. Benton (1979) cited Muth's work to argue that the core-fluid velocity u just below ?K can be estimated from measurements on ?V of the magnetic field B and its time derivative ?tB. A simple model for westward drift is discussed which shows that Benton's conclusion is also wrong.In case (a), it is shown that knowledge of σ in M is unnecessary for estimating ES on ?K with a relative error |r?r 1nσ|?1from measurements of ES on ?V and knowledge of ?tB in M (calculable from ?tB on ?V if σ is small). Then, in case (a), u just below ?K can be estimated as ?r×ES/Br. The method is impractical unless the contribution to ES on ?V from ocean currents can be removed.The perturbation theory appropriate when σ in M is small is considered briefly; smallness of σ and of |r?r ln σ|?1 a independent questions. It is found that as σ → 0, B approaches the vacuum field in M but E does not; the explanation lies in the hydromagnetic approximation, which is certainly valid in M but fails as σ → 0. It is also found that the singular perturbation theory for |r?r ln σ|?1 is a useful tool in the perturbation calculations for σ when both σ and |r?r ln σ|?1 are small.  相似文献   

4.
Abstract

We discuss the steady states of the αω-dynamo in a thin disc which arise due to α-quenching. Two asymptotic regimes are considered, one for the dynamo numberD near the generation thresholdD 0, and the other for |D| ? 1. Asymptotic solutions for |D—D 0| ? |D 0| have a rather universal character provided only that the bifurcation is supercritical. For |D| ? 1 the asymptotic solution crucially depends on whether or not the mean helicity α, as a function ofB, has a positive root (hereB is the mean magnetic field). When such a root exists, the field value in the major portion of the disc is O(l), while near the disc surface thin boundary layers appear where the field rapidly decreases to zero (if the disc is surrounded by vacuum). Otherwise, when α = O(|B|?s) for |B| → ∞, we demonstrate that |B| = O(|D|1/s ) and the solution is free of boundary layers. The results obtained here admit direct comparison with observations of magnetic fields in spiral galaxies, so that an appropriate model of nonlinear galactic dynamos hopefully could be specified.  相似文献   

5.
The mean tangential stresses at a corrugated interface between a solid, electrically insulating mantle and a liquid core of magnetic diffusivity λ are calculated for uniform rotation of both mantle and core at an angular velocity Ω in the presence of a corotating magnetic field B. The core and mantle are assumed to extend indefinitely in the horizontal plane. The interface has the form z = η(x, y), where z is the upward vertical distance and x, y are the zonal and latitudinal distances respectively. The function η(x, y) has a planetary horizontal length scale (i.e. of the order of the radius of the Earth) and small amplitude and vertical gradient. The liquid core flows with uniform mean zonal velocity U0 relative to the mantle. Ω and B possess vertical and horizontal components.The vertical (poloidal) component Bp is uniform and has a value of 5 G while the horizontal (toroidal) field BT = Bpαz, where α is a constant. When |α| ? 1, the mean horizontal stresses are found to have the same order of magnitude (10?2 N m?2) as those inferred from variations in the decade fluctuations in the length of the day, although the exact numerical values depend on the orientation of Ω as well as on the wavenumbers in the zonal and latitudinal directions.The influence of the steepness (as measured by α) of the toroidal field on the stresses is investigated to examine whether the constraint that the mean horizontal stresses at the core-mantle interface be of the order of 10?2 N m?2 might provide a selection mechanism for the behaviour of the toroidal field in the upper reaches of the outer core of the Earth. The results indicate that the restriction imposed on α is related to the value assigned to the toroidal field deep into the core. For example, if |α| ? 1 then the tangential stresses are of the right order of magnitude only if the toroidal field is comparable with the poloidal field deep in the core.  相似文献   

6.
The dependence of bulk sound speed Vφ upon mean atomic weight m and density ρ can be expressed in a single equation:
Vφ=Bρλ(m0m[12+λ(1?c)](km/sec)
Here B is an empirically determined “universal” parameter equal to 1.42, m0 = 20.2, a reference mean atomic weight for which well-determined elastic properties exist, and λ = 1.25 is a semi empirical parameter equal to γ ? 13 where γ is a Grüneisen parameter. The constant c = (? ln VM/? ln m)X, where VM is molar volume, is in general different for different crystal structure series and different cation substitutions. However, it is possible to use cFe = 0.14 for Fe2+Mg2+ and GeSi substitutions and cCa ? 1.3 for CaMg substitutional series. With these values it is pos to deduce from the above equation Birch's law, its modifications introduced by Simmons to account for Ca-bearing minerals, variations in the seismic equation of state observed by D.L. Anderson, and the apparent proportionality of bulk modulus K to VM?4.  相似文献   

7.
Lowes (1966, 1974) has introduced the function Rn defined by Rn =(n + 1) m=0 [(gmn)2 + (hmn)2] where gnmand hnm are the coefficients of a spherical harmonic expansion of the scalar potential of the geomagnetic field at the Earth's surface. The mean squared value of the magnetic field B = ??V on a sphere of radius r > α is given by B ·〉 =n=1 Rn(a/r)2n=4where a is the Earth's radius. We refer to Rn as the spherical harmonic spatial power spectrum of the geomagnetic field.In this paper it is shown that Rn = RMn = RCn where the components RnM due to the main (or core) field and RnC due to the crustal field are given approximately by RMn = [(n =1)/(n + 2)](1.142 × 109)(0.288n Λ2 RCn = [(n =1){[1 — exp(-n/290)]/(n/290)} 0.52 Λ2where Iγ = 1 nT. The two components are approximately equal for n = 15.Lowes has given equations for the core and crustal field spectra. His equation for the crustal field spectrum is significantly different from the one given here. The equation given in this paper is in better agreement with data obtained on the POGO spacecraft and with data for the crustal field given by Alldredge et al. (1963).The equations for the main and crustal geomagnetic field spectra are consistent with data for the core field given by Peddie and Fabiano (1976) and data for the crustal field given by Alldredge et al. The equations are based on a statistical model that makes use of the principle of equipartition of energy and predicts the shape of both the crustal and core spectra. The model also predicts the core radius accurately. The numerical values given by the equations are not strongly dependent on the model.Equations relating average great circle power spectra of the geomagnetic field components to Rn are derived. The three field components are in the radial direction, along the great circle track, and perpendicular to the first two. These equations can, in principle, be inverted to compute the Rn for celestial bodies from average great circle power spectra of the magnetic field components.  相似文献   

8.
Virtually all dynamo models may be expected to give rise to a permanent differential rotation between mantle and core. Weak conductivity in the mantle permits small leakage currents which couple to the radial component of the magnetic field, producing a Lorentz torque. Mechanical equilibrium is achieved when a zero net torque is established at a critical rotation rate. An estimate of the drift is determined easily given the magnetic field structure predicted by any dynamo model. The result for the drift rate at the core-mantle interface along the equator is given by the product of three factors
Uφ1=UφRλ*L*
The first of these is a geometrical factor which depends only on the structural character of the field. For a variety of model fields, this factor ranges from 16 to 35. The second factor is the ratio of r.m.s. toroidal to poloidal field. This ratio is an (implicitly) adjustable parameter of both α2 and α-ω dynamos, and is a measure of the relative efficiency of the generation process for each component. The third (dimensional) term is the ratio of core magnetic diffusivity to core radius, 10?4 cm s?1.The result is essentially independent of the value of mantle diffusivity and its effective depth. The sign of the result may be positive or negative. For α2 dynamos a westward drift is produced by choosing α > 0 in the Northern Hemisphere, which constitutes a dynamical assertion about the dynamo process. For an r.m.s. toroidal field of the order of 15 Gs, based on fairly general considerations, a drift rate comparable to observation is expected.  相似文献   

9.
Investigation of magnetic field generation by convective flows is carried out for three values of kinematic Prandtl number: P = 0.3, 1 and 6.8. We consider Rayleigh–Bénard convection in Boussinesq approximation assuming stress-free boundary conditions on horizontal boundaries and periodicity with the same period in the x and y directions. Convective attractors are modelled for increasing Rayleigh numbers for each value of the kinematic Prandtl number. Linear and non-linear dynamo action of these attractors is studied for magnetic Prandtl numbers P m ≤ 100. Flows, which can act as magnetic dynamos, have been found for all the three considered values of P, if the Rayleigh number R is large enough. The minimal R, for which of magnetic field generation occurs, increases with P. The minimum (over R) of critical Pm for magnetic field generation in the kinematic regime is admitted for P = 0.3. Thus, our study indicates that smaller values of P are beneficial for magnetic field generation.  相似文献   

10.
11.
12.
A crucial step in the investigation of the energetics of motions in the Earth's core and the generation of the geomagnetic field by the hydromagnetic dynamo process is the estimation of the average strength B of the magnetic field B = Bp + BT in the core. Owing to the probability that the toroidal field BT in the core, which has no radial component, is a good deal stronger than the poloidal field Bp, direct downward extrapolation of the surface field to the core-mantle interface gives no more than an extreme lower limit to B. This paper outlines the indirect methods by which B can be estimated, arguing that B is probably about 10?2 T (100 Γ) but might be as low as 10?3 T (10 Γ) or as high as 5 × 10?2 T (500 Γ).  相似文献   

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

14.
The “anomalous” layer in the lowermost mantle, identified as D″ in the notation of K.E. Bullen, appears in the PREM Earth model as a 150 km-thick zone in which the gradient of incompressibility with pressure, dKdP, is almost 1.6, instead of 3.2 as in the overlying mantle. Since PREM shows no accompanying change in density or density gradient, we identify D″ as a thermal boundary layer and not as a chemically distinct zone. The anomaly in dKdP is related to the temperature gradient by the temperature dependence of Ks, for which we present a thermodynamic identity in terms of accessible quantities. This gives the numerical result (?Ks/?T)P=?1.6×107 Pa K?1 for D″ material. The corresponding temperature increment over the D″ range is 840 K. Such a layer cannot be a static feature, but must be maintained by a downward motion of the lower mantle toward the core-mantle boundary with a strong horizontal flow near the base of D″. Assuming a core heat flux of 1.6 × 1012 W, the downward speed is 0.07 mm y?1 and the temperature profile in D″, scaled to match PREM data, is approximately exponential with a scale height of 73 km. The inferred thermal conductivity is 1.2 W m?1 K?1. Using these values we develop a new analytical model of D″ which is dynamically and thermally consistent. In this model, the lower-mantle material is heated and softened as it moves down into D″ where the strong temperature dependence of viscosity concentrates the horizontal flow in a layer ~ 12 km thick and similarly ensures its removal via narrow plumes.  相似文献   

15.
Abstract

Convection in a rotating spherical shell has wide application for understanding the dynamics of the atmospheres and interiors of many celestial bodies. In this paper we review linear results for convection in a shell of finite depth at substantial but not asymptotically large Taylor numbers, present nonlinear multimode calculations for similar conditions, and discuss the model and results in the context of the problem of solar convection and differential rotation. Detailed nonlinear calculations are presented for Taylor number T = 105, Prandtl number P = 1, and Rayleigh number R between 1 |MX 104 and 4 |MX 104 (which is between about 4 and 16 times critical) for a shell of depth 20% of the outer radius. Sixteen longitudinal wave numbers are usually included (all even wave numbers m between 0 and 30) the amplitudes of which are computed on a staggered grid in the meridian plane.

The kinetic energy spectrum shows a peak in the wave number range m = 12–18 at R = 104, which straddles the critical wave number m = 14 predicted by linear theory. These are modes which peak near the equator. The spectrum shows a second strong peak at m = 0, which represents the differential rotation driven by the peak convective modes. As R is increased, the amplitude of low wave numbers increases relative to high wave numbers as convection fills in in high and middle latitudes, and as the longitudinal scale of equatorial convection grows. By R = 3 |MX 104, m = 8 is the peak convective mode. There is a clear minimum in the total kinetic energy at middle latitudes relative to low and high, well into the nonlinear regime, representing the continued dominance of equatorial and polar modes found in the linear case. The kinetic energy spectrum for m > 0 is maintained primarily by buoyancy work in each mode, but with substantial nonlinear transfer of kinetic energy from the peak modes to both lower and higher wave numbers.

For R = 1 to 2 |MX 104, the differential rotation takes the form of an equatorial acceleration, with angular velocity generally decreasing with latitude away from the equator (as on the sun) and decreasing inwards. By R = 4 |MX 104, this equatorial profile has completely reversed, with angular velocity increasing with depth and latitude. Also, a polar vortex which has positive rotation relative to the reference frame (no evidence of which has been seen on the sun) builds up as soon as polar modes become important. Meridional circulation is quite weak relative to differential rotation at R = 104, but grows relative to it as R is increased. This circulation takes the farm of a single cell of large latitudinal extent in equatorial regions, with upward flow near the equator, together with a series of narrower cells in high latitudes. It is maintained primarily by axisymmetric buoyancy forces. The differential rotation is maintained at all R primarily by Reynolds stresses, rather than meridional circulation. Angular momentum transport toward the equator for R = 1–2 |MX 104 maintains the equatorial acceleration while radially inward transport maintains the opposite profile at R = 4 |MX 104.

The total heat flux out the top of the convective shell always shows two peaks for the range of R studied, one at the equator and the other near the poles (no significant variation with latitude is seen on the sun), while heat flux in at the bottom shows only a polar peak at large R. The meridional circulation and convective cells transport heat toward the equator to maintain this difference.

The helicity of the convection plus the differential rotation produced by it suggest the system may be capable of driving a field reversing dynamo, but the toroidal field may migrate with lime in each cycle toward the poles and equator, rather than just toward the equator as apparently occurs on the sun.

We finally outline additions to the physics of the model to make it more realistic for solar application.  相似文献   

16.
The radiation belt dynamics during the extreme solar events in November 2004 and January 2005 is studied based on the measurements of relativistic electrons (with energies of 0.8–8 MeV) on the Express-A2 geostationary satellite and Meteor-3M polar satellite. New radiation belts of relativistic electrons in the space (L ~ 3) between the stationary outer and inner belts were formed as a result of either superstorm (|Dst|max = 373 and 289 nT). The position of the maximums of these belts (L max = 2.9 and 3.1) coincides with the known dependence of L max on the magnetic storm Dst variation amplitude: |Dst|max = 2.75 × 104/L max 4 . In November–December the new belt very slowly (ΔL ~ 0.1 per month) shifted toward the Earth. During the series of moderate (~100 nT) magnetic storms that developed as a result of the extreme solar events in January 2005, the belt in the space shifted toward deeper L shells (L ~ 2.5). The moderate January storms produced new belts with L max ≥ 4.  相似文献   

17.
The dependence of the maximal values of the |Dst| and AE geomagnetic indices observed during magnetic storms on the value of the interplanetary electric field (E y ) was studied based on the catalog of the large-scale solar wind types created using the OMNI database for 1976–2000 [Yermolaev et al., 2009]. An analysis was performed for eight categories of magnetic storms caused by different types of solar wind streams: corotating interaction regions (CIR, 86 storms); magnetic clouds (MC, 43); Sheath before MCs (ShMC, 8); Ejecta (95); Sheath (ShE, 56); all ICME events (MC + Ejecta, 138); all compression regions Sheaths before MCs and Ejecta (ShMC + ShE, 64); and an indeterminate type of storm (IND, 75). It was shown that the |Dst| index value increases with increasing electric field E y for all eight types of streams. When electric fields are strong (E y > 11 mV m−1), the |Dst| index value becomes saturated within magnetic clouds MCs and possibly within all ICMEs (MC + Ejecta). The AE index value during magnetic storms is independent of the electric field value E y for almost all streams except magnetic clouds MCs and possibly the compressed (Sheath) region before them (ShMC). The AE index linearly increases within MC at small values of the electric field (E y < 11 mV m−1) and decrease when these fields are strong (E y > 11 mV m−1). Since the dynamic pressure (Pd) and IMF fluctuations (σB) correlate with the E y value in all solar wind types, both geomagnetic indices (|Dst| and AE) do not show an additional dependence on Pd and IMF δB. The nonlinear relationship between the intensities of the |Dst| and AE indices and the electric field E y component, observed within MCs and possibly all ICMEs during strong electric fields E y , agrees with modeling the magnetospheric-ionospheric current system of zone 1 under the conditions of the polar cap potential saturation.  相似文献   

18.
In the present study, the deterministic chaotic behaviour of interplanetary magnetic field (IMF) under various geomagnetic conditions of low and high solar active periods was analyzed, using the time series of IMF |B| and Bz, by employing chaotic quantifiers like, Lyapunov exponent, Tsallis entropy, correlation dimension, and non-linear prediction error. We have investigated whether the chaotic behaviour of interplanetary magnetic field would modify, when it produces major geomagnetic storms, and how it depends on the phase of solar activity. The yearly average values of Lyapunov exponent for the time series of IMF |B| and Bz, show solar flux dependence, whereas those values of entropy, correlation dimension and non-linear prediction error had no significant solar flux dependence. The yearly average values of entropy for quiet periods are higher compared to those values for major storm periods belonging to low/high solar active conditions, for both the time series |B| and Bz.  相似文献   

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

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