Analysis of seabed instability using element free Galerkin method

Wave-induced seabed instability, either momentary liquefaction or shear failure, is an important topic in ocean and coastal engineering. Many factors, such as seabed properties and wave parameters, affect the seabed instability. A non-dimensional parameter is proposed in this paper to evaluate the occurrence of momentary liquefaction. This parameter includes the properties of the soil and the wave. The determination of the wave-induced liquefaction depth is also suggested based on this non-dimensional parameter. As an example, a two-dimensional seabed with finite thickness is numerically treated with the EFGM meshless method developed early for wave-induced seabed responses. Parametric study is carried out to investigate the effect of wavelength, compressibility of pore fluid, permeability and stiffness of porous media, and variable stiffness with depth on the seabed response with three criteria for liquefaction. It is found that this non-dimensional parameter is a good index for identifying the momentary liquefaction qualitatively, and the criterion of liquefaction with seepage force can be used to predict the deepest liquefaction depth.     

Estimation of Land Subsidence Based on Groundwater Flow Model

This article presents an approach for estimating land subsidence due to withdrawal of groundwater. The proposed method calculates the groundwater seepage in 3-D-condition and calculates the land subsidence one-dimensionally. The governing equation on groundwater seepage is based on the three-dimensional mass conservation law and the principle of effective stress. The land subsidence calculation method is derived based on the following assumptions: (1) displacements occur only in the vertical direction, and (2) in vertical direction the total stresses do not change. The governing equation is solved by numerical method, i.e., finite element method (FEM) in spatial discretization and finite difference method (FDM) in time series discretization. In FEM Galerkin method is adopted and in FDM, lumped matrix method is employed. The proposed method is calibrated via analyzing 1-D consolidation problem and the results are compared with those from Terzaghi's one-dimensional consolidation theory and oedometer test. The proposed method is employed to analyze the consolidation of a soft layer due to withdrawal of groundwater from an aquifer under it. Moreover, this method is also applied to a field case of land subsidence due to groundwater pumping in a gas production field in Japan. The analytical results are compared with the field observed data. The results show that this approach simulates the field case well.     

Topographic Effect of Izu Ridge on the Distribution of the North Pacific Intermediate Water South of Japan

The topographic effect of the Izu Ridge on the horizontal distribution of the North Pacific Intermediate Water (NPIW) south of Japan has been studied using observational data obtained by the Seisui-Maru of Mie University (Mie Univ. data) and those compiled by Japan Oceanographic Data Center (JODC data). Both data sets show that water of salinity less than 34.1 psu on potential density () surface of 26.8 is confined to the eastern side of the Izu Ridge, while water of salinity less than 34.2 psu is confined to the southern area over the Izu Ridge at a depth greater than 2000 m and to the southeastern area in the Shikoku Basin. It is also shown by T-S analysis of Mie Univ. data over the Izu Ridge that water of salinity less than 34.2 psu dominates south of 30°N, where the depth of the Izu Ridge is deeper than 2000 m and NPIW can intrude westward over the Izu Ridge. JODC data reveal that relatively large standard deviations of the salinity on surface of 26.7, 26.8 and 26.9 are detected along the mean current path of the Kuroshio and the Kuroshio Extension. Almost all of the standard deviations are less than 0.05 psu in other area with the NPIW, which shows that the time variation in the salinity can be neglected. This observational evidence shows that the topographic effect of the Izu Ridge on the horizontal distribution of the NPIW, which is formed east of 145°E by the mixing of the Kuroshio water and the Oyashio water, is prominent north of 30°N with a depth shallower than 2000 m.     

A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

A coupling of mixed and discontinuous Galerkin finite-element methods for poroelasticity

In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c _o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.     

变差函数参数的计算

为了求出地质统计学中变差函数理论模型的参数，首先把相应模型化成线性函数的形式，多面手应用线性方程组非负解的理论计算出它的参数。     

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.     

局部区域模型重力异常快速算法

地球重力场位系数模型可以用于计算局部重力扰动场元。然而随着地球重力场模型阶次的提高、局域重力场计算范围的增大,其计算速度往往不能满足工程需求。针对这一问题,在对位系数模型泰勒级数展开的基础上提出了采用向量运算、混合编程的方法,同时对连带勒让德函数Belikov递推方法中与经纬度无关的量进行了预先计算,有效提高了计算速度。提出的方法对于利用超高阶次重力场模型快速解算大范围、高分辨率重力场元数据以及累加求和计算具有一定的参考与借鉴意义。     

改进的无网格法及其在冻土区桩基温度场中的应用

采用带权正交基函数对传统的无网格法中的基函数进行了改进,避免了计算过程中可能出现的矩阵不可逆情况,而且编程容易实现且计算效率高。并将其应用到冻土区桩基包括了热传导和相变潜热的温度场中,同时考虑了混凝土水化热释放对桩周冻土的影响,将计算结果和有限元计算结果及现场实测数据进行了比较分析,从趋势和最大数值看都反映了实测曲线的趋势,又把不同深度处桩侧和不同桩径处温度随时间的变化规律计算结果和有限元计算结果做了对比,都验证了该方法的可行性和优越性。     

FINITE ELEMENT MODELLING FOR LAND DISPLACEMENTS DUE TO PUMPING

Equations of equilibrium (force balance) and flow in multidimensions were coupled in this paper to describe land displacements due to pressure decline in aquifers. A Galerkin finite element model based on these equations was developed. The saturated/unsaturated behaviour and the isotropic/anisotropic properties of permeability and elasticity were considered when the model was formulated. This model was verified by comparing its simulation results with those of known analytical solutions for simplified cases. The simulation of displacements due to pressure decline in unsaturated media was also performed. Those results demonstrated that the choice of boundary ranges for an aquifer with infinite domain may significantly affect the estimated horizontal and vertical displacements. To obtain a good estimation of land displacements, the boundary ranges should be carefully chosen. The displacements occurring in unconfined aquifers are caused by the drop of the water table and the change in body force in the dewatering zone. Simulation results also indicated that the change in body force should be considered once an unconfined aquifer has been pumped. Otherwise, the horizontal and vertical displacements in unconfined aquifers would be overestimated and underestimated, respectively. The behaviour of land displacements due to pumping was shown to be affected by changes in the total stresses in aquifers.