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

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell heated from below and within have been carried out with a nonlinear, three-dimensional, time-dependent pseudospectral code. The investigated phenomena include the sequence of transitions to chaos and the differential mean zonal rotation. At the fixed Taylor number T a =106 and Prandtl number Pr=1 and with increasing Rayleigh number R, convection undergoes a series of bifurcations from onset of steadily propagating motions SP at R=R c = 13050, to a periodic state P, and thence to a quasi-periodic state QP and a non-periodic or chaotic state NP. Examples of SP, P, QP, and NP solutions are obtained at R = 1.3R c , R = 1.7 R c , R = 2R c , and R = 5 R c , respectively. In the SP state, convection rolls propagate at a constant longitudinal phase velocity that is slower than that obtained from the linear calculation at the onset of instability. The P state, characterized by a single frequency and its harmonics, has a two-layer cellular structure in radius. Convection rolls near the upper and lower surfaces of the spherical shell both propagate in a prograde sense with respect to the rotation of the reference frame. The outer convection rolls propagate faster than those near the inner shell. The physical mechanism responsible for the time-periodic oscillations is the differential shear of the convection cells due to the mean zonal flow. Meridional transport of zonal momentum by the convection cells in turn supports the mean zonal differential rotation. In the QP state, the longitudinal wave number m of the convection pattern oscillates among m = 3,4,5, and 6; the convection pattern near the outer shell has larger m than that near the inner shell. Radial motions are very weak in the polar regions. The convection pattern also shifts in m for the NP state at R = 5R c , whose power spectrum is characterized by broadened peaks and broadband background noise. The convection pattern near the outer shell propagates prograde, while the pattern near the inner shell propagates retrograde with respect to the basic rotation. Convection cells exist in polar regions. There is a large variation in the vigor of individual convection cells. An example of a more vigorously convecting chaotic state is obtained at R = 50R c . At this Rayleigh number some of the convection rolls have axes perpendicular to the axis of the basic rotation, indicating a partial relaxation of the rotational constraint. There are strong convective motions in the polar regions. The longitudinally averaged mean zonal flow has an equatorial superrotation and a high latitude subrotation for all cases except R = 50R c , at this highest Rayleigh number, the mean zonal flow pattern is completely reversed, opposite to the solar differential rotation pattern.  相似文献   

2.
Abstract

Measurements have been made of the net horizontal force F acting on a sphere moving with horizontal velocity U (Reynolds numbers in the range 102-104) through a stratified fluid rotating about a vertical axis with uniform angular velocity Ω. In both homogeneous and stratified rotating fluids with small Rossby number R(R = U/Ωa ? 1 where a is the radius of the sphere) the force F is of magnitude 2ΩρUV (where ρ is the density of the fluid and V is the volume of the sphere). In a homogeneous fluid the relative directions of F and U were found to depend on the quantity F = 8Ωa 2/UD (where D is the depth of the fluid in which the object is placed (Mason, 1975)). In a rotating stratified fluid the relative directions of F and U are found to depend on the inverse Froude number k(k = Na/U where N 2 = (g/δ)?ρ/?z) provided D > 4aΩ/N. In a homogeneous fluid with F ? 1 the force F is mainly in the U direction (a drag force due to inertial wave radiation) and is ～ ?0.4 |MX 2ΩρUV For F ? 1 a “Taylor column” occurs and the force, in correspondence with theoretical expectations, is ～ - 2Ω |MX UρV In a rotating stratified fluid with N ～2Ω and k ? 1 the force F is mainly in the U direction but is roughly one half of that occurring in the homogeneous situation with F ? 1 (tentatively explained as due to the evanescence of inertia-gravity disturbances). In a rotating stratified fluid with k ? 1 the flow should have no vertical motion (as with F ? 1) and again in correspondence with theoretical expectations the drag is ～ ?2 Ω |MX UρV. In a non-rotating stratified fluid the drag coefficient C D(C D = F U/½?ρU 2) was measured in the range k = 0.1 to 10 and had a maximum value ～ 1.2 for k ～ 3.  相似文献   

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

