首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Results are presented from both linear stability analysis and numerical simulations of three-dimensional nonlinear convection in a Boussinesq fluid in an annular channel, under experimental boundary conditions, rotating about a vertical axis uniformly heated from below. The focus is placed on the Prandtl number Pr = 7.0, representing liquid water at room temperature. The linear analysis shows that, when the aspect ratio is sufficiently small, there exists only one stationary mode that occupies the whole fluid container. When the aspect ratio is moderate or large, however, there exist three different linear solutions: (i) the outer sidewall-localized traveling wave propagating against the sense of rotation; (ii) the inner sidewall-localized traveling wave propagating in the same sense as rotation; and (iii) both the counter-traveling waves occurring simultaneously. Guided by the result of the linear stability analysis, fully three-dimensional simulations are then performed for a channel with a moderate aspect ratio. It is found that neither the prograde nor the retrograde mode is physically realizable near threshold and beyond. The dynamics of nonlinear convection in a rotating channel are chiefly characterized by the interaction between the sidewall-localized waves and the interior convection cells/rolls, producing an interesting and unusual nonlinear phenomenon. In order to compare with the classical Rayleigh–Bénard problem without vertical sidewalls, we also study linear and nonlinear convection at exactly the same parameters but in an infinitely extended layer with periodic horizontal conditions. This reveals that both the linear instability and nonlinear convection in a rotating channel are characteristically different from those in a rotating layer with periodic horizontal conditions.  相似文献   

2.
Abstract

This paper explores the properties of a two-dimensional, Boussinesq convection model with an ad hoc term in the buoyancy tendency equation that represents a positive external feedback process acting on the buoyancy fluctuations. Linear stability analyses and nonlinear integrations are presented for the case of constant heat flux boundary conditions. Although the large wavenumber modes grow the fastest from a state of rest, the nonlinear solutions progressively evolve to cells of small wavenumber. Applications to mesoscale cellular convection in the atmosphere are discussed.  相似文献   

3.
Abstract

Models of a differentially rotating compressible convection zone are calculated, considering the inertial forces in the poloidal components of the equations of motion. Two driving mechanisms have been considered: latitude dependent heat transport and anisotropic viscosity. In the former case a meridional circulation is induced initially which in turn generates differential rotation, whereas in the latter case differential rotation is directly driven by the anisotropic viscosity, and the meridional circulation is a secondary effect.

In the case of anisotropic viscosity the choice of boundary conditions has a big influence on the results: depending on whether or not the conditions of vanishing pressure perturbation are imposed at the bottom of the convection zone, one obtains differential rotation with a fast (≥ 10 ms?1) or a slow (～ 1 ms?1) circulation. In the latter case the rotation law is mainly a function of radius and the rotation rate increases inwards if the viscosity is larger in radial direction than in the horizontal directions.

The models with latitude dependent heat transport exhibit a strong dependence on the Prandtl number. For values of the Prandtl number less than 0.2 the pole-equator temperature difference and the surface velocity of the meridional circulation are compatible with observations. For sufficiently small values of the Prandtl number the convection zone becomes globally unstable like a layer of fluid for which the critical Rayleigh number is exceeded.  相似文献   

4.
Abstract

We study the bifurcation to steady two-dimensional convection with the heat flux prescribed on the fluid boundaries. The fluid is weakly non-Boussinesq on account of a slight temperature dependence of its material properties. Using expansions in the spirit of shallow water theory based on the preference for large horizontal scales in fixed flux convection, we derive an evolution equation for the horizontal structure of convective cells. In the steady state, this reduces to a simple nonlinear ordinary differential equation. When the horizontal scales of the cells exceed a certain critical size, the bifurcation to steady convection is subcritical and the degree of subcriticality increases with increasing cell size.  相似文献   

5.
Abstract

In the present paper the behavior of (internal) magneto-acoustic-gravity waves near the cusp resonance in a lossless, compressible, isothermal, stratified, electrically conducting atmosphere that is permeated by a uniform, nearly horizontal magnetic field is re-addressed (Kamp, 1989). The previously analyzed linear conversion of long acoustic-gravity waves into short magneto-acoustic waves that carry off the energy from the resonance region along the magnetic field, is re-analyzed with boundary layer techniques that are based on the smallness of the vertical component of the magnetic field. More specifically the existence of the so-called valve effect for the generated magneto-acoustic mode near the critical level is explicitly demonstrated and shown to be governed by two rivalling effects.  相似文献   

6.
Abstract

