Dynamics of unsaturated poroelastic solids at finite strain

We derive the governing equations for the dynamic response of unsaturated poroelastic solids at finite strain. We obtain simplified governing equations from the complete coupled formulation by neglecting the material time derivative of the relative velocities and the advection terms of the pore fluids relative to the solid skeleton, leading to a so‐called u^s ? p^w ? p^a formulation. We impose the weak forms of the momentum and mass balance equations at the current configuration and implement the framework numerically using a mixed finite element formulation. We verify the proposed method through comparison with analytical solutions and experiments of quasi‐static processes. We use a neo‐Hookean hyperelastic constitutive model for the solid matrix and demonstrate, through numerical examples, the impact of large deformation on the dynamic response of unsaturated poroelastic solids under a variety of loading conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite element analysis of the effect of corrosion on the behavior of reinforced earth walls

A simple model for the corrosion‐induced loss of stiffness and strength of the steel strips of earth‐reinforced walls was introduced in a finite element simulation of the long‐term behavior of the wall, in which the backfill‐strips interactions are taken into account by means of a generalized homogenization procedure (called a multiphase model). The results show an initial phase of slow displacements induced by the loss of stiffness, followed after a few decades by a steep acceleration of the displacements, leading to wall failure. The influences of the parameter controlling corrosion, the backfill cohesion and the heterogeneity of the corrosion process are discussed. Results are used to discuss a strategy for reinforced earth wall surveillance. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling of groundwater flow in heterogeneous porous media by finite element method

The HySuf‐FEM code (Hydrodynamic of Subsurface Flow by Finite Element Method) is proposed in this article in order to estimate the spatial variability of the transmissivity values of the Berrechid aquifer (Morocco). The calibration of the model is based on the hydraulic head, hydraulic conductivity and recharge. Three numerical tests are used to validate the model and verify its convergence. The first test case consists in using the steady analytical solution of the Poisson equation. In the second, the model has been compared with the hydrogeological system which is characterized by an unconfined monolayer (isotropic layer) and computed by using PMWIN‐MODFLOW software. The third test case is based on the comparison between the results of HySuf‐FEM and the multiple cell balance method in the aquifer system with natural boundaries case. Good agreement between the Hydrodynamic of Subsurface Flow, the numerical tests and the spatial distribution of the thickening of the hydrogeological system is deduced from the analysis and the interpretations of hydrogeological wells. Copyright © 2011 John Wiley & Sons, Ltd.     

A simplified computational approach for prediction of transient climatic conditions in underground mine airways

Summary A simplified approach for prediction of transient climatic conditions in an underground airway os presented. The dry surface temperature of the airway is calculated using a finite difference scheme, whereas the wet surface temperature is determined following the concept of pseudo-base-temperature of air and using the value of thermal gradient at the dry surface as obtained from the finite difference solution of the temperature profile in the dry rock. The air-current temperatures are regularly evaluated after a suitable time span so as to simulate a real-case situation, and the re-evaluated air-current temperatures are used for further calculations. For increasing age of the airway the wet surface temperature obtained by the finite difference scheme, with an appropriate value of thermal diffusivity, is compared with its value determined by the pseudo-base-temperature approach.     

佳木斯市地下水水量水质模型

佳木斯市是一个以开采地下水为主要供水水源的城市.本文根据当地的水文地质条件,建立了佳木斯市的地下水的水量水质模型.用有限元法求解水量模型进行地下水资源评价.用特征有限元法求解水质模型进行地下水污染预测,该方法在求解对流一弥散方程时能有效地消除数值弥散和数值振荡,精度较高.为使该方法付诸实施,文中也提出了当利用特征有限元法求解时确定运动水质点所属单元的方法.     

用解析有限元法嵌入地下水管理模型求最优解

为了弥补当前地下水资源优化管理模型中通用的两种构模方法的缺点,即嵌入法的代数方程组过于庞大、上机困难,和响应矩阵法不保证全局最优,提出了用解析有限元法建立整体非稳定(即带有时间变量)抽水响应矩阵作为约束条件嵌入线性规划模型.这一方法由于具备了空间上表达任一节点在任一时刻的地下水运动规律的能力,克服了上述构模的两大缺点,其推广应用将有利于提高地下水资源的优化开发和管理.     

无结构网格上二维浅水流动的数值模拟

提出在无结构网格上建立有限体积高性能格式族的统一框架,通过引入跨单元界面法向数值通量的逆风分解,将一维Osher、TVD两种通量分裂格式自然地推广至二维浅水方程组。给出了各种情况下浅水方程组有限体积法边界处理的计算公式。最后,利用该格式对穹包溢流、陆地动边界和河口潮流三个问题进行计算。结果表明,这类守恒格式具有高精度、无振荡性以及处理复杂流态过渡、自动捕俘间断和模拟陆地动边界的功能。     

拟合法15°波动方程有限差分串联偏移

本文把最佳拟合算法（三点和五点差分格式）及串联偏移方法应用在15°波动方程偏移中。由于拟合法的差分因子随着采样间隔、步长、速度等参数的变化而变化,因此与原有的波动方程差分偏移法相比,串联偏移方法能有效的使大倾角归位,克服频散,速度及步长可以灵活变化,运算效率高等优点。