FLOW FIELD IN SCOURED ZONE OF CHANNEL CONTRACTIONS

Experiments were conducted in a laboratory flume to measure the two-dimensional turbulent flow field in the scoured zone of channel contractions under a clear-water scour condition. The Acoustic Doppler Velocimeter (ADV) was used to detect the flow field at different vertical lines along the centerline of uncontracted (main channel) and contracted zones of the channel. The distributions of time-averaged velocity components, turbulent intensity, turbulent kinetic energy, and Reynolds stresses are presented in nondimensional graphical form. The bed shear stresses are computed from the measured Reynolds stresses being in threshold condition within the zone of contraction where bed was scoured. The data presented in this paper would be useful to the investigators for the development of kinematic flow model and morphological model of scour at a channel or river contraction.     

重金属迁移转化模型研究

大量重金属污染物迁移转化的现象和研究成果都表明重金属以泥沙颗粒为载体迁移转化,描述重金属在天然水体中的迁移转化必须紧紧抓住泥沙颗粒运动及重金属与泥沙之间的转化关系进行.为此在对现有的重金属迁移转化数学模型进行概括分类的基础上,根据水沙运动与污染物相互作用关系,分析了泥沙颗粒运动及重金属吸附解吸不平衡过程,并结合水沙数学模型,建立了重金属迁移转化的耦合模型.同时在模型合理性分析的基础上,对模型进行了计算分析,分析表明模型能够合理地反映重金属污染物在水体中的迁移转化过程.     

Hydro-thermal performance of multilayer capillary barriers in arid lands

In this study, a capillary barrier system was designed and tested for an arid land environment. To simulate arid land conditions of high temperature and sub-irrigation systems, the barrier was subjected to thermal and hydraulic gradients in opposite directions; to test the barrier system under these severe conditions, an experimental apparatus was designed and fabricated. The multilayer capillary barrier consisted of three layers made of silica sand, a mixture of sand and bentonite in equal portions, and a mixture of clay (25%) and aggregate (75%). Several one dimensional coupled heat and moisture tests were performed. Temperature variations along the thickness of the barrier were recorded as a function of time, and at the end of each test, the barrier was sliced into small sections, for the determination of volumetric water content as a function of distance from the heat source. The experimental results were discussed in view of the barrier's intended purpose of its ability to store moisture for long time durations. Coupled heat and moisture flow equations were developed and solved numerically via a finite difference method. Diffusivity parameters were calculated by using experimental results, a numerical model, and Powell's conjugate directions method of nonlinear optimization. The model was calibrated and the results were discussed. Good agreement between calculated and experimental results was obtained. This revised version was published online in July 2006 with corrections to the Cover Date.     

Transport of dissolved methoxy-nonafluorobutane through a saturated column

Solute transport experiments were conducted in a one-dimensional saturated column using dissolved methoxy-nonafluorobutane (HFE-7100), a Novec engineered fluid developed by the 3M Corporation, as the solute. Novec engineered fluids are considered dense non-aqueous phase liquids (DNAPLs) because they are immiscible with water and have a specific gravity greater than one. The HFE-7100 fluid is safer and environmentally friendlier than common DNAPL contaminants such as tetrachloroethylene (PCE) or trichloroethylene (TCE); thus, it is an ideal substitute DNAPL for laboratory groundwater contamination research. Three sets of solute transport experiments were conducted. The first set of experiments was conducted in a glass-bead-packed column using dissolved HFE-7100 as the solute. The second set of experiments was conducted in a sand-packed column using dissolved HFE-7100 as the solute. The third set of experiments was conducted in a sand-packed column using dissolved PCE as the solute. The dissolved HFE-7100 column breakthrough concentrations were compared with dissolved PCE breakthrough concentrations. Results show that the one-dimensional solute transport equation was successful in describing the transport behavior of dissolved HFE-7100. This study demonstrates that the HFE-7100 fluid can be used as a safer substitute DNAPL for groundwater contaminant dissolution and transport research.     

Effects of temperature and humidity upon the transport of sedimentary particles by wind

Wind tunnel simulations of aeolian transport carried out over a range in mean temperature between 32 °C and ?9 °C suggest that cold airflows support higher mass transport rates (Q) than very warm air. The magnitude of this increase is larger than expected, so that analytical and semi‐empirical models underestimate Q. Extrapolation of the results suggests that, at ?40 °C, as for example in the dry valleys of Antarctica in winter, Q may be as much as 70% higher than for the equivalent wind speed in hot deserts at air temperatures of 40 °C. Temperature‐dependent changes in air density and turbulence contribute to this result. The decreased tension of water adsorbed onto particle surfaces at low temperatures is postulated to reduce interparticle cohesion and, thus, to increase the elasticity of particle impacts on cold beds. Definition of the roles that temperature and humidity play in aeolian transport is relevant to studies of palaeoenvironmental reconstruction and extraterrestrial (or planetary) geology. Investigation of present‐day, cold climate features and of climate change effects also requires knowledge of these fundamental relations.     

Nuclear Waste Disposal Simulations: Couplex Test Cases

