Delineation of Microscale Regimes of Fully-Developed Drainage and Implications for Continuum Models     

Yanis C. Yortsos  Baomin Xu  Dominique Salin 《Computational Geosciences》2001,5(3):257-278

Despite significant progress made in recent years, a fundamental understanding of immiscible displacements at the macroscale is lacking. In this paper we use a version of percolation theory, based on invasion percolation in a gradient, to connect drainage processes at the pore-network scale with the displacement at the macroscale. When the mobility ratio M is sufficiently small, the displacement is stabilized and can be described by invasion percolation in a stabilizing gradient. In the opposite case, it has common features with invasion percolation in a destabilizing gradient. A diagram delineating the regimes of fully developed drainage is developed. The transition between stabilized displacement and fingering is controlled by the change of the sign of the gradient of the percolation probability, and the transition boundary is described by a scaling law involving the capillary number and the viscosity ratio. We review recent work for random networks and extend the method to correlated pore networks. As the regimes of stabilized displacement are also those for which conventional theories (such as the Buckley–Leverett equation) are expected to apply, the phase diagram helps to delineate their validity.  相似文献   

Analysis of Rayleigh waves in saturated porous elastic media by finite element method     

H. Hirai 《Soil Dynamics and Earthquake Engineering》1992,11(6)

A numerical procedure for the analysis of Rayleigh waves in saturated porous elastic media is proposed by use of the finite element method. The layer stiffness matrix, the layer mass matrix and the layer damping matrix in a layered system are presented for the discretized form of the solid-fluid equilibrium equation proposed by Biot. In order to consider the influence of the permeability coefficient on the behavior of Rayleigh waves, attention is focused on the following states: ‘drained’ state, ‘undrained’ state and the states between two extremes of ‘drained’ and ‘undrained’ states. It is found from computed results that the permeability coefficient exerts a significant effect on dispersion curves and displacement distributions of Rayleigh waves in saturated porous media.  相似文献   

Flow through anisotropic networks     

L. Rothenburg  A. Bogobowicz 《Stochastic Environmental Research and Risk Assessment (SERRA)》1990,4(3):175-191

This paper presents an analysis of Hagen-Poiseulle flow through plane random anisotropic networks of interconnected channels. Macroscopic permeability tensor of the network is expressed in terms of statistico-geometrical characteristics like the degree of anisotropy in channel orientations, average co-ordination number of the network and first two moments of channel length distribution. Analytical results are illustrated and verified using numerical analysis of flow in a simulated random network. The emphasis of the paper is on the effects of anisotropy on distributions of flow rates in channels. It is shown that, due to anisotropy the maximum flow rate generally occurs in channels that are not aligned along the direction of the macroscopic pressure gradient.  相似文献   

Stochastic analysis of multiphase flow in porous media: 1. Spectral/perturbation approach     

C. Chang  M. W. Kemblowski  J. Kaluarachchi  A. Abdin 《Stochastic Environmental Research and Risk Assessment (SERRA)》1995,9(3):239-267

Stochastic analysis of steady-state multiphase (water, oil, and air) flow in heterogeneous porous media was performed using the perturbation theory and spectral representation techniques. The gas phase is assumed to have constant pressure. The governing equations describing the flow of oil and water are coupled and nonlinear. The key stochastic input variables are intrinsic permeability,k, and the soil grain size distribution index, . Three different stochastic combinations of these two input parameters were considered. The perturbation/spectral analysis was used to develop closed-form expressions that describe stochastic variability of key output processes, such as capillary and individual phase pressures and specific discharges. The analysis also included the derivation of the mean flow equations and estimation of the effective flow properties. The impact of the spatial variability ofk and on the effective conductivities and the variances of pressures and specific discharges was examined.  相似文献   

A semi-analytical solution for the transient response of one-dimensional unsaturated single-layer poroviscoelastic media     

《Computers and Geotechnics》2016

This paper develops a semi-analytical solution for the transient response of an unsaturated single-layer poroviscoelastic medium with two immiscible fluids by using the Laplace transformation and the state-space method. Using the elastic–viscoelastic correspondence principle, we first introduce the Kelvin–Voigt model into Zienkiewicz’s unsaturated poroelastic model. The vibrational response for unsaturated porous material can be obtained by combining these two models and assuming that the wetting and non-wetting fluids are compressible, the solid skeleton and solid particles are viscoelastic, and the inertial and mechanical couplings are taken into account. The Laplace transformation and state-space method are used to solve the basic equations with the associated initial and boundary conditions, and the analytical solution in the Laplace domain is developed. To evaluate the responses in the time domain, Durbin’s numerical inverse Laplace transform method is used to obtain the semi-analytical solution. There are three compressional waves in porous media with two immiscible fluids. Moreover, to observe the three compressional waves clearly, we assume the two immiscible fluids are water and oil. Finally, several examples are provided to show the validity of the semi-analytical solution and to assess the influences of the viscosity coefficients and dynamic permeability coefficients on the behavior of the three compressional waves.  相似文献   

