Effective Diffusivity Tensors of Point-Like Molecules in Anisotropic Porous Media by Monte Carlo Simulation

Monte Carlo simulations of random walks in anisotropic structured media are performed to determine the dependence of effective diffusivities on geometrical properties. The anisotropic media used in this study are periodic systems, which are generated by extending primitive, face-centered, and body-centered unit cells indefinitely in all axial directions. Results of simulations compare well with published experimental data and the calculations by the volume averaging method. In addition, these results suggest that if the 2D media with percolation thresholds subtantially differ from those of 3D, 2D approximations of 3D media are not satisfactory. When percolation thresholds are the same, the effective diffusivity tensors depend solely on the porosity. This fact has been suggested for isotropic media and it seems to hold for anisotropic media.     

Tortuosity based on Anisotropic Diffusion Process in Structured Plate-like Obstacles by Monte Carlo Simulation

The diffusion of fluids in porous media, composed of regularly aligned plate-like obstacles, was studied by Monte Carlo simulation. The diffusion coefficients and all diagonal components of the diffusion tensor were estimated for these media. The calculated tortuosities were modeled as a function of porosity by using the Koponen’s equation related to percolation threshold. These results indicated that a media with a homogeneous porosity has a heterogeneous tortuosity, is affected by the alignments of the plate-like obstacles. Furthermore, the calculation results were compared with the experimental results for fixed perpendicular plates of Comiti and Renaud as a function of porosity. The results for tortuosity compared well for porosity larger than 0.86.     

Monte Carlo Simulations of Observation Time-Dependent Self-Diffusion in Porous Media Models

Observation time-dependent self-diffusion coefficients can be used to obtain microstructural information of porous media. This paper presents two different kinds of Monte Carlo simulations of the self diffusion process of fluids like water in porous systems, a lattice-free method and a lattice-based method. The results for simple porous media model geometries agree well with each other and with published analytical as well as semi-analytical equations. The use of these equations, which are important for the interpretation of Pulsed Field Gradient-Nuclear Magnetic Resonance (PFG-NMR) time-dependent diffusion data with respect to properties of porous media, is discussed.     

Thermodynamic Automaton Simulations of Fluid Flow and Diffusion in Porous Media

A thermodynamic automaton model of fluid flow in porous media is presented. The model is a nonrelativistic version of a Lorentz invariant lattice gas model constructed by Udey et al. (1998). In the previous model it was shown that the energy momentum tensor and the relativistic Boltzman equation can be rigorously derived from the collision and propagation rules. In the present paper we demonstrate that this nonrelativistic model can be used to accurately simulate well known results involving single phase flow and diffusion in porous media. The simulation results show that (1) one-phase flow simulations in porous media are consistent with Darcy's law; (2) the apparent diffusion coefficient decreases with a decrease in permeability; (3) small scale heterogeneity does not affect diffusion significantly in the cases considered.     

膨胀性非牛顿流体多孔介质流动的模拟分析

在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程．采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间．对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效．分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度．当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑．压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显．     

Acoustics of Rotating Deformable Saturated Porous Media

We investigate wave propagation in elastic porous media which are saturated by incompressible viscous Newtonian fluids when the porous media are in rotation with respect to a Galilean frame. The model is obtained by upscaling the flow at the pore scale. We use the method of multiple scale expansions which gives rigorously the macroscopic behaviour without any prerequisite on the form of the macroscopic equations. For Kibel numbers A A(1), the acoustic filtration law resembles a Darcys law, but with a conductivity which depends on the wave frequency and on the angular velocity. The bulk momentum balance shows new inertial terms which account for the convective and Coriolis accelerations. Three dispersive waves are pointed out. An investigation in the inertial flow regime shows that the two pseudo-dilatational waves have a cut-off frequency.     

Review on Scale Dependent Characterization of the Microstructure of Porous Media

