Superconvergence of a combined mixed finite element and discontinuous Galerkin method for a compressible miscible displacement problem

A combined mixed finite element and discontinuous Galerkin method for a compressible miscible displacement problem which includes molecular diffusion and dispersion in porous media is investigated. That is to say, the mixed finite element method with Raviart-Thomas space is applied to the flow equation, and the transport one is solved by the symmetric interior penalty discontinuous Galerkin (SIPG) approximation. Based on projection interpolations and induction hypotheses, a superconvergence estimate is obtained. During the analysis, an extension of the Darcy velocity along the Gauss line is also used in the evaluation of the coefficients in the Galerkin procedure for the concentration.     

Superconvergence of a combined mixed finite element and discontinuous Galerkin approximation for an incompressible miscible displacement problem

A combined mixed finite element and discontinuous Galerkin approximation for an incompressible miscible displacement problem which includes molecular diffusion and dispersion in porous media is studied. That is to say, the mixed finite element method is applied to the flow equation, and the transport equation is solved by an interior penalty discontinuous Galerkin method. Convolution of the Darcy velocity approximation with the Bramble-Schatz kernel function and averaging are applied in the evaluation of the coefficients in the Galerkin procedure for the concentration. A superconvergence estimate is obtained. Numerical experimental results are presented to verify the theoretical analysis.     

A combined mixed and discontinuous Galerkin method for compressible miscible displacement problem in porous media

In this paper, we present a numerical scheme for solving the coupled system of compressible miscible displacement problem in porous media. The flow equation is solved by the mixed finite element method, and the transport equation is approximated by a discontinuous Galerkin method. The scheme is continuous in time and a priori hp error estimates is presented.     

非牛顿幂律流体球向不定常渗流

本文研究了弱压缩非牛顿幂律流体球向不定常渗流,导出了抛物型偏微分非线性方程.球向扩散方程是其特殊情况.用Laplace变换的方法,找到了线性化后方程的解析解和渐近解.用影响半径的概念和平均值方法求得了近似解.渐近解和近似解的结构是相似的,从而丰富了非牛顿流体一维不定常渗流的理论.     

Front tracking for two-phase flow in reservoir simulation by adaptive mesh

In this article, an algorithm for the numerical approximation of two-phase flow in porous media by adaptive mesh is presented. A convergent and conservative finite volume scheme for an elliptic equation is proposed, together with the finite difference schemes, upwind and MUSCL, for a hyperbolic equation on grids with local refinement. Hence, an IMPES method is applied in an adaptive composite grid to track the front of a moving solution. An object-oriented programmation technique is used. The computational results for different examples illustrate the efficiency of the proposed algorithm. © 1997 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 13: 673–697, 1997     

多孔介质中三维混溶驱动问题在三角剖分复合网格上的有限差分格式

针对多孔介质中不可压缩流体的混溶驱动问题,基于平衡方程,利用有限体积方法建立了其三维问题在三角剖分单元中心空间局部加密复合网格上的有限差分格式,分析了差分格式的稳定性和收敛性,得到了关于饱和度的能量模误差估计,最后给出了数值算例.     

Degenerate two-phase incompressible flow III. Sharp error estimates

Summary. This is the third paper of a series in which we analyze mathematical properties and develop numerical methods for a degenerate elliptic-parabolic partial differential system which describes the flow of two incompressible, immiscible fluids in porous media. In this paper we consider a finite element approximation for this system. The elliptic equation for the pressure and velocity is approximated by a mixed finite element method, while the degenerate parabolic equation for the saturation is approximated by a Galerkin finite element method. A fully discrete approximation is analyzed. Sharp error estimates in energy norms are obtained for this approximation. The error analysis does not use any regularization of the saturation equation; the error estimates are derived directly from the degenerate equation. Also, the analysis does not impose any restriction on the nature of degeneracy. Finally, it respects the minimal regularity on the solution of the differential system. Received March 9, 1998 / Revised version received July 17, 2000 / Published online May 30, 2001     

Numerical simulations of water–gas flow in heterogeneous porous media with discontinuous capillary pressures by the concept of global pressure

