多孔材料中声波的传播与演化

采用两相多孔介质的拉格朗日模型来描述一种理论流体充填的多孔弹性固体材料,其中孔隙度的变化满足一个附加的平衡方程。     

含气多孔介质卸压层裂的间隔特征——突出的前兆

含气多孔介质在卸压下的拉伸破坏有两种模式－层裂和突出，达到或超过临界破坏条件，引发骨架层裂破坏，达到或超过更强一些的临界突出条件，则引发突出，从采煤实践中逐次掘进的长过程看，层裂是已达破坏条件而未到突出条件的情况，进一步发展可能引发突出。实验表明层裂显微镜规律的间隔特征，因而与层裂相关的瓦斯浓度变化和地音信息可以认作是突出的前兆，如果在内部形成破坏区，渗流场与应力场的耦合有可以导致延时突出，石门自     

液饱和多孔介质中三维应力波的传播

为了深入研究液饱和多孔介质中应力波的传播，提出了三维两相细观计算模型．基于此模型。应用Galerkin余量法并计及流-固耦合界面的耦合效应，利用直接耦合的技术，开发了三维流-固混合显式动力有限元计算程序．在此基础上对冲击载荷作用下液饱和多孔介质中三维应力波的传播现象进行了数值模拟，并详细讨论了孔隙率，孔隙形状等因素对应力波传播主导波形的影响．     

内源瞬态荷载作用下弹性介质与饱和多孔介质中衬砌隧道的动力响应对比分析

衬砌隧道周围介质不同,其内源瞬态荷载作用下的动力响应不同,该文对弹性介质与饱和多孔介质中作用有内源荷载的衬砌隧道的动力响应进行对比分析。首先,令饱和多孔介质孔隙率n=0,将衬砌周围饱和介质退化为弹性介质,动力响应解答相应退化为弹性介质中的解答,退化解与已有的隧道周围为弹性介质时的动力响应解答完全一致。其次,对不同孔隙率的饱和多孔介质和理想弹性介质中衬砌隧道的动力响应进行计算,结果表明：与弹性介质相比,饱和多孔介质中衬砌内表面的动力响应较大,轴向衰减速度较快;随着孔隙率的减小,饱和土中衬砌内表面的动力响应逐渐减小,位移和应力时程曲线均趋近于弹性介质。     

An Analytical Solution to the One-Dimensional Solute Advection-Dispersion Equation in Multi-Layer Porous Media

An analytical solution to the one-dimensional solute advection-dispersion equation in multi-layer porous media is derived using a generalized integral transform method. The solution was derived under conditions of steady-state flow and arbitrary initial and inlet boundary conditions. The results obtained by this solution agree well with the results obtained by numerically inverting Laplace transform-generated solutions previously published in the literature. The analytical solution presented in this paper provides more flexibility with regard to the inlet conditions. The numerical evaluation of eigenvalues and matrix exponentials required in this solution technique can be accurately and efficiently computed using the sign-count method and eigenvalue evaluation methods commonly available. The illustrative calculations presented herein have shown how an analytical solution can provide insight into contaminant distribution and breakthrough in transport through well defined layered column systems. We also note that the method described here is readily adaptable to two and three-dimensional transport problems.     

Asymptotic Analysis of Solute Transport with Linear Nonequilibrium Sorption in Porous Media

We analysed the asymptotic behaviour of breakthrough curves (BTCs) obtained after a single pulse injection in a 1D flow domain. Five different types of solute transport with nonequilibrium sorption were considered. The properties of the porous medium were assumed to be spatially constant. For long times, the concentration at a fixed position in time was found to decay like exp(–t) where depends on both the transport parameters and the parameters describing the nonequilibrium process. The results from the asymptotic analysis were compared with 1D numerical transport calculations. For all cases examined a good agreement between numerical calculations and the asymptotic analysis was found. The results from the asymptotic analysis provide an alternative way to determine transport and sorption related parameters from BTCs. The derived relationships between and the model parameters are however only valid for large times. This requires that the very low concentrations need to be measured and not only the bulk mass, too. By either increasing or decreasing the velocity during BTC experiments additional asymptotic equations are obtained which can be used to determine the value of the model parameters. The results from the asymptotic analysis can also be used in standard inverse modelling techniques to either obtain good initial guesses or to reduce the parameter space. The fact that linear nonequilibrium processes decay like exp(–t) can be used to qualitatively evaluate observed BTCs. The asymptotic analysis can also be easily extended to a larger class of transport problems (e.g. transport of solutes with microbial decay) provided that the overall transport problem remains linear in the concentration.     

线性高温壁面多孔介质传热研究

以多孔介质为研究对象,采用整场求解法分析了多孔介质温度场和流场随Ra数的变化.研究结果表明:随着Ra增加,由开始流线均匀分布为一个顺时针大窝,逐渐变化为一个顺时针大窝,一个逆时针小窝;流线出现分叉现象;流线分布在两个窝中间和外侧,即流场在小窝中间剧烈流动而在大窝中央基本上保持静止状态;当Ra很小时,热的传输主要依靠传导作用,随着Ra数逐渐增大,对流作用成为热传输的主要动力;当Ra数小于105时,高温壁面Nu数基本在-110～-17间变化,最大值发生在底部,最大值为17,最小值发生在上部,最小值为-110.     

Experiments on Solute Transport in Aggregated Porous Media: Are Diffusions Within Aggregates and Hydrodynamic Dispersion Independent?