The paper discusses local porosity theory and its relation with other geometric characterization methods for porous media such as correlation functions and contact distributions. Special emphasis is placed on the charcterization of geometric observables through Hadwigers theorem in stochastic geometry. The four basic Minkowski functionals are introduced into local porosity theory, and for the first time a relationship is established between the Euler characteristic and the local percolation probabilities. Local porosity distributions and local percolation probabilities provide a scale dependent characterization of the microstructure of porous media that can be used in an effective medium approach to predict transport.     

On the Nield-Kuznetsov Theory for Convection in Bidispersive Porous Media

We revisit the problem of thermal convection in a bidispersive porous medium, first addressed by Nield and Kuznetsov (Int. J. Heat Mass Transfer, 49: 3068–3074, 2006). We investigate the possibility of oscillatory convection by using a highly accurate Chebyshev tau numerical method. We also develop a nonlinear energy stability theory for the same problem. This yields a global stability threshold below which instabilities cannot arise. These thresholds together with the linear instability boundaries yield a zone where thermal instability may be found. The results and theory of Nield and Kuznetsov (Int. J. Heat Mass Transfer, 49: 3068–3074, 2006) are thus proven to be a highly important development in the modern theory of designer porous materials, cf. Nield and Bejan (Convection in Porous Media, Springer, New York, 2006), pp. 94–97. This work was supported in part by a Research Project Grant of the Leverhulme Trust—Grant Number F/00128/AK.     

Three Pressures in Porous Media

In a thermodynamic setting for a single phase (usually fluid), the thermodynamically defined pressure, involving the change in energy with respect to volume, is often assumed to be equal to the physically measurable pressure, related to the trace of the stress tensor. This assumption holds under certain conditions such as a small rate of deformation tensor for a fluid. For a two-phase porous medium, an additional thermodynamic pressure has been previously defined for each phase, relating the change in energy with respect to volume fraction. Within the framework of Hybrid Mixture Theory and hence the Coleman and Noll technique of exploiting the entropy inequality, we show how these three macroscopic pressures (the two thermodynamically defined pressures and the pressure relating to the trace of the stress tensor) are related and discuss the physical interpretation of each of them. In the process, we show how one can convert directly between different combinations of independent variables without re-exploiting the entropy inequality. The physical interpretation of these three pressures is investigated by examining four media: a single solid phase, a porous solid saturated with a fluid which has negligible physico-chemical interaction with the solid phase, a swelling porous medium with a non-interacting solid phase, such as well-layered clay, and a swelling porous medium with an interacting solid phase such as swelling polymers.     

Surfactant Concentration and End Effects on Foam Flow in Porous Media

Foaming injected gas is a useful and promising technique for achieving mobility control in porous media. Typically, such foams are aqueous. In the presence of foam, gas and liquid flow behavior is determined by bubble size or foam texture. The thin-liquid films that separate foam into bubbles must be relatively stable for a foam to be finely textured and thereby be effective as a displacing or blocking agent. Film stability is a strong function of surfactant concentration and type. This work studies foam flow behavior at a variety of surfactant concentrations using experiments and a numerical model. Thus, the foam behavior examined spans from strong to weak.Specifically, a suite of foam displacements over a range of surfactant concentrations in a roughly 7m², one-dimensional sandpack are monitored using X-ray computed tomography (CT). Sequential pressure taps are employed to measure flow resistance. Nitrogen is the gas and an alpha olefin sulfonate (AOS 1416) in brine is the foamer. Surfactant concentrations studied vary from 0.005 to 1wt%. Because foam mobility depends strongly upon its texture, a bubble population balance model is both useful and necessary to describe the experimental results thoroughly and self consistently. Excellent agreement is found between experiment and theory.     

Macroscopic Conductivity of Vugular Porous Media

A systematic numerical study of the macroscopic electrical conductivity of vugular porous media has been conducted by solving the Laplace equation. The structure of these bimodal media is characterized by the micro and macroporosities and by the micro and macro correlation lengths. The overall correlation function is analyzed. It is shown that mainly depends on the total porosity, and very little on the other parameters. Moreover, similar results are obtained when the microporous medium is considered as a continuum. Finally, these predictions are compared to experimental data for the formation factor and the thermal conductivity of limestones; the agreement is good in the first case, and excellent in the second.     