We present an approach and numerical results for a new formulation modeling immiscible compressible two-phase flow in heterogeneous porous media with discontinuous capillary pressures. The main feature of this model is the introduction of a new global pressure, and it is fully equivalent to the original equations. The resulting equations are written in a fractional flow formulation and lead to a coupled degenerate system which consists of a nonlinear parabolic (the global pressure) equation and a nonlinear diffusion–convection one (the saturation equation) with nonlinear transmission conditions at the interfaces that separate different media. The resulting system is discretized using a vertex-centred finite volume method combined with pressure and flux interface conditions for the treatment of heterogeneities. An implicit Euler approach is used for time discretization. A Godunov-type method is used to treat the convection terms, and the diffusion terms are discretized by piecewise linear conforming finite elements. We present numerical simulations for three one-dimensional benchmark tests to demonstrate the ability of the method to approximate solutions of water–gas equations efficiently and accurately in nuclear underground waste disposal situations.     

多孔介质中可压缩可混溶驱动问题的Potempa混合元方法

Potempa-mixed finite element methods are considered in this paper for computing the compressible miscible displacement in porous media. The concentration equation is treated by Potempa‘s scheme,while the pressure equation is treated by a mixed finite element procedure. A H^1 error estimate for concentration with L2 for velocity is derived.The methods can provide the mass conservation and possesses minimal grid orientation.     

Taylor dispersion in a packed tube

Presented in this paper is a theoretical analysis for longitudinal scalar spread of mean concentration under a fully developed flow in a tube packed with porous media. A general form of momentum equation for superficial flow in porous media is introduced as a combination of the Navier–Stokes equation and Darcy’s law plus a superficial dispersion term due to phase discontinuity between the fluid flow and solid frame. The analytical solution presented for the fully developed superficial flow includes that for the Poiseulle flow in an evacuated tube as a limiting case. As an extension of Taylor’s classical work on dispersion of soluble matter in solvent flowing slowly through an evacuated tube, a one-dimensional dispersion equation valid for overall environmental assessment of contaminant is rigorously derived by cross-sectionally averaging the superficial mass equation and introducing a closure relation for a new unknown out of the averaging procedure, and corresponding Taylor dispersivity determined is shown to be a generalization of Taylor’s well-known result for the Poiseulle flow.     

不可压缩核废料污染问题的变网格特征有限元法

1 引言 多孔介质中的核废料污染问题是环境保护领域的重要课题。对于不可压缩二维模型,它是地层中迁移型耦合抛物型方程组的初边值问题:     

Three-dimensional finite/infinite elements analysis of fluid flow in porous media

Fluid flow in naturally fractured porous media can always be regarded as an unbounded domain problem and be better solved by finite/infinite elements. In this paper, a three-dimensional two-direction mapped infinite element is generated and combined with conventional finite elements and one direction infinite element to simulate poroelasticity. Therefore, the entire semi-infinite domain can be included in the numerical analysis. Both single- and dual-porosity porous media are considered. For the purpose of validation, we compare the results of finite/infinite elements with those of finite elements under two extreme boundary conditions. The comparison indicated that mapped infinite element is an appropriate approach to model fluid flow in porous media and provides an intermediate solution.     

Error estimates for an Euler implicit-mixed finite element scheme for reactive transport in saturated/unsaturated soil

We present convergence results for a fully discrete scheme based on the mixed finite element (MFE) method and an one-step Euler implicit (EI) method for simulating reactive solute transport in saturated/unsaturated soil. The results considered the low regularity of the solution of the degenerate parabolic equation describing the water flow in porous media. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Steady flow in a deformable porous medium

A nonlinear model for a steady flow in a deformable porous medium is considered. The flow is governed by the poroelasticity system consisting of an elasticity equation for the displacement of the porous medium and Darcy's equation for the pressure in the fluid. This poroelasticity system is nonlinear when the permeability in Darcy's equation is assumed to depend on the dilatation of the porous medium. Existence and uniqueness of a weak solution of this poroelasticity system is established under rather weak assumptions on the regularity of the data. Convergence of a finite element approximation is proved and verified through numerical experiments. Copyright © 2013 John Wiley & Sons, Ltd.     

Convergence analysis of an approximation to miscible fluid flows in porous media by combining mixed finite element and finite volume methods