Two region models for solute transport in porous media assume that hydrodynamic dispersion in mobile water and solute diffusion within immobile water regions are independent. Experimental and theoretical results for transport through a macropore indicate that hydrodynamic dispersion and solute exchange are interdependent. Experiments were carried out to investigate this problem for a column packed with spherical porous aggregates. The effective diffusion coefficient of a tracer within the agreggates was determined from specific experiments. The dispersivity of the bed was determined from experiments carried out with a column filled with nonporous beads. We took advantage of the dependence of hydrodynamic dispersion on density ratios between the invading and displaced solutions to obtain a set of breakthrough curves corresponding to situations where the diffusion coefficient remains constant, whereas the dispersivity varies. Simulations reproduce correctly the experiments. Small discrepancies are noted that can be corrected either by increasing the dispersion coefficient or by fitting the external mass transfer coefficient. Increased dispersion coefficients probably reveal a modification of Taylor dispersion due to solute exchange. The fitted external mass transfer coefficients are close to the values obtained with classical correlations of the chemical engineering literature.     

Smoothed Particle Hydrodynamics Model for Diffusion through Porous Media

A smoothed particle hydrodynamics (SPH) model is presented for the study of diffusion in spatially periodic porous media. The method of SPH is formulated to solve the convection–diffusion equation for tracer diffusion under steady state and transient conditions. Solutions obtained using SPH are compared with other available solutions and the model is used to calculate diffusion coefficients of spatially periodic porous media for the steady state diffusion problem. Diffusion coefficients are then used to calculate nondimensional diffusivities of the media. The effects of media properties on the values of nondimensional diffusivity are also presented.     

Removal of Particle Bridges From a Porous Material by Ultrasonic Irradiation

Particle bridge formation during the flow of a liquid with particles through a porous material is a fouling mechanism that can block the pores and, hence, decrease the permeability of the material. Ultrasonic irradiation of the material is a cleaning method that can restore the permeability. We make a numerical study of this cleaning method using the lattice-Boltzmann method. We start from a pore blocked by two spherical particles attached to the pore wall by colloidal adhesion forces, thus forming a particle bridge. Next we calculate the hydrodynamic force exerted by a high-frequency acoustic wave on the two particles. By comparing the hydrodynamic force and the adhesion force we investigate, whether the particle bridge will be removed by the ultrasonic irradiation. A sensitivity study is carried out to investigate the influence of some relevant parameters, such as the acoustic wave amplitude, the acoustic frequency, the fluid flow velocity and the ratio of particle diameter and pore diameter. An upscaling procedure is applied to translate the microscopic results for the removal of the particles at the pore level to the permeability improvement of the material at the macroscopic level. A comparison is made between numerical results and experimental data. The agreement is reasonable.     

Predicting the Capillary Imbibition of Porous Rocks from Microstructure

The kinetics of capillary imbibition into porous rocks is studied experimentally and theoretically. The Washburn law is modified by introducing various corrections relating to the microstructure of the rocks, such as tortuosity, pore shape (obtained experimentally), and applying the effective medium approximation (EMA) in order to calculate the effective radius that defines the hydraulic conductance and the topology of the capillary imbibition. The application of the EMA shows that capillary imbibition is mainly produced in 1-D, and the pore structure is constituted by different pore throats in series, linked by chamber pores. The capillary process has been discussed as a function of their petrography and pore structure. Our study of the Washburn equation and the addition of correction factors for the pore structure allows a very accurate prediction of the weight rate.     

Finite Element Analysis of Convection Flow Through a Porous Medium in a Horizontal Channel

We analyse the convection flow of a viscous fluid through a horizontal channel enclosing a fully saturated porous medium. The Galerkin finite element analysis is used to discuss the flow and heat transfer through the porous medium using serendipity elements. The velocity, the temperature distributions and the rate of heat transfer are analysed for variations in the governing parameters. The profiles at different vertical levels are asymmetric curves, exhibiting reversal flow everywhere except on the midplane. In a given porous medium, for fixed G or N, the temperature in the fluid region at any position in fluids with a higher Prandtl number, is much higher than in fluids with a lower Prandtl number. Likewise, other parameters being fixed, lesser the permeability of the medium, lower the temperature in the flow field. Nu reduces across the flow at all axial positions, while it enhances along the axial direction of the channel. Nu reduces with decrease in the Darcy parameter D, and thus lesser the permeability of the medium, lesser the rate of heat transfer across the boundary at any axial position of the channel.     

The Peculiarities of Linear Wave Propagation in Double Porous Media

The features of propagation of longitudinal and transverse waves (LW and TW) in fractured porous medium (FPM) saturated with liquid are investigated by methods of multiphase mechanics. The mathematical model of FPM accounting for inequality of velocities and pressures of liquid in pores and fractures, liquid mass exchange and nonstationary interaction forces is developed. Processes of monochromatic wave propagation are studied. The dispersion relation is obtained and the effect of model parameters on wave propagation is analysed. It is established that one transverse and three longitudinal waves propagate in FPM saturated with liquid. The fastest LW is a deformational wave and the two others are filtrational. Filtrational waves attenuate much stronger than deformational and transverse waves. Distinction of velocities and pressures in liquid in various pore systems provides an explanation for the existence of the two filtrational waves in porous medium with two different characteristic sizes of pores.