透射边界高频耦合失稳机理及稳定实现——SH波动

就大型近场波动的高效数值模拟而言,稳定实现高阶人工边界是一个尚未圆满解决的问题.本文针对使用多次透射公式的SH波动集中质量有限元模拟,依据GKS定理的群速度解释,进一步阐明了人工边界与内域离散格式耦合所导致高频失稳的机理,即两者支持群速度指向内域的外行高频平面谐波,波动能量自发地从人工边界进入內域,从而导致失稳,而这类谐波是由集中质量有限元离散引入的.本文提出了消除此种耦合失稳的一种方法:通过修改有限元刚度阵来改变内域离散格式,并保证修改格式的精度不低于原有格式的精度.理论分析和数值实验表明此法能稳定实现透射边界.本文研究结果具有推广应用前景.     

Instability of flows near the onset of convection in a rotating layer with stress-free horizontal boundaries

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 ε (ε² 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 α².     

Reducing the non-axisymmetry of a planetary dynamo and an application to saturn

Abstract

If a conducting fluid shell is undergoing spin-axisymmetric differential rotation and overlies the dynamo generating region of a planet then it is capable of greatly reducing the non-spin-axisymmetric components of the generated field, provided the appropriate magnetic Reynolds number is large. The model, closely related to the electromagnetic skin effect, is quantified and applied to Saturn. The observed small dipole tilt (～ 1°) of Saturn's magnetic field can be explained because of the presence of a stably stratified conducting layer overlying the dynamo region. This layer is a predicted consequence of the thermal evolution, arises because of the limited solubility of helium in metallic hydrogen (Stevenson, 1980), and appears to be required by the Voyager infrared observations indicating depletion of helium from Saturn's atmosphere. The much larger dipole tilt angles of Jupiter and the Earth indicate the absence of any such stable, differentially rotating layer with a large magnetic Reynolds number.     

FRICTIONAL SLIDING OF PLAGIOCLASE GOUGE UNDER LOWER-CRUST TEMPERATURE AND RELATIVELY LOW EFFECTIVE NORMAL STRESS

The discovery of tremors on the lower crust portion of the San Andreas Fault has attracted more attention on the mechanical properties of the lower crust in recent years, and some experimental studies have been carried out to understand the mechanical behavior. Previous experiments under effective normal stresses of 200MPa have shown that pyroxene and plagioclase mineral separated from the gabbro and their mixtures all show velocity weakening in the lower-crust temperature range, which results in unstable slip when frictional sliding is the dominant deformation mechanism. This work is to examine whether the velocity-weakening behavior of plagioclase gouge also applies to relatively lower effective normal stress. Our experiments were performed under effective normal stress of about 100MPa, with a constant confining pressure control, with pore pressure of 30MPa and temperature of 100℃ to 600℃. We found that the frictional sliding of plagioclase are basically the same with the previous results obtained under effective normal stress of 200MPa, both of which show velocity weakening over the entire temperature range. The only difference is the out-of-trend drop of constitutive parameter a at 600℃ for the lower effective normal stress of 100MPa. It is thus concluded that reducing the effective normal stress has little effect on the sliding stability of plagioclase, and the previous conclusion made for mechanical behavior of the lower crust that unstable slips are possible therein also applies to the lower effective normal stress of 100MPa.     

Delayed baroclinic response of the Antarctic circumpolar current to surface wind stress

The Antarctic Circumpolar Current (ACC) responds to the surface windstress via two processes, i.e., the instant barotropic process and the delayed baroclinic process. This study focuses on the baroclinic instability mechanism in ACC, which was less reported in the literatures. Results show that the strengthening of surface zonal windstress causes the enhanced tilting of the isopycnal surface, leading to more intense baroclinic instability. Simultaneously, the mesoscale eddies resulting from the baro- clinic instability facilitate the transformation of mean potential energy to eddy energy, which causes the remarkable decrease of the ACC volume transport with the 2-year lag time. This delayed negative cor- relation between the ACC transport and the zonal windstress may account for the steadiness of the ACC transport during last two decades.     

大功率高频电波加热电离层的非线性物理过程

地基大功率电波加热电离层是通过地基大功率短波发射机向电离层发射无线电波,通过波-粒和波-波的相互作用将无线电波的能量注入电离层.通过这种有目的可操控的方式改变电离层电子密度和温度的分布,可以深入研究电离层中等离子体能量和物质的非线性演化过程,特别是电离层电子的非平衡态分布和加速问题.本文通过对电离层加热中几个比较重要物理过程的评述,对过去20年来我国研究学者在这一研究方向上取得的重要进展进行了介绍.

     

地球流体力学的研究与进展

简要介绍中国科学院大气物理研究所七十多年来在理论与计算地球流体力学方面的若干研究及其新的进展.在理论地球流体力学方面,介绍了长波动力学及线性稳定性问题、弱非线性理论及行星波动力学以及用Arnold方法(能量-Casimir方法)研究大气和海洋中各种流体运动的非线性稳定性问题的成果.此外,对扰动演变、扰动和基流相互作用及热带大气动力学中的第二类不稳定条件(CISK)也作了简要的介绍.在计算地球流体力学方面,主要内容包括:用物理观点和数学分析相结合的方法阐述了造成计算紊乱和计算不稳定的机理,论证计算稳定性、算     

高地应力围岩分区碎裂化的时间效应分析和相关参数研究

本文阐述了高地应力围岩分区碎裂化现象的形成需要经过一定的时间,但这段时间不会很长;采用蠕变理论对岩石分区碎裂化的时间效应进行了分析;推导了能够描述加速蠕变阶段的流变模型-改进的西原模型的本构方程和蠕变方程,在此基础上求出了岩石分区碎裂化发生时破裂带半径(破碎带距巷道中心的距离)的公式。     

Baroclinic energetics and zonal plane distribution of monsoon disturbances

The adiabatic, quasi-geostraphic, 25-layer, numerical, linear model with Ekman boundary layer friction is utilised to perform the baroclinic stability analysis of the mean monsoon zonal wind profile. It is shown thec _i is a function of the resultant wavenumber alone. This relation is able to explain the effects of the lateral walls on the unstable waves.The energetics and zonal plane distribution of the short and long preferred viscous waves are computed. The upward motion of the short wave together with the warm (cold) core lies to the west of the surface trough position above (below) 850 mb. Further, it is shown that the main source of kinetic energy for the wave lies in the middle layer (850–700 mb) which is transported to the lower and upper layers. Computed is found to be in good agreement with observed values.     

On the physical model of earthquake precursor fields and the mechanism of precursors＇timespace distribution──origin and evidences of the strong body earthquake──generating model

Ｏｎｔｈｅｐｈｙｓｉｃａｌｍｏｄｅｌｏｆｅａｒｔｈｑｕａｋｅｐｒｅｃｕｒｓｏｒｆｉｅｌｄｓａｎｄｔｈｅｍｅｃｈａｎｉｓｍｏｆｐｒｅｃｕｒｓｏｒｓ＇ｔｉｍｅｓｐａｃｅｄｉｓｔｒｉｂｕｔｉｏｎ──ｏｒｉｇｉｎａｎｄｅｖｉｄｅｎｃｅｓｏｆ　ｔｈｅｓｔｒｏｎ...