A micropolar mixture theory of multi-component porous media

A mixture theory is developed for multi-component micropolar porous media with a combination of the hybrid mixture theory and the micropolar continuum theory. The system is modeled as multi-component micropolar elastic solids saturated with multi- component micropolar viscous fluids. Balance equations are given through the mixture theory. Constitutive equations are developed based on the second law of thermodynamics and constitutive assumptions. Taking account of compressibility of solid phases, the volume fraction of fluid as an independent state variable is introduced in the free energy function, and the dynamic compatibility condition is obtained to restrict the change of pressure difference on the solid-fluid interface. The constructed constitutive equations are used to close the field equations. The linear field equations are obtained using a linearization procedure, and the micropolar thermo-hydro-mechanical component transport model is established. This model can be applied to practical problems, such as contaminant, drug, and pesticide transport. When the proposed model is supposed to be porous media, and both fluid and solid are single-component, it will almost agree with Eringen＇s model.

Determination of a Pore Structure Parameter of Porous Media by Analysis of Percolation Network Model

According to the simulation of nitrogen sorption process in porous media with three-dimensional network model, and the analysis for such a process with percolation theory, a new method is proposed to determine a pore structure parameter--mean coordination number of pore network, which represents the connectivity among a great number of pores. Here the 'chamber-throat' model and the Weibull distribution are used to describe the pore geometry and the pore size distribution respectively. This method is based on the scaling law of percolation theory after both effects of sorption thermodynamics and pore size on the sorption hysteresis loops are considered. The results show that it is an effective procedure to calculate the mean coordination number for micro- and meso-porous media.

The Kinematics and Field Equations for Porous Media

With a porous medium regarded as an immiscible mixture of multiphase and each phase as a miscible mixture of multi constituent, a systematical research on the kinematics and field equations for porous media is carried out from the point of view of mixture theory. It is shown that the motion of each phase is the mathematical average of the motions of all constituents in the phase, and that the motion of porous media may be described as the motion of the skeleton and the relative motion of each phase with respect to the skeleton. The influence of mass exchange between different constituents in each phase and the influence of mass exchange of same constituent between different phases in porous media are considered in field equations which are self consistent in theory. All the field equations in the references are special cases of the equations proposed in this paper.

GURTIN-TYPE VARIATIONAL PRINCIPLES FOR DYNAMICS OF A NON-LOCAL THERMAL EQUILIBRIUM SATURATED POROUS MEDIUM

Based on the porous media theory and by taking into account the efects of the pore fuid viscidity, energy exchanges due to the additional thermal conduction and convection between solid and fuid phases, a mathematical model for the dynamic-thermo-hydro-mechanical coupling of a non-local thermal equilibrium fuid-saturated porous medium, in which the two constituents are assumed to be incompressible and immiscible, is established under the assumption of small de- formation of the solid phase, small velocity of the fuid phase and small temperature changes of the two constituents. The mathematical model of a local thermal equilibrium fuid-saturated porous medium can be obtained directly from the above one. Several Gurtin-type variational principles, especially Hu-Washizu type variational principles, for the initial boundary value problems of dy- namic and quasi-static responses are presented. It should be pointed out that these variational principles can be degenerated easily into the case of isothermal incompressible fuid-saturated elastic porous media, which have been discussed previously.

Stress dependence of elastic wave dispersion and attenuation in fluid-saturated porous layered media

The fluid-saturated porous layered(FSPL)media widely exist in the Earth's subsurface and their overall mechanical properties,microscopic pore structure and wave propagation characteristics are highly relevant to the in-situ stress.However,the effect of in-situ stress on wave propagation in FSPL media cannot be well explained with the existing theories.To fill this gap,we propose the dynamic equations for FSPL media under the effect of in-situ stress based on the theories of poroacoustoelasticity and anisotropic elasticity.Biot loss mechanism is considered to account for the stress-dependent wave dispersion and attenuation induced by global wave-induced fluid flow.Thomsen's elastic anisotropy parameters are used to represent the anisotropy of the skeleton.A plane-wave analysis is implemented on dynamic equations yields the analytic solutions for fast and slow P waves and two S waves.Modelling results show that the elastic anisotropy parameters significantly determine the stress dependence of wave velocities.Vertical tortuosity and permeability have remarkable effects on fast and slow P-wave velocity curves and the corresponding attenuation peaks but have little effect on S-wave velocity.The difference in velocities of two S waves occurs when the FSPL medium is subjected to horizontal uniaxial stress,and the S wave along the stress direction has a larger velocity,which implies that the additional anisotropy other than that induced by the beddings appears due to horizontal stress.Besides,the predicted velocity results have the reasonable agreement with laboratory measurements.Our equations and results are relevant to a better understanding of wave propagation in deep strata,which provide some new theoretical insights in the rock physics,hydrocarbon exploration and stress detection in deep-strata shale reservoirs.

