含随机裂纹网络孔隙材料渗透率的逾渗模型研究

采用逾渗理论对含随机裂纹网络的孔隙材料渗透性进行研究. 开裂孔隙材料渗透率的影响因素包括裂纹网络的几何特征、孔隙材料本体渗透率以及裂纹开度, 本文使用连续区逾渗理论模型建立了渗透率的标度律. 对于裂纹网络的几何特征, 本文基于连续区逾渗理论并考虑裂纹网络的分形特征提出了有限区域内二维随机裂纹网络的连通度定义; 对随机裂纹网络的几何分析表明, 随机裂纹局部团簇效应会降低裂纹网络的整体连通性, 随机裂纹网络的标度指数并非经典逾渗理论给出的固定值, 而是随着网络的分形维数的减小而增大. 本文在网络连通度和主裂纹团的曲折度的基础上, 提出了开裂孔隙材料渗透率标度律的解析表达, K=K₀(K_m,b)(ρ-ρ_c)^μ, 分别考虑了裂纹网络的几何逾渗特征 (ρ-ρ_c)^μ、孔隙材料渗透率K_m 以及裂纹开度比b; 对有限区域含有随机裂纹网络的孔隙材料渗透过程的有限元模拟表明, K₀ 在裂纹逾渗阈值附近与b呈指数关系, 但当裂纹的局部渗透率与K_m比值高于10⁶ 后, 开度比b对渗透率不再有影响.     

无规分形衬底上银膜逾渗系统的交流特性

通过对沉积在具有无规自相似结构的α-Al₂O₃断裂面上的银薄膜逾渗系统交流特性的测量,证明了在逾渗阈值p_c附近这一系统的交流电导率和介电常数满足指数临界规律和普适标度关系。还给出了这种逾渗系统的电容与直流电阻的关系,并对其意义作了讨论。 关键词：     

随机多孔介质的蒙特卡罗重构与渗流特性模拟

根据多孔介质微观结构的分形尺度标度特征,采用蒙特卡罗方法分别重构随机多孔介质的微观颗粒和孔隙结构,并基于分形毛管束模型研究多尺度多孔介质的气体渗流特性,建立多孔介质微观结构和宏观渗流特性的定量关系。结果表明：分形蒙特卡罗重构的多孔介质微细结构接近真实介质结构,气体渗流特性的计算结果与格子玻尔兹曼模拟数据较为吻合; 多孔介质气体渗透率随着克努森数的增加而增大,孔隙分形维数对于气体渗流的微尺度效应具有显著影响,而迂曲度分形维数对于表观渗透率和固有渗透率的比值影响可以忽略。提出的分形蒙特卡罗方法具有收敛速度快且计算误差与维数无关的优点,有利于深入理解多尺度多孔介质的渗流机理。     

二维孔隙裂隙双重介质逾渗规律研究

在孔隙介质逾渗理论的基础上，将另外一个非常重要的渗透通道——裂隙引入到介质的逾渗研究中，提出了更为普遍的孔隙裂隙双重介质的逾渗研究方法.通过对二维平面孔隙裂隙双重介质的数值计算，得到了孔隙裂隙双重介质三个重要参数：孔隙率，裂隙分形维数、裂隙数量分布初值与逾渗概率的关系，给出了孔隙裂隙双重介质逾渗阈值的数学表达式，揭示了孔隙裂隙双重介质的逾渗规律. 关键词： 孔隙 裂隙 双重介质 逾渗 逾渗阈值     

多孔介质中卡森流体的分形分析

研究了非牛顿流体中的卡森流体在多孔介质中的流动特性.基于服从分形分布的弯曲毛细管束模型,运用分形几何理论推导出了该流体在多孔介质中流动的流量、流速、启动压力梯度和有效渗透率的分形解析解.模型中的每一个参数都有明确的物理意义,它将卡森流体在多孔介质中的流动特性与多孔介质的微结构参数有机联系起来.文中给出了卡森流体的流速、启动压力梯度和有效渗透率随着各影响因素的变化趋势,并进行了讨论.所得分形模型可以更深刻地理解卡森流体在多孔介质中流动的内在物理机理. 关键词： 多孔介质 卡森流体 分形     

