Ageing of random porous media following fluid deterministic displacement, freezing, thawing

This paper introduces and investigates a simple model of random porous media degradation via several fluid displacing, freezing, and thawing cycles. The fluid transport is based on the deterministic method. The result shows that the topology and the geometry of porous media have a strong effect on displacement processes. The cluster size of the viscous fingering (VF) pattern in the percolation cluster increases with the increase of iteration parameter n. When iteration parameter , the VF pattern does not change with n. When and , the peak value of the distribution increases as n increases; is the normalized distribution of throat sizes after different displacement-damage but before the freezing. The distribution of throat size N(r) after displacement but before freezing damage, shows that the major change, after successive cycles, happens at r>0.9. The peak value of the distribution reaches a maximum when and r=1, where is the normalized distribution of the size of invaded throats for different iterations. This result is different from invasion percolation. The distribution of velocities normal to the interface of VF in the percolation cluster is also studied. When , the scaling function distribution is very sharp. The sweep efficiency E increases along with the increasing of iteration parameter n and decreases with the network size L. And E has a minimum as L increases to the maximum size of the lattice. The VF pattern in the percolation cluster has one frozen zone and one active zone. Received 30 March 1999 and Received in final form 8 August 1999     

Two—phase flow in correlated pore—throat random porous media

We have constructed a porous media model in which there are percolation clusters with varying percolation probability P and correlated site-bonds. Taking into account both the pore and the throat geometry, the viscous fingering (VF) in porous media has been investigated by using the standard over-relaxed Gauss-Seidel scheme. The simulation results show that the VF structure varies with the correlation parameter ε, the viscosity ratio M and the percolation probability P. The smaller the correlation parameter ε, the greater the deviation of the normalized size distribution of the invaded throat N_inv(r) from the truncated Rayleigh distribution. For a larger viscosity ratio M, the VF pattern looks like a diffusion-limited-aggregation structure in percolation clusters. The fractal dimension D increases with the increase of the percolation probability P and the correlation parameter ε. The velocity distribution f(α) of VF in percolation clusters is of a parabola-like curve. The tail of the distribution (large α) is longer for a larger correlation parameter ε. For a smaller ε, the distribution is very sharp. The sweep efficiency E decreases along with the decrease of the correlation parameter ε and the increase of the network size L_nz. E has a minimum as L_nz increases up to the maximum no matter what the values of P, M and ε. The E～ L_nz curve has a frozen zone and an active zone. The geometry and the topology of the porous media have strong effects on the displacement processes and the structure of VF.     

Computer simulations of two-phase flow with surface tension in the percolation cluster

Percolation clusters with varying occupation probability were constructed. Viscous fingering (VF) in the percolation cluster, based on the assumption that throat radii are Rayleigh distributed, is investigated by means of a successive over-relaxation technique. The fractal dimension and the sweep efficiency of VF in the percolation cluster when surface tension is considered are larger than when surface tension is neglected. The fractal dimension of VF will increase as the percolation probability increases or the viscous ratio decreases. VF's fractal dimension of porous media in the limit viscous ratio → ∞ is found to be identical with the DLA. The topology and the geometry of the porous medium have a strong effect on the displacement processes and the structure of the VF.     

FRACTAL VISCOUS FINGERING AND ITS SCALING STRUCTURE IN RANDOM SIERPINSKI CARPET

Viscous fingering (VF) in random Sierpinski carpet is investigated by means of successive over-relaxation technique and under the assumption that bond radii are of Rayleigh distribution. In the random Sierpinski network, the VF pattern of porous media in the limit M→∞ (M is the viscosity ratio and equals to η₂/η₁ where η₁ and η₂ are the viscosities of the injected and displaced fluids, respectively) is found to be similar to the diffusion-limited aggregation (DLA) pattern. The interior of the cluster of the displacing fluid is compact on long length scales when M=1, and the pores in the interior of the cluster have been completely swept by the displacing fluid. For finite values of M such as M≥10, the pores in the interior of the cluster have been only partly swept by the displacing fluid on short length scales. But for values of M in 1f(α) sites have velocites scaling as L^-α; and the scaling function f(α) is measured and its variation with M is found.     

Simulations of migration, fragmentation and coalescence of non-wetting fluids in porous media

