MODEL FOR PREDICTING TRANSPORT OF PESTICIDES AND THEIR METABOLITES IN THE UNSATURATED ZONE1

ABSTRACT: A finite element model based on Galerkin's upstream weighted residual technique was developed to predict the simultaneous convective-dispersion transport and transformations of pesticides and their metabolites in the unsaturated zone. Transformations of the parent compound and its metabolites were assumed to be first-order reactions for oxidation and hydrolysis, while adsorption of the pesticide species (parent compound and metabolites) to the soil components was assumed to be represented by a linear equilibrium (Freundlich type) isotherm. Volatilization and plant root uptake of pesticides in the solution phase were neglected in the analysis. The proposed model was used to simulate the transport and transformation of aldicarb and its metabolites, aldicarb sulfoxide and aldicarb sulfone, in the soil profile. Several examples are used to demonstrate the accuracy, validity, and applicability of the proposed model. Simulated results indicate that the proposed model can potentially be used to estimate the mass flux of water, and pesticide and pesticide metabolite concentrations in the subsurface environment. However, further verification of the model against actual field data is needed to fully demonstrate the model's potential.     

黄土滑坡稳定性评价的反演分析

本文根据黄土滑坡的变形监测结果，同时考虑了边坡的初期应变分布状况，通过有限元反分析法，反推出边坡地基的力学参数，为黄土滑坡的稳定性评价和预报提供了可靠的依据。     

NUMERICAL MODELING OF SALTWATER INTRUSION IN THE NORTHERN GUAM LENS1

ABSTRACT This paper deals with the application of a two-dimensional, saltwater intrusion model to the aquifer in Northern Guam. The model used finite element theory and the Galerkin, weighted-residual technique as its basis. The Northern Guam lens was discretized into 299 linear, triangular elements and 189 nodes. The model was calibrated using 1978 hydrologic data. The output of the model was compared with measured water levels in six observation wells. The calibrated values of permeability and porosity were then used to verify the model using 1979 data. A calibrated and verified model can be used to make an infinite variety of management and planning studies. In this study, three applications are provided that would be considered typical management runs. Steady state runs were made to compare the four conditions of no pumping, 1978 pumping levels, twice 1978 pumping levels, and five times 1978 pumping levels. The water levels due to these conditions are shown in plan and in cross sections of the aquifer. The effect of zero recharge to the aquifer is next demonstrated for the pumping levels existing during 1978. The final run shows how long the aquifer takes to reach steady state when the pumping rate is increased from the 1978 pumping level to twice that value. The program can be used for numerous other studies for management and planning purposes.     

EFFECTS OF TIME STEP SIZE IN IMPLICIT DYNAMIC ROUTING1

ABSTRACT The effects of the size of the Δt time step used in the integration of the implicit difference equations of unsteady open-channel flow are determined for numerous typical hydrographs with durations in the order of days or even weeks. Truncation errors related to the size of the Δt time step cause a numerical distortion (dispersion and attenuation) of the computed transient. The magnitude of the distortion is related directly to the size of the time step, the length of channel reach, and the channel resistance and inversely to the time of rise of the hydrograph. The type of finite difference expression which replaces spatial derivatives and non-derivative terms in the partial differential equations of unsteady flow has an important influence on the magnitude of the numerical distortion, as well as the numerical stability of the implicit difference equations. Time step sizes in the range of 3 to 6 hrs generally tend to minimize the combination of required computation time and numerical distortion of transients having a time of rise of the order of several days.     

考虑纵横弯曲的空间杆系有限元计算

在普通刚架有限元计算模型的基础上 ,推导出了考虑空间纵横弯曲有限元的单元刚度矩阵和形函数 ,单元刚度矩阵和形函数考虑了轴力的影响。并以此编制了“kjgj”有限元软件 ,软件考虑了轴力引起的二阶效应。用平面和空间算例进行验证 ,单根杆件只剖分一个单元计算得到的精度可以达到商业软件剖分多个单元达到的精度。     

基于优势隐含周期分析的气象灾害预报

建立了优势因子和优势隐含周期的概念,并提出了多因子逐步回归周期分析的新方法。该方法既能增长预报期,又能充分利用物理因子的重要信息,同时在计算过程中还能挑选出起显著作用的优势因子和优势隐含周期,对预报气象要素的变化趋势效果较好,有一定的实用价值,是气象灾害预报的一条比较有效的途径。     

A FINITE ELEMENT MODEL OF CONTAMINANT MOVEMENT IN GROUNDWATER1

