Thermomechanics of Swelling Biopolymeric Systems

A two-scale theory for the swelling biopolymeric media is developed. At the microscale, the solid polymeric matrix interacts with the solvent through surface contact. The relaxation processes within the polymeric matrix are incorporated by modeling the solid phase as viscoelastic and the solvent phase as viscous at the mesoscale. We obtain novel equations for the total stress tensor, chemical potential of the solid phase, heat flux and the generalized Darcy's law all at the mesoscale. The constitutive relations are more general than those previously developed for the swelling colloids. The generalized Darcy's law could be used for modeling non-Fickian fluid transport over a wide range of liquid contents. The form of the generalized Fick's law is similar to that obtained in earlier works involving colloids. Using two-variable expansions, thermal gradients are coupled with the strain rate tensor for the solid phase and the deformation rate tensor for the liquid phase. This makes the experimental determination of the material coefficients easier and less ambiguous.     

Thermo-Chemo-Electro-Mechanical Formulation of Saturated Charged Porous Solids

A thermo-chemo-electro-mechanical formulation of quasi-static finite deformation of swelling incompressible porous media is derived from a mixture theory including the volume fraction concept. The model consists of an electrically charged porous solid saturated with an ionic solution. Incompressible deformation is assumed. The mixture as a whole is assumed locally electroneutral. Different constituents following different kinematic paths are defined: solid, fluid, anions, cations and neutral solutes. Balance laws are derived for each constituent and for the mixture as a whole. A Lagrangian form of the second law of thermodynamics for incompressible porous media is used to derive the constitutive restrictions of the medium. The material properties are shown to be contained in one strain energy function and a matrix of frictional tensors. A principle of reversibility results from the constitutive restrictions. Existing theories of swelling media should be evaluated with respect to this principle.     

Mixture theory and unsaturated flow in swelling soils

The current state of knowledge on various aspects of mixture theory applied to unsaturated/saturated swelling soils is discussed. Two and three phase problems are studied. On the smallest scale (micro) the individual platelets and adsorbed (vicinal) water exist as separate phases. On the intermediate scale (meso) the platelets and vicinal water are homogenized to form a saturated particle where vicinal water and solid are overlaying continua. On the macroscale, bulk water, water vapor, and the mesoscale particles are homogenized resulting in three overlaying continua for solid, bulk water, vicinal water, and water vapor. Stress tensor formulations and Darcy’s laws are presented at the mesoscale and macroscales. A theoretical formulation for surface crusting is presented at the mesoscale. General viscoelastic mesoscale and macroscale models are discussed and related to classical approaches.     

煤层气在非饱和水流阶段的非定常渗流摄动解

煤层甲烷由煤层的割理裂隙系统流入生产井一般经历：单相水流、非饱和流和气、水两相饱和流三个阶段，在非饱和流阶段，储层压力降至临界解吸压力之后，储存在煤基质中的吸附气体少量被解吸出来形成互不连续的气泡并阻止水的流动，含气量尚未达到饱和程度。同时煤层甲烷运移包含渗流场、变形场和应力场的动态耦合过程。本文考虑渗流过程中水-气两相不溶混流体与固体耦合作用，建立了非饱和水流阶段非定常渗流问题的流固耦合数学模型，对该强非线性一维数学模型采用摄动法和积分变换法进行解析求解，并讨论了其压力动态特性，分析了压力随饱和度S及时间t变化的规律和气相及耦合作用的影响，这些研究对煤层气、石油和天然气的开采等地下工程领域具有一定的指导意义。     

Derivation of a Darcy’s Law for a Porous Medium Composed of Two Solid Phases Saturated by a Single- Phase Fluid: A Homogenization Approach

The objective of this article is to derive a macroscopic Darcy’s law for a fluid-saturated moving porous medium whose matrix is composed of two solid phases which are not in direct contact with each other (weakly coupled solid phases). An example of this composite medium is the case of a solid matrix, unfrozen water, and an ice matrix within the pore space. The macroscopic equations for this type of saturated porous material are obtained using two-space homogenization techniques from microscopic periodic structures. The pore size is assumed to be small compared to the macroscopic scale under consideration. At the microscopic scale the two weakly coupled solids are described by the linear elastic equations, and the fluid by the linearized Navier–Stokes equations with appropriate boundary conditions at the solid–fluid interfaces. The derived Darcy’s law contains three permeability tensors whose properties are analyzed. Also, a formal relation with a previous macroscopic fluid flow equation obtained using a phenomenological approach is given. Moreover, a constructive proof of the existence of the three permeability tensors allows for their explicit computation employing finite elements or analogous numerical procedures.     