Permeability analysis of fractured vuggy porous media based on homogenization theory

Based on the characteristics of fractured vuggy porous media,a novel mathematical model was proposed to model fluid flow in such media on fine scale,i.e.,the discrete fracture-vug network model.The new model consists of three systems:porous rock system,fracture system,and vug system.The fractures and vugs are embedded in porous rock,and the isolated vugs could be connected via the discrete fracture network.The flow in porous rock and fractures follows Darcy's law,and the vugs system is free fluid region.Using a two-scale homogenization limit theory,we obtained a macroscopic Darcy's law governing the media on coarse scale.The theoretical formula of the equivalent permeability of the fractured vuggy porous media was derived.The model and method of this paper were verified by some numerical examples.At the end the permeability of some fractured vuggy porous media with typical fracture-vug structures was analyzed.

Dynamic response of bilayered saturated porous media based on fractional thermoelastic theory

Considering the thermal contact resistance and elastic wave impedance at the interface,in this paper we theoretically investigate the thermo-hydro-mechanical(THM)coupling dynamic response of bilayered saturated porous media.Fractional thermoelastic theory is applied to porous media with imperfect thermal and mechanical contact.The analytical solutions of the dynamic response of the bilayered saturated porous media are obtained in frequency domain.Furthermore,the effects of fractional derivative parameters and thermal contact resistance on the dynamic response of such media are systematically discussed.Results show that the effects of fractional derivative parameters on the dynamic response of bilayered saturated porous media are related to the thermal contact resistance at the interface.With increasing thermal contact resistance,the displacement,pore water pressure,and stress decrease gradually.

流体饱和多孔弹性柱体动力响应的微分求积法

基于多孔介质混合物理论,在小变形的假设下,建立了两相不可压流体饱和多孔介质弹性空间轴对称问题的控制方程。在空间域和时间域内分别采用微分求积方法和二阶向后差分格式来离散控制方程,给出了处理对称轴处奇异性条件的方法,并在空间离散后采用消元法来缩减未知量,提高计算速度。作为应用,分析了流体饱和弹性多孔介质圆柱体的动力响应,考察了所布节点数对数值结果的影响。

热局部非平衡流体多孔弹性半平面问题的非线性研究

在几何非线性和热局部非平衡条件下,对不可压流体饱和多孔热弹性半平面在表面温度载荷作用下的热动力学特性进行研究。首先基于多孔介质混合物理论给出了问题的动力学数学模型。然后提出了一种综合的数值计算方法来求解问题,该方法通过微分求积法和二阶后向差分格式在空间域和时间域离散数学模型,利用Newton-Raphson法求解非线性代数方程组,从而可得到问题的数值结果。研究表明,求解方法是有效可靠的,且具有计算量小、精度高等优点。最后,研究了流体饱和多孔热弹性半平面在表面温度载荷下的热力学特性,考察了材料参数和几何非线性对半平面动力学行为的影响。

Bifurcation of Nonlinear Problems Modeling Flows Through Porous Media

This paper deals with a class of nonlinear boundary value problems which appears in the study of models of flows through porous media. Existence results of asymptotic bifurcation and continua are reported both for operator equations and for boundary value problems.

Seismic responses of a hemispherical alluvial valley to SV waves: a three-dimensional analytical approximation

