Kinematic wave solutions for pollutant transport over an infiltrating plane with finite‐period mixing and mixing zone

Kinematic wave solutions are derived for transport of a conservative non‐point‐source pollutant during a rainfall‐runoff event over an infiltrating plane for two cases: (i) finite‐period mixing and (ii) soil‐mixing zone. Rainfall is assumed to be steady, uniform and finite in duration, and it is assumed to have zero concentration of pollutants. Infiltration is assumed constant in time and space. Prior to the start of rainfall, the pollutant is distributed uniformly over the plane. In the first case, when rainfall occurs, the mixing of pollutant in the runoff water occurs in a finite period of time. In the second case, the chemical concentration is assumed to be a linearly decreasing function of rainfall intensity and overland flow. The solute concentration and discharge are found to depend on the flow characteristics as well as the solute concentration parameters. The characteristics of solute concentration and discharge graphs seem to be similar to those reported in the literature and observed in laboratory experiments. Copyright © 2002 John Wiley & Sons, Ltd.     

A Generalized Backward Euler algorithm for the numerical integration of an isotropic hardening elastoplastic model for mechanical and chemical degradation of bonded geomaterials

An extended version of the classical Generalized Backward Euler (GBE) algorithm is proposed for the numerical integration of a three‐invariant isotropic‐hardening elastoplastic model for cemented soils or weak rocks undergoing mechanical and non‐mechanical degradation processes. The restriction to isotropy allows to formulate the return mapping algorithm in the space of principal elastic strains. In this way, an efficient and robust integration scheme is developed which can be applied to relatively complex yield surface and plastic potential functions. Moreover, the proposed algorithm can be linearized in closed form, thus allowing for quadratic convergence in the global Newton iteration. A series of numerical experiments are performed to illustrate the accuracy and convergence properties of the algorithm. Selected results from a finite element analysis of a circular footing on a soft rock layer undergoing chemical weathering are then presented to illustrate the algorithm performance at the boundary value problem level. Copyright © 2002 John Wiley & Sons, Ltd.     

A relook at NEH‐4 curve number data and antecedent moisture condition criteria

This paper investigates the variation of the popular curve number (CN) values given in the National Engineering Hand Book–Section 4 (NEH‐4) of the Soil Conservation Service (SCS) with antecedent moisture condition (AMC) and soil type. Using the volumetric concept, involving soil, water, and air, a significant condensation of the NEH‐4 tables is achieved. This leads to a procedure for determination of CN for gauged as well as ungauged watersheds. The rainfall‐runoff events derived from daily data of four Indian watersheds exhibited a power relation between the potential maximum retention or CN and the 5‐day antecedent rainfall amount. Including this power relation, the SCS‐CN method was modified. This modification also eliminates the problem of sudden jumps from one AMC level to the other. The runoff values predicted using the modified method and the existing method utilizing the NEH‐4 AMC criteria yielded similar results. Copyright © 2006 John Wiley & Sons, Ltd.     

A bimodal four‐parameter lognormal linear model of soil water repellency persistence

Soil water repellency may be characterized in terms of the delayed infiltration time of a water droplet resting on the soil surface, which is, water drop penetration time (WDPT), or repellency persistence. Such repellency persistence varies nonlinearly with soil water content (θ_g), although no models have been proposed to reproduce the variation of WDPT with θ_g in soils. Dynamic factor analysis (DFA) is used to identify two common patterns of unexplained variability in a scattered dataset of WDPT versus θ_g measurements. A four‐parameter lognormal distribution was fitted to both common patterns obtained by DFA, and these were combined additively in a weighted multiple linear bimodal model. We show how such an empirical model is capable of reproducing a large variety of WDPT versus θ_g curve shapes (N = 80) both within a wide range of measured WDPTs (0–17 000 s) and for samples with organic matter content ranging from 21·7 to 80·6 g (100 g)^?1. Copyright © 2009 John Wiley & Sons, Ltd.     

基岩弯曲河段洪水水流结构的试验研究

低坡度基岩弯曲河流在地质构造控制区域广泛存在（床面坡度小于5‰）,洪水对基岩弯曲河流的河床淤积与侵蚀具有较大影响,但以往研究对基岩弯曲河流的洪水动力结构认识不足。通过几何概化基岩弯曲河段,考虑Froude相似与边壁粗糙,建立基岩弯曲河道概化模型,分析洪水下的弯道水面线分布、时均流场与湍流结构特性。结果表明:在洪水流量下,弯顶上游出现最小水面横比降、凸岸水流分离、凹岸双环流发育、流速下潜且二次环流强度达到最大,在弯顶下游水面横比降达到最大并出现环流分裂;床面切应力分布于凸岸水流分离以及弯顶上游中心区域,横向动量输移集中于弯顶上游。试验结果为基岩弯曲河道中的床面侵蚀与沿程淤积提供了水动力方面的解释。     

沉积盆地走滑构造的力学机制和变形特征

走滑构造是地壳上最重要也是最复杂的构造之一,尽管经历了长期研究历史,但在其某些特征方面还存在着不确定认识。笔者等基于力学和脆性破裂理论以及安德森断层模式,分析了走滑构造的成因类型及相关构造的几何学特点。在总结走滑构造基本特征的基础上,说明了沉积盆地中大型走滑构造是基底先存断层或薄弱面复活形成,局部有新生型的,提出了准摩擦滑动的概念。明确了直扭、压扭和张扭的含义,并将亚简单剪切应用于分析断层活动。细化了各种构造的发育顺序,讨论了不同性质构造的形成条件,强调了构造几何关系与走滑带应力性质的关系。归纳了花状构造形成的构造部位,明确了走滑构造带结构及其主控因素,并用于地表和地下走滑构造的识别。     

中国地质学会公布第三届科普奖评选结果

为推动中国地质学会地学科普研学工作的开展,根据 《中国地质学会精品地学研学路线、课程评选办法(试行)》, 学会组织开展了第二批精品地学研学路线、第一批精品地学 研学课程的申报、推荐和评审工作。     

Inverse Analysis of Deep Excavation Using Differential Evolution Algorithm

This paper presents the applications of the differential evolution (DE) algorithm in back analysis of soil parameters for deep excavation problems. A computer code, named Python‐based DE, is developed and incorporated into the commercial finite element software ABAQUS, with a parallel computing technique to run an FE analysis for all trail vectors of one generation in DE in multiple cores of a cluster, which dramatically reduces the computational time. A synthetic case and a well‐instrumented real case, that is, the Taipei National Enterprise Center (TNEC) project, are used to demonstrate the capability of the proposed back‐analysis procedure. Results show that multiple soil parameters are well identified by back analysis using a DE optimization algorithm for highly nonlinear problems. For the synthetic excavation case, the back‐analyzed parameters are basically identical to the input parameters that are used to generate synthetic response of wall deflection. For the TNEC case with a total of nine parameters to be back analyzed, the relative errors of wall deflection for the last three stages are 2.2, 1.1, and 1.0%, respectively. Robustness of the back‐estimated parameters is further illustrated by a forward prediction. The wall deflection in the subsequent stages can be satisfactorily predicted using the back‐analyzed soil parameters at early stages. Copyright © 2014 John Wiley & Sons, Ltd.