不可压饱和粘弹性多孔介质的变分原理及其无网格模拟

基于饱和多孔介质理论,在固相和液相微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,建立了流体饱和粘弹性多孔介质动力响应的若干Gurtin型变分原理,包括Hu-Washizu变分原理.利用所建立的变分原理,导出了流体饱和粘弹性多孔介质动力响应无网格数值模拟的离散控制方程,此方程是一个关于时间的对称微分方程组,便于分析计算.作为数值例子,研究了流体饱和粘弹性多孔柱体的一维动力响应,数值结果揭示了流体饱和粘弹性多孔柱体中波的传播特性以及固相粘性的影响.     

不可压饱和粘弹性多孔介质固结问题的Gurtin型变分原理

基于多孔介质理论,在固相骨架和孔隙流体微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,利用卷积积分的性质,本文首先建立了以固相骨架位移、孔隙流体相对速度和孔隙流体压力为宗量的流体饱和粘弹性多孔介质固结问题的一个Gurtin型变分原理.其次,利用Lagrange乘子法解除相关的变分约束条件,建立了流体饱和粘弹性多孔介质固结问题的若干广义Gurtin型变分原理,包括第三类的Hu-Washizu型变分原理.最后,简单讨论了等价初边值问题的相应变分原理.这些Gurtin型变分原理的建立不仅丰富了饱和粘弹性多孔介质的相关理论,而且为相关数值模拟方法,如有限元法、无网格法等的建立奠定了理论基础.     

On the effective stress in unsaturated porous continua with double porosity

Using mixture theory we formulate the balance laws for unsaturated porous media composed of a double-porosity solid matrix infiltrated by liquid and gas. In this context, the term ‘double porosity’ pertains to the microstructural characteristic that allows the pore spaces in a continuum to be classified into two pore subspaces. We use the first law of thermodynamics to identify energy-conjugate variables and derive an expression for the ‘effective’, or constitutive, stress that is energy-conjugate to the rate of deformation of the solid matrix. The effective stress has the form , where σ is the total Cauchy stress tensor, B is the Biot coefficient, and is the mean fluid pressure weighted according to the local degrees of saturation and pore fractions. We identify other emerging energy-conjugate pairs relevant for constitutive modeling of double-porosity unsaturated continua, including the local suction versus degree of saturation pair and the pore volume fraction versus weighted pore pressure difference pair. Finally, we use the second law of thermodynamics to determine conditions for maximum plastic dissipation in the regime of inelastic deformation for the unsaturated two-porosity mixture.     

Convective Instability of Oldroyd-B Fluid Saturated Porous Layer Heated from Below using a Thermal Non-equilibrium Model

The linear stability of a viscoelastic fluid saturated densely packed horizontal porous layer heated from below and cooled from above is investigated by considering the Oldroyd-B type fluid. A generalized Darcy model, which takes into account the viscoelastic properties, is employed as momentum equation and a two-field model is used for energy equation each representing solid and fluid phases separately. Linear stability analysis suggests that, there is a competition between the processes of viscous relaxation and thermal diffusion that causes the first convective instability to be oscillatory rather than stationary. Analytical expression for the occurrence of oscillatory onset is obtained, and it is found that the necessary condition for the existence of the same is Λ < 1. Besides, the effect of viscoelastic parameters and the thermal non-equilibrium on the stability of the system is analyzed.     

Transport of Viruses Through Saturated and Unsaturated Columns Packed with Sand

Laboratory-scale virus transport experiments were conducted in columns packed with sand under saturated and unsaturated conditions. The viruses employed were the male-specific RNA coliphage, MS2, and the Salmonella typhimurium phage, PRD1. The mathematical model developed by Sim and Chrysikopoulos (Water Resour Res 36:173–179, 2000) that accounts for processes responsible for removal of viruses during vertical transport in one-dimensional, unsaturated porous media was used to fit the data collected from the laboratory experiments. The liquid to liquid–solid and liquid to air–liquid interface mass transfer rate coefficients were shown to increase for both bacteriophage as saturation levels were reduced. The experimental results indicate that even for unfavorable attachment conditions within a sand column (e.g., phosphate-buffered saline solution; pH = 7.5; ionic strength = 2 mM), saturation levels can affect virus transport through porous media.     

Modeling Large-deformation-induced Microflow in Soft Biological Tissues

