共查询到20条相似文献，搜索用时 46 毫秒
1.
In this paper, we present a numerical model for simulating two-phase (oil–water and air–water) incompressible and immiscible flow in porous media. The mathematical model which is based on a fractional flow formulation is formed of two nonlinear partial differential equations: a mean pressure equation and a water saturation equation. These two equations can be solved in a sequential manner. Two numerical methods are used to discretize the equations of the two-phase flow model: mixed hybrid finite elements are used to treat the pressure equation, h-based Richards' equation and the diffusion term in the saturation equation, the advection term in the saturation equation is treated with the discontinuous finite elements. We propose a better way to calculate the nonlinear coefficients contained in our equations on each element of the discretized domain. In heterogeneous porous media, the saturation becomes discontinuous at the interface between two porous media. We show in this paper how to use the capillary pressure–saturation relationship in order to handle the saturation jump in the mixed hybrid finite element method. The two-phase flow simulator is verified against analytical solutions for some flow problems treated by other authors. 相似文献
2.
B‐bar based algorithm applied to meshfree numerical schemes to solve unconfined seepage problems through porous media 
下载免费PDF全文
下载免费PDF全文 The object of this work is to establish a meshfree framework for solving coupled, steady and transient problems for unconfined seepage through porous media. The Biot's equations are formulated in displacements (or u − w) assuming an elastic solid skeleton. The free surface location and its evolution in time are obtained by interpolation of pore water pressures throughout the domain. Shape functions based on the principle of local maximum entropy are chosen for the meshfree approximation schemes. In order to avoid the locking involved in the fluid phase of the porous media, a B‐bar based algorithm is devised to compute the average volumetric strain in a patch composed of various integration points. The efficiency of such an implementation for one phase problems is shown through the Benchmark problem, Cook's membrane loaded by a distributive shear load. The proposed methodology is firstly applied to various classical examples in unconfined steady seepage problems through earth dams, then to the dynamic consolidation of a soil column. The results obtained for both problems are quite satisfactory and demonstrate the feasibility of the proposed method in solving coupled problems in porous media. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
3.
Smoothed particle hydrodynamics (SPH) is a versatile technique which can be applied to single and multiphase flow through porous media. The versatility of SPH is offset by its computational expense which limits the practicability of SPH for large problems involving low Reynolds number flow. A parallel pore-scale numerical model based on SPH is described for modeling flow phenomena in porous media. Aspects of SPH which complicate parallelization are emphasized. The speed of the method is demonstrated to be proportional to the number of processors for test cases where load balance was achieved. The parallel algorithm permits the application of SPH to more complicated porous media problems than previously considered. For such problems, best performance is achieved when several soil grains are simulated by each processor. Finally, future applications of the method and possible extensions are discussed. 相似文献
4.
Mohammad Reza Maleki Javan 《国际地质力学数值与分析法杂志》2015,39(3):229-250
It is well known that for a sufficiently high seepage velocity, the governing flow law of porous media is nonlinear (J. Computers & Fluids 2010; 39 : 2069–2077). However, this fact has not been considered in the studies of soil‐pore fluid interaction and in conventional soil mechanics. In the present paper, a fully explicit dynamic finite element method is developed for nonlinear Darcy law. The governing equations are expressed for saturated porous media based on the extension of the Biot (J. Appl. Phys. 1941; 12 : 155–164) formulation. The elastoplastic behavior of soil under earthquake loading is simulated using a generalized plasticity theory that is composed of a yield surface along with non‐associated flow rule. Numerical simulations of porous media subjected to horizontal and vertical components of ground motion excitations with different permeability coefficients are carried out; while computed maximum pore water pressure is specially taken into consideration to make the difference between Darcy and non‐Darcy flow regimes tangible. Finally, the effect of non‐Darcy flow on the evaluated liquefaction potential of sand in comparison to conventional Darcy law is examined. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
5.
A numerical approach is proposed to model the flow in porous media using homogenization theory. The proposed concept involves the analyses of micro‐true flow at pore‐level and macro‐seepage flow at macro‐level. Macro‐seepage and microscopic characteristic flow equations are first derived from the Navier–Stokes equation at low Reynolds number through a two‐scale homogenization method. This homogenization method adopts an asymptotic expansion of velocity and pressure through the micro‐structures of porous media. A slightly compressible condition is introduced to express the characteristic flow through only characteristic velocity. This characteristic flow is then numerically solved using a penalty FEM scheme. Reduced integration technique is introduced for the volumetric term to avoid mesh locking. Finally, the numerical model is examined using two sets of permeability test data on clay and one set of permeability test data on sand. The numerical predictions agree well with the experimental data if constraint water film is considered for clay and two‐dimensional cross‐connection effect is included for sand. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
6.
A primal C0-conforming virtual element discretization for the approximation of the bidimensional two-phase flow of immiscible fluids in porous media using general polygonal meshes is discussed. This work investigates the potentialities of the Virtual Element Method (VEM) in solving this specific problem of immiscible fluids in porous media involving a time-dependent coupled system of non-linear partial differential equations. The performance of the fully discrete scheme is thoroughly analysed testing it on general meshes considering both a regular problem and more realistic benchmark problems that are of interest for physical and engineering applications.
相似文献7.
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. 相似文献
8.
George Zyvoloski 《国际地质力学数值与分析法杂志》1983,7(1):75-86
Two finite element algorithms suitable for long term simulation of geothermal reservoirs are presented. Both methods use a diagonal mass matrix and a Newton iteration scheme. The first scheme solves the 2N unsymmetric algebraic equations resulting from the finite element discretization of the equations governing the flow of heat and mass in porous media by using a banded equation solver. The second method, suitable for problems in which the transmissibility terms are small compared to the accumulation terms, reduces the set of N equations for the Newton corrections to a symmetric system. Comparison with finite difference schemes indicates that the proposed algorithms are competitive with existing methods. 相似文献
9.
Variationally consistent homogenization is exploited for the analysis of transient uncoupled consolidation in micro‐heterogeneous porous solids, whereby the classical approach of first‐order homogenization for stationary problems is extended to transient problems. Homogenization is then carried out in the spatial domain on representative volume elements (RVE), which are introduced in quadrature points in standard fashion. Along with the classical averages, a higher‐order conservation quantity is obtained. An iterative FE2‐algorithm is devised for the case of nonlinear permeability and storage coefficients, and it is applied to pore pressure changes in asphalt‐concrete (particle composite). Various parametric studies are carried out, in particular, with respect to the influence of the ‘substructure length scale’ that is represented by the size of the RVE's. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
10.
A novel two-dimensional mixed fracture–pore seepage model for fluid flow in fractured porous media is presented based on the computational framework of finite-discrete element method (FDEM). The model consists of a porous seepage model in triangular elements bonded by unbroken joint elements, as well as a fracture seepage model in broken joint elements. The principle for determining the fluid exchange coefficient of the unbroken joint element is provided to ensure numerical accuracy and efficiency. The mixed fracture–pore seepage model provides a simple but effective tool for solving fluid flow in fractured porous media. In this paper, examples of 1D and 2D seepage flow in porous media and porous media with a single fracture or multiple fractures are studied. The simulation results of the model match well with theoretical solutions or results obtained by commercial software, which verifies the correctness of the mixed fracture–pore seepage model. Furthermore, combining FDEM mechanical calculation and the mixed fracture–pore seepage model, a coupled hydromechanical model is built to simulate fluid-driven dynamic propagation of cracks in the porous media, as well as its influence on pore seepage and fracture seepage.
相似文献11.
This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek−1 is 𝒪(1). These variations are estimated to be less than 5% if Ek−1<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
12.
多尺度有限单元法求解非均质多孔介质中的三维地下水流问题 总被引：7，自引：0，他引：7
应用多尺度有限单元法模拟非均质多孔介质中的三维地下水流问题。与传统有限单元法相比，多尺度有限单元法的基函数具有能反映单元内参数变化的优点，所以这种方法能在大尺度上抓住解的小尺度特征获得较精确的解。在介绍多尺度有限单元法求解非均质多孔介质中三维地下水流问题的基本原理之后，对参数水平方向渐变垂直方向突变的非均质多孔介质中的三维地下水流和Borden实验场的三维地下水流分别用多尺度有限单元法和传统等参有限单元法进行了计算，结果表明在模拟高度非均质多孔介质中的三维地下水流问题时，多尺度有限单元法比传统有限单元法有效，既节省计算量又有较高的精度；在模拟非均质性弱的多孔介质中的三维地下水流问题时，多尺度有限单元法虽然也能在大尺度上获得较为精确的解，但效果不明显。 相似文献
13.
A finite element method for modeling coupled flow and deformation in porous fractured media 下载免费PDF全文
下载免费PDF全文 Ahmad Pouya 《国际地质力学数值与分析法杂志》2015,39(16):1836-1852
Modeling the flow in highly fractured porous media by finite element method (FEM) has met two difficulties: mesh generation for fractured domains and a rigorous formulation of the flow problem accounting for fracture/matrix, fracture/fracture, and fracture/boundary fluid mass exchanges. Based on the recent theoretical progress for mass balance conditions in multifractured porous bodies, the governing equations for coupled flow and deformation in these bodies are first established in this paper. A weak formulation for this problem is then established allowing to build a FEM. Taking benefit from recent development of mesh‐generating tools for fractured media, this weak formulation has been implemented in a numerical code and applied to some typical problems of hydromechanical coupling in fractured porous media. It is shown that in this way, the FEM that has proved its efficiency to model hydromechanical phenomena in porous media is extended with all its performances (calculation time, couplings, and nonlinearities) to fractured porous media. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
14.
In this paper, a three‐dimensional smooth particle hydrodynamics (SPH) simulator for modeling grain scale fluid flow in porous media is presented. The versatility of the SPH method has driven its use in increasingly complex areas of flow analysis, including the characterization of flow through permeable rock for both groundwater and petroleum reservoir research. SPH provides the means to model complex multi‐phase flows through such media; however, acceptance of the methodology has been hampered by the apparent lack of actual verification within the literature, particulary in the three‐dimensional case. In this paper, the accuracy of SPH is addressed via a comparison to the previously recognized benchmarks of authors such as Sangani and Acrivos (Int. J. Multiphase Flow 1982; 8 (4): 343–360), Zick and Homsy (J. Fluid Mech. 1982; 115 :13–26) and Larson and Higdon (Phys. Fluids A 1989; 1 (1):38–46) for the well‐defined classical problems of flow through idealized two‐ and three‐dimensional porous media. The accuracy of results for such low Reynolds number flows is highly dependent on the implementation of no‐slip boundary conditions. A new, robust and numerically efficient, method for implementing such boundaries in SPH is presented. Simulation results for friction coefficient and permeability are shown to agree well with the available benchmarks. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
15.
This work investigates the profiles of surface water flow over a pervious pavement, for example, highway, during a uniform rainfall. The pavement, such as porous asphalt or open‐graded asphalt friction course, is regarded as a porous medium, and the flow inside the layer is porous media flow. At first, the velocity distributions are solved based on the simplified Navier–Stokes equations and Biot's theory of poroelasticity for the water layer and the permeable layer, respectively. Then, the flow profiles can be found via the continuity equation. Because the closed form of the water depth cannot be obtained, a numerical technique is employed to find the flow profiles on the pavement surface. A critical permeability factor related to the material and structure of pavement is illustrated to be η = 1.36 × 10? 4m as the half road width L = 20 m, the cross slope S = 0.02,the pavement thickness H = 0.05 m, and rainfall intensity i = 100 mm/h. The water on the road surface can be drained very fast, whereas the parameter is greater than the critical value. Another illustration of allowable/safety water depth is made for the design of the pavement thickness. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
16.
Two formulations for calculating dynamic response of a cylindrical cavity in cross‐anisotropic porous media based on complex functions theory are presented. The basis of the method is the solution of Biot's consolidation equations in the complex plane. Employing two groups of potential functions for solid skeleton and pore fluid (each group includes three functions), the u–w formulation of Biot's equations are solved. Difference of these two solutions refers to use of two various potential functions. Equations for calculating stress, displacement and pore pressure fields of the medium are mentioned based on each two formulations. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
17.
Joachim Berdal Haga Harald Osnes Hans Petter Langtangen 《Computational Geosciences》2012,16(3):723-734
Large-scale simulations of coupled flow in deformable porous media require iterative methods for solving the systems of linear
algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps
in material properties, which is often the case in realistic geological applications, such as basin evolution at regional
scales. The success of iterative methods for such problems depends strongly on finding effective preconditioners with good
parallel scaling properties, which is the topic of the present paper. We present a parallel preconditioner for Biot’s equations
of coupled elasticity and fluid flow in porous media. The preconditioner is based on an approximation of the exact inverse
of the two-by-two block system arising from a finite element discretisation. The approximation relies on a highly scalable
approximation of the global Schur complement of the coefficient matrix, combined with generally available state-of-the-art
multilevel preconditioners for the individual blocks. This preconditioner is shown to be robust on problems with highly heterogeneous
material parameters. We investigate the weak and strong parallel scaling of this preconditioner on up to 512 processors and
demonstrate its ability on a realistic basin-scale problem in poroelasticity with over eight million tetrahedral elements. 相似文献
18.
Jeremy Edward Kozdon Bradley T. Mallison Margot G. Gerritsen 《Computational Geosciences》2011,15(3):399-419
Truly multidimensional methods for hyperbolic equations use flow-based information to determine the computational stencil,
as opposed to applying one-dimensional methods dimension by dimension. By doing this, the numerical errors are less correlated
with the underlying computational grid. This can be important for reducing bias in flow problems that are inherently unstable
at simulation scale, such as in certain porous media problems. In this work, a monotone, multi-D framework for multiphase
flow and transport in porous media is developed. A local coupling of the fluxes is introduced through the use of interaction
regions, resulting in a compact stencil. A relaxed volume formulation of the coupled hyperbolic–elliptic system is used that
allows for nonzero residuals in the pressure equation to be handled robustly. This formulation ensures nonnegative masses
and saturations (volume fractions) that sum to one (Acs et al., SPE J 25(4):543–553, 1985). Though the focus of the paper is on immiscible flow, an extension of the methods to a class of more general scalar hyperbolic
equations is also presented. Several test problems demonstrate that the truly multi-D schemes reduce biasing due to the computational
grid. 相似文献
19.
Large‐scale simulations of flow in deformable porous media require efficient iterative methods for solving the involved systems of linear algebraic equations. Construction of efficient iterative methods is particularly challenging in problems with large jumps in material properties, which is often the case in geological applications, such as basin evolution at regional scales. The success of iterative methods for this type of problems depends strongly on finding effective preconditioners. This paper investigates how the block‐structured matrix system arising from single‐phase flow in elastic porous media should be preconditioned, in particular for highly discontinuous permeability and significant jumps in elastic properties. The most promising preconditioner combines algebraic multigrid with a Schur complement‐based exact block decomposition. The paper compares numerous block preconditioners with the aim of providing guidelines on how to formulate efficient preconditioners. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
20.
Based on the Biot theory, the exact solutions for one‐dimensional transient response of single layer of fluid‐saturated porous media and semi‐infinite media are developed, in which the fluid and solid particles are assumed to be compressible and the inertial, viscous and mechanical couplings are taken into account. First, the control equations in terms of the solid displacement u and a relative displacement w are expressed in matrix form. For problems of single layer under homogeneous boundary conditions, the eigen‐values and the eigen‐functions are obtained by means of the variable separation method, and the displacement vector u is put forward using the searching method. In the case of nonhomogeneous boundary conditions, the boundary conditions are first homogenized, and the displacement field is constructed basing upon the eigen‐functions. Making use of the orthogonality of eigen‐functions, a series of ordinary differential equations with respect to dimensionless time and their corresponding initial conditions are obtained. Those differential equations are solved by the state‐space method, and the series solutions for three typical nonhomogeneous boundary conditions are developed. For semi‐infinite media, the exact solutions in integral form for two kinds of nonhomogeneous boundary conditions are presented by applying the cosine and sine transforms to the basic equations. Finally, three examples are studied to illustrate the validity of the solutions, and to assess the influence of the dynamic permeability coefficient and the fluid inertia to the transient response of porous media. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献