共查询到20条相似文献，搜索用时 31 毫秒
1.
A space-time discontinuous Galerkin finite element method is proposed and applied to a convection-dominant single-phase flow
problem in porous media. The numerical scheme is based on a coupled space-time finite element discretization allowing for
discontinuous approximations in space and in time. The continuities on the element interfaces are weakly enforced by the flux
treatments, so that no extra penalty factor has to be determined. The resulting space-time formulation possesses the advantage
of capturing the steep concentration front with sharp gradients efficiently. The stability and reliability of the proposed
approach is demonstrated by numerical experiments.
The author is grateful to the DFG (German Science Foundation—Deutsche Forschungsgemeinschaft) for the financial support under
the grant number Di 430/4-2. 相似文献
2.
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. 相似文献
3.
This paper presents a space-time adaptive framework for solving porous media flow problems, with specific application to reservoir simulation. A fully unstructured mesh discretization of space and time is used instead of a conventional time-marching approach. A space-time discontinuous Galerkin finite element method is employed to achieve a high-order discretization on the anisotropic, unstructured meshes. Anisotropic mesh adaptation is performed to reduce the error of a specified output of interest, by using a posteriori error estimates from the dual-weighted residual method to drive a metric-based mesh optimization algorithm. The space-time adaptive method is tested on a one-dimensional two-phase flow problem, and is found to be more efficient in terms of computational cost (degrees-of-freedom and total runtime) required to achieve a specified output error level, when compared to a conventional first-order time-marching finite volume method and the space-time discontinuous Galerkin method on structured meshes. 相似文献
4.
In this paper, we present a fast streamline-based numerical method for the two-phase flow equations in high-rate flooding scenarios for incompressible fluids in heterogeneous and anisotropic porous media. A fractional flow formulation is adopted and a discontinuous Galerkin method (DG) is employed to solve the pressure equation. Capillary effects can be neglected in high-rate flooding scenarios. This allows us to present an improved streamline approach in combination with the one-dimensional front tracking method to solve the transport equation. To handle the high computational costs of the DG approximation, domain decomposition is applied combined with an algebraic multigrid preconditioner to solve the linear system. Special care at the interior interfaces is required and the streamline tracer has to include a dynamic communication strategy. The method is validated in various two- and three-dimensional tests, where comparisons of the solutions in terms of approximation of flow front propagation with standard fully implicit finite-volume methods are provided. 相似文献
5.
Primary drainage in a water-wet saturated medium in the absence of capillarity is typically a combination of shock (discontinuous) and rarefaction (continuous) waves. Using nonlinear relative permeability functions for the host fluid and the invading fluid leads to the existence of a shock wave front, and the degree of nonlinearity of the relative permeability functions has an inverse relationship with the size of the shock wave (i.e., difference of saturation between upstream and downstream of the shock wave), whereas for linear relative permeability functions, the shock wave size approaches 0. Injection of a lower-viscosity immiscible phase such as gas or solvent into a water-wet porous medium in the presence of large capillary pressure leads to development of an extended and growing saturation transition zone that follows the discontinuous shock wave front. In this article, a semianalytical solution for the position of equisaturation contours (isosats) in the transition zone in the presence of gravity is obtained for a set of linearized relative permeability functions. The capillary (diffusive) and buoyancy terms are neglected, and the generalized convective equation for mass conservation is obtained. The set of equations is then reduced to a one-dimensional steady-state differential equation through forcing the isosat formulation to obey mass conservation. This scheme allows the isosat distribution to be solved, and the case of injection into an axisymmetric geometry for a confined planar configuration is solved and presented. A finite element model was developed to demonstrate the reasonable agreement between analytical and numerical solutions. 相似文献
6.
We study and compare five different combinations of finite element spaces for approximating the coupled flow and solid deformation system, so-called Biot’s equations. The permeability and porosity fields are heterogeneous and depend on solid displacement and fluid pressure. We provide detailed comparisons among the continuous Galerkin, discontinuous Galerkin, enriched Galerkin, and two types of mixed finite element methods. Several advantages and disadvantages for each of the above techniques are investigated by comparing local mass conservation properties, the accuracy of the flux approximation, number of degrees of freedom (DOF), and wall and CPU times. Three-field formulation methods with fluid velocity as an additional primary variable generally require a larger number of DOF, longer wall and CPU times, and a greater number of iterations in the linear solver in order to converge. The two-field formulation, a combination of continuous and enriched Galerkin function space, requires the fewest DOF among the methods that conserve local mass. Moreover, our results illustrate that three out of the five methods conserve local mass and produce similar flux approximations when conductivity alteration is included. These comparisons of the key performance indicators of different combinations of finite element methods can be utilized to choose the preferred method based on the required accuracy and the available computational resources. 相似文献
7.
We consider the slightly compressible two-phase flow problem in a porous medium with capillary pressure. The problem is solved using the implicit pressure, explicit saturation (IMPES) method, and the convergence is accelerated with iterative coupling of the equations. We use discontinuous Galerkin to discretize both the pressure and saturation equations. We apply two improvements, which are projecting the flux to the mass conservative H(div)-space and penalizing the jump in capillary pressure in the saturation equation. We also discuss the need and use of slope limiters and the choice of primary variables in discretization. The methods are verified with two- and three-dimensional numerical examples. The results show that the modifications stabilize the method and improve the solution. 相似文献
8.
本文建立了在考虑黄土中有节理存在时,水分场数值计算的有限元方法。采用质量守恒的观点推导了质点元中饱和度的变化与流速的关系。进而利用达西定律得到以水头为变量的渗流基本方程。针对黄土垂直节理的渗流特点,确定了节理渗流基本方程中的参数。采用四边形等参元,利用Galerkin加权余量法建立考虑节理影响的黄土非饱和渗流的有限元形式。对局部水头边界条件下的黄土节理二维水分入渗问题进行了数值分析。结果表明,节理对黄土场地湿润峰的迁移有很大影响。 相似文献
9.
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. 相似文献
10.
In this paper, we show how to couple the local discontinuous Galerkin method and the Raviart–Thomas mixed finite element method for elliptic equations modeling flow problems. We then show that the approximation of the velocity converges with the optimal order of k when we take the local discontinuous Galerkin that uses polynomials of degree k and the Raviart–Thomas space of polynomials of degree k?1. 相似文献
11.
We consider the finite element (FE) approximation of the two dimensional shallow water equations (SWE) by considering discretizations in which both space and time are established using a stable FE method. Particularly, we consider the automatic variationally stable FE (AVS-FE) method, a type of discontinuous Petrov-Galerkin (DPG) method. The philosophy of the DPG method allows us to establish stable FE approximations as well as accurate a posteriori error estimators upon solution of a saddle point system of equations. The resulting error indicators allow us to employ mesh adaptive strategies and perform space-time mesh refinements, i.e., local time stepping. We establish a priori error estimates for the AVS-FE method and linearized SWE and perform numerical verifications to confirm corresponding asymptotic convergence behavior. In an effort to keep the computational cost low, we consider an alternative space-time approach in which the space-time domain is partitioned into finite sized space-time slices. Hence, we can perform adaptive mesh refinements on each individual slice to preset error tolerances as needed for a particular application. Numerical verifications comparing the two alternatives indicate the space-time slices are superior for simulations over long times, whereas the solutions are indistinguishable for short times. Multiple numerical verifications show the adaptive mesh refinement capabilities of the AVS-FE method, as well the application of the method to some commonly applied benchmarks for the SWE. 相似文献
12.
This paper reviews and extends the theory and application of mimetic finite difference methods for the solution of diffusion problems in strongly heterogeneous anisotropic materials. These difference operators satisfy the fundamental identities, conservation laws and theorems of vector and tensor calculus on nonorthogonal, nonsmooth, structured and unstructured computational grids. We provide explicit approximations for equations in two dimensions with discontinuous anisotropic diffusion tensors. We mention the similarities and differences between the new methods and mixed finite element or hybrid mixed finite element methods. 相似文献
13.
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. 相似文献
14.
In this paper, we formulate and test numerically a fully-coupled discontinuous Galerkin (DG) method for incompressible two-phase flow with discontinuous capillary pressure. The spatial discretization uses the symmetric interior penalty DG formulation with weighted averages and is based on a wetting-phase potential/capillary potential formulation of the two-phase flow system. After discretizing in time with diagonally implicit Runge-Kutta schemes, the resulting systems of nonlinear algebraic equations are solved with Newton’s method and the arising systems of linear equations are solved efficiently and in parallel with an algebraic multigrid method. The new scheme is investigated for various test problems from the literature and is also compared to a cell-centered finite volume scheme in terms of accuracy and time to solution. We find that the method is accurate, robust, and efficient. In particular, no postprocessing of the DG velocity field is necessary in contrast to results reported by several authors for decoupled schemes. Moreover, the solver scales well in parallel and three-dimensional problems with up to nearly 100 million degrees of freedom per time step have been computed on 1,000 processors. 相似文献
15.
In this paper, we study the properties of approximate solutions to a doubly nonlinear and degenerate diffusion equation, known
in the literature as the diffusive wave approximation of the shallow water equations (DSW), using a numerical approach based
on the Galerkin finite element method. This equation arises in shallow water flow models when special assumptions are used
to simplify the shallow water equations and contains as particular cases the porous medium equation and the p-Laplacian. Diverse
numerical schemes have been implemented to approximately solve the DSW equation and have been successfully applied as suitable
models to simulate overland flow and water flow in vegetated areas such as wetlands; yet, no formal mathematical analysis
has been carried out in order to study the properties of approximate solutions. In this study, we propose a numerical approach
as a means to understand some properties of solutions to the DSW equation and, thus, to provide conditions for which the use
of the DSW equation may be inappropriate from both the physical and the mathematical points of view, within the context of
shallow water modeling. For analysis purposes, we propose a numerical method based on the Galerkin method and we obtain a
priori error estimates between the approximate solutions and weak solutions to the DSW equation under physically consistent
assumptions. We also present some numerical experiments that provide relevant information about the accuracy of the proposed
numerical method to solve the DSW equation and the applicability of the DSW equation as a model to simulate observed quantities
in an experimental setting. 相似文献
16.
We present the results of a study on a posteriori error control strategies for finite volume element approximations of second order elliptic differential equations. Finite volume methods ensure local mass conservation and, combined with some up-wind strategies, give monotone solutions. We adapt the local refinement techniques known from the finite element method to the finite volume discretizations of various boundary value problems for steady-state convection–diffusion–reaction equations. In this paper we derive and study a residual type error estimator and illustrate its practical performance on a series of computational tests in 2 and 3 dimensions. Our tests show that the discussed locally conservative approximation methods with a posteriori error control can be used successfully in numerical simulation of fluid flow and transport in porous media. 相似文献
17.
A key ingredient in simulation of flow in porous media is accurate determination of the velocities that drive the flow. Large‐scale irregularities of the geology (faults, fractures, and layers) suggest the use of irregular grids in simulation. This paper presents a control‐volume mixed finite element method that provides a simple, systematic, easily implemented procedure for obtaining accurate velocity approximations on irregular (i.e., distorted logically rectangular) block‐centered quadrilateral grids. The control‐volume formulation of Darcy’s law can be viewed as a discretization into element‐sized “tanks” with imposed pressures at the ends, giving a local discrete Darcy law analogous to the block‐by‐block conservation in the usual mixed discretization of the mass‐conservation equation. Numerical results in two dimensions show second‐order convergence in the velocity, even with discontinuous anisotropic permeability on an irregular grid. The method extends readily to three dimensions. 相似文献
18.
We study the gravity driven flow of two fluid phases in a one dimensional homogeneous porous column when history dependence of the pressure difference between the phases (capillary pressure) is taken into account. In the hyperbolic limit, solutions of such systems satisfy the Buckley-Leverett equation with a non-monotone flux function. However, solutions for the hysteretic case do not converge to the classical solutions in the hyperbolic limit in a wide range of situations. In particular, with Riemann data as initial condition, stationary shocks become possible in addition to classical components such as shocks, rarefaction waves and constant states. We derive an admissibility criterion for the stationary shocks and outline all admissible shocks. Depending on the capillary pressure functions, flux function and the Riemann data, two cases are identified a priori for which the solution consists of a stationary shock. In the first case, the shock remains at the point where the initial condition is discontinuous. In the second case, the solution is frozen in time in at least one semi-infinite half. The predictions are verified using numerical results. 相似文献
19.
构造了非线性守恒型浅水方程的伴随方程,并给出无结构三角网格的开边界校正表达式,利用Godunov格式和Riemann间断解Roe通量格式的设计思想,建立了无结构三角网格有限体积法的伴随开边界反演模型。进行了概化涌潮河口和钱塘江涌潮河口M2分潮开边界反演的数值实验,经过同化得到开边界上振幅的平均误差分别为0.000 6 m和0.053 3 m,证实了本文构造有限体积伴随模型的可行性,也表明本模型能够适应间断解问题。 相似文献
20.
We consider the modeling and simulation of compositional two-phase flow in a porous medium, where one phase is allowed to
vanish or appear. The modeling of Marchand et al. (in review) leads to a nonlinear system of two conservation equations. Each
conservation equation contains several nonlinear diffusion terms, which in general cannot be written as a function of the
gradients of the two principal unknowns. Also the diffusion coefficients are not necessarily explicit local functions of them.
For the generalised mixed finite elements approximation, Lagrange multipliers associated to each principal unknown are introduced,
the sum of the diffusive fluxes of each component is explicitly eliminated and the static condensation leads to a “global”
nonlinear system of equations only in the Lagrange multipliers also including complementarity conditions to cope with vanishing
or appearing phases. After time discretisation, this system can be solved at each time step using a semi-smooth Newton method.
The static condensation involves “local” nonlinear systems of equations associated to each element, solved also by a semismooth
Newton method. The algorithm is successfully applied to 1D and 2D examples of water–hydrogen flow involving gas phase appearance
and disappearance. 相似文献
|