大孔穴多孔介质内湍动特性研究

建立了用于过滤燃烧数值模拟的大孔穴多孔介质3-D均质网络模型。利用该模型研究了四组不同参数多孔介质内压力梯度与渗流速度之间的变化规律。数值模拟与实验结果吻合得很好。数值计算表明:相对于孔隙率,孔径亦即渗透率对压力梯度的影响更为显著;入口效应仅局限在距入口界面2～3个单元长度内,几何参数和渗流速度对其影响可以忽略;多孔介质的存在抑制了湍动能的发展,在受入口效应影响的区域,湍动能衰减得比较剧烈,之后变化趋缓,最终与耗散率形成动态平衡。     

随机四参数法生成多孔介质及渗流模拟

本文基于多孔介质内流动的广泛应用背景,模拟了多孔介质渗流。采用随机四参数生长法生成了多孔介质模型,并对模型中的微小固相进行了过滤处理。将随机四参数生长法生成的多孔介质结合格子玻尔兹曼方法D2Q9模型模拟了压力驱动的多孔介质中的渗流。通过将模拟结果与著名的CK理论公式的计算结果进行比较,发现在孔隙率高时误差较小(φ=0.9时,误差约为6%)。研究表明在高孔隙率时采用文中的方法得到的多孔介质渗透率及孔隙尺度的流动结果较为准确。     

气体和气溶胶在多孔介质中迁移的数值模拟

为了研究多孔介质对气体和气溶胶的封闭能力,用数值方法模拟了气体和气溶胶在多孔介质中的迁移.考虑了渗透、吸附过滤、扩散、热传导等因素,在柱对称坐标下用双孔双渗模型进行了数值模拟,并与实验结果进行了比较.计算结果显示,多孔介质对气溶胶有较强的阻挡作用,模拟实验结果也证明了这一点.     

具有树形流道双孔隙率多孔介质渗流特性研究

设计具有合适孔道结构的吸附剂能有效提高吸附过程的效率。本文根据构形理论将树形流道嵌入到多孔介质基质中,构建了高渗透率的双孔隙率多孔介质。计算了树形流道内的无量纲压降分布、流道体积分数和渗透率随分叉级数的变化趋势。并利用双孔隙率多孔介质渗流模型计算了该多孔介质的有效渗透率。与并列直管道相比,树形流道在较大的分叉级数下具有更小的压降;树形流道的嵌入能够有效提高多孔介质的有效渗透率。     

多孔介质模型在管壳式换热器数值模拟中的应用

本文详细讨论了多孔介质模型在管壳式换热器数值模拟中的应用,开发了一套能自动生成多孔介质特性参数的通用程序。该程序主要基于三维交错网格及SIMPLE算法,然后运用该模型,采用改进的κ-ε模型和壁面函数法,对换热器壳侧的湍流流动进行了数值模拟。计算结果与换热器冷态实验结果符合良好,表明该模型和计算方法是切实可行的。     

Bulk transport properties and exponent inequalities for random resistor and flow networks

The properties of random resistor and flow networks are studied as a function of the density,p, of bonds which permit transport. It is shown that percolation is sufficient for bulk transport, in the sense that the conductivity and flow capacity are bounded away from zero wheneverp exceeds an appropriately defined percolation threshold. Relations between the transport coefficients and quantities in ordinary percolation are also derived. Assuming critical scaling, these relations imply upper and lower bounds on the conductivity and flow exponents in terms of percolation exponents. The conductivity exponent upper bound so derived saturates in mean field theory.Research supported by the NSF under Grant No. DMR-8314625Research supported by the DOE under Grant No. DE-AC02-83ER13044     

Generic sandpile models have directed percolation exponents

We study sandpile models with stochastic toppling rules and having sticky grains so that with a nonzero probability no toppling occurs, even if the local height of pile exceeds the threshold value. Dissipation is introduced by adding a small probability of particle loss at each toppling. Generically for the models with a preferred direction, the avalanche exponents are those of critical directed percolation clusters. For undirected models, avalanche exponents are those of directed percolation clusters in one higher dimension.     

