Material instabilities in inelastic saturated porous media under dynamic loadings

Instabilities in inelastic saturated porous media are investigated here for general three-dimensional states under dynamic loadings using a perturbation approach.Under quasi-static conditions, unbounded growth of perturbations is related to the emergence of stationary discontinuities under drained or undrained conditions, while under dynamic conditions, unbounded growth is related either to the emergence of stationary discontinuities (and these are set by drained conditions) or to the appearance of the flutter phenomenon (acceleration waves).For associative behaviour unbounded growth always corresponds to localization under drained conditions and the onset of growth of perturbations occurs here only through divergence growth. It is only for non-associative flow that unbounded growth may correspond to undrained localization in quasi-static conditions and to flutter under dynamic conditions.     

Material coefficients of multiphase thermoporoelasticity for anisotropic micro-heterogeneous porous media

Constitutive equations of thermoporoelasticity for two-phase fluid system derived by Coussy (2007) are extended to those for multiphase system with anisotropic micro-heterogeneous porous media. All the material coefficients are found to be expressed by coefficients common to saturated thermoporoelasticity and saturation-dependent effective stress coefficients. Experimental conditions to evaluate thermal pore expansivity and saturation dependent effective stress coefficients are also presented. Material coefficients for orthotropy and transverse isotropy, which are often seen in actual porous media, are explicitly listed to facilitate application.     

Thawing in saturated porous media

A mathematical model for thawing in a saturated porous medium is considered. The well-posedness of the corresponding mathematical problem is proved and similarity solutions are found.

---

Sommario Si considera un modello matematico per to scongelamento in un mezzo poroso saturo. Viene dimostrata la buona posizione del corrispondente problema matematico e si trovano soluzioni di similarità.

     

Simulation of multiphase flows in porous media

The ability to numerically simulate single phase and multiphase flow of fluids in porous media is extremely important in developing an understanding of the complex phenomena governing the flow. The flow is complicated by the presence of heterogeneities in the reservoir and by phenomena such as diffusion, dispersion, and viscous fingering. These effects must be modeled by terms in coupled systems of nonlinear partial differential equations which form the basis of the simulator. The simulator must be able to handle both single and multiphase flows and the transition regimes between the two. A discussion of some of the aspects of modeling dispersion and viscous fingering is presented along with directions for future work.The partial differential equation models are convection-dominated and contain important local effects. An operator-splitting technique is used to address these different effects accurately. Convection is treated by time stepping along the characteristics of the associated pure convection problem, and diffusion is modeled via a Galerkin method for single phase flow and a Petrov-Galerkin technique for multiphase regimes. ELLAM (Eulerian-Lagrangian Localized Adjoint Methods) are discussed to effectively treat the advection-dominated processes. Accurate approximations of the fluid velocities needed in the Eulerian-Lagrangian time-stepping procedure are obtained by mixed finite element methods. Adaptive local grid refinement techniques are then indicated to resolve important local phenomena around wells and large heterogeneities or to resolve the moving internal boundary layers which often govern the mass transfer between phases.     

Diffusion in anisotropic porous media

An experimental system was constructed in order to measure the two distinct components of the effective diffusivity tensor in transversely isotropic, unconsolidated porous media. Measurements were made for porous media consisting of glass spheres, mica particles, and disks made from mylar sheets. Both the particle geometry and the void fraction of the porous media were determined experimentally, and theoretical calculations for the two components of the effective diffusivity tensor were carried out. The comparison between theory and experiment clearly indicates that the void fraction and particle geometry are insufficient to characterize the process of diffusion in anisotropic porous media. Roman Letters A interfacial area between - and -phases for the macroscopic system, m² - A _e area of entrances and exits of the -phase for the macroscopic system, m² - A interfacial area contained within the averaging volume, m² - a characteristic length of a particle, m - b average thickness of a particle, m - c _A concentration of species A, moles/m³ - c _o reference concentration of species A, moles/m³ - c _A intrinsic phase average concentration of species A, moles/m³ - c _a c _A–c _A, spatial deviation concentration of species A, moles/m³ - C c _A/c ₀, dimensionless concentration of species A - binary molecular diffusion coefficient, m²/s - D _eff effective diffusivity tensor, m²/s - D _xx component of the effective diffusivity tensor associated with diffusion parallel to the bedding plane, m²/s - D _yy component of the effective diffusivity tensor associated with diffusion perpendicular to the bedding plane, m²/s - D _eff effective diffusivity for isotropic systems, m₂/s - f vector field that maps c _A on to c _a , m - h depth of the mixing chamber, m     