The homogenization approach to multiscale modeling of soft biological tissues is presented. The homogenized model describes the relationship between the macroscopic hereditary creep behavior and the microflow in a fluid-saturated dual-porous medium at the microscopic level. The micromodel is based on Biot’s system for quasistatic deformation processes, modified for the updated Lagrangian formulation to account for coupling the fluid diffusion through a porous solid undergoing large deformation. Its microstructure is constituted by fluid-filled inclusions embedded in the porous matrix. The tangential stiffness coefficients and the retardation stress for the macromodel are derived for a time-stepping algorithm. Numerical examples are discussed, showing the strong potential of the model for simulations of deformation-driven physiological processes at the microscopic scale.     

Revisiting poromechanics in the context of microporous materials

Poromechanics offers a consistent theoretical framework for describing the mechanical response of porous solids, fully or partially saturated with a fluid phase. When dealing with fully saturated microporous materials, which exhibit pores of the nanometre size, aside from the fluid pressure acting on the pore walls additional effects due to adsorption and confinement of the fluid molecules in the smallest pores must be accounted for. From the mechanical point of view, these phenomena result into volumetric deformations of the porous solid: the so-called “swelling” phenomenon. The present work investigates how the poromechanical theory should be refined in order to describe adsorption and confinement induced swelling in microporous solids. Firstly, we report molecular simulation results that show that the pressure and density of the fluid in the smallest pores are responsible for the volumetric deformation of the material. Secondly, poromechanics is revisited in the context of a microporous material with a continuous pore size distribution. Accounting for the thermodynamic equilibrium of the fluid phase in the overall pore space, the new formulation introduces an apparent porosity and an interaction free energy. We use a prototype constitutive relation relating these two quantities to the Gibbs adsorption isotherm, and then calculate the induced deformation of the solid matrix. Agreement with experimental data found in the literature is observed. As an illustrating example, we show the predicted strains in the case of adsorption of methane on activated carbon.     

A comparison of boundary element and finite element methods for modeling axisymmetric polymeric drop deformation

A modified boundary element method (BEM) and the DEVSS‐G finite element method (FEM) are applied to model the deformation of a polymeric drop suspended in another fluid subjected to start‐up uniaxial extensional flow. The effects of viscoelasticity, via the Oldroyd‐B differential model, are considered for the drop phase using both FEM and BEM and for both the drop and matrix phases using FEM. Where possible, results are compared with the linear deformation theory. Consistent predictions are obtained among the BEM, FEM, and linear theory for purely Newtonian systems and between FEM and linear theory for fully viscoelastic systems. FEM and BEM predictions for viscoelastic drops in a Newtonian matrix agree very well at short times but differ at longer times, with worst agreement occurring as critical flow strength is approached. This suggests that the dominant computational advantages held by the BEM over the FEM for this and similar problems may diminish or even disappear when the issue of accuracy is appropriately considered. Fully viscoelastic problems, which are only feasible using the FEM formulation, shed new insight on the role of viscoelasticity of the matrix fluid in drop deformation. Copyright © 2001 John Wiley & Sons, Ltd.     

基于介观结构的饱和与非饱和多孔介质有效应力

基于描述含液颗粒材料介观结构的Voronoi 胞元模型和离散颗粒集合体与多孔连续体间的介-宏观均匀化过程, 定义饱和与非饱和多孔介质有效应力. 导出了计及孔隙液压引起之颗粒体积变形的饱和多孔介质广义有效应力. 用以定义广义有效应力的Biot 系数不仅依赖于颗粒材料的多孔连续体固体骨架及单个固体颗粒的体积模量(材料参数),同时与固体骨架当前平均广义有效应力及单个固体颗粒的体积应变(状态量) 有关. 提出了描述非饱和多孔介质中非混和固体颗粒、孔隙液体和气体等三相相互作用的具介观结构的Voronoi 胞元模型.具体考虑在低饱和度下双联(binary bond) 模式的摆动(pendular) 液桥系统介观结构. 导出了基于介观水力-力学模型的非饱和多孔介质的各向异性有效应力张量与有效压力张量. 考虑非饱和多孔介质Voronoi 胞元模型介观结构的各向同性情况,得到了与非饱和多孔连续体理论中唯象地假定的标量有效压力相同的有效压力形式.但本文定义的与确定非饱和多孔介质有效应力和有效压力相关联的Bishop 参数由基于三相介观水力-力学模型, 作为饱和度、孔隙度和介观结构参数的函数导出,而非唯象假定.      

Coupled Solvent and Heat Transport of a Mixture of Swelling Porous Particles and Fluids: Single Time-Scale Problem

