Reflection and Transmission of Elastic Waves from the Interface of a Fluid-saturated Porous Solid and a Double Porosity Solid

The reflection and transmission characteristics of an incident plane P1 wave from the interface of a fluid-saturated single porous solid and a fluid-saturated double porosity solid are investigated. The fluid-saturated porous solid is modeled with the classic Biot’s theory and the double porosity medium is described by an extended Biot’s theory. In a double-porosity model with dual-permeability there exist three compressional waves and a shear wave. The effects of the incident angle and frequency on amplitude ratios of the reflected and transmitted waves to the incident wave are discussed. Two boundary conditions are discussed in detail: (a) Open-pore boundary and (b) Sealed-pore boundary. Numerical results reveal that the characteristics of the reflection and transmission coefficients to the incident angle and the frequency are quite different for the two cases of boundary conditions. Properties of the bulk waves existing in the fluid-saturated porous solid and the double porosity medium are also studied.     

Reflection and Transmission of Elastic Waves at the Interface Between Water and a Double Porosity Solid

The reflection and transmission of an incident plane wave at an interface between water and a fluid-saturated double porosity solid are investigated. The properties of the three compressional waves and one shear wave in double porosity media are discussed in detail. The effect of the viscosity, permeability, and porosity on the phase velocity and attenuation of four bulk waves is presented. Comparison between the double porosity theory and the Biot theory reveals that there exists obvious difference in the phase velocity, attenuation and the reflection and transmission coefficients. Two cases of boundary conditions, i.e., the sealed-pore and the open-pore boundaries, are introduced in the numerical calculation. At last, the effect of the fracture permeability and porosity on the reflection and transmission coefficients considered. It is revealed that the amplitudes of the reflected and transmitted waves strongly depend the boundary condition, fracture permeability, and porosity.     

Modeling Fluid Flow in Fractured Porous Media with the Interfacial Conditions Between Porous Medium and Fracture

One of the most popular models that has been applied to predict the fluid velocity inside the fracture with impermeable walls is the cubic law. It highlights that the mean flux along the fracture is proportional to the cubic of fracture aperture. However, for a fractured porous medium, the normal and tangential interface conditions between the fracture and porous matrix can change the velocity profile inside the fracture. In this paper, a correction factor is introduced for flow equation along the fracture by imposing the continuity of normal and tangential components of velocity at the interface between the fracture and porous matrix. As a result, the mean velocity inside the fracture depends not only on the fracture aperture, but also on a set of non-dimensional numbers, including the matrix porosity, the ratio of intrinsic permeability of fracture to that of matrix, the wall Reynolds number, and the ratio of normal velocity on one wall to the other. Finally, the introduced correction factor is employed within the extended finite element method, which is widely used for numerical simulation of fluid flow within the fractured porous media. Several numerical results are presented for the fluid flow through a specimen containing single fracture, in order to investigate the deviation from the cubic law in different case studies.

     

Modeling Effective Interface Laws for Transport Phenomena Between an Unconfined Fluid and a Porous Medium Using Homogenization

In this chapter of the special issue of the journal “Transport in Porous Media,” on the topic “Flow and transport above permeable domains,” we present modeling of flow and transport above permeable domains using the homogenization method. Our goal is to develop a heuristic approach which can be used by the engineering community for treating this type of problems and which has a solid mathematical background. The rigorous mathematical justification of the presented results is given in the corresponding articles of the authors. The plan is as follows: We start with the section “Introduction” where we give an overview and comparison with interface conditions obtained using other approaches. In Sect. 2, we give a very short derivation of the Darcy law by homogenization, using the two-scale expansion in the typical pore size parameter ε. It gives us the definition of various auxiliary functions and typical effective properties as permeability. In Sect. 3, we introduce our approach to the effective interface laws on a simple 1D example. The approximation is obtained heuristically using the two steps strategy. For the 1D problem we calculate the approximation and the effective interface law explicitly and show that it is valid at order O(ε ²). Next, in Sect. 4 we give a derivation of the Beavers–Joseph–Saffman interface condition and of the pressure jump condition, using homogenization. We construct the corresponding boundary layer and present a heuristic calculation, leading to the interface law and being based on the rigorous mathematical result. In addition, we show the invariance of the law with respect to the small variations in the choice of the interface position. Finally, there is a short concluding section. The research of A.M. was partially supported by the GDR MOMAS (Modélisation Mathématique et Simulations numériques liées aux problèmes de gestion des déchets nucléaires) (PACEN/CNRS, ANDRA, BRGM, CEA, EDF, IRSN).     

