西昆仑山南侧甜水海湖岩芯铁变化的环境记录

阐述了青藏高原西昆仑山南侧甜水海湖岩芯沉积物中氧化铁的丰度、频率分布和富集系数的变化规律。其丰度值范围FeO0．39％-3.95％，平均值为1．25％；Fe2O30．44％－8．09％，平均值为2．60％；Fe2O3／FeO0．21－3．60，平均比值1．12。说明该彻岩芯沉积物中铁平均丰度值较低，但其变化幅度较宽且清楚，在一定程度上反映了古气候环境具有较大幅度的波动。变异系数值表明，该湖岩芯沉积物中FeO的离散程度较Fe2O3大。大约距今23万年来氧化铁的沉积地球化学特征可划分出7个环境演化阶段，并基本上可与深海氧同位索阶段对比。揭示了甜水海湖演化与青藏高原强烈隆升引起的气候环境改变的密切关系。在距今23一15万年，沉积以河流相或冲积、洪积相为主，沉积环境以氧化环境为主导，气候以温干型波动；在距今15－7．4万年，沉积以湖相为主，沉积环境以还原环境为主导，气候以温湿型为主；在距今74万年以来，湖盆沉积环境多样化，沉积物以湖相和近源冲积、洪积相交替出现。气候出现冷湿、冷干、温湿和温干多种波动模式。这些资料分析初步认为甜水海古湖大约形成于距今15万年，即青藏高原第三成湖期，这也与青藏高原在距今15万年左右的再次构造隆升有关。     

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.     

Asymptotic Analysis and Layer Decomposition for the Couplex Exercise

This document explores the possibilities of multiscale expansions and domain decomposition to solve part of the Couplex 1 exercise. We concentrate on the hydrostatic pressure and show that the slenderness of the domain and the large variations of the Darcy constants allows an analytical approximation which our test reveals to be true to relative errors smaller than 10^–3. The numerical tests are done in 2D with freefem+ and in 3D with freefem3D. Some considerations are also given for Iodine transport.     

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).     

Scattering of plane SV waves by cylindrical canyons in saturated porous medium

On the basis of Biot dynamic theory, an analytic solution of two-dimensional scattering and diffraction of plane SV waves by circular cylindrical canyons in a half space of saturated porous media is presented in this paper for the first time. The solution is obtained by employing the Fourier–Bessel series expansion technique. Parametric studies had been carried out, which includes: the angle of incidence, the frequency of the incident SV wave, the porosity of saturated porous medium and the stiffness and Poisson's ratio of the solid-skeleton. All the outcomes are useful for the seismic analysis of the surface topography conditions.     

SiO2-TiO2薄膜制备及其光催化性能研究

Sol-gel法制备SiO2-TiO2复合薄膜，进行甲醛的光催化降解实验，研究了SiO2掺杂量、薄膜焙烧温度、甲醛溶液pH值及初始浓度对光催化性能的影响。结果表明，500℃焙烧1h，掺杂质量SiO2：TiO2=12：100时，SiO2-TiO2薄膜的光催化活性最高，150min后甲醛降解率达76％，是单纯TiO2薄膜的1．3倍，该薄膜稳定性较好，经重复8次使用，最大降解率只下降4％。     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

A model for moisture and heat flow in fractured unsaturated porous media

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.     

The effective stress concept and the evaluation of changes in pore air pressure under jacketed isotropic compression tests for dry sand based on the 2‐phase porous theory

The effective stress concept for solid‐fluid 2‐phase media was revisited in this work. In particular, the effects of the compressibility of both the pore fluid and the soil particles were studied under 3 different conditions, i.e., undrained, drained, and unjacketed conditions based on a Biot‐type theory for 2‐phase porous media. It was confirmed that Terzaghi effective stress holds at the moment when soil grains are assumed to be incompressible and when the compressibility of the pore fluid is small enough compared to that of the soil skeleton. Then, isotropic compression tests for dry sand under undrained conditions were conducted within the triaxial apparatus in which the changes in the pore air pressure could be measured. The ratio of the increment in the cell pressure to the increment in the pore air pressure, m, corresponds to the inverse of the B value by Bishop and was obtained during the step loading of the cell pressure. In addition, the m values were evaluated by comparing them with theoretically obtained values based on the solid‐fluid 2‐phase mixture theory. The experimental m values were close to the theoretical values, as they were in the range of approximately 40 to 185, depending on the cell pressure. Finally, it was found that the soil material with a highly compressible pore fluid, such as air, must be analyzed with the multi‐phase porous mixture theory. However, Terzaghi effective stress is practically applicable when the compressibilities of both the soil particles and the pore fluid are small enough compared to that of the soil skeleton.