On a fully nonlinear degenerate parabolic system modeling immiscible gas–water displacement in porous media

In this paper, we are interested in the simultaneous flow of two immiscible fluid phases within a porous medium. We consider a two-phase flow model where the fluids are immiscible and there is no mass transfer between the phases. The medium is saturated by compressible/incompressible phase flows. We study the gas–water displacement without simplified assumptions on the state law of gas density. We establish an existence result for the nonlinear degenerate parabolic system based on new energy estimate on pressures.     

Stability of magnetic fluid motions in a saturated porous medium

We examine the asymptotic stability of both equilibrium and arbitrary basic flows of a magnetic fluid saturated in a porous medium. In both cases, we determine the stability bounds and determine the conditions when these flows are asymptotically stable. We also establish the uniqueness for an initial boundary value problem of magnetic fluids in the porous medium.     

Solvability of problems with degeneration: Imbibition of fluid in a porous medium

For a degenerate system of equations such as the equations of motion of immiscible fluids in porous media, we study the solvability of an initial–boundary value problem. Using the process of capillary imbibition of a wetting fluid as an example, we study a class of self-similar solutions with degeneration on the movable boundary and on the entry into the porous layer. The considered problem can be reduced to the analysis of properties of a nonlinear operator equation. For the classical solution of the original problem, we prove existence and uniqueness theorems.     

Simulation of uncompressible fluid flow through a porous media

Recently, a great interest has been focused for investigations about transport phenomena in disordered systems. One of the most treated topics is fluid flow through anisotropic materials due to the importance in many industrial processes like fluid flow in filters, membranes, walls, oil reservoirs, etc. In this work is described the formulation of a 2D mathematical model to simulate the fluid flow behavior through a porous media (PM) based on the solution of the continuity equation as a function of the Darcy’s law for a percolation system; which was reproduced using computational techniques reproduced using a random distribution of the porous media properties (porosity, permeability and saturation). The model displays the filling of a partially saturated porous media with a new injected fluid showing the non-defined advance front and dispersion of fluids phenomena.     

Slightly compressible and immiscible two-phase flow in porous media

We describe the flow of two compressible phases in a porous medium. We consider the case of slightly compressible phases for which the density of each phase follows an exponential law with a small compressibility factor. A nonlinear parabolic system including quadratic velocity terms is derived to describe compressible and immiscible two-phase flow in porous media. In one-dimensional space, we establish the existence and uniqueness of a local strong solution for the regularized system. We show also that the saturation is physically admissible. We describe the asymptotic behavior of the solutions when the compressibility factor goes to zero.     

弹性固体和充满两种互不相溶粘性流体的多孔固体之间界面的反射和折射及其衰减波

在充满两种互不相溶粘性流体的多孔固体中,研究弹性波的传播.用3个数性的势函数描述3个纵波的传播,用1个矢性的势函数单独描述横波的传播.根据这些势函数,在不同的组合相中,定义出质点的位移.可以看出,可能存在3个纵波和1个横波.在一个弹性固体半空间与一个充满两种互不相溶粘性流体的多孔固体半空间之间,研究其界面上入射纵波和横波所引起的反射和折射现象.由于孔隙流体中有粘性,折射到多孔介质中的波,朝垂直界面方向偏离.将入射波引起的反射波和折射波的波幅比,作为非奇异的线性代数方程组计算.进一步通过这些波幅比,计算出各个被离散波在入射波能量中所占的份额.通过一个特殊的数值模型,计算出波幅比和能量比系数随入射角的变化.超过SV波的临界入射角,反射波P将不再出现.越过界面的能量守恒原理得到了验证.绘出了图形并对不同孔隙饱和度以及频率的变化,讨论它们对能量分配的影响.     

不可压饱和多孔弹性杆动力响应的多辛方法