Waves in partially saturated porous media

The propagation of compressional waves in a porous medium is investigated in case the pore liquid contains a small volume fraction of gas. The effect of oscillating gas bubbles is taken into account by introducing a frequency-dependent fluid bulk modulus, which is incorporated in the Biot theory. Using a shock tube technique, new experimental data are obtained for a porous column subjected to a pressure step wave. An oscillatory behaviour is observed, consisting of two distinct frequency bands, which is predicted by the theoretical analysis.     

Linear waves in saturated porous media

The peculiarity of linear and nonlinear wave propagation in porous media saturated with liquid or gas has been investigated by the methods of multiphase media mechanics. It has been shown that for the analysis and interpretation of experimental data, it is expedient to build models taking into account the nonstationary powers of interaction between the solid and liquid phases and the viscouselastic behaviour of the porous media skeleton. Inertia and inertia-viscous powers principally influence wave attentuation in porous media. Two interphase mechanisms of momentum transfer (two stress tensors — in the solid phase and liquid) lead to two types of waves. Attenuation is determined not only by interphase friction, but also by dissipation resulting from intergrain friction in the solid phase, the influence of which multiplexly exceeds the liquid viscosity influence. The real decrement of attenuation may exceed the sphere restricted by the limiting curves corresponding to the frozen and equilibrium schemes of intergrain deformation. The attenuation of momentum perturbation has been studied. The method of discrete Fourier transform has been used. The analysis of experimental data contained in the literature and their comparison with the results of calculations has been carried out.     

Fundamentals of thermomechanics of porous saturated media

Institute of Mechanics, Academy, of Sciences of the Ukrainian SSR, Kiev. Translated from Prikladnaya Mekhanika, Vol. 24, No. 4, pp. 3–13, April, 1988.     

Models and problems for saturated porous media

A multiphase model of a liquid- and gas-saturated porous medium is proposed. The model takes into account the finite deformations of a skeleton, arbitrary flows of liquids and gases, and the phase mass transfer between the skeleton and the liquids. Some basic relations are given for the corresponding boundary value problems in the cases of arbitrary and small motions. The constitutive relations describing the properties of the skeleton resistance to deformation (skeleton stresses) and the mutual resistance of components (internal interaction) are considered in detail. The internal actions exerted on the skeleton by a moving liquid (gas) are discussed when these actions may take the form of drag forces, lift (displacing) forces, overturning moments, and rotational (screw) moments.     

Shock waves in saturated thermoelastic porous media

The objective of this paper is to develop and present the macroscopic motion equations for the fluid and the solid matrix, in the case of a saturated porous medium, in the form of coupled, nonlinear wave equations for the fluid and solid velocities. The nonlinearity in the equations enables the generation of shock waves. The complete set of equations required for determining phase velocities in the case of a thermoelastic solid matrix, includes also the energy balance equation for the porous medium as a whole, as well as mass balance equations for the two phase. In the special case of a rigid solid matrix, the wave after an abrupt change in pressure propagates only through the fluid.     

Nonstationary filtration in partially saturated porous media

From the mathematical formulation of a one-dimensional flow through a partially saturated porous medium, we arrive at a nonlinear free boundary problem, the boundary being between the saturated and the unsaturated regions in the medium. In particular we obtain an equation which is parabolic in the unsaturated part of the domain and elliptic in the saturated part.Existence, uniqueness, a maximum principle and regularity properties are proved for weak solutions of a Cauchy-Dirichlet problem in the cylinder {(x,t): 0x1, t0} and the nature, in particular the regularity, of the free boundary is discussed.Finally, it is shown that solutions of a large class of Cauchy-Dirichlet problems converge towards a stationary solution as t and estimates are given for the rate of convergence.     