A three-spatial scale, single time-scale model for both moisture and heat transport is developed for an unsaturated swelling porous media from first principles within a mixture theoretic framework. On the smallest (micro) scale, the system consists of macromolecules (clay particles, polymers, etc.) and a solvating liquid (vicinal fluid), each of which are viewed as individual phases or nonoverlapping continua occupying distinct regions of space and satisfying the classical field equations. These equations are homogenized forming overlaying continua on the intermediate (meso) scale via hybrid mixture theory (HMT). On the mesoscale the homogenized swelling particles consisting of the homogenized vicinal fluid and colloid are then mixed with two bulk phase fluids: the bulk solvent and its vapor. At this scale, there exists three nonoverlapping continua occupying distinct regions of space. On the largest (macro) scale the saturated homogenized particles, bulk liquid and vapor solvent, are again homogenized forming four overlaying continua: doubly homogenized vicinal fluid, doubly homogenized macromolecules, and singly homogenized bulk liquid and vapor phases. Two constitutive theories are developed, one at the mesoscale and the other at the macroscale. Both are developed via the Coleman and Noll method of exploiting the entropy inequality coupled with linearization about equilibrium. The macroscale constitutive theory does not rely upon the mesoscale theory as is common in other upscaling methods. The energy equation on either the mesoscale or macroscale generalizes de Vries classical theory of heat and moisture transport. The momentum balance allows for flow of fluid via volume fraction gradients, pressure gradients, external force fields, and temperature gradients.     

Flow in a Porous Matrix with Anisotropic Structure

The flow of a viscous fluid through a porous matrix undergoing only infinitesimal deformation is described in terms of intrinsic variables, namely, the density, velocity and stress occurring in coherent elements of each material. This formulation arises naturally when macroscopic interfaces are conceptually partitioned into area fractions of fluid–fluid, fluid–solid, and solid–solid contact. Such theory has been shown to yield consistent jump conditions of mass, momentum and energy across discontinuities, either internal or an external boundary, unlike the standard mixture theory jump conditions. In the previous formulation, the matrix structure has been considered isotropic; that is, the area fractions are independent of the interface orientation. Here, that is not assumed, so in particular, the cross-section area of a continuous fluid tube depends on its orientation, which influences the directional fluxes, and in turn the directional permeability, anisotropy of the structure. The simplifications for slow viscous flow are examined, and particularly for an isotropic linearly elastic matrix in which area partitioning induces anisotropic elastic response of the mixture. A final specialization to an incompressible fluid and stationary matrix leads to potential flow, and a simple plane flow solution is presented to illustrate the effects of anisotropic permeability.     

Deformation and stress from in-pore drying-induced crystallization of salt

The deformation and the fracture of porous solids from internal crystallization of salt is explored in the framework of the thermodynamics of unsaturated brittle poroelasticity. In the first place the usual theory of crystal growth in confined conditions is further developed in order to include both the deformation and the drying of the porous solid. The thermodynamics reveals the existence of a dilation coefficient associated with the crystallization process, and provides a solute-crystal equilibrium condition which involves the relative humidity, the supersaturation, and the salt characteristics. This thermodynamic condition and the mechanical equilibrium of the solution-crystal interface combine to give the current crystallization pore radius. Upscaling this information at the macroscopic scale, and taking into account the salt mass supplied by the invading solution, the approach leads to a quantitative analysis of the role of the pore size distribution on the crystal growth under repeated imbibition-drying cycles. The deformation and the fracture of the porous solid from drying-induced crystallization are then considered in the context of brittle poroelasticity. The current unsaturated macroscopic poroelastic properties are upscaled from the microscopic elastic properties of the solid matrix and from the current liquid, crystal and gas saturations. The adoption of a fracture criterion based on the elastic energy that the solid matrix can ultimately store finally leads to the determination of how long a stone can resist repeated cycles of drying-induced crystallization of salt.     

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.     

含液颗粒材料液固耦合分析的离散颗粒模型及特征线SPH法

对含液颗粒材料流固耦合分析建议了一个基于离散颗粒模型与特征线SPH法的显式拉格朗日-欧拉无网格方案。在已有的用以模拟固体颗粒集合体的离散颗粒模型[1]基础上,将颗粒间间隙内的流体模型化为连续介质,对其提出并推导了基于特征线的SPH法。数值例题显示了所建议方案在模拟颗粒材料与间隙流相互作用的能力和性能以及间隙流体对颗粒结构承载能力及变形的影响。     

An elastoplastic model for unsaturated rocks and concrete

This paper presents an elastoplastic model for unsaturated rocks and concrete. The formulation is based on the poroplastic theory of porous media. The effect of water saturation on plastic deformation is described through a generalised effective stress. Plastic shrinkage and swelling due to suction change are taken into account. All model’s parameters can be determined from specific laboratory tests. Comparisons between numerical simulations and experimental data are provided for a claystone.