Percolation of diffusional creep: a new universality class

We study the percolation aspects of diffusional "Coble" creep on heterogeneous grain boundary networks, assuming free grain boundary sliding. A novel percolation threshold is obtained for the honeycomb lattice when two representative types of grain boundaries are randomly distributed, p(cc)=0.5416+/-0.0036. The creep viscosity diverges near the percolation threshold with power-law exponents t=1.69+/-0.09 and s=1.88+/-0.12, different from the standard conduction and rigidity percolation exponents. The moments of both the force and flux distributions all conform to finite-size scaling at p(cc), but with new exponents. These new scaling behaviors seen in the creeping system are proposed to arise from the unique coupling of both force and flux balances in the network.     

Efficient method of finding scaling exponents from finite-size Monte-Carlo simulations

Monte-Carlo simulations are routinely used for estimating the scaling exponents of complex systems. However, due to finite-size effects, determining the exponent values is often difficult and not reliable. Here, we present a novel technique of dealing with the problem of finite-size scaling. This new method allows not only to decrease the uncertainties of the scaling exponents, but makes it also possible to determine the exponents of the asymptotic corrections to the scaling laws. The efficiency of the technique is demonstrated by finding the scaling exponent of uncorrelated percolation cluster hulls.     

Random walk on topologically biased anisotropic percolation clusters and transport properties

Unbiased random walks are performed on topologically biased anisotropic percolation clusters (APC). Topologically biased APCs are generated using suitable anisotropic percolation models. New walk dimensions are found to characterize the anisotropic behaviour of the unbiased random walk on the biased topology. Critical properties of electro and magneto conductivities are characterized estimating respective dynamical critical exponents. A dynamical scaling theory relating dynamical and static critical exponents has been developed. The dynamical critical exponents satisfy the scaling relations within error bar.     

Can surfactant be present at pinch-off of a liquid filament?

Surfactants lower surface tension and are used to facilitate breakup and spreading. How much surfactant remains where a filament of initial radius R breaks is set by the ratio of convection, which sweeps surfactant away, to diffusion, which replenishes it, or Peclet number Pe proportional, variantR. Thus, as is well known, surfactant concentration Gamma-->0 when a macroscale filament breaks. Here theory and simulation are used to investigate pinch-off of microscopic filaments. At breakup, Gamma is shown to be nonzero but uniform on a filament of negligible Pe. Since R must be finite, the zero-Pe limit is transitory and yields to a final regime. Two such regimes with distinct dynamics characterized by different scaling exponents are reported.     

Critical Behaviors and Universality Classes of Percolation Phase Transitions on Two-Dimensional Square Lattice

We have investigated both site and bond percolation on two-dimensional lattice under the random rule and the product rule respectively. With the random rule, sites or bonds are added randomly into the lattice. From two candidates picked randomly, the site or bond with the smaller size product of two connected clusters is added when the product rule is taken. Not only the size of the largest cluster but also its size jump are studied to characterize the universality class of percolation. The finite-size scaling forms of giant cluster size and size jump are proposed and used to determine the critical exponents of percolation from Monte Carlo data. It is found that the critical exponents of both size and size jump in random site percolation are equal to that in random bond percolation. With the random rule, site and bond percolation belong to the same universality class. We obtain the critical exponents of the site percolation under the product rule, which are different from that of both random percolation and the bond percolation under the product rule. The universality class of site percolation differs different from that of bond percolation when the product rule is used.     

Quantum versus geometric disorder in a two-dimensional Heisenberg antiferromagnet

We present a numerical study of the spin-1/2 bilayer Heisenberg antiferromagnet with random interlayer dimer dilution. From the temperature dependence of the uniform susceptibility and a scaling analysis of the spin correlation length we deduce the ground state phase diagram as a function of nonmagnetic impurity concentration p and bilayer coupling g. At the site percolation threshold, there exists a multicritical point at small but nonzero bilayer coupling g(m)=0.15(3). The magnetic properties of the single-layer material La(2)Cu(1-p)(Zn,Mg)(p)O4 near the percolation threshold appear to be controlled by the proximity to this new quantum critical point.