Advective fluxes in multiphase porous media under nonisothermal conditions

We consider the case in which more than one fluid phase occupies the void space of a porous medium. The advective flux law is formulated for a fluid phase, under nonisothermal conditions and with the presence of solutes in the fluid phases. The derivation of the flux laws is based on an approximated version of the averaged balance equation for linear momentum. Taking into account momentum transfer through the interface between the fluid phases, leads to coupling between the flow in adjacent phases. Fluxes are also shown to depend on the surface tension at the interface between the adjacent fluid phases. Since the latter depends on temperature and solute concentration in the two phases, the advective flux is shown to depend on both temperature and solute concentration gradients in the two phases. A preliminary order of magnitude analysis gives conditions under which the coupling phenomena are not negligible. The approach is applied to the unsaturated zone, as a typical example of a multiphase porous medium.The main conclusion is that the well known Darcy law for single phase flow, may have to be modified for a multi fluid phase system, especially when temperature and solute concentration are not uniform.     

流体饱和多孔介质黏弹性动力人工边界

基于Biot流体饱和多孔介质本构方程,分别考察具有辐射阻尼性质的外行柱面波和球 面波在圆柱面和球面人工边界上引起的法向、切向应力的表达式. 在应力表达形式上,固相 介质和孔隙流体的法向和切向应力都是由两项组成,它们分别与质点的位移和速度成正比, 因此,可在人工边界的法向和切向设置连续分布的并联弹簧------黏滞阻尼器,用来模拟人工边 界以外的无限域介质对来自有限计算域的外行波动的能量吸收作用,从而形成了流体饱和多 孔介质的黏弹性动力人工边界. 流体饱和多孔介质的黏弹性动力人工边界可方便地与大型通 用软件结合,用于分析饱和土中复杂的结构-地基动力相互作用问题. 算例表明流体饱和多 孔介质黏弹性动力人工边界具有较好的精度和稳定性.     

Rapoport-leas two-phase flow model for anisotropic porous media

The Rapoport-Leas mathematical model of two-phase flow is generalized to include the case of anisotropic porous media. The formula for the capillary pressure, which specifies the relationship between the phase pressures, contains a scalar function of a vector argument. In order to determine the scalar function, the capillary pressure tensor and the tensor inverse to the tensor of characteristic linear dimensions are introduced. The capillary pressure is determined by the contraction of the second-rank tensors with a unit vector collinear to the phase pressure gradients, also assumed to be collinear. It is shown that the saturation function introduced for isotropic porous media (Leverett function) can be generalized to include anisotropic media and is now determined by a fourth-rank tensor. Generalized expressions for the Leverett and relative phase permeability functions are given for orthotropic and transversely isotropic media with account for the hysteresis of the phase permeabilities and capillary pressure.     

An improved numerical simulator for multiphase flow in porous media

In basin modelling the thermodynamics of a multicomponent multiphase fluid flux are computationally too expensive when derived from an equation of state and the Gibbs equality constraints. In this article we present a novel implicit molar mass formulation technique using binary mixture thermodynamics. The two proposed solution methods, with and without cross derivative terms between components, are based on a preconditioned Newton‐GMRES scheme for each time‐step with analytical computation of the derivatives. These new algorithms reduce significantly the numerical effort for the computation of the molar masses, and we illustrate the behavior of these methods with numerical computations. Copyright © 2004 John Wiley & Sons Ltd.     

Determining equations of two-phase flows through anisotropic porous media

A new interpretation of the concept of relative phase permeability is given. Relative phase permeabilities are represented in the form of fourth-rank tensors. It is shown that in the case of anisotropic porous media functions depending not only on the saturation but also on the anisotropy parameters represented in the form of ratios of the principal values of the absolute permeability coefficient tensor correspond to the classical representation of the relative phase permeabilities. For a two-phase flow in anisotropic porous media with orthotropic and transversely-isotropic symmetry a generalized two-term Darcy’s law is analyzed. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 87–94, March–April, 1998. The work was carried out with support from the Russian Foundation for Fundamental Research (project No. 96-01-00623).