A model based on invasion percolation was used to simulate the migration on a non-wetting fluid through a porous medium filled with an immiscible wetting fluid under the influence of a gradient such as that provided by gravity. The migrating fluid clusters undergo both fragmentation and coalescence. The fragment size distribution obtained from two-dimensional simulations in which the gradient g is slowly increased from 0 can be represented by the scaling form N_s(g)s^-2ƒ(s|g|^-z where z=1+(D−1)ν(ν+1). Here D is the fractal dimensionality of invasion percolation, with trapping, and ν is the ordinary percolation correlation length exponent.     

Aging kinetics of porous media due to freezing-thawing cycles

The two-dimensional Ausloos et al. model of fluid invasion, freezing and thawing in a porous medium is elaborated upon and investigated in order to take into account the pore volume redistribution and conservation during freezing. The results are qualitatively different from previous work, since the damaged pore sizes are found to be much less than the possible maximum value and is reached after a large number of invasion-freezing-thawing cycles, e.g. the material is “slowly damaged”. The pore size distribution is thus found in better agreement with expected practical findings. The successive invasion percolation clusters are still found to be self-avoiding with aging. The cluster size decreases with a power law as a function of invasion-frost-thaw iterations. The aging kinetics is also discussed through the normalized totally invaded pore volume. Received 24 September 1999 and Received in final form 5 January 2000     

自仿射Sierpinski地毯中的粘滞指进

构造了自仿射Sierpinski地毯.在自仿射Sierpinski地毯中,认为喉管半径服从截断瑞利分布,采用逐次超松弛技术,模拟了自仿射Sierpinski地毯中的粘滞指进.计算了粘滞指进分维,结果表明,当粘滞比M→∞时,指进图样与在自仿射Sierpinski地毯中DLA模拟的结果类似：当M=1时,驱替流体在长标度范围内具有紧凑的结构,且具有稳定的位移. 关键词：     

沉降分布孔隙率多孔介质模型及数值研究

研究沉降分布孔隙率多孔介质流动和传热,根据"O"形圈理论和现场测定确定孔隙率系数,建立坐标方向孔隙率分布函数;考虑流体密度变化,并引入Brinkman-Forchheimer的扩展Darcy模型,能量方程采用界面连续条件,建立沉降分布孔隙率多孔介质流动和传热求解模型.采用差分法对模型进行离散化,应用高斯-赛德尔方法迭代求解.数值分析表明:沉降分布孔隙率条件下多孔介质内流体流动速度在壁面附近较大,中心部位较小,壁面附近孔隙率的增大使得低流速区域减小,较高流速区域增大;当孔隙率小值时,温度按线性减小;当孔隙率大值时,温度在高低温壁面附近迅速减小,在中部减小较缓,热量按导热和对流共同传递;孔隙率增大能使平均怒谢尔数增大,对流换热作用增强.     

Tricritical point in explosive percolation

The suitable interpolation between classical percolation and a special variant of explosive percolation enables the explicit realization of a tricritical percolation point. With high-precision simulations of the order parameter and the second moment of the cluster size distribution a fully consistent tricritical scaling scenario emerges yielding the tricritical crossover exponent 1/φ(t)=1.8 ± 0.1.     

随机多孔介质逾渗模型渗透率的临界标度性质

研究了一类非零键渗透率满足均匀分布的随机多孔介质逾渗模型-数值计算了该模型系统渗透率在临界点处的标度指数-结果表明该指数并不能看作是普适常数,而与均匀分布的参数有关-这意味着即使非零键渗透率值的概率密度函数满足负一阶矩存在条件,系统渗透率在逾渗临界点处的标度指数仍然依赖于分布函数的具体参数,并不是常数-这一数值结果与Sahimi对此问题的结论不同- 关键词： 逾渗 随机多孔介质 标度指数 渗透率     

In situ X-ray diffraction investigation of porous silicon strains induced by the freezing of a confined organic fluid

High resolution X-ray diffraction is used to perform an in situ measurement of the variations of the lattice parameter of the nanometer size crystallites of porous silicon, induced by the freezing of a confined organic fluid, dodecane. Two p⁺ type PS layers of 60 and 70% porosity are investigated, and the variations of their lattice parameter with the temperature (in the range 150–300 K) are measured. The experimental curves are discussed in relation with the results of a previous calorimetric study of the freezing of confined dodecane. We explain the observed strains by the presence of capillary stresses, that appear in the layer due to the formation of internal liquid-vapour meniscus during the freezing process of the confined fluid. Received 22 October 1999 and Received in final form 25 March 2000     

Numerical simulation of heat transfer and fluid flow of Water-CuO Nanofluid in a sinusoidal channel with a porous medium

This study has compared the convection heat transfer of Water-based fluid flow with that of Water-Copper oxide (CuO) nanofluid in a sinusoidal channel with a porous medium. The heat flux in the lower and upper walls has been assumed constant, and the flow has been assumed to be two-dimensional, steady, laminar, and incompressible. The governing equations include equations of continuity, momentum, and energy. The assumption of thermal equilibrium has been considered between the porous medium and the fluid. The effects of the parameters, Reynolds number and Darcy number on the thermal performance of the channel, have been investigated. The results of this study show that the presence of a porous medium in a channel, as well as adding nanoparticles to the base fluid, increases the Nusselt number and the convection heat transfer coefficient. Also the results show that As the Reynolds number increases, the temperature gradient increases. In addition, changes in this parameter are greater in the throat of the flow than in convex regions due to changes in the channel geometry. In addition, porous regions reduce the temperature difference, which in turn increases the convective heat transfer coefficient.     

Spanning lengths of percolation clusters

The spanning length of a percolation cluster is defined as the difference between the maximum and minimum coordinates of the cluster with respect to some chosen direction. It is statistically related to the number size of the cluster by an exponent that differs from the iriverse dimension that would characterize a compact cluster. This exponent for large percolation clusters in simple cubic lattice sites was studied by the Monte Carlo technique, and results are presented. Previous theoretical treatments of this exponent and its relationship with other critical exponents are discussed.In the present paper we shall refer exclusively to the site percolation problem, and all our definitions will be within that context.     

Stability analysis of a Maxwell fluid in a porous medium heated from below

Based on a modified-Darcy–Brinkman–Maxwell model, stability analysis of a horizontal layer of Maxwell fluid in a porous medium heated from below is performed. By solving the eigenvalue problems, the critical Rayleigh number, wave number and frequency for overstability are determined. It is found that the critical Rayleigh number for overstability decreases as the relaxation time increases and the elasticity of a Maxwell fluid has a destabilizing effect on the fluid layer in porous media. On the other hand, the critical Rayleigh number for overstability increases by increasing the porous parameter which acts to stabilize the system. In limiting cases, some previous results for viscoelastic fluids in nonporous media are recovered from our results.     

Molecular dynamics simulation of decomposition and thermal conductivity of methane hydrate in porous media

The hydrate has characteristics of low thermal conductivity and temperature sensitivity. To further analysis the mechanism of thermal conductivity and provide method for the exploitation, transportation and utilization of hydrate, the effect of decomposition and thermal conductivity of methane hydrate in porous media has been studied by using the molecular dynamics simulation. In this study, the simulation is carried out under the condition of temperature 253.15 K-273.15 K and pressure 1 MPa. The results show that the thermal conductivity of methane hydrate increases with the increase of temperature and has a faster growth near freezing. With the addition of porous media, the thermal conductivity of the methane hydrate improves significantly. The methane hydrate-porous media system also has the characteristics of vitreous body.With the decrease of the pore size of the porous media, thermal conductivity of the system increases gradually at the same temperature. It can be ascertained that the porous media of different pore sizes have strengthened the role of the thermal conductivity of hydrates.     

Influence of the size distribution of granules and of their attractive interaction on the percolation threshold in granulated alloys

Numerical simulation was applied to study the influence of the size distribution of granules and the interaction between them on the percolation threshold in granulated metal-insulator alloys. An alloy model was considered in which metal granules have two characteristic sizes, l and L (with L>l), and the size distribution of granules of greater size L having an average value of approximately L ₀ is described by a normal distribution with a standard deviation d, by a step function with a halfwidth d, or by a delta function. A model with attraction between granules and mechanisms of trapping of an additional granule by an already developed cluster with a characteristic trapping range R was also considered. The percolation threshold significantly grows with the ratio L ₀/l and with R for both two-and three-dimensional cases and tends to flattening at large L/l or R. The calculated results make it possible to explain the high percolation threshold observed for the majority of granulated alloys.