Modeling of Heat and Mass Transfer in a Fractured Porous Medium

While fractured formations are possibly the most important contributors to the production of oil worldwide, modeling fractured formations with rigorous treatments has eluded reservoir engineers in the past. To date, one of the most commonly used fractured reservoir models remains the one that was suggested by Warren and Root nearly four decades ago. In this paper, a new model for fractures embedded in a porous medium is proposed. The model considers the Navier-Stokes equation in the fracture (channel flow) while using the Brinkman equation for the porous medium. Unlike the previous approach, the proposed model does not require the assumption of orthogonality of the fractures (sugar cube assumption) nor does it impose incorrect boundary conditions for the interface between the fracture and the porous medium. Also, the transfer coefficient between the fracture and matrix interface does not need to be specified, unlike the cases for which Darcy's law is used. In order to demonstrate the usefulness of the approach, a two-dimensional model of a fractured formation is developed and numerical simulation runs conducted.

The proposed model is derived through a series of finite element modeling runs for various cases using the Navier-Stokes equation in the channel while maintaining the Brinkman equation in the porous medium. Various cases studied include different fracture orientations, fracture frequencies, and thermal and solutal constraints. The usefulness of the proposed model in modeling complex formations is discussed. Finally, a series of numerical runs also provided validity of the proposed model for the cases in which thermal and solutal effects are important. Such a study of double diffusive phenomena, coupled with forced convection, in the context of fractured formations has not been reported before.     

Effects of Pressure Oscillations on Drainage in an Elastic Porous Medium

The effects of seismic stimulation on the flow of two immiscible fluids in an elastic synthetic porous medium is experimentally investigated. A wetting fluid is slowly evacuated from the medium, while a pressure oscillation is applied on the injected non- wetting fluid. The amplitude and frequency of the pressure oscillations as well as the evacuation speed are kept constant throughout an experiment. The resulting morphology of the invading structure is found to be strongly dependent on the interplay between the amplitude and the frequency of the applied pressure oscillations and the elasticity of the porous medium. Different combinations of these properties yield morphologically similar structures, allowing a classification of structures that is found to depend on a proposed dimensionless number.     

A Numerical Study of the Critical Gas Saturation in a Porous Medium

We use porenetwork simulations to study the dependence of the critical gas saturation in solutiongas drive processes on the geometric parameters of the porous medium. We show that for a variety of growth regimes (including global and local percolation, instantaneous and sequential nucleation, and masstransfer driven processes), the critical gas saturation, Sgc, follows a powerlaw scaling with the final nucleation fraction (fraction of sites activated), fq. For 3D processes, this relation reads Sgcfq0.16, indicating a sensitive dependence of Sgc to fq at very small values of fq.     

Reactive Solute Transport in a Nonstationary,Fractured Porous Medium: A Dual-Porosity Approach

