排序方式: 共有110条查询结果,搜索用时 15 毫秒
1.
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. 相似文献
2.
Estimation of Land Subsidence Based on Groundwater Flow Model 总被引:8,自引:0,他引:8
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. 相似文献
3.
A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow
Clint Dawson 《Computational Geosciences》2008,12(4):451-472
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. 相似文献
4.
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. 相似文献
5.
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. 相似文献
6.
7.
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. 相似文献
8.
In this paper, the seismic response of ‘infinitely’ long slopes is numerically analysed via the formulation of a 1D analytical/numerical model, in which the soil mechanical behaviour is assumed to be elasto‐perfectly viscoplastic and simple shear (SS) kinematical constraints are imposed. In order to simplify the problem, a theoretically based procedure to set up a fully 1D shear constitutive model is defined, within which the mechanical response of a multiaxial relationship is condensed. The use of a 1D shear constitutive model is aimed at reducing the number of unknowns and, therefore, the computational costs. In particular, the case of the Mohr–Coulomb yield criterion is considered, while an enhanced Taylor–Galerkin finite element algorithm is employed to simulate the seismic wave propagation within the soil stratum. The proposed ‘condensation’/calibration procedure captures both the ‘pseudo’‐hardening pre‐failure behaviour and the influence of dilation on the occurrence of strain‐localization, which characterize, under SS conditions, the static response of virgin perfectly plastic soils. The effectiveness of the conceived method is shown with reference to freshly deposited deposits, while, in the case of highly overconsolidated strata, some difficulties arise because of the brittle behaviour induced both by unloading and non‐associativeness. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
9.
无网格伽辽金法应用的参数选择及内部边界处理 总被引:8,自引:0,他引:8
无网格辽金法作为-种新的岩土工程数值计算方法, 该法其只需节点信息的无单元特性, 使其具有计算优势。本文结合固结EFGM刚度矩阵公式, 对不同的计算参数进行计算分析, 找出其影响规律。并采用跳跃函数处理内部边界条件, 计算结果表明, EFGM处理内部场函数不连续是准确的。 相似文献
10.
In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation
model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin
method. Theoretical convergence error estimates are derived in a discrete-in-time setting. Of particular interest is the case
when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation
and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently
applied to existing reservoir engineering simulators. 相似文献