Abstract An analytical solution to the three-dimen-sional scattering and diffraction of plane SV-waves by a saturated hemispherical alluvial valley in elastic half-space is obtained by using Fourier-Bessel series expan-sion technique. The hemispherical alluvial valley with saturated soil deposits is simulated with Biot＇s dynamic theory for saturated porous media. The following conclusions based on numerical results can be drawn： （1） there are a significant differences in the seismic response simulation between the previous single-phase models and the present two-phase model; （2） the nor-malized displacements on the free surface of the alluvial valley depend mainly on the incident wave angles, the dimensionless frequency of the incident SV waves and the porosity of sediments; （3） with the increase of the incident angle, the displacement distributions become more complicated; and the displacements on the free surface of the alluvial valley increase as the porosity of sediments increases.

求解粘弹性饱和土隧道频域响应的DQEM

用微分求积单元法研究了粘弹性饱和土隧道的频域响应问题.首先,基于饱和多孔介质混合物理论,分别建立了频域内隧道耦合系统中粘弹性饱和土和粘弹性衬砌的控制方程,给出了耦合系统的边界条件以及饱和土和衬砌连接面上满足的连接条件.其次,发展了求解该耦合系统的微分求积单元法,利用微分求积单元法,在频率域内对耦合系统的数学模型进行了离散,并利用Newton-Raphson迭代法得到了系统的数值解.最后,对耦合系统的频域响应进行了分析,考察了参数的影响,同时也验证了数值方法的有效性.

粘弹性饱和多孔介质中圆柱孔洞的瞬态响应

饱和多孔介质流-固耦合系统的动力学特性研究具有重要的学术理论价值和潜在的应用前景.研究了含衬砌的无限长粘弹性饱和多孔介质中圆柱孔洞的瞬态响应问题,相关结论可用于饱和土隧道分析、核废料等污染物处置、地热资源开发等基础工程设计.首先,基于多孔介质混合物理论,建立了耦合系统动力学分析的控制方程、边界条件、连接条件和初始条件.其次,采用微分求积单元法和2阶向后差分格式分别在空间和时间域内离散了系统的控制方程.最后,利用Newton-Raphson迭代方法在初始条件下求得了系统的数值解,分析了耦合系统的瞬态响应,并考察了参数的影响,同时也验证了数值方法的有效性.

基于DQEM的分层流体饱和热弹性多孔介质轴对称问题的动力响应分析

研究了分层不可压流体饱和热弹性多孔介质轴对称问题的动力响应问题,基于多孔介质理论(PMT),给出了该问题的数学模型.在空间域内采用微分求积单元法(DQEM)设置离散的控制微分方程、边界条件和连接条件,在时间域内采用二阶向后差分格式处理时间导数.在离散化的初始条件下,运用Newton-Raphson法进行迭代求解,得到各离散点处未知物理量的数值结果.研究表明:该方法有效、可靠,且具备精度较高、计算量较小、数值稳定等优点.

Nonlinear deformation of a cantilever incompressible poroelastic beam under a uniform load

Nonlinear bending of cantilever incompressible poroelastic beams subjected to a uniform load is investigated with the constraint that fluid flow is only in the axial direction. The governing equations for large deflection of the poroelastic beam are derived from theory of incompressible saturated porous media. Then, nonlinear responses of a cantilever beam with impermeable fixed end and permeable free end are examined with the Galerkin truncation method. The deflections and bending moments of the poroelastic beam and the equivalent couples of the pore fluid pressures are shown in figures. The differences of the results between the large deflection and the small deflection theories are analyzed. It is shown that the results of the large deflection theory are smaller than those of the small deflection theory, and the time needed to approach their stationary states for the large deflection theory is shorter than that for the small deflection theory.

A nonlinear mathematical model for large deflection of incompressible saturated poroelastic beams

Nonlinear governing equations are established for large deflection of incompressible fluid saturated poroelastic beams under constraint that diffusion of the pore fluid is only in the axial direction of the deformed beams. Then, the nonlinear bending of a saturated poroelastic cantilever beam with fixed end impermeable and flee end permeable, subjected to a suddenly applied constant concentrated transverse load at its free end, is examined with the Gaierkin truncation method. The curves of deflections and bending moments of the beam skeleton and the equivalent couples of the pore fluid pressure are shown in figures. The results of the large deflection and the small deflection theories of the cantilever poroelastic beam are compared, and the differences between them are revealed. It is shown that the results of the large deflection theory are less than those of the corresponding small deflection theory, and the times needed to approach its stationary states for the large deflection theory are much less than those of the small deflection theory.