We give some results obtained for the Couplex test cases proposed by the ANDRA. In this paper our aim is twofold. Firstly, to compute the release of nuclides out of the repository by concentrating on the 3D near field (Couplex 2). The simulation of the transport phenomena takes into account the dissolution of the glass containers and congruent emissions of the radio-nuclides including filiation chains and some simplified chemistry. Secondly, it is to use the near field computations in order to simulate the nuclide migrations in a 2D far field (Couplex 3). Coupling in between the two simulations takes into consideration the periodicity of the disposal modules and the geometry of the repository described in Couplex 1. The mixed finite element and discontinuous Galerkin methods are used to solve the convection–diffusion equations. In order to handle the nonlinear precipitation/dissolution term, we developed a new iterative technique that combines Picard and Newton–Raphson methods.     

Modélisation du transport réactif en milieux poreux non saturés avec une méthode ELLAM en maillage variable

Modelling contaminant transfer with biological/chemical/radioactive processes needs appropriate numerical methods able to reproduce sharp concentration fronts. In this work, we develop a new Eulerian–Lagrangian Localized Adjoint Method (ELLAM) for solving the reactive transport equation with non-constant coefficients. To avoid interpolation (leading to errors), we use a moving grid to define the solution and test functions. The method is used to simulate first the infiltration of solute into a column of unsaturated porous medium and second the multispecies transport. The developed ELLAM gives accurate results without non-physical oscillations or numerical diffusion, even when using large time steps. To cite this article: A. Younes, C. R. Geoscience 336 (2004).     

Long-term bed load transport rate based on aerial-photo and ground penetrating radar surveys of fan-delta growth, Coast Mountains, British Columbia

Sequential aerial photography, sonar bathymetry, ground-penetrating radar (GPR), and sediment sampling and analysis provide the basis for calculating the volumetric and mass rate of progradation of the delta of Fitzsimmons Creek, a steep, high-energy, debris-flow-dominated channel draining about 100 km² of the southern Coast Mountains of British Columbia. Fitzsimmons Creek is typical of small mountain rivers in the region. GPR imaging is used to define the pre-depositional morphology of the receiving basin, a technique that improves the accuracy of the volumetric survey. The 52-year record (1947–1999) of progradation yielded an average annual volumetric transport rate of 1.00±0.16×10⁴ m³ year⁻¹ for bed load, corresponding to a mass transport rate of 1.60±0.28×10⁴ Mg year⁻¹. Bed load yields are consistent with those obtained in hydrogeomorphically similar basins in the region and elsewhere. Decade-based annual rates, which vary from 0.64±0.11×10⁴ to 2.85±0.38×10⁴ Mg year⁻¹, provide poor estimates of the 52-year average. Indeed, the 52-year record may also not be long enough to fully integrate the significant fluctuations in the sediment efflux from Fitzsimmons Creek. The methodology proposed in this paper can be transferred to other comparable mountain environments worldwide.     

A parameterisation for the threshold shear velocity to initiate deflation of dry and wet sediment

A new parameterisation for the threshold shear velocity to initiate deflation of dry and wet particles is presented. It is based on the balance of moments acting on particles at the instant of particle motion. The model hence includes a term for the aerodynamic forces, including the drag force, the lift force and the aerodynamic-moment force, and a term for the interparticle forces. The effect of gravitation is incorporated in both terms. Rather than using an implicit function for the effect of the aerodynamic forces as reported earlier in literature, a constant aerodynamic coefficient was introduced. From consideration of the van der Waals force between two particles, it was further shown that the effect of the interparticle cohesion force between two dry particles on the deflation threshold should be inversely proportional to the particle diameter squared. The interparticle force was further extended to include wet bonding forces. The latter were considered as the sum of capillary forces and adhesive forces. A model that expresses the capillary force as a function of particle diameter squared and the inverse of capillary potential was deduced from consideration of the well-known model of Fisher and the Young–Laplace equation. The adhesive force was assumed to be equal to tensile strength, and a function which is proportional to particle diameter squared and the inverse of the potential due to adhesive forces was derived. By combining the capillary-force model and the adhesive force model, the interparticle force due to wet bonding was simplified and written as a function of particle diameter squared and the inverse of matric potential. The latter was loglinearly related to the gravimetric moisture content, a relationship that is valid in the low-moisture content range that is important in the light of deflation of sediment by wind. By introducing a correction to force the relationship to converge to zero moisture content at oven dryness, the matric potential–moisture content relationship contained only one unknown model parameter, viz. moisture content at −1.5 MPa. Working out the model led to a rather simple parameterisation containing only three coefficients. Two parameters were incorporated in the term that applies to dry sediment and were determined by using experimental data as reported by Iversen and White [Sedimentology 29 (1982) 111]. The third parameter for the wet-sediment part of the model was determined from wind-tunnel experiments on prewetted sand and sandy loam aggregates. The model was validated using data from wind-tunnel experiments on the same but dry sediment, and on data obtained from simulations with the model of Chepil [Soil Sci. Soc. Am. Proc. 20 (1956) 288]. The experiments showed that soil aggregates should be treated as individual particles with a density equal to their bulk density. Furthermore, it was shown that the surface had to dry to a moisture content of about 75% of the moisture content at −1.5 MPa before deflation became sustained. The threshold shear velocities simulated with our model were found to be in good agreement with own observations and with simulations using Chepil's model.