渗流方程自适应非均匀网格Dagan粗化算法

在粗网格内先统计渗透率在粗网格中的概率分布，利用Dagan渗透率粗化积分方程通过渗透率概率分布计算粗化网格的等效渗透率，并由等效渗透率计算了粗化网格的压强分布，计算压强时还将渗透率自适应网格技术应用于三维渗流方程的网格粗化算法中，在渗透率或孔隙度变化异常区域自动采用精细网格，用直接解法求解渗透率或孔隙度变化异常区域的压强分布。整个求解区采用不均匀网格粗化，在流体流速高的区域采用精细网格。利用本文方法计算了三维渗流方程的压强分布，结果表明这种算法的解在渗透率或孔隙度异常区的压强分布规律非常逼近精细网格的解，在其他区域压强分布规律非常逼近粗化算法的解，计算速度比采用精细网格提高了约100倍。     

Steady-State Source Flow in Heterogeneous Porous Media

The flow of fluids in heterogeneous porous media is modelled by regarding the hydraulic conductivity as a stationary random space function. The flow variables, the pressure head and velocity field are random functions as well and we are interested primarily in calculating their mean values. The latter had been intensively studied in the past for flows uniform in the average. It has been shown that the average Darcy's law, which relates the mean pressure head gradient to the mean velocity, is given by a local linear relationship. As a result, the mean head and velocity satisfy the local flow equations in a fictitious homogeneous medium of effective conductivity. However, recent analysis has shown that for nonuniform flows the effective Darcy's law is determined by a nonlocal relationship of a convolution type. Hence, the average flow equations for the mean head are expressed as a linear integro-differential operator. Due to the linearity of the problem, it is useful to derive the mean head distribution for a flow by a source of unit discharge. This distribution represents a fundamental solution of the average flow equations and is called the mean Green function G _d (x). The mean head G _d(x) is derived here at first order in the logconductivity variance for an arbitrary correlation function (x) and for any dimensionality d of the flow. It is obtained as a product of the solution G _d ⁽⁰⁾(x) for source flow in unbounded domain of the mean conductivity K _A and the correction _d (x) which depends on the medium heterogeneous structure. The correction _d is evaluated for a few cases of interest.Simple one-quadrature expressions of _d are derived for isotropic two- and three-dimensional media. The quadratures can be calculated analytically after specifying (x) and closed form expressions are derived for exponential and Gaussian correlations. The flow toward a source in a three-dimensional heterogeneous medium of axisymmetric anisotropy is studied in detail by deriving ₃ as function of the distance from the source x and of the azimuthal angle . Its dependence on x, on the particular (x) and on the anisotropy ratio is illustrated in the plane of isotropy (=0) and along the anisotropy axis ( = /2).The head factor k ^* is defined as a ratio of the head in the homogeneous medium to the mean head, k ^*=G _d ⁽⁰⁾/G _d= _d ^–1. It is shown that for isotropic conductivity and for any dimensionality of the flow the medium behaves as a one-dimensional and as an effective one close and far from the source, respectively, that is, lim _x0 k ^*(x) = K _H/K _A and lim _x k ^*(x) = K ^efu/K _A, where K _A and K _H are the arithmetic and harmonic conductivity means and K ^efu is the effective conductivity for uniform flow. For axisymmetric heterogeneity the far-distance limit depends on the direction. Thus, in the coordinate system of (x) principal directions the limit values of k ^* are obtained as . These values differ from the corresponding components of the effective conductivities tensor for uniform flow for = 0 and /2, respectively. The results of the study are applied to solving the problem of the dipole well flow. The dependence of the mean head drop between the injection and production chambers on the anisotropy of the conductivity and the distance between the chambers is analyzed.     

Coriolis Effects on Filtration Law in Rotating Porous Media