Quasi-static and dynamical bending of a cantilever poroelastic beam

Based on the theory of porous media, the quasi-static and dynamical bending of a cantilever poroelastic beam subjected to a step load at its free end is investigated, and the influences of its permeability on bending deformation is examined. The initial boundary value problems for dynamical and quasi-static responses are solved with the Laplace transform technique, and the deflections, the bending moments of the solid skeleton and the equivalent couples of the pore fluid pressure are shown in figures. It is shown that the dynamical and quasi-static behavior of the saturated poroelastic beam depends closely on the permeability conditions at the beam ends. Under the different permeability conditions, the deflections of the beam may oscillate or not. The Mandel-Cryer effect also exists in liquid-saturated poroelastic beams.

Bending of simply supported incompressible saturated poroelastic beams with axial diffusion

Based on the mathematical model of the bending of the incompressible saturated poroelastic beam with axial diffusion, the qUasi-static bendings of the simply supported poroelastic beam subjected to a suddenly applied constant load were investigated, and the analytical solutions were obtained for different diffusion conditions of the pore fluid at the beam ends. The deflections, the bending moments of the solid skeleton and the equivalent couples of the pore pressures were presented in figures. It is also shown that the behavior of the saturated poroelastic beams depends closely on the diffusion conditions at the beam ends, especially for the equivalent couples of the pore pressures. It is found that the Mandel-Cryer effect also exists in the bending of the saturated poroelastic beams under specific diffusion conditions at the beam ends.

Finite element implementation of coupled hydro-mechanical modeling of transversely isotropic porous media in DEAL.II

In geotechnical engineering,modeling geo-structures is challenging,particularly in cases where the interaction between the structures and soil or rock is complex.Most wellknown commercial modeling software is based on homogenous and isotropic materials.However,soil and rock are often modeled in heterogeneous and anisotropic media because of the inherent anisotropy of sedimentary rock masses and their stratified structure.In recent decades,coupled hydro-mechanical(HM)interactions in isotropic porous media have been studied;however,the behavior of transversely isotropic porous media is rarely considered.In addition,it is difficult for commercial software such as Plaxis and Flac3D to express complex rock formation where the anisotropy of the material and the associated cracks and fractures could be assembled into a single model.In this study,a finite element implementation using Differential Equation Analysis Library(DEAL.II),an open-source library of finite element codes,was developed to model the fully coupled HM behavior of transversely isotropic porous media.The proposed implementation can be applied to both isotropic and transversely isotropic porous media based on Biot’s theory.The developed code can be used to model poroelastic media with(1)equations of linear elasticity for the solid matrix and(2)diffusion equations for fluid flow based on mass and linear-momentum conservation laws.We verified the performance and accuracy of the code through two examples,i.e.,Mandel’s problem with a compared analytical solution and a tunnel excavation process with the Flac3D software.On the basis of these numerical applications,we present the code to model the behavior of various geo-structures such as tunnels and pile–soil interactions with anisotropic materials.

Modelling the formation and dissolution behavior of alumina agglomerate in the cryolite

The presence of alumina agglomerates seriously affects the current efficiency of the aluminum electrolysis process.Clarify the dynamic dissolution process of agglomerates is essential to improve the current efficiency of aluminum electrolysis.A mathematical model is proposed to describe the different phenomena from the formation until complete dissolution of agglomerates.Considering permeation and solidification processes of cryolite,a semi-analytical mathematical model is developed to formulate the formation,melting and dissolution processes of agglomerates,and the time duration for each stage is deduced.Porosity and heat mass transfer of agglomerates are explored based on the packing theory and mechanism of heat mass transfer in wet porous media.Dimensionless approach is applied to investigate the main factors affecting the dissolution stages and porosity of agglomerates.The results show that the superheat has a great influence on the formation and melting stage,the diameter of agglomerates can reach 14.93 mm for 200 particles agglomerated.The density decreases with the increase of agglomerated particle number,which varies in the range of 2.27-2.28 g/cm^(3).The average dissolution rate of agglomerates is about 1.83×10^(-5)-2.95×10^(-5) kg/s within the range of alumina concentration in this study.