A positivity preserving method for simulation of steam injection for NAPL site remediation

Numerical modelling of steam injection methods for cleanup of non-aqueous phase liquid (NAPL) contamination of groundwater requires consideration of multiphase, multicomponent convective and dispersive transport. Standard techniques do not ensure that the solution of the discrete equations has positive mole fractions, for finite mesh sizes. Negative mole fractions may cause the simulation to abort due to failure of the Newton iteration. A method for alleviating this problem is described. This method ensures that the mole fractions are positive, and results in an error that is the same size as the usual finite element discretization error. Example computations are presented for cartesian and axisymmetric two-dimensional geometries.     

CLIMATIC ANALYSIS OF IMPACTS OF THE TROPICAL CONVECTION IN WESTERN PACIFIC ON WESTERN PACIFIC SUBTROPICAL HIGH’S SHORT-TERM SCALES DURING BOREAL SUMMER*

Based on the multi-year average NCEP/NCAR reanalysis data and NOAA's OLR data,the climatic characteristics of the tropical convection in tropical western Pacific and Indian Oceans as well as its relationship with western Pacific subtropical high (WPSH) is shown as follows:on short-term scales,the tropical convection that has significant influence on western Pacific high's latitudinal movement is located in the area of the South China Sea to the Philippines,which is the 2-day precursor prior to WPSH's latitudinal fluctuation,that is,WPSH is shifting to north 2 days after the tropical convection becomes more active,and vice versa.Moreover,the tropical convection has less effect on WPSH's longitudinal movement.     

半拉格朗日、半隐式欧拉方程组大气数值模式研究

介绍了一个具有较高时间积分效率的三维弹性大气数值模式。其中,控制方程为能够描述大气非静力及可压缩性的欧拉方程组,时间积分采用了高效率的“半拉格朗日、半隐式”方案。通过引用一个简化的“云物理过程”参数化方案,针对大气中深厚湿对流过程的若干问题进行了数值试验,结果表明,模式能够对该过程中的一些基本和复杂现象进行有效地模拟。试验结果还表明,动力学框架具有很好的稳定性,能够实现高效率的时间积分,它还具有较好的频散特性,保证了空间计算的精度,从而完好地描述了风暴的形态。     

Salt transport on islands in the Okavango Delta: Numerical investigations

This study uses a numerical model to investigate the groundwater flow and salt transport mechanisms below islands in the Okavango Delta. Continuous evapotranspiration on the islands results in accumulation of solutes and the formation of a saline boundary layer, which may eventually become unstable. A novel Lagrangian method is employed in this study and compared to other numerical methods. The numerical results support the geophysical observations of density fingering on Thata Island. However, the process is slow and it takes some hundreds of years until density fingering is triggered. The results are sensitive to changes of the hydraulic gradient and the evapotranspiration rate. Small changes may lead to different plume developments. Results further demonstrate that density effects may be entirely overridden by lateral flow on islands embedded in a sufficiently high regional hydraulic gradient.     

Upper mantle convection beneath northwest China and its adjacent region driven by density anomaly

Introduction The northwest of China includes Tarim, Junggar and Qaidam basins, and Kunlun, Tianshan, Altun and Qilian mountains, as well as the north part of the Tibetan Plateau. For a long time, the study of lithosphere structures and dynamics in this area has been a popular topic in geoscience, and has yielded many results. For example, TANG (1994) and LI et al (1998) suppose that the tectonic structure of Tarim Basin is various with geological periods, which changes many times betwee…     

Tangling Turbulence and Semi-Organized Structures in Convective Boundary Layers

A new mean-field theory of turbulent convection is developed based on the idea that only the small-scale region of the spectrum is considered as turbulence, whereas its large-scale part, including both regular and semi-organized motions, is treated as the mean flow. In the shear-free regime, this theory predicts the convective wind instability, which causes the formation of large-scale semi-organized motions in the form of cells. In the presence of wind shear, the theory predicts another type of instability, which causes the formation of large-scale semi-organized structures in the form of rolls and the generation of convective-shear waves propagating perpendicular to the convective rolls. The spatial characteristics of these structures, such as the minimum size of the growing perturbations and the size of perturbations with the maximum growth rate, are determined. This theory might be useful for understanding the origin of large-scale cells and rolls observed in the convective boundary layer and laboratory turbulent convection     

库尔勒地区两次强对流天气过程的形势和雷达回波对比分析

通过对2006年5月3日和7月17日库尔勒地区两次强对流天气的形势和物理量对比分析,以及库尔勒新一代天气雷达(CINRAD-CC)产品中的组合反射率(CR)、风暴相对平均径向速度(V)等产品的分析,得出库尔勒地区夏季强对流天气过程发生时的有利天气形势、物理量场分布及多谱勒天气雷达产品特征。     

从大洋中脊反思大陆漂移说——兼论地球膨裂说

本文通过讨论地幔对流的可能性以及大洋中脊的运动方式,证明了大陆漂移说对大陆漂移的源动力的解释仍存在漏洞,并提出大陆漂移源动力与地球膨裂的关系。对大陆漂移说具有补充性。     

The Role of History-Dependent Rheology in Plate Boundary Lubrication for Generating One-Sided Subduction

We have developed a two-dimensional dynamical model of asymmetric subduction integrated into the mantle convection without imposed plate velocities. In this model we consider that weak oceanic crust behaves as a lubricator on the thrust fault at the plate boundary. We introduce a rheological layer that depends on the history of the past fracture to simulate the effect of the oceanic crust. The thickness of this layer is set to be as thin as the Earth's oceanic crust. To treat 1-kilometer scale structure at the plate boundary in the 1000-kilometer scale mantle convection calculation, we introduce a new numerical method to solve the hydrodynamic equations using a couple of uniform and nonuniform grids of control volumes. Using our developed models, we have systematically investigated effects of basic rheological parameters that determine the deformation strength of the lithosphere and the oceanic crust on the development of the subducted slab, with a focus on the plate motion controlling mechanism. In our model the plate subduction is produced when the friction coefficient (0.004–0.008) of the modeled oceanic crust and the maximum strength (400 MPa) of the lithosphere are in plausible range inferred from the observations on the plate driving forces and the plate deformation, and the rheology experiments. In this range of the plate strength, yielding induces the plate bending. In this case the speed of plate motion is controlled more by viscosity layering of the underlying mantle than by the plate strength. To examine the setting of the overriding plate, we also consider the two end-member cases in which the overriding plate is fixed or freely-movable. In the case of the freely-movable overriding plate, the trench motion considerably changes the dip angle of the deep slab. Especially in the case with a shallow-angle plate boundary, retrograde slab motion occurs to generate a shallow-angle deep slab.