We investigate the filtration law of incompressible viscous Newtonian fluids in rigid non-inertial porous media, for example, rotating porous media. The filtration law is obtained by upscaling the flow at the pore scale. We use the method of multiple scale expansions which gives rigorously the macroscopic behaviour without any prerequisite on the form of the macroscopic equations. For finite Ekman numbers the filtration law is shown to resemble a Darcy's law, but with a non-symmetric permeability tensor which depends on the angular velocity of the porous matrix. We obtain the filtration analog of the Hall effect. For large Ekman numbers the filtration law is a small correction to the classical Darcy's law. The corrector is antisymmetric. In this case we recover a structure of law which is similar to phenomenological laws introduced in the literature, but with a dissimilar effective coefficient.     

Contaminant Transport in Fractured Media with Sources in the Porous Domain

We study contaminant flow with sources in a fractured porous mediumconsisting of a single fracture bounded by a porous matrix. In the fracturewe assume convection, decay, surface adsorption to the interface, and lossto the porous matrix; in the porous matrix we include diffusion, decay,adsorption, and contaminant sources. The model leads to a nonhomogeneous,linear parabolic equation in a quarter-space with a differential equationfor an oblique boundary condition. Ultimately, we study the problemu _t = u _yy – u + f(x,y,t),x,y>0, t>0, u _t = –u _x + u _y – u on y = 0; u(0,0,t) =u₀(t), t>0,with zero initial data. Using Laplace transforms we obtain the Green'sfunction for the problem, and we determine how contaminant sources in theporous media are propagated in time.     

Asymptotic Research of Nonlinear Wave Processes in Saturated Porous Media

The problem of nonlinear wave dynamics of a fluid-saturated porous medium is investigated. The mathematical model proposed is based on the classical Frenkel--Biot--Nikolaevskiy theory concerning elastic wave propagation and includes mass, momentum, energy conservation laws, as well as rheological and thermodynamic relations. The model describes nonlinear, dispersive, and dissipative medium. To solve the system of differential equations, an asymptotic modified two-scales method is developed and a Cauchy problem for initial equations system is transformed to a Cauchy problem for nonlinear generalized Korteweg--de Vries--Burgers equation for modulated quick wave amplitudes and an inhomogeneous set of equations for slow background motion. Stationary solutions of the derived evolutionary equation that have been constructed numerically reflect different regimes of elastic wave attenuation: diffusive, oscillating, and soliton-like.     

Modeling of Two-Phase Flow in Fractured Porous Media on Unstructured Non-Uniformly Coarsened Grids

In this article, we investigate two strategies for coarsening fractured geological models. The first approach, which generates grids that resolve the fractures, is referred to as explicit fracture-matrix separation (EFMS). The second approach is based on a non-uniform coarsening strategy introduced in Aarnes et al. (Adv Water Resour 30(11):2177–2193, 2007a). A series of two-phase flow simulations where the saturation is modeled on the respective coarse grids are performed. The accuracy of the resulting solutions is examined, and the robustness of the two strategies is assessed with respect to number of fractures, degree of coarsening, well locations, phase viscosities, and fracture permeability. The numerical results show that saturation solutions obtained on the non-uniform coarse grids are consistently more accurate than the corresponding saturation solutions obtained on the EFMS grids. The numerical results also reveal that it is much easier to tune the upscaling factor with the non-uniform coarsening approach.     

Convective Flows on Reactive Surfaces in Porous Media

A model for the convective flow in a fluidsaturated porous medium containing a reactive component is considered. This component undergoes an exothermic reaction (modelled by a first order mechanism) on an impermeable bounding surface, the resulting heat released driving the convective flow. Large Rayleigh number flow near a stagnation point is treated in detail by first considering the steady states. Multiple solution branches and critical points arising from a hysteresis bifurcation are identified. The form that these solution branches take depends on whether or not the effects of reactant consumption are included. An initialvalue problem is then discussed. This shows that both the lower (slow reaction) and upper (fast reaction) solution branches are stable (and the ultimate state of the system). When the parameter values are such that there is no steady state, the solution develops a finitetime singularity, the nature of which is analysed.