This article deals with development and analysis of a numerical method for a coupled system describing miscible displacement of one incompressible fluid by another through heterogeneous porous media. A mixed finite element (MFE) method is employed to discretize the Darcy flow equation combined with a conservative finite volume (FV) method on unstructured grids for the concentration equation. It is shown that the FV scheme satisfies a discrete maximum principle. We derive L^∞ and BV estimates under an appropriate CFL condition. Then we prove convergence of the approximate solutions to a weak solution of the coupled system. Numerical results are presented to see the performance of the method in two space dimensions. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A numerical scheme for multi-point special boundary-value problems and application to fluid flow through porous layers

In this paper we propose a numerical scheme based on finite differences for the numerical solution of nonlinear multi-point boundary-value problems over adjacent domains. In each subdomain the solution is governed by a different equation. The solutions are required to be smooth across the interface nodes. The approach is based on using finite difference approximation of the derivatives at the interface nodes. Smoothness across the interface nodes is imposed to produce an algebraic system of nonlinear equations. A modified multi-dimensional Newton’s method is proposed for solving the nonlinear system. The accuracy of the proposed scheme is validated by examples whose exact solutions are known. The proposed scheme is applied to solve for the velocity profile of fluid flow through multilayer porous media.     

Modeling non-Darcian flow and solute transport in porous media with the Caputo–Fabrizio derivative

In this study, the non-Darcian flow and solute transport in porous media are modeled with a revised Caputo derivative called the Caputo–Fabrizio fractional derivative. The fractional Swartzendruber model is proposed for the non-Darcian flow in porous media. Furthermore, the normal diffusion equation is converted into a fractional diffusion equation in order to describe the diffusive transport in porous media. The proposed Caputo–Fabrizio fractional derivative models are addressed analytically by applying the Laplace transform method. Sensitivity analyses were performed for the proposed models according to the fractional derivative order. The fractional Swartzendruber model was validated based on experimental data for water flows in soil–rock mixtures. In addition , the fractional diffusion model was illustrated by fitting experimental data obtained for fluid flows and chloride transport in porous media. Both of the proposed fractional derivative models were highly consistent with the experimental results.     

Existence of solutions to various rock types odel model of two-phase flow in porous media

We study the flow of two immiscible and incompressible fluids through a porous media c,onsisting of different rock types: capillary pressure and relative permeablities curves are different in each type of porous media. This process can be formulated as a coupled system of partial differential equations which includes an elliptic pressurevelocity equation and a nonlinear degenerated parabolic saturation equation. Moreover the transmission conditions are nonlinear and the saturation is discontinuous at interfaces separating different media. A change of unknown leads to a new formulation of this problem. We derive a weak form for this new problem, which is a combination of a mixed formulation for the elliptic pressure-velocity equation and a standard variational formulation for the new parabolic equation. Under some realistic assumptions, we prove the existence of weak solutions to the implicit system given by time discretization.     

Solution of double nonlinear problems in porous media by a combined finite volume–finite element algorithm

The combined finite volume–finite element scheme for a double nonlinear parabolic convection-dominated diffusion equation which models the variably saturated flow and contaminant transport problems in porous media is extended. Whereas the convection is approximated by a finite volume method (Multi-Point Flux Approximation), the diffusion is approximated by a finite element method. The scheme is fully implicit and involves a relaxation-regularized algorithm. Due to monotonicity and conservation properties of the approximated scheme and in view of the compactness theorem we show the convergence of the numerical scheme to the weak solution. Our scheme is applied for computing two dimensional examples with different degrees of complexity. The numerical results demonstrate that the proposed scheme gives good performance in convergence and accuracy.     

不可压混流驱动问题的最小二乘混合元方法

1引言考虑多孔介质中两相不可压缩可混溶渗流驱动问题,它是由一组非线性耦合的椭园型压力方程和抛物型浓度方程组成：dVV。—一山人V什）gVV却）一q,VEn,（．1）＆,,。＿．、。。—一。x）＿＋u·grade－dlv（D（u）grade）一（1－c）q－,xEn,tEJ,（1．2）＆”－－’”””‘”－”””——－’——,、—’一其中a（）一a（x,c）一是（x）／卢（c）,J一［0,Ti,DcyR‘为水平油藏区域．方程式（1．l）一（1．2）中各物理量的意义如下：广为流体压力,c为流体的浓度,u为流体的Darer速度,叶为源汇项,／一—。x（q,O）,…