研究了不可压饱和多孔弹性杆的一维动力响应问题.基于多孔介质理论,在流相和固相微观不可压、固相骨架小变形的假定下,建立了不可压流体饱和多孔弹性杆一维轴向动力响应的数学模型.利用Hamilton空间体系的多辛理论,构造了不可压饱和多孔弹性杆轴向振动方程的多辛形式及其多种局部守恒律.采用中点Box离散方法得到轴向振动方程的多辛离散格式和局部能量守恒律以及局部动量守恒律的离散格式;数值模拟了不可压饱和多孔弹性杆的轴向振动过程,记录了每一时间步的局部能量数值误差和局部动量数值误差.结果表明,已构造的多辛离散格式具有很高的精确性和较长时间的数值稳定性,这为解决饱和多孔介质的动力响应问题提供了新的途径.     

A non-equilibrium model of a porous medium saturated with immiscible fluids

A model of a thermoelastic porous medium, saturated with two immiscible fluids, is considered. It is assumed that there are no phase transitions, the contribution of pulsations to the stress and kinetic energy is small, and that the components of the medium are in thermodynamic equilibrium. The non-equilibrium of the state, related to the finite time of redistribution of the fluids among the pores of the channels due to the presence of surface forces, is taken into account. A general form of the governing relations, necessary and sufficient to satisfy the principles of thermodynamic compatibility and independence of the choice of system of coordinates, is obtained. It is shown that the establishment of equilibrium is accompanied by dissipation due to capillary forces, which does not lead to seepage dissipation or thermal dissipation. For the case when the deformation of the skeleton and the deviation of the mean porous pressure and the temperature from the initial values are small, while the saturation and the non-equilibrium parameter undergo finite changes, an approximation of the potential of the skeleton is proposed in the form of a quadratic expansion in small parameters. A feature of the expansion is the presence of an initial value of the potential, which depends on the saturation and non-equilibrium. The relationship between the thermodynamic potential and the non-equilibrium kinetics, related to the requirement that the dissipation by the capillary forces should be non-negative, is determined. A generalized Darcy's law is formulated, which takes cross terms into account. It is shown that the proposed approximations enable key effects, which accompany the motion of immiscible fluids in a porous medium, to be described.     

Control of displacement front in a model of immiscible two-phase flow in porous media

For the Buckley–Leverett equation describing the flow of two immiscible fluids in porous media, an exact parametric representation of the solution is constructed with the help of the Bäcklund transformation. As a result, the advance of the displacement front can be controlled to a high degree of accuracy. The method is illustrated using an example of a typical oil well with actual parameters.     

On a front evolution in porous media with a source–analysis and numerics

We analyze evolution of the interface between immiscible fluids of different densities in porous media. The fluids can be compressible (CO₂ or natural gases) or incompressible (oil, water). We rigorously prove that, if the heavier fluid is on the top and there are no sink or source, a tip of the interface will move in the direction of the gravity (if the tip is directed toward the bottom) or the buoyancy (if the tip is directed toward the top).We also showhowthe sink/source influence propagation of the interface and provide numerical examples.     

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.     

A hybrid numerical model for multiphase fluid flow in a deformable porous medium

In this paper, a fully coupled finite volume-finite element model for a deforming porous medium interacting with the flow of two immiscible pore fluids is presented. The basic equations describing the system are derived based on the averaging theory. Applying the standard Galerkin finite element method to solve this system of partial differential equations does not conserve mass locally. A non-conservative method may cause some accuracy and stability problems. The control volume based finite element technique that satisfies local mass conservation of the flow equations can be an appropriate alternative. Full coupling of control volume based finite element and the standard finite element techniques to solve the multiphase flow and geomechanical equilibrium equations is the main goal of this paper. The accuracy and efficiency of the method are verified by studying several examples for which analytical or numerical solutions are available. The effect of mesh orientation is investigated by simulating a benchmark water-flooding problem. A representative example is also presented to demonstrate the capability of the model to simulate the behavior in heterogeneous porous media.     

Flow of oil and water in a porous medium

The flow of two immiscible and incompressible fluids in a porous medium is described by a system of quasilinear degenerate partial differential equations. In this paper the existence of a weak solution by regularization is shown.     