A numerical method of moment is developed for solute flux through a nonstationary, fractured porous medium. Solute flux is described as a space-time process where time refers to the solute flux breakthrough and space refers to the transverse displacement distribution at a control plane. A first-order mass diffusion model is applied to describe interregional mass diffusion between fracture (advection) and matrix (nonadvection) regions. The chemical is under linear equilibrium sorption in both fracture and matrix regions. Hydraulic conductivity in the fracture region is assumed to be a spatial random variable. In this study, the general framework of Zhang et al.(2000) is adopted for solute flux in a nonstationary flow field. A time retention function related to physical and chemical sorption in the dual-porosity medium is developed and coupled with solute advection along random trajectories. The mean and variance of total solute flux are expressed in terms of the probability density function of the parcel travel time and transverse displacement. The influences of various factors on solute transport are investigated. These factors include the interregional mass diffusion rate between fracture and matrix regions, chemical sorption coefficients in both regions, water contents in both regions, and location of the solute source. In comparison with solute transport in a one-region medium, breakthrough curves of the mean and variance of the total solute flux in a two-region medium have lower peaks and longer tails. As compared with the classical stochastic studies on solute transport in fractured media, the numerical method of moment provides an approach for applying the stochastic method to study solute transport in more complicated fractured media.     

Experimental and Numerical Study of the Onset of Transient Natural Convection in a Fractured Porous Medium

Transport in Porous Media - In this study, analysis of transient natural convection in a porous medium with a vertical fracture across it in a rectangular enclosure was performed experimentally and...     

Properties of Elastic Waves in a Non-Newtonian (Maxwell) Fluid-Saturated Porous Medium

The present study investigates novelties brought into the classic Biot's theory of propagation of elastic waves in a fluid-saturated porous solid by inclusion of non-Newtonian effects that are important, for example, for hydrocarbons. Based on our previous results (Tsiklauri and Beresnev, 2001), we investigated the propagation of rotational and dilatational elastic waves by calculating their phase velocities and attenuation coefficients as a function of frequency. We found that the replacement of an ordinary Newtonian fluid by a Maxwell fluid in the fluid-saturated porous solid results in: (a) an overall increase of the phase velocities of both the rotational and dilatational waves. With the increase of frequency these quantities tend to a fixed, higher level, as compared to the Newtonian limiting case, which does not change with the decrease of the Deborah number . (b) The overall decrease of the attenuation coefficients of both the rotational and dilatational waves. With the increase of frequency these quantities tend to a progressively lower level, as compared to the Newtonian limiting case, as decreases. (c) Appearance of oscillations in all physical quantities in the deeply non-Newtonian regime.     

Fully Developed Mixed Convection Flow Between Inclined Parallel Plates Filled with a Porous Medium

A fully developed mixed convection flow between inclined parallel flat plates filled with a porous medium is considered through which there is a constant flow rate and with heat being supplied to the fluid by the same uniform heat flux on each plate. The equations governing this flow are made non-dimensional and are seen to depend on two dimensionless parameters, a mixed convection parameter λ and the Péclet number Pe, as well as the inclination γ of the plates to the horizontal. The velocity and temperature profiles are obtained in terms of λ, Pe and γ when the channel is inclined in an upwards direction as well as for horizontal channels. The limiting cases of small and large λ and small Pe are considered with boundary-layer structures being seen to develop on the plates for large values of λ.     

Flow from the Soil Surface into a Fractured Porous Aeration Zone

The one-dimensional problem of the contamination of a fractured porous aeration zone as a result of a fast spill of fluid over the soil surface is investigated. The block capillary imbibition rate is approximated with allowance for the experimental data. An analytic dependence describing the trajectory of the leading contamination front is obtained and the depth of penetration of the spill into the soil is found. The block contamination profile is determined.     

Static Deformations of a Linear Elastic Porous Body Filled with an Inviscid Fluid

We study infinitesimal deformations of a porous linear elastic body saturated with an inviscid fluid and subjected to conservative surface tractions. The gradient of the mass density of the solid phase is also taken as an independent kinematic variable and the corresponding higher-order stresses are considered. Balance laws and constitutive relations for finite deformations are reduced to those for infinitesimal deformations, and expressions for partial surface tractions acting on the solid and the fluid phases are derived. A boundary-value problem for a long hollow porous solid cylinder filled with an ideal fluid is solved, and the stability of the stressed reference configuration with respect to variations in the values of the coefficient coupling deformations of the two phases is investigated. An example of the problem studied is a cylindrical cavity leached out in salt formations for storing hydrocarbons.     