4.
Linear α2Ω-dynamo waves are investigated in a thin turbulent, differentially rotating convective stellar shell. A simplified one-dimensional model is considered and an asymptotic solution constructed based on the small aspect ratio of the shell. In a previous paper Griffiths et al. (Griffiths, G.L., Bassom, A.P., Soward, A.M. and Kuzanyan, K.M., Nonlinear α2Ω-dynamo waves in stellar shells, Geophys. Astrophys. Fluid Dynam., 2001, 94, 85–133) considered the modulation of dynamo waves, linked to a latitudinal-dependent local α-effect and radial gradient of the zonal shear flow. These effects are measured at latitude θ by the magnetic Reynolds numbers R α f(θ) and R Ω g(θ). The modulated Parker wave, which propagates towards the equator, is localised at some mid-latitude θp under a Gaussian envelope. In this article, we include the influence of a latitudinal-dependent zonal flow possessing angular velocity Ω*(θ) and consider the possibility of non-axisymmetric dynamo waves with azimuthal wave number m. We find that the critical dynamo number D c?=?R α R Ω is minimised by axisymmetric modes in the αΩ-limit (Rα→0). On the other hand, when Rα?≠?0 there may exist a band of wave numbers 0?m?m ? for which the non-axisymmetric modes have a smaller D c than in the axisymmetric case. Here m ? is regarded as a continuous function of R α with the property m?→0 as R α→0 and the band is only non-empty when m??>1, which happens for sufficiently large R α. The preference for non-axisymmetric modes is possible because the wind-up of the non-axisymmetric structures can be compensated by phase mixing inherent to the α2Ω-dynamo. For parameter values resembling solar conditions, the Parker wave of maximum dynamo activity at latitude θp not only propagates equatorwards but also westwards relative to the local angular velocity Ω*(θ p ). Since the critical dynamo number D c?=?R α R Ω is O (1) for small R α, the condition m ??>?1 for non-axisymmetric mode preference imposes an upper limit on the size of |dΩ*/dθ|.  相似文献   

5.
The effects of rotation and a toroidal magnetic field on the preferred pattern of small amplitude convection in spherical fluid shells are considered. The convective motions are described in terms of associated Legendre functions Pl|m| (cos θ). For a given pair of Prandtl number P and magnetic Prandtl number Pm the physically realized solution is represented either by m = 0 or |m| = l depending on the ratio of the rotation rate Λ to the magnetic field amplitude H. The case of m = 0 is preferred if this ratio ranges below a critical value, which is a function of the shell thickness, and |m| = l otherwise.  相似文献   

6.
In this article we study the linear instability of the two-dimensional strongly stratified model for global MHD in the diffusive solar tachocline. Gilman and Fox [Gilman, P.A. and Fox, P., Joint instability of the latitudinal differential rotation and toroidal magnetic fields below the solar convection zone. Astrophys. J., 1997, 484, 439–454] showed that for ideal MHD, the observed surface differential rotation becomes more unstable than is predicted by Watson's [Watson, M., Shear instability of differential rotation in stars. Geophys. Astrophys. Fluid Dyn., 1981, 16, 285–298] nonmagnetic analysis. They showed that the solar differential rotation is unstable for essentially all reasonable values of the differential rotation in the presence of an antisymmetric toroidal field. They found that for the broad field case B φ～sinθcosθ, θ being the co-latitude, instability occurs only for the azimuthal m?=?1 mode, and concluded that modes which are symmetric (meridional flow in the same direction) about the equator onset at lower field strengths than the antisymmetric modes. We study the effect of viscosity and magnetic diffusivity in the strongly stably stratified case where diffusion is primarily along the level surfaces. We show that antisymmetric modes are now strongly preferred over symmetric modes, and that diffusion can sometimes be destabilising. Even solid body rotation can be destabilised through the action of magnetic field. In addition, we show that when diffusion is present, instability can occur when the longitudinal wavenumber m?=?2.  相似文献   

7.
Abstract

A linear analysis is used to study the stability of a rapidly rotating, electrically-conducting, self-gravitating fluid sphere of radius r 0, containing a uniform distribution of heat sources and under the influence of an azimuthal magnetic field whose strength is proportional to the distance from the rotation axis. The Lorentz force is of a magnitude comparable with that of the Coriolis force and so convective motions are fully three-dimensional, filling the entire sphere. We are primarily interested in the limit where the ratio q of the thermal diffusivity κ to the magnetic diffusivity η is much smaller than unity since this is possibly of the greatest geophysical relevance.

Thermal convection sets in when the temperature gradient exceeds some critical value as measured by the modified Rayleigh number Rc. The critical temperature gradient is smallest (Rc reaches a minimum) when the magnetic field strength parameter Λ ? 1. [Rc and Λ are defined in (2.3).] The instability takes the form of a very slow wave with frequency of order κ/r 2 0 and its direction of propagation changes from eastward to westward as Λ increases through Λ c ? 4.