ABSTRACT. Transient, two-dimensional solutions are developed which describe the movement and distribution of a conservative substance in a stream-aquifer system. The solutions are obtained by solving sequentially the groundwater flow and mass transport equations. A variational approach in conjunction with the finite element method is used to solve the groundwater flow equation. Galerkin's approach coupled with the finite element method is used to solve the mass transport equation. Linear approximated triangular elements and a centered scheme of numerical integration are employed to calculate the hydraulic head distribution and the concentration of solute in the flow region. The linear approximation used to define the concentration function within each element is not appropriate for cases involving steep concentration gradients. For such cases, higher order approximations are necessary to assure the continuity of gradients across interelemental boundaries. Numerical examples that illustrate the applicability of the model are presented.     

天津市夏季PM10浓度及组成元素垂直分布特征

通过天津市气象局大气边界层观测站高255 m的气象铁塔,获得了10m、40m、120m和220 m 4个高度共80个PM10样品,分析了PM10的日均质量浓度及其元素的垂直分布特征.结果表明,4个高度处PM10质量浓度时序变化基本一致;PM10质量浓度及其组成元素总量随高度增加递减;10 m、40m和120m3个高度处的PM10质量浓度时间序列分布的相关性较好,组分分歧系数较小,而220 m处的PM10质量浓度时间序列分布与其他3个高度处的相关性较差,组分分歧系数较大.样品中元素质量浓度由高到低依次是Si、Al、Ca、Fe、Na、Mg、K、Zn、Ti、Cu、Mn、Pb、Cr、Ni、Co、V、As、Cd和Hg.其中,Si、Al、Na、K质量分数随高度变化不大;Ca、Fe、Mg质量分数随高度增加递减明显;Zn、Cu、Mn、Pb、Cr在10m、40m和120 m3个高度处质量分数接近,而在220m处明显降低;Ti元素在4个高度处的质量浓度接近,质量分数总体上随高度增加递增;Ni、Co、V、As、Cd、Hg质量分数随高度的变化不同,其中与燃煤有密切关系的Hg在120 m处质量分数明显高于其他高度处,而在220m处最低.Al、Fe、Na、K等元素相对富集因子小于或接近1,其质量浓度与PM10质量浓度时序分布的相关性较好,其中Al、Na相关系数随高度增加显著递减;Zn,Cu、Pb、Cr、As、Ni的相对富集因子显著大于1,其质量浓度与PM10质量浓度时序分布的相关性较差,相关系数随高度增加变化不大.     

Development and Validation of the Total HUman Model for Safety (THUMS) Toward Further Understanding of Occupant Injury Mechanisms in Precrash and During Crash

Objective: Active safety devices such as automatic emergency brake (AEB) and precrash seat belt have the potential to accomplish further reduction in the number of the fatalities due to automotive accidents. However, their effectiveness should be investigated by more accurate estimations of their interaction with human bodies. Computational human body models are suitable for investigation, especially considering muscular tone effects on occupant motions and injury outcomes. However, the conventional modeling approaches such as multibody models and detailed finite element (FE) models have advantages and disadvantages in computational costs and injury predictions considering muscular tone effects. The objective of this study is to develop and validate a human body FE model with whole body muscles, which can be used for the detailed investigation of interaction between human bodies and vehicular structures including some safety devices precrash and during a crash with relatively low computational costs.

Methods: In this study, we developed a human body FE model called THUMS (Total HUman Model for Safety) with a body size of 50th percentile adult male (AM50) and a sitting posture. The model has anatomical structures of bones, ligaments, muscles, brain, and internal organs. The total number of elements is 281,260, which would realize relatively low computational costs. Deformable material models were assigned to all body parts. The muscle–tendon complexes were modeled by truss elements with Hill-type muscle material and seat belt elements with tension-only material. The THUMS was validated against 35 series of cadaver or volunteer test data on frontal, lateral, and rear impacts. Model validations for 15 series of cadaver test data associated with frontal impacts are presented in this article. The THUMS with a vehicle sled model was applied to investigate effects of muscle activations on occupant kinematics and injury outcomes in specific frontal impact situations with AEB.

Results and Conclusions: In the validations using 5 series of cadaver test data, force–time curves predicted by the THUMS were quantitatively evaluated using correlation and analysis (CORA), which showed good or acceptable agreement with cadaver test data in most cases. The investigation of muscular effects showed that muscle activation levels and timing had significant effects on occupant kinematics and injury outcomes. Although further studies on accident injury reconstruction are needed, the THUMS has the potential for predictions of occupant kinematics and injury outcomes considering muscular tone effects with relatively low computational costs.