Wave reflection from vertical breakwater with porous structure     

Hong-Bin Chen  Ching-Piao Tsai 《Ocean Engineering》2006,33(13):1705-1717

This paper presents numerical solutions for the wave reflection from submerged porous structures in front of the impermeable vertical breakwater. A new time-dependent mild-slope equation involves the parameters of the porous medium including the porosity, the friction factor and the inertia coefficient, etc. is derived for solving the boundary value problem. A comprehensive comparison between the present model and the existing analytical solution for the case of simple rectangular geometries of the submerged structure is performed first. Then, more complicated cases such as the inclined and trapezoidal submerged porous structures in front of the vertical breakwater with sloping bottom are considered. This study also examines the effects of the permeable properties and the geometric configurations of the porous structure to the wave reflection. It is found that the submerged porous structure with trapezoidal shape has more efficiency to reduce the wave reflection than that of triangular shape. The numerical results show that the minimum wave reflection is occurred when the breakwater is located at the intermediate water depth.  相似文献   

Wave characteristics past a flexible fishnet   总被引：1，自引：0，他引：1  

A. T. Chan  S. W. C. Lee 《Ocean Engineering》2001,28(11)

The scattering of surface waves by a flexible fishnet is studied analytically. The fishnet is modelled as a porous flexible barrier displaced solely by hydrodynamic force like a catenary. The objective is to investigate how a flexible permeable barrier affects the passing waves in the way they are transmitted and reflected, as observed by the fact that the water inside a fishfarm surrounded by fishnets is significantly calmer than that outside. The boundary value problems are solved by defining the reflection coefficient in terms of velocity potential and then the full solutions are obtained by suitable application of the eigenfunction expansion method and the least squares approximation method. The variations of the reflection coefficient, hydrodynamic pressure, barrier deformation and surface wave elevation are determined with respect to the barrier length, porosity and stiffness. It is observed that as the fishnet gets more flexible, its deformation increases and the reflection coefficient decreases, whereas as the fishnet gets more porous, more water can pass through it and thus the reflection coefficient, barrier deformation and the hydrodynamic force are reduced. The flexibility of the barrier behaves like its porosity by allowing more wave energy to act on it through its deformation and hence reduce the reflection and hydrodynamic force of the incident waves acting on the barrier.  相似文献   

Wave Action on Structures of Combined Cylinders     

孙路 李玉成腾斌  张宁川 《中国海洋工程》2005,19(3):375-384

1.Introduction Owingtoitssignificanteffectofreductionofwaveloadsandwaverun up,theperforatedwallon anarrayofcylindricalstructureshasreceivedconsiderableattentioninrecentyears.Manyresearches havebeencarriedoutinthisfield.Anexactsolutionforthediffractionoflinearwaterwaveswithan arrayofimpermeablecylinderswasfirstgivenbySpringandMonkmeyer(1974)usinganeigenfunction expansionapproach.Subsequently,LintonandEvans(1990)madeamajorsimplificationtothetheo ry,whichallowedthenear fieldquantitiessuchasload…  相似文献   

Displacement ring load Green's functions for saturated porous transversely isotropic tri‐material full‐space          下载免费PDF全文

K. Sahebkar  M. Eskandari‐Ghadi 《国际地质力学数值与分析法杂志》2017,41(3):359-381

Based on the Biot's poroelastic theory and using scalar potential functions both the ring load and point load displacement Green's functions for a transversely isotropic saturated porous full‐space composed of an upper half‐space, a finite thickness middle layer and a lower half‐space is analytically presented for the first time. It is assumed that each region consists of a different transversely isotropic material. The equations of poroelastodymanics in terms of the solid displacements and the pore fluid pressure are uncoupled with the help of two scalar potential functions, so that the governing equations for the potential functions are either a second order wave equation or a repeated wave‐heat transfer equation of sixth order. With the aid of Fourier expansion with respect to circumferential direction and Hankel integral transforms with respect to the radial direction in cylindrical coordinate system, the response is determined in the form of line integrals in the real space, followed by theorem of inverse Hankel integral transforms. The solutions degenerate to a single phase elastic material, and the results are compared with previous studies, where an excellent agreement may be observed with the results provided in the literature. Some examples of displacement Green's functions are finally given to illustrate the solution. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

Patterned ground formation and penetrative convection in porous media     

J. H. George  R. D. Gunn  B. Straughan 《地球物理与天体物理流体动力学》2013,107(3):135-158

Abstract

A theoretical explanation is advanced consisting of a five stage process for the formation of polygonal ground which consists of stone borders forming regular hexagons and soil centres. One of these stages, namely the onset of convection in a porous soil between temperatures of 0°C and approximately 4-6°C, is studied analytically. Darcy's law is employed but variable permeability is allowed for and a parabolic density dependence on temperature is assumed. It is found that the theoretical predictions of the aspect ratio agree very well with field studies when a constant upper surface heat flux condition is imposed and an upwardly stratified permeability is chosen. Field study data, which agree very well with the theory, are reported in detail.  相似文献