Fluid flow,relative permeabilities and capillary pressure curves through heterogeneous porous media

There is a growing interest in fluid flow through rocks which show a large range, or heterogeneity, in pore sizes (from order tens of nanometres up to a hundred microns) which lead to small porosities and permeabilities. Digital imagery of these rocks yield under-resolved pore structures or a gray-scale sample. In this work we use the gray-scale lattice Boltzmann (LB) method to determine flow of two, immiscible fluids through such three-dimensional porous media. The method yields transient and equilibrium flow fields from which we determine relative permeability and capillary pressure curves. The method is applied to real digital rock data to extract flow properties through the sample. Both single phase and multiphase LB simulations are carried out. We find low permeabilities (of the order of milli-Darcy) together with large residual saturations (between 50–75%), indicative of a rock sample which has void regions remote from each other, i.e. either disconnected or connected by very narrow pore throats.     

On the Two‐Dimensional Muskat Problem with Monotone Large Initial Data

We consider the evolution of two incompressible, immiscible fluids with different densities in porous media, known as the Muskat problem [21], which in two dimensions is analogous to the Hele‐Shaw cell [24]. We establish, for a class of large and monotone initial data, the global existence of weak solutions. The proof is based on a local well‐posedness result for the initial data with certain specific asymptotics at spatial infinity and a new maximum principle for the first derivative of the graph function.© 2016 Wiley Periodicals, Inc.     

Viscous-inviscid model for the linear stability of core-annular flow

In this paper the linear stability of two immiscible fluids of widely different kinematic viscosity and equal density flowing through a circular pipe is analysed. A viscous-inviscid model is used which offers a consistent zeroth-order approximation to the stability problem as long as the thickness of the ring flow, where the inviscid fluid is located, is large enough. In this way the laminar sublayer at the pipe wall does not interact with the fluid interface. A closed form expression for the complex dispersion relation due to an arbitrary wavelength perturbation is derived, which determines a stability criterion, and then simplified for large and short wavelength values.     

Uncoupling evolutionary groundwater‐surface water flows using the Crank–Nicolson Leapfrog method

Consider an incompressible fluid in a region Ω_f flowing both ways across an interface into a porous media domain Ω_p saturated with the same fluid. The physical processes in each domain have been well studied and are described by the Stokes equations in the fluid region and the Darcy equations in the porous media region. Taking the interfacial conditions into account produces a system with an exactly skew symmetric coupling. Spatial discretization by finite element method and time discretization by Crank–Nicolson LeapFrog give a second‐order partitioned method requiring only one Stokes and one Darcy subphysics and subdomain solver per time step for the fully evolutionary Stokes‐Darcy problem. Analysis of this method leads to a time step condition sufficient for stability and convergence. Numerical tests verify predicted rates of convergence; however, stability tests reveal the problem of growth of numerical noise in unstable modes in some cases. In such instances, the addition of time filters adds stability. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

On the Rayleigh-Taylor Instability for Two Uniform Viscous Incompressible Flows

The authors study the Rayleigh-Taylor instability for two incompressible immis- cible fluids with or without surface tension, evolving with a free interface in the presence of a uniform gravitational field in Eulerian coordinates. To deal with the free surface, instead of using the transformation to Lagrangian coordinates, the perturbed equations in Eule- rian coordinates are transformed to an integral form and the two-fluid flow is formulated as a single-fluid flow in a fixed domain, thus offering an alternative approach to deal with the jump conditions at the free interface. First, the linearized problem around the steady state which describes a denser immiscible fluid lying above a light one with a free interface separating the two fluids, both fluids being in （unstable） equilibrium is analyzed. By a general method of studying a family of modes, the smooth （when restricted to each fluid domain） solutions to the linearized problem that grow exponentially fast in time in Sobolev spaces are constructed, thus leading to a global instability result for the linearized problem. Then, by using these pathological solutions, the global instability for the corresponding nonlinear problem in an appropriate sense is demonstrated.