When the fluid is sufficiently stably stratified and when Λ > Λm ? 22 a new mode of instability sets in. It is magnetically driven but requires some stratification before the energy stored in the magnetic field can be released. The instability takes the form of an eastward propagating wave with azimuthal wavenumber m = 1.  相似文献   

8.
An axisymmetric model of convection in a rotating cylinder in an external uniform magnetic field has been considered. In the considered model, the meridional circulation is created by a nonuniform rotation of the lower boundary relative to the other boundaries. In the considered model, the time of formation of the stationary regime in the magnetic field considerably increases if the vertical density (compressibility) inhomogeneity is taken into account for Ekman numbers of E = E M = 3 × 10−3. This example shows that the compressibility of a medium should be taken into account in the convection and dynamics of the magnetic field when the magnetohydrodynamics of the Earth is analyzed.  相似文献   

9.
Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell with and without inhomogeneous temperature anomalies on the top boundary have been carried out using a three-dimensional, time-dependent, spectral-transform code. The spherical shell of Boussinesq fluid has inner and outer radii the same as those of the Earth's liquid outer core. The Taylor number is 107, the Prandtl number is 1, and the Rayleigh number R is 5Rc (Rc is the critical value of R for the onset of convection when the top boundary is isothermal and R is based on the spherically averaged temperature difference across the shell). The shell is heated from below and cooled from above; there is no internal heating. The lower boundary of the shell is isothermal and both boundaries are rigid and impermeable. Three cases are considered. In one, the upper boundary is isothermal while in the others, temperature anomalies with (l,m) = (3,2) and (6,4) are imposed on the top boundary. The spherically averaged temperature difference across the shell is the same in all three cases. The amplitudes of the imposed temperature anomalies are equal to one-half of the spherically averaged temperature difference across the shell. Convective structures are strongly controlled by both rotation and the imposed temperature anomalies suggesting that thermal inhomogeneities imposed by the mantle on the core have a significant influence on the motions inside the core. The imposed temperature anomaly locks the thermal perturbation structure in the outer part of the spherical shell onto the upper boundary and significantly modifies the velocity structure in the same region. However, the radial velocity structure in the outer part of the shell is different from the temperature perturbation structure. The influence of the imposed temperature anomaly decreases with depth in the shell. Thermal structure and velocity structure are similar and convective rolls are more columnar in the inner part of the shell where the effects of rotation are most dominant.  相似文献   

10.
Summary The problem of the gross nature of the Jovian atmospheric circulation is examined from the viewpoint of the following previous findings of the writer and others. 1) The equatorial acceleration cannot be accounted for by axisymmetric motions. 2) The departures from symmetry in a rotating system having an equatorial acceleration must impart angular momentum selectively to those particles moving toward the jet maximum and abstract it from those moving away. 3) These selective (pressure) torques and associated sorting processes arise spontaneously in the presence of a vertical convection mode involving motions not independent of longitude, if the cell sizes and other conditions are right.Since there is evidence that Jovian dark spots have statistical maxima of occurrence along the tropical shear lines flanking the equator, these are assumed to be vertical convective systems forming, in effect,convective vortex sheets which generate the high angular momentum of the equatorial zone. Various additional concepts are discussed, and many comparisons with conditions in the sun and in the earth's atmosphere are made.  相似文献   

11.
We investigate instability of convective flows of simple structure (rolls, standing and travelling waves) in a rotating layer with stress-free horizontal boundaries near the onset of convection. We show that the flows are always unstable to perturbations, which are linear combinations of large-scale modes and short-scale modes, whose wave numbers are close to those of the perturbed flows. Depending on asymptotic relations of small parameters α (the difference between the wave number of perturbed flows and the critical wave number for the onset of convection) and ε (ε2 being the overcriticality and the perturbed flow amplitude being O(ε)), either small-angle or Eckhaus instability is prevailing. In the case of small-angle instability for rolls the largest growth rate scales as ε8/5, in agreement with results of Cox and Matthews (Cox, S.M. and Matthews, P.C., Instability of rotating convection. J. Fluid. Mech., 2000, 403, 153–172) obtained for rolls with k = k c . For waves, the largest growth rate is of the order ε4/3. In the case of Eckhaus instability the growth rate is of the order of α2.  相似文献   

