Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

A THREE-DIMENSIONAL ELASTIC NESTED-GRID MESO-(β-γ)SCALE ATMOSPHERIC MODEL——PARTI:MODEL DESCRIPTION

A three-dimensional elastic nonhydrostatic mesoscale(β-γ)model with nested-grid is presented.It uses a set of fullequations in terrain-following coordinates as its basic dynamic frame,which is solved with a time-splitting algorithmfor acoustic and gravity waves.The model physical parameterization includes a K-theory subgrid eddy mixing for cloudand free atmosphere,a bulk planetary boundary layer parameterization,and three types of sofisticated cloudmicrophysics schemes with double-parameters for hail-bearing clouds,warm clouds and snowing clouds respectively.The model is designed to be used flexibly for simulations of a variety of meso-and small-scale atmospheric processes,and can be improved as a regional and local operational NWP system in future.     

黄金的成色

黄金的成色并不是一成不变的，它成色的高低是与成矿温度、成矿阶段、深度、时代、变质程度、水溶液及风化作用等有一定关系，我们从这些关系中可以寻找到我们找矿的捷径和依据。     

Uncoupling of coupled flows in soil—a finite element method

Coupled flow of water, chemicals, heat and electrical potential in soil are of significance in a variety of circumstances. The problem is characterized by the coupling between different flows, i.e. a flow of one type driven by gradients of other types, and by the dual nature of certain flows, i.e. combined convection and conduction. Effective numerical solutions to the problem are challenged due to the coupling and the dual nature. In this paper, we first present a general expression that can be used to represent various types of coupled flows in soil. A finite element method is then proposed to solve the generalized coupled flows of convection-conduction pattern. The unknown vector is first decomposed into two parts, a convective part forming a hyperbolic system and a conductive part forming a parabolic system. At each time step, the hyperbolic system is solved analytically to give an initial solution. To solve the multi-dimensional hyperbolic system, we assume that a common eigenspace exists for the coefficient matrices, so that the system can be uncoupled by transforming the unknown vector to the common eigenspace. The uncoupled system is solved by the method of characteristics. Using the solution of the hyperbolic system as the initial condition, we then solve the parabolic system by a Galerkin finite element method for space discretization and a finite difference scheme for time stepping. The proposed technique can be used for solving multi-dimensional, transient, coupled or simultaneous flows of convection-conduction type. Application to a flow example shows that the technique indeed exhibits optimality in convergence and in stability.     

扬子地台钾镁煌斑岩 Nd、Sr、Pb 同位素特征

扬子地台周缘钾镁煌斑岩 Nd、Sr同位素特征介于南非Ⅰ、Ⅱ型金伯利岩和西澳与北美钾镁煌斑岩两种变化趋势之间，不同岩区间同位素组成差异反映地台周边构造性质不同，Pb同位素特征表明晚太古代地幔岩石圈具U亏损特征。     

草簇状石墨烯片限域Pd纳米粒子及其催化应用(英文)

通过氨基离子液体改性石墨烯,并将其固载于堇青石表面,作为负载型Pd催化剂的载体.所制备的Pd催化剂经加氢老化后,表面石墨烯呈草簇状结构,将Pd纳米粒子限域于片层内,有效防止了Pd的流失和团聚.在重要的工业反应对羧基苯甲醛(4-CBA)加氢中,此结构催化剂与传统的钯碳催化剂相比,表现出很好的稳定性     

Typical basin-fill sequences and basin migration in Yanshan, North China--Response to Mesozoic tectonic transition

Basin-fill sequences of Mesozoic typical basins in the Yanshan area, North China may be divided into four phases, reflecting lithosphere tectonic evolution from flexure (T3), flexure with weak rifting (J1+2), tectonic transition (J3), and rifting (K). Except the first phase, the other three phases all start with lava and volcaniclastic rocks, and end with thick coarse clastic rocks and/or conglomerates, showing cyclic basin development rather than simple cyclic rift mechanism and disciplinary basin-stress change from extension to compression in each phase. Prototype basin analysis, based on basin-fill sequences, paleocurrent distribution and depositional systems, shows that single basin-strike and structural-line direction controlling basin development had evidently changed from east-west to northeast in Late Jurassic in the Yanshan area, although basin group still occurred in east-west zonal distribution. Till Early Cretaceous, main structural-line strike controlling basins just turned to northeast by north in the studied area.     

Soil respiration in tropical seasonal rain forest in Xishuangbanna, SW China

With the static opaque chamber and gas chromatography technique, from January 2003 to January 2004 soil respiration was investigated in a tropical seasonal rain forest in Xishuangbanna, SW China. In this study three treatments were applied, each with three replicates: A (bare soil), B (soil+litter), and C (soil+litter+seedling). The results showed that soil respiration varied seasonally, low from December 2003 to February 2004, and high from June to July 2004. The annual average values of CO2 efflux from soil respiration differed among the treatments at 1% level, with the rank of C (14642 mgCO2· m-2. h-1)＞B (12807 mgCO2· m-2. h-1)＞A (9532 mgCO2· m-2. h-1). Diurnal variation in soil respiration was not apparent due to little diurnal temperate change in Xishuangbanna. There was a parabola relationship between soil respiration and soil moisture at 1% level. Soil respiration rates were higher when soil moisture ranged from 35% to 45%. There was an exponential relationship between soil respiration and soil temperature (at a depth of 5cm in mineral soil) at 1% level. The calculated Q1o values in this study,ranging from 2.03 to 2.36, were very near to those of tropical soil reported. The CO2 efflux in 2003was 5.34 kgCO2· m-2. a-1 from soil plus litter plus seedling, of them 3.48 kgCO2· m-2. a-1 from soil (accounting for 62.5%), 1.19 kgCO2· m-2. a-1 from litter (22.3%) and 0.67 kgCO2·m-2. a-1 from seedling (12.5%).     

Delayed geochemical hazard: Concept, digital model and case study

Delayed Geochemical Hazard (DGH briefly) presents the whole process of a kind of serious ecological and environmental hazard caused by sudden reactivation and sharp release of long-term accumulated pollutant from stable species to active ones in soil or sediment system due to the change of physical-chemical conditions (such as temperature, pH, Eh, moisture, the concentrations of organic matters, etc.) or the decrease of environment capacity. The characteristics of DGH are discussed. The process of a typical DGH can be expressed as a nonlinear polynomial. The points where the derivative functions of the first and second orders of the polynomial reach zero, minimum and maximum are keys for risk assessment and harzard pridication.The process and mechanism of the hazard is due to the transform of pollutant among different species principally. The concepts of "total releasable content of pollutant", TRCP, and "total concentration of active specie", TCAS, are necessarily defined to describe the mechanism of DGH. The possibility of the temporal and spatial propagation is discussed. Case study shows that there exists a transform mechanism of "gradual release" and "chain reaction" among the species of the exchangeable and the bounds to carbonate, iron and manganese oxides and organic matter, thus causing the delayed geochemical hazard.