A system is considered in which electrically conducting fluid is contained between two rigid horizontal planes and bounded laterally by a circular cylinder. The fluid is permeated by a strong azimuthal magnetic field. The strength of the field increases linearly with distance from the vertical axis of the cylinder, about which the entire system rotates rapidly. An unstable temperature gradient is maintained by heating the fluid from below and cooling from above. When viscosity and inertia are neglected, an arbitrary geostrophic velocity, which is aligned with the applied azimuthal magnetic field and independent of the axial coordinate, can be superimposed on the basic axisymmetric state. In this inviscid limit, the geostrophic velocity which occurs at the onset of convection is such that the net torque on geostrophic cylinders vanishes (Taylor's condition). The mathematical problem which describes the ensuing marginal convection is nonlinear, and was discussed previously for the planar case by Soward (1986). Here new features are isolated which result from the cylindrical geometry. New asymptotic solutions are derived valid when Taylor's condition is relaxed to include viscous effects.  相似文献   

7.
Abstract

As the sun evolves, a sharp compositional peak of He 3 builds up in the core. Nuclear reactions involving He 3 are very temperature sensitive, as a result, this He 3 layer is susceptible to thermal instability. The small horizontal wavenumber g-modes have large time scales, comparable to the thermal time scale. Using a two-layer model, we find that such “shellular modes” are the most unstable. As a result of nuclear heating, these modes may be excited in the solar core in a shallow layer confined to the He 3 zone. A possible effect of such shellular convection on the solar neutrino problem is discussed. In this paper we discuss the linear theory; the nonlinear effects will be treated in a subsequent paper.  相似文献   

8.
Abstract

Experiments have been performed in a cylindrical annulus with horizontal temperature gradients imposed upon the horizontal boundaries and in which the vertical depth was smaller than the width of the annulus. Qualitative observations were made by the use of small, suspended, reflective flakes in the liquid (water).

Four basic regimes of flow were observed: (1) axisymmetric flow, (2) deep cellular convection, (3) boundary layer convective rolls, and (4) baroclinic waves. In some cases there was a mix of baroclinic and convective instabilities present. As a “mean” interior Richardson number was decreased from a value greater than unity to one less than zero, axisymmetric baroclinic instability of the Solberg type was never observed. Rather, the transition was from non-axisymmetric baroclinic waves, to a mix of baroclinic and convective instability, to irregular cellular convection.  相似文献   

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

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

11.
Abstract

The linear problem of the onset of convection in rotating spherical shells is analysed numerically in dependence on the Prandtl number. The radius ratio η=r i/r o of the inner and outer radii is generally assumed to be 0.4. But other values of η are also considered. The goal of the analysis has been the clarification of the transition between modes drifting in the retrograde azimuthal direction in the low Taylor number regime and modes traveling in the prograde direction at high Taylor numbers. It is shown that for a given value m of the azimuthal wavenumber a single mode describes the onset of convection of fluids of moderate or high Prandtl number. At low Prandtl numbers, however, three different modes for a given m may describe the onset of convection in dependence on the Taylor number. The characteristic properties of the modes are described and the singularities leading to the separation with decreasing Prandtl number are elucidated. Related results for the problem of finite amplitude convection are also reported.  相似文献   

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

13.
Abstract

The onset of convection in a cylindrical fluid annulus is analyzed in the case when the cylindrical walls are rotating differentially, a temperature gradient in the radial direction is applied, and the centrifugal force dominates over gravity. The small gap approximation is used and no-slip conditions on the cylindrical walls are assumed. It is found that over a considerable range of the parameter space either convection rolls aligned with the axis of rotation or rolls in the perpendicular (azimuthal) direction are preferred. It is shown that by a suitable redefinition of parameters, results for finite amplitude Taylor vortices and for convection rolls in the presence of shear can be applied to the present problem. Weakly nonlinear results for transverse rolls in a Couette flow indicate the possibility of subcritical bifurcation for Prandtl numbers P less than 0.82. Heat and momentum transports are derived as functions of P and the problem of interaction between transverse and longitudinal rolls is considered. The relevance of the analysis for problems of convection in planetary and stellar atmospheres is briefly discussed.  相似文献   

14.
ABSTRACT

The present study aims to link the dynamics of geophysical fluid flows with their vortical structures in physical space and to study the transition of these structures due to the control parameters. The simulations are carried in a rectangular box filled with liquid gallium for three different cases, namely, Rayleigh–Bénard convection (RBC), magnetoconvection (MC) and rotating magnetoconvection (RMC). The physical setup and material properties are similar to those considered by Aurnou and Olson in their experimental work. The simulated results are validated with theoretical results of Chandrasekhar and experimental results of Aurnou and Olson. The results are also topologically verified with the help of Euler number given by Ma and Wang. For RBC, the onset is obtained at Ra greater than 1708 and at this Ra, the symmetric rolls are orientated in/along a horizontal axis. As the value of Ra increases further, the width of the horizontal rolls starts to amplify. It is observed that these two-dimensional rolls are nothing but the cross-sections of three-dimensional (3D) cylindrical rolls with wave structures. When the vertically imposed magnetic field is added to RBC, the onset of convection is delayed due to the effect of Lorentz force on the thermal buoyancy force. The presence of 3D rectangular structures is highlighted and analysed. When the magnetically influenced rectangular box rotates about vertical axis at low rotation rates in magnetoconvection model, the onset of convection gets further delayed by magnetic field, which is in general agreement with the theoretical predictions. The critical Ra increases linearly with magnetic field intensity. Coherent thermal oscillations are detected near the onset of convection, at moderate rotation rates.  相似文献   

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

16.
Abstract

Drift rates and amplitudes of convection columns driven by centrifugal bouyancy in a cylindrical fluid annulus rotating about a vertical axis have been measured by thermistor probes. Conical top and bottom boundaries of the annular fluid region are responsible for the prograde Rossby wave like dynamics of the convection columns. A constant positive temperature difference between the outer and the inner cylindrical boundaries is generated by the circulation of thermostatically controled water. Mercury and water have been used as converting fluids. The measurements extend the earlier visual observations of Busse and Carrigan (1974) and provide quantitative data for an eventual comparison with nonlinear theories of thermal Rossby waves. The measured drift frequencies are in general agreement with linear theory. Of particular interest is the decline of the amplitude of convection with increasing Rayleigh number in a region beyond the onset of convection.  相似文献   

17.
Abstract

A series of experiments has been undertaken to investigate the onset of convection in a rapidly rotating fluid filled sphere. The boundary is subjected to a time varying temperature allowing the simulation of radial temperature profiles associated with internal heating. The system is similar to that treated theoretically by Roberts (1968), Busse (1970) and Soward (1977). It is found that Busse's modification of Roberts' linear analysis, taking into account velocity perturbations which are antisymmetric about the equatorial plane, provides a good estimate of the temperature gradient required to initiate convection. As observed in the experiments of Carrigan and Busse (1983) and predicted by linear theory, convection appears in the form of rolls or columns, aligned parallel to the rotation axis. As in earlier experiments, observed azimuthal wavenumbers are consistently smaller than predicted which we postulate to be a consequence of nonlinear effects. Owing to the presence of a centrifugally driven thermal wind, the predicted azimuthal drift of the rolls has not been observed.  相似文献   

18.
  相似文献   

19.
Non-linear Rayleigh-Bénard convection in a fluid layer is considered as a model of convection in the Earth's upper mantle. Previous studies have shown that when the temperature is held fixed at one of the boundaries of the layer, convection takes place in cells of width of the order of the layer depth or less. We investigate the effects of a different thermal boundary condition, in which the flux of heat is held fixed on both layer boundaries; then if this flux is just greater than that required for the onset of convection, motion takes place on horizontal scales much greater than the layer depth. An analytical treatment of the equations, based on an expansion in the depth-to-width ratio of the cells, shows that cells of a definite horizontal scale are the fastest growing according to linearised theory, but that these cells are unstable to ones of larger wavelength than themselves. Thus the dominant wavelength lengthens with time. The results hold whether the heat flux is generated internally of comes from beneath the layer. These results produce flow patterns similar to those found when the heat flux is much greater than the critical value. The results have important consequences for the understanding of mantle convection.  相似文献   

20.
Abstract

Chandrasekhar (1961) has summarized the stability results of Bénard convection in a rotating fluid for the cases where the boundary surfaces are both rigid and free, and for both exchange of stabilities and overstability. His analysis provides very accurate results for a limited range of Taylor number J. Bisshopp and Niiler (1965) presented an asymptotic analysis of the rigid boundary problem for exchange of stabilities which is valid for very large Taylor number. The present paper makes use of modern rotating fluid theory to develop an approximate scheme for evaluating the Rayleigh number and other parameters and variables. Known asymptotic results for the free boundary problem at large J are used and an expansion in powers of E1/6 (the Ekman number, E = 2J ?½) yields a sequence of equations and appropriate boundary conditions for the rigid boundary problem. After the algorithm for the calculation is developed, results are given for the problem to second order in the expansion parameter for the case of exchange of stabilities and to first order in the expansion parameters for the overstable case. Ekman boundary layers are important in the development as one might expect. However, an additional, diffusive boundary layer of thickness E? is necessary to provide the details of the temperature field. This boundary layer is the thermal response in the vertical direction to the horizontal spacing of the cells which is also order E?. The horizontal spacing of the cells is essentially a series of contiguous, Stewartson (1957) layers of thickness E?.  相似文献   

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

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