12.
朱涛  马宗晋  冯锐 《地球物理学报》2006,49(5):1347-1358
建立了三维黏度扰动下的变黏度地幔对流模型,并提供了在引入地幔的三维地震波速度结构下相应的求解方法. 依此反演了瑞利数Ra = 106时,两种不同边界条件下的极、环型场对流图像,这有助于深化对地幔物质流动和大地构造运动的深部动力学过程的认识和理解. 研究结果表明,不但地幔浅部的极型场对流图像显示出了与大地构造运动的相关性并揭示了其深部动力学过程,更重要的是,地幔浅部的环型场对流图像首次为我们认识和理解板块构造的水平与旋转运动提供了重要的信息：环型场速度剖面中在赤道附近存在一条大致南东东—北西西向的强对流条带,可能与环赤道附近大型剪切带的形成相关,进而表明可能是该带强震发生的深部动力学背景;在南北半球存在的旋转方向相反的对流环表明它们整体上可能存在差异旋转.  相似文献   

13.
Regional seismic characteristics in northeastern and southwestern China are quite different. Earthquakes are more frequent in SW China than in NE China and intensity attenuates faster in SW China than in NE China. Through regression analyses of data from these two regions we have found the attenuation of peak horizontal accelerations also to be noticeably different. For northeastern China the regression is log10(a) = ?0.474 + 0.613M ? 0.873 log10 (R) ?0.00206R and for southwestern China, the regression is log10 (a) = 0.437 + 0.454M ? 0.739 log10 (R) ?0.00279R where a is acceleration in cm/sec2, M is the Chinese surface wave magnitude (Institute of Geophysics, State Seismological Bureau), and R is the distance in kilometres. The formula for northeastern China has distance coefficients similar to those for Central United States as derived by Nuttli and Herrmann14.  相似文献   

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

15.
Measurements with a HF Doppler sounder at Kodaikanal (10.2°N, 77.5°E, geomagnetic latitude 0.8°N) showed conspicuous quasi-periodic fluctuations (period 25/35 min) in F region vertical plasma drift, Vz in the interval 0047/0210 IST on the night of 23/24 December, 1991 (Ap = 14, Kp < 4–). The fluctuations in F region vertical drift are found to be coherent with variations in Bz (north-south) component of interplanetary magnetic field (IMF), in geomagnetic H/X components at high-mid latitude locations both in the sunlit and dark hemispheres and near the dayside dip equator, suggestive of DP2 origin. But the polarity of the electric field fluctuations at the midnight dip equator (eastward) is the same as the dayside equator inferred from magnetic variations, contrary to what is expected of equatorial DP2. The origin of the coherent occurrence of equatorial electric field fluctuations in the DP2 range of the same sign in the day and night hemispheres is unclear and merits further investigations.  相似文献   

16.
  相似文献   

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

18.
Abstract

To model penetrative convection at the base of a stellar convection zone we consider two plane parallel, co-rotating Boussinesq layers coupled at their fluid interface. The system is such that the upper layer is unstable to convection while the lower is stable. Following the method of Kondo and Unno (1982, 1983) we calculate critical Rayleigh numbers Rc for a wide class of parameters. Here, Rc is typically much less than in the case of a single layer, although the scaling Rc～T2/3 as T → ∞ still holds, where T is the usual Taylor number. With parameters relevant to the Sun the helicity profile is discontinuous at the interface, and dominated by a large peak in a thin boundary layer beneath the convecting region. In reality the distribution is continuous, but the sharp transition associated with a rapid decline in the effective viscosity in the overshoot region is approximated by a discontinuity here. This source of helicity and its relation to an alpha effect in a mean-field dynamo is especially relevant since it is a generally held view that the overshoot region is the location of magnetic field generation in the Sun.  相似文献   

19.
Summary Using the geocentric constant GM=398 601.3 × 10 9 m 3s –2 , the known value of the angular velocity of the Earth's rotation , Stokes' constants J n (k) and S n (k) upto n=21 (zonal), n=16 (tesseral and sectorial) [2], the geocentric co-ordinates and heights above sea-level of SAO satellite stations [2], the following will be derived: the potential on the geoid Wo, the scale factor for lengths Ro=GM/Wo, the radius-vector of the surface W=Wo, the parameters of the best-fitting Earth tri-axial ellipsoid, and the components of the deflections of the vertical with respect to the geocentric rotational IAG ellipsoid (Lucerne 1967), as well as to the best-fitting geocentric tri-axial ellipsoid. Some of the differences in the structure of the gravity field over the Northern and Southern Hemispheres will be given, and the mean values of gravity over the equatorial zone, determined from the dynamics of satellite orbits, on the one hand, and from terrestrial gravity data, on the other, will be compared.Presented at the Fifteenth IUGG General Assembly, Moscow, July 30 — August 14, 1971.  相似文献   

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

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

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