共查询到20条相似文献，搜索用时 312 毫秒
1.
The evolution of a gravity-driven free-surface flow of varying horizontal extent which couples with a field evolving within the flow is solved using a finite difference discretization of a mapping of the problem onto the unit square. Since the size of the solution domain may show several orders of magnitude of variation, while the normalized geometry of the domain and the internal field may not vary significantly, this procedure avoids excessively fine or coarse discretizations, as well as interpolations at the boundary. The parabolic and hyperbolic evolution equations for the internal field are considered. The evolution of the coupled system is solved by an implicit marching scheme. The discretizations in space and in time are accurate to second order. Multipoint upwinding is used to avoid an instability arising advective terms are large. The evolution equations are nonlinear, and are solved using a nested Newton–Raphson procedure. The nesting is achieved by using successively better approximations to the ture evolution equations. The matrix equation that arises is solved by a conjugate-gradient-like (ORTHOMIN) iteration procedure with an incomplete Cholesky factorization preconditioning. The method has a wide variety of potential applications in the earth sciences, with the ability to describe glacier flow, lava flow, avalanching and landslides. Some calculations of the thermomechanical evolution of ice-sheets are given as illustrations, and the possible existence of thermally induced instabilities is considered. 相似文献
2.
Uncoupled analysis of stabilizing piles in weathered slopes 总被引：15，自引：0，他引：15
This paper describes a simplified numerical approach for analyzing the slope/pile system subjected to lateral soil movements. The lateral one-row pile response above and below the critical surface is computed by using load transfer approach. The response of groups was analyzed by developing interaction factors obtained from a three-dimensional nonlinear finite element study. An uncoupled analysis was performed for stabilizing piles in slope in which the pile response and slope stability are considered separately. The non-linear characteristics of the soil–pile interaction in the stabilizing piles are modeled by hyperbolic load transfer curves. The Bishop's simplified method of slope stability analysis is extended to incorporate the soil-pile interaction and evaluate the safety factor of the reinforced slope. Numerical study is performed to illustrate the major influencing parameters on the pile-slope stability problem. Through comparative studies, it has been found that the factor of safety in slope is much more conservative for an uncoupled analysis than for a coupled analysis based on three-dimensional finite element analysis. 相似文献
3.
4.
Philip Cardiff 《国际地质力学数值与分析法杂志》2015,39(13):1410-1430
Accurate prediction of the interactions between the nonlinear soil skeleton and the pore fluid under loading plays a vital role in many geotechnical applications. It is therefore important to develop a numerical method that can effectively capture this nonlinear soil‐pore fluid coupling effect. This paper presents the implementation of a new finite volume method code of poro‐elasto‐plasticity soil model. The model is formulated on the basis of Biot's consolidation theory and combined with a perfect plasticity Mohr‐Coulomb constitutive relation. The governing equation system is discretized in a segregated manner, namely, those conventional linear and uncoupled terms are treated implicitly, while those nonlinear and coupled terms are treated explicitly by using any available values from previous time or iteration step. The implicit–explicit discretization leads to a linearized and decoupled algebraic system, which is solved using the fixed‐point iteration method. Upon the convergence of the iterative method, fully nonlinear coupled solutions are obtained. Also explored in this paper is the special way of treating traction boundary in finite volume method compared with FEM. Finally, three numerical test cases are simulated to verify the implementation procedure. It is shown in the simulation results that the implemented solver is capable of and efficient at predicting reasonable soil responses with pore pressure coupling under different loading situations. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
5.
In this work, we measure the performance of the fixed stress split algorithm for the immiscible water-oil flow coupled with linear poromechanics. The two-phase flow equations are solved on general hexahedral elements using the multipoint flux mixed finite element method whereas the poromechanics equations are discretized using the conforming Galerkin method. We introduce a rigorous calculation of the update in poroelastic properties during the iterative solution of the coupled system equations. The effects of the coupling parameter on the performance of the fixed stress algorithm is demonstrated in two field studies: the Frio oil reservoir and the Cranfield injection site.
相似文献6.
E. Bruce Pitman 《国际地质力学数值与分析法杂志》1993,17(6):385-400
The equations governing the elastic-plastic deformation of granular materials are typically hyperbolic, or contain small-magnitude damping or rate effects. A finite element algorithm is the standard method for the numerical integration of these systems. In particular, finite elements allow great flexibility in the design of grid geometry. However, modern finite difference methods for hyperbolic systems have been successful in aerodynamics computations, resolving wave structures more sharply than finite element schemes. In this paper we develop a finite difference scheme for granular flow problems. We report on a second-order Godunov-type scheme for the integration of hyperbolic equations for the elastoplastic deformation of a simple model of granular flow. The Godunov method includes a characteristic tracing step in the integration, providing minimal wave dispersion, and a slope limiting step, preventing unphysical oscillations. The granular flow model we consider is hyperbolic, but hyperbolicity is lost at a large value of accumulated plastic strain. This loss of hyperbolicity is a tell-tale signal for the formation of a shear band within the sample. Typically, when systems lose hyperbolicity a regularization mechanism is added to the model equations in order to maintain the well posedness of the system. These regularizations include viscosity, viscoplasticity, higher-order gradient effects or stress coupling. Here we appeal to a very different kind of regularization. When the system loses hyperbolicity and a shear band forms, we treat the band as an internal boundary, and impose jump conditions at this boundary. Away from the band, the system remains hyperbolic and the integration step proceeds as usual. 相似文献
7.
不连续面在双重介质热-水-力三维耦合分析中的有限元数值实现 总被引：1，自引：0，他引：1
在建立双重介质热-水-力耦合微分控制方程的基础上,提出了裂隙岩体热-水-力耦合的三维力学模型,对不同介质分别建立以节点位移、水压力和温度为求解量的三维有限元格式,开发了双重介质热-水-力耦合分析的的三维有限元计算程序,在有限元数值分析中不连续面应力计算采用等厚度空间8节点节理单元进行离散,而不连续面渗流和热能计算时采用平面4节点等参单元进行离散,这样保证了不同介质之间的水量、热量交换和两类模型接触处节点水头、温度和位移相等。通过高温岩体地热开发算例,揭示了在热-水-力耦合作用下不连续面处于低应力区,其张开度随运行时间的延长呈非线性增加,非稳定渗流阶段不连续面显著地控制着渗流场的整体分布,它的水头远高于拟连续岩体介质的水头,而进入稳定渗流阶段不连续面的控渗作用不明显,由于高温岩体地热开发系统中存在大规模的热量补给,不连续面对岩体温度场分布的影响并不显著。 相似文献
8.
One of the major causes of instability in geotechnical structures such as dikes or earth dams is the phenomenon of suffusion including detachment, transport and filtration of fine particles by water flow. Current methods fail to capture all these aspects. This paper suggests a new modeling approach under the framework of the porous continuous medium theory. The detachment and transport of the fine particles are described by a mass exchange model between the solid and the fluid phases. The filtration is incorporated to simulate the filling of the inter-grain voids created by the migration of the fluidized fine particles with the seepage flow, and thus, the self-filtration is coupled with the erosion process. The model is solved numerically using a finite difference method restricted to one-dimensional (1-D) flows normal to the free surface. The applicability of the model to capture the main features of both erosion and filtration during the suffusion process has been validated by simulating 1-D internal erosion tests and by comparing the numerical with the experimental results. Furthermore, the influence of the coupling between erosion and filtration has been highlighted, including the development of material heterogeneity induced by the combination of erosion and filtration.
相似文献9.
Clarisse Alboin Jérôme Jaffré Patrick Joly Jean E. Roberts Christophe Serres 《Computational Geosciences》2002,6(3-4):523-543
Contaminant transport in a fractured porous medium can be modeled, under appropriate conditions, with a double porosity model. Such a model consists of a parabolic equation with a coupling term describing contaminant exchange between the fractures, which have high permeability, and the matrix block, which has low permeability. A locally conservative method based on mixed finite elements is used to solve the parabolic problem, and the calculation of the coupling term, which involves the solution of diffusion equations in the matrix blocks, is based on an analytic expression. Numerical experiments show that this semi-analytic method for the coupling term is accurate and faster than several other methods but at a small expense of computer memory. 相似文献
10.
11.
12.
The present paper proposes a new family of multiscale finite volume methods. These methods usually deal with a dual mesh resolution, where the pressure field is solved on a coarse mesh, while the saturation fields, which may have discontinuities, are solved on a finer reservoir grid, on which petrophysical heterogeneities are defined. Unfortunately, the efficiency of dual mesh methods is strongly related to the definition of up-gridding and down-gridding steps, allowing defining accurately pressure and saturation fields on both fine and coarse meshes and the ability of the approach to be parallelized. In the new dual mesh formulation we developed, the pressure is solved on a coarse grid using a new hybrid formulation of the parabolic problem. This type of multiscale method for pressure equation called multiscale hybrid-mixed method (MHMM) has been recently proposed for finite elements and mixed-finite element approach (Harder et al. 2013). We extend here the MH-mixed method to a finite volume discretization, in order to deal with large multiphase reservoir models. The pressure solution is obtained by solving a hybrid form of the pressure problem on the coarse mesh, for which unknowns are fluxes defined on the coarse mesh faces. Basis flux functions are defined through the resolution of a local finite volume problem, which accounts for local heterogeneity, whereas pressure continuity between cells is weakly imposed through flux basis functions, regarded as Lagrange multipliers. Such an approach is conservative both on the coarse and local scales and can be easily parallelized, which is an advantage compared to other existing finite volume multiscale approaches. It has also a high flexibility to refine the coarse discretization just by refinement of the lagrange multiplier space defined on the coarse faces without changing nor the coarse nor the fine meshes. This refinement can also be done adaptively w.r.t. a posteriori error estimators. The method is applied to single phase (well-testing) and multiphase flow in heterogeneous porous media. 相似文献
13.
A model of chemo-osmotic consolidation of clays in multi-dimensional domains is presented, with volumetric strains induced by both changes in the chemistry and osmotically driven pore water flow considered. Three fully coupled governing equations considering force equilibrium, pore water transport and chemical transport are presented and solved using the transient finite element method. The proposed approach is verified via consideration of chemo-osmotic consolidation of a compacted clay landfill liner and then applied to investigation of a hypothetical case with a local leachate leak in the compacted clay liner. An assessment of the impact of the two-dimensional nature of the system is made. The consolidation process is found to be dominated by osmotic consolidation in the early stages and subsequently by chemical consolidation. It is found that the surface settlement and the leachate concentration in the compacted clay liner may be highly overestimated by a one-dimensional analysis. Moreover, the peak negative excess pore pressure predicted by the two-dimensional solution remains in a shallow region under the leak while in the one-dimensional solution it moves progressively downwards. 相似文献
14.
Benjamin Ganis Mark E. Mear A. Sakhaee-Pour Mary F. Wheeler Thomas Wick 《Computational Geosciences》2014,18(5):613-624
We describe an algorithm for modeling saturated fractures in a poroelastic domain in which the reservoir simulator is coupled with a boundary element method. A fixed stress splitting is used on the underlying fractured Biot system to iteratively couple fluid and solid mechanics systems. The fluid system consists of Darcy’s law in the reservoir and is computed with a multipoint flux mixed finite element method, and a Reynolds’ lubrication equation in the fracture solved with a mimetic finite difference method. The mechanics system consists of linear elasticity in the reservoir and is computed with a continuous Galerkin method, and linear elasticity in the fracture is solved with a weakly singular symmetric Galerkin boundary element method. This algorithm is able to compute both unknown fracture width and unknown fluid leakage rate. An interesting numerical example is presented with an injection well inside of a circular fracture. 相似文献
15.
Zhang Hong-Wu 《国际地质力学数值与分析法杂志》1995,19(12):851-867
Finite element procedures for numerical solution of various engineering problems are often based on variational formulations. In this paper, a parametric variational principle applicable to elastic-plastic coupled field problems in consolidation analysis of saturated porous media is presented. This principle can be used to solve problems where materials are inconsistent with Drucker's postulate of stability, such as in non-associated plasticity flow or softening problems. The finite element formulation was given, and it can be solved by either the conventional method or a parametric quadratic programming method. 相似文献
16.
构造-流体-成矿体系的反应-输运-力学耦合模型和动力学模拟 总被引：12，自引：0，他引：12
建立了一个综合的构造流体成矿体系的反应输运力学耦合动力学模型。利用有限元方法求解岩石变形、断裂作用和断裂网络统计动力学、流体流动、有机和无机地球化学反应及成岩成矿作用、压力溶液和其它压实力学、热迁移的方程组 ,可以对构造流体成矿体系的动力学演化过程进行 1～ 3维数值模拟。模拟的主要内容是在各种过程耦合作用下描述构造流体成矿体系的主要变量的时空演化 :( 1)与成矿流体的形成和性质有关的变量 ,如地层中矿物 (包括成矿物质 )的溶解速率、流体中各组分的浓度与饱和度、流体温度、压力、离子强度等 ;( 2 )与构造变形和流体运移有关的各变量 ,如应力与变形速率、岩石孔隙度、构造 (断裂 )渗透率等 ;( 3 )与沉淀成矿有关的变量 ,如矿物 (金属矿物和脉石矿物 )的成核速率、各矿物的沉淀量等 ;( 4 )上述各有关变量间的时空耦合关系 ,如断裂渗透率时空演化与流体流动、汇聚和成矿的耦合关系等。以湖南沃溪金锑钨矿床为例 ,应用该模型和方法对成矿动力学过程和动力学机制进行了初步的模拟与分析。 相似文献
17.
热液系统流体输运－化学反应耦合动力学综述 总被引：2，自引：0，他引：2
输运－反应模式通过连续质量守恒方程进行描述，它表示为主要组分的Ｎ＋Ｍ个耦合的非线性偏微分方程。这些方程可以通过有限差分或有限元等数值方法进行求解。利用准静态近似法，涉及平流、分子扩散和水动力弥散作用的地球化学系统的演化在时间上通过一系列静止状态模拟出来。模拟实际地质系统的动力学实验对输运－反应耦合动力学的应用具有十分重要的影响。 相似文献
18.
Eirik Keilegavlen Jeremy E. Kozdon Bradley T. Mallison 《Computational Geosciences》2012,16(4):1021-1042
The governing equations for multiphase flow in porous media have a mixed character, with both nearly elliptic and nearly hyperbolic variables. The flux for each phase can be decomposed into two parts: (1) a geometry- and rock-dependent term that resembles a single-phase flux; and (2) a mobility term representing fluid properties and rock–fluid interactions. The first term is commonly discretized by two- or multipoint flux approximations (TPFA and MPFA, respectively). The mobility is usually treated with single-point upstream weighting (SPU), also known as dimensional or donor cell upstream weighting. It is well known that when simulating processes with adverse mobility ratios, SPU suffers from grid orientation effects. An important example of this, which will be considered in this work, is the displacement of a heavy oil by water. For these adverse mobility ratio flows, the governing equations are unstable at the modeling scale, rendering a challenging numerical problem. These challenges must be addressed in order to avoid systematic biasing of simulation results. In this work, we present a framework for multidimensional upstream weighting for multiphase flow with buoyancy on general two-dimensional grids. The methodology is based on a dual grid, and the resulting transport methods are provably monotone. The multidimensional transport methods are coupled with MPFA methods to solve the pressure equation. Both explicit and fully implicit approaches are considered for time integration of the transport equations. The results show considerable reduction of grid orientation effects compared to SPU, and the explicit multidimensional approach allows larger time steps. For the implicit method, the total number of non-linear iterations is also reduced when multidimensional upstream weighting is used. 相似文献
19.
Abstract The calculations of unsteady flow to a multiple well system with the application of boundary element method (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem of parabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time variable. The image function of water head H can be solved by BEM. We derived the boundary integral equation of the transformed variable H and the discretization form of it, so that there is no need to discretize the boundaries of well walls and it becomes easier to solve the groundwater head H by numerical inversion. 相似文献