各向异性弹性体光滑接触界面上的滑移脉冲波

利用Stroh公式,Fourier分析和奇异积分方程技术研究了两各向异性弹性半空间光滑接触可分离界面上滑移脉冲波的存在及其传播特性。结果表明,如果至少能在一种介质中存在Rayleigh波,且其波速小于两种介质中的最小极限速度,则滑移脉冲波就可以存在。这种脉冲波传播速度不确定,可在最小极限波速与较低的Rayleigh波速之间取值,而该取值范围又取决于无界面分离情况下的第一、第二滑移波的解。分离区大小取决于扰动的强度,界面法向力和质点速度在分离区两端有 1 /2奇异性。     

Natural Convection Inside an Inclined Wavy Enclosure Filled with a Porous Medium

The Darcy flow model with the Boussinesq approximation is used to investigate numerically the natural convection inside an inclined wavy cavity filled with a porous medium. Finite Element Method is used to discretize the governing differential equations with non-staggered variable arrangement. Results are presented for and , where ϕ, Ra, A and λ correspond to the cavity inclination angle, Rayleigh number, aspect ratio and surface waviness parameter, respectively. Stream and isotherm lines representing the corresponding flow and thermal fields, and local and average Nusselt numbers distribution expressing the rate of heat transfer are determined and shown on graphs and tables. A good agreement is observed between the present results and those known from the open literature. The flow and thermal structures found to be highly dependent on surface waviness for inclination angles less than 45°, especially for high Rayleigh numbers.     

On the Elastic Field of a Shpherical Inhomogeneity with an Imperfectly Bonded Interface

This study is devoted to the development of a unified and explicit elastic solution to the problem of a spherical inhomogeneity with an imperfectly bonded interface. Both tangential and normal displacement discontinuities at the interface are considered and a linear interfacial condition, which assumes that the tangential and the normal displacement jumps are proportional to the associated tractions, is adopted. The elastic disturbance due to the presence of an imperfectly bonded inhomogeneity is decomposed into two parts: the first is formulated in terms of an equivalent nonuniform eigenstrain distributed over a perfectly bonded spherical inclusion, while the second is formulated in terms of an imaginary Somigliana dislocation field which models the interfacial sliding and normal separation. The exact form of the equivalent nonuniform eigenstrain and the imaginary Somigliana dislocation are fully determined in this paper.     

Numerical Study of a Magneto-fluid Motion Through a Porous Medium Between Two Wavy Plates with Constant Suction

In this note, the problem of an incompressible viscous fluid moving through a porous medium (Brinkman model) between two wavy plates under the effects of a constant inclined magnetic field that makes an angle with the vertical axis and constant suction, are studied numerically by a method related to the method of Takabatake and Ayukawa in 1982. The present approach is not restricted by any of the parameters appearing in the problem such as Reynolds number, magnetic parameter, suction parameter, the wave number and amplitude ratio. The variations in velocity, flow rate and pressure gradient with the above governing parameters are presented. Moreover, the effect of varying the porous medium and the inclined angle is also studied.     

Stress Concentration Near an Elliptic Crack in the Interface Between Elastic Bodies under Steady-State Vibrations

The paper addresses the three-dimensional problem on steady-state vibrations of an elastic body consisting of two perfectly joined dissimilar half-spaces with an elliptic mode I crack located in one of the half-spaces normally to the interface. The problem is reduced to a boundary integral equation for the crack opening function. The integration domain of the equation is bounded by the crack domain, and the interaction between the crack and the interface is described by a regular kernel. The equation is solved using the mapping method. Numerical results are obtained for the case where the surfaces of the elliptic crack are subjected to harmonic loading with constant amplitude. The dependences of the stress intensity factors on the wave number are presented for various relationships among the mechanical constants that ensure the absence of near-surface waves