Numerical models of Rayleigh-Taylor instabilities superimposed upon convection
Numerical models of Rayleigh-Taylor instabilities superimposed upon convection
复制标题
对流叠加的瑞利-泰勒不稳定性数值模型
DOI:
--
复制
发表时间:
1988
期刊:
影响因子:
--
通讯作者:
H. Schmeling
中科院分区:
文献类型:
--
作者:
H. Schmeling
Gravity driven instabilities in the earth's crust are caused by inverted density distributions which may be of combined compositional and thermal (convective) origin. An order of magnitude estimate suggests that such instabilities may develop in the lower crust of ther mally anomalous regions if the effective viscosity is less than 101H-1019 Pa s. To compare the vigour of compositionally driven gravity overturns with convection instabilities, and to describe their effect on the temperature field, a nondimensional number Rt ("Rayleigh Taylor number.") is defined, which is analogous to the Rayleigh number, Ra, used for con· vection. Rt, multiplied by a factor depending on the geometrical arrangement, C, describes the diapiric heat advection relative to conduction (i.e. Rt C is a diapiric Peclet no.). Nu merical models are presented in which convective instabilities are superimposed upon com positional (Rayleigh-Taylor) instabilities. Ra-Rt-space is explored up to Ra = 30000 and Rt = 10000. Several f10w regimes are identified, in general agreement with previous results based on linear stability analysis; after one Rayleigh-Taylor overturn the f10w may become stable, it may oscillate or it may continue to convect in either one or two layers separated by an irregular oblique interface. In the case of compositional density contrasts only, Rt gives a good measure for the distortion of the temperature field or the anomalous heat f1ow: diapiric advection of heat becomes significant for Rt of a few hundred or more. Heat f10w anomalies associated with diapirism may be interpreted in terms of Rt which contains infor mation about the effective viscosity.