Crystallization,pore relaxation and micro-cryosuction in cohesive porous materials

A poro-elastic analysis is undertaken to account for the pressure time history of water-infiltrated pores within a material subjected to freezing. The thermodynamic-mechanical equilibrium of undercooled water and ice crystal, and Poiseuille-like flow through the connection channels, combine to reveal three successive mechanisms: in-pore crystallization, in-pore partial melting and a micro-cryosuction process, driving liquid water from the yet unfrozen pores to the frozen sites. The model turns out to be apt to predict the macroscopic relaxation process observed at the onset of crystallization as reported in the literature for cement-based materials. To cite this article: O. Coussy, T. Fen-Chong, C. R. Mecanique 333 (2005).     

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.     

Elastic behavior of saturated porous materials under undrained freezing

The elastic behavior of saturated porous materi- als under undrained freezing is investigated by using a poro- mechanical approach. Thermodynamic equilibria are used to describe the crystallization process of the partially frozen solution in bulk state and confined state in pores. By phase transition at freezing, fusion energy, thermal contraction of solid, solution and ice crystals, volume changes of crystallization build up remarkable pore pressure that induces expansion or shrinkage of solid matrix. Owing to the lower chemical potential when pore water mixes with salts, fewer ice forms in pores. Penetration of ice into the porous materials increases the capillary pressure, but limits effect on the pore liquid pressure and the strain of solid matrix. On the contrary, the pore pressure induced by solution density rises as salt concentration increases and causes significant shrinkage of solid matrix.     

Freezing and thawing porous media: experimental study with a dielectric capacitive method

A capacitive sensor-based apparatus has been used to study the ice/water phase change in consolidated porous media subjected to freezing and thawing. This technique relies on the dielectric properties of water, ice, air, and the mineral substrate in the radio-frequency range. It gives directly the freezing and thawing temperature depressions and indirectly provides an estimation of pore size distribution through the Gibbs–Thomson relation. It also holds good promise for evaluating the amount of liquid water in frozen porous media by combining drying and freezing tests. To cite this article: T. Fen-Chong, A. Fabbri, C. R. Mecanique 333 (2005).     

Heat and mass transfer in unsaturated porous media with solid–liquid change

 A model on heat and mass transfer in unsaturated porous media with solid/liquid phase change was developed with extending the three-variable model previously proposed. The movement of air phase and its effect on the motion of water is considered. The model has been checked with comparison of the experimental results of the temperature distribution for two dimensional freezing process. The evolution of air pressure, water and ice saturation were predicted by solving the governing equations. The ice segregation and moisture movement toward the front of freezing were numerically simulated. Received on 8 December 1999     

Dependency of Tortuosity and Permeability of Porous Media on Directional Distribution of Pore Voids

Single-phase fluid flow in porous media is usually direction dependent owing to the tortuosity associated with the internal structures of materials that exhibit inherent anisotropy. This article presents an approach to determine the tortuosity and permeability of porous materials using a structural measure quantifying the anisotropic distribution of pore voids. The approach uses a volume averaging method through which the macroscopic tortuosity tensor is related to both the average porosity and the directional distribution of pore spaces. The permeability tensor is derived from the macroscopic momentum balance equation of fluid in a porous medium and expressed as a function of the tortuosity tensor and the internal structure of the material. The analytical results generally agree with experimental data in the literature.     

Stress-induced crystal preferred orientation in the poromechanics of in-pore crystallization

A non-hydrostatic stress field affects the orientation of crystals growing in the pore network of an elastic porous medium. The hypothesis of a hydrostatic state of stress within the crystal has been implicitly made in the recent extension of poromechanics to in-pore crystalization (Coussy, 2006). This underlying hypothesis is revisited on a small-scale conceptual model based on Eshelby's problem and shows that chemo-mechanical equilibrium requires that the crystal adapts its shape and orientation to the far-field stress, therefore resulting at equilibrium in a hydrostatic state of stress within the crystal. The optimum crystal shape as a function of the far-field stress is consistently investigated, highlighting limiting cases. The small scale model allows to understand the macroscopic effects associated with deviatoric stresses in the poromechanics of in-pore crystallization. Moreover, it provides the building block for an up-scaling of the macroscopic tangent poroelastic properties, which depend on both the current crystal saturation and the state of stress. A dilute micromechanical scheme illustrates the variation of the macroscopic Biot's coefficient tensor as a function of deviatoric stresses. A simple configuration akin to a potential laboratory experiment finally illustrates the strong induced anisotropy of the crystallization induced macroscopic strain when deviatoric stresses are applied to the material prior to crystallization.     

多孔介质中稳定泡沫的流动计算及实验研究

将多孔介质简化为一簇变截面毛管束,根据多孔介质的颗粒直径、颗粒排列方式、孔喉尺度比及束缚水饱和度,计算出变截面毛细管的喉道半径和孔隙半径. 在考虑多孔介质喉道和孔隙中单个气泡的受力和变形基础上,利用动量守恒定理,推导出单个孔隙单元内液相的压力分布和孔隙单元两端的压差计算公式,最终得到多孔介质的压力分布计算公式. 利用长U型填砂管对稳定泡沫的流动特性进行了实验研究. 研究结果表明:稳定泡沫流动时多孔介质中的压力分布呈线性下降,影响泡沫在多孔介质中流动特性的因素包括:多孔介质的孔喉结构、泡沫流体的流量和干度、气液界面张力、气泡尺寸,其中孔喉结构和泡沫干度是影响泡沫封堵能力的主要因素.关键词: 稳定泡沫;多孔介质;变截面毛管;流动;表观粘度;压力分布;实验研究      

Analysis of the Propagation of Sound Waves in Partially Saturated Soils by Means of a Macroscopic Linear Poroelastic Model

In this paper, the propagation of sound waves in partially saturated soils is investigated. A macroscopic linear model that is based on the two-component model of Biot and on the Simple Mixture Model by Wilmanski is used. For the construction of the model by a micro-macro transition, see Albers, Géotechnique, 2007. We investigate a porous medium consisting of a deformable skeleton and two compressible, chemically non-reacting, pore fluids (liquid and gas). The wave analysis of the poroelastic model reveals the number of acoustic waves and the dependence of velocities and attenuations of these waves on the initial saturation and frequency. There appear four body waves: three longitudinal waves, P1, P2, P3, and one shear wave, S. The P2-wave shows a similar feature as in air–water mixtures: from the theory of suspensions, it is well known that the existence of air bubbles in water reveals a minimum in the sonic velocity. This is also the case for the P2 -speed in the unsaturated porous medium. The P1-velocity increases very abruptly for a certain degree of saturation. This provides the hope for the development of a nondestructive testing method.     

Plastic potential and yield function of porous materials with aligned and randomly oriented spheroidal voids

Based on an energy approach, the plastic potential and yield function of a porous material containing either aligned or randomly oriented spheroidal voids are developed at a given porosity and pore shape. The theory is applicable to both elastically compressible and incompressible matrix and, it is proved that, in the incompressible case, the theory with spherical and aligned spheroidal voids also coincides with Ponte Castaneda's bounds of the Hashin-Shtrikman and Willis types, respectively. Comparison is also made between the present theory and those of Gurson and Tvergaard, with a result giving strong overall support of this new development. For the influence of pore shape, the yield function and therefore the stress-strain curve of the isotropic porous material are found to be stiffest when the voids are spherical, and those associated with other pore shapes all fall below these values, the weakest one being caused by the disc-shaped voids. The transversely isotropic nature of the yield function and stress-strain curves of a porous material containing aligned pores are also demonstrated as a function of porosity and pore shape, and it is further substantiated with a comparison with an exact, local analysis when the void shape becomes cylindrical.     

轴向扩散下简支饱和多孔弹性梁的大挠度分析

基于多孔介质理论和弹性梁的大挠度理论,并考虑轴向变形,在孔隙流体仅沿轴向扩散的假设下,建立了微观不可压饱和多孔弹性梁大挠度弯曲变形的一维非线性数学模型.在此基础上,忽略饱和多孔弹性梁的轴向应变,并利用Galerkin截断法,研究了两端可渗透的简支饱和多孔弹性梁在突加横向均布载荷作用下的拟静态弯曲,给出了饱和多孔梁弯曲时挠度、弯矩和轴力以及孔隙流体压力等效力偶等沿轴线的分布曲线.揭示了大挠度非线性和小挠度线性模型的结果差异,指出大挠度非线性模型的结果小于相应小挠度线性模型的结果,并且这种差异随着载荷的增大而增大.计算表明:当无量纲载荷参数q>5时,应该采用大挠度非线性数学模型进行研究.     

Forced convection freezing characteristics along the convex wall of a return bend with rectangular cross section

An experimental study was made of the forced convection freezing characteristics on the convex wall of a return bend with a rectangular cross section. Observations were carried out for duct heights of 17 and 30 mm, a duct width of 300 mm, and a radius of curvature of 159 mm. The convex wall temperature was uniformly maintained below the freezing temperature of water, and the concave wall was insulated. It was found that a stepwise ice layer forms on the convex wall of a return bend and that the step position at the steady-state condition is closely dependent on both the water flow velocity and the cooling temperature ratio.     

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.     

简支饱和多孔弹性梁的非线性弯曲

基于饱和多孔介质理论和弹性梁的大挠度弯曲假设,在多孔弹性梁轴线不可伸长,孔隙流体仅沿轴向方向扩散的限制下,建立了微观不可压饱和多孔弹性梁大挠度拟静态响应的一维非线性数学模型.在此基础上,利用Galerkin截断法,分析了两端可渗透的简支多孔弹性梁在突加横向均布载荷作用下的非线性弯曲,给出了梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶随时间的响应曲线.数值结果表明:当载荷较小时,大挠度非线性与小挠度线性理论的结果相差很小,而当载荷较大时,非线性大挠度理论的结果小于相应线性小挠度理论的结果,并且这种差异随着载荷的增大而增大.同时,在载荷突加于梁上时,多孔弹性梁骨架起初不变形,孔隙流体压力等效力偶由零突增为非零,其值与外载荷保持平衡.随着时间的增加,梁的挠度增加,等效力偶逐渐减小为零,最终多孔梁骨架承担全部的外载荷.     

饱和多孔弹性Timoshenko梁的大挠度分析

基于微观不可压饱和多孔介质理论和弹性梁的大挠度变形假设,考虑梁剪切变形效应,在梁轴线不可伸长和孔隙流体仅沿轴向扩散的限定下,建立了饱和多孔弹性Timoshenko梁大挠度弯曲变形的非线性数学模型.在此基础上,利用Galerkin截断法,研究了两端可渗透简支饱和多孔Timoshenko梁在突加均布横向载荷作用下的拟静态弯曲,给出了饱和多孔 Timoshenko梁弯曲变形时固相挠度、弯矩和孔隙流体压力等效力偶等随时间的响应.比较了饱和多孔Timoshenko梁非线性大挠度和线性小挠度理论以及饱和多孔 Euler-Bernoulli梁非线性大挠度理论的结果,揭示了他们间的差异,指出当无量纲载荷参数q＞l0时,应采用饱和多孔Timoshenko梁或Euler-Bernoulli梁的大挠度数学模型进行分析,特别的,当梁长细比λ＜30时,应采用饱和多孔Timoshenko梁大挠度数学模型进行分析.     

Dynamic void nucleation and growth in solids: A self-consistent statistical theory

We present a framework for a self-consistent theory of spall fracture in ductile materials, based on the dynamics of void nucleation and growth. The constitutive model for the material is divided into elastic and “plastic” parts, where the elastic part represents the volumetric response of a porous elastic material, and the “plastic” part is generated by a collection of representative volume elements (RVEs) of incompressible material. Each RVE is a thick-walled spherical shell, whose average porosity is the same as that of the surrounding porous continuum, thus simulating void interaction through the resulting lowered resistance to further void growth. All voids nucleate and grow according to the appropriate dynamics for a thick-walled sphere made of incompressible material. The macroscopic spherical stress in the material drives the response in all volume elements, which have a distribution of critical stresses for void nucleation, and the statistically weighted sum of the void volumes of all RVEs generates the global porosity. Thus, macroscopic pressure, porosity, and a distribution of growing microscopic voids are fully coupled dynamically. An example is given for a rate-independent, perfectly plastic material. The dynamics of void growth gives rise to a rate effect in the macroscopic material even though the parent material is rate independent.     

饱和多孔弹性Timoshenko悬臂梁的拟静力弯曲

在经典单相Timoshenko梁变形和孔隙流体仅沿多孔梁轴向运动的假定下,基于不可压饱和多孔介质的三维理论,本文首先建立了横观各向同性饱和多孔弹性Timoshenko悬臂梁拟静力弯曲的一维数学模型,并给出了相应的边界条件。其次,利用Laplace变换及其数值逆变换,分析了端部不同渗透条件下,饱和多孔弹性Timoshenko悬臂梁在端部梯载荷作用下的拟静力响应,给出了饱和多孔Timoshenko悬臂梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶等随时间的响应曲线,并与饱和多孔Euler-Bernoulli悬臂梁的响应进行了比较,考察了梁长细比对弯曲的影响。数值结果表明：固相骨架与孔隙流体的相互作用具有粘性效应,梁弯曲的拟静态挠度具有蠕变行为,端部渗透条件对梁的弯曲响应有显著的影响,并且,饱和多孔弹性Timoshenko悬臂梁的拟静态响应亦存在Mandel-Cryer现象。     

不可压饱和多孔弹性简支梁的动力响应

在杆件弯曲小变形的假定下,考虑杆件的侧向变形因素,根据多孔介质理论,本文首先建立了不可压饱和多孔弹性梁弯曲变形时动力响应的控制方程。其次,基于所建立的控制微分方程,利用变量分离法,研究了两端可渗透的饱和多孔弹性简支梁在梁中间集中载荷作用下的动力响应,得到了不同物性参数下简支梁动态弯曲时挠度和孔隙流体压力等效力偶等随时间的响应曲线。研究发现由于孔隙流体和固相骨架的相互作用,不可压饱和多孔弹性梁挠度的动力响应具有粘性特征,同时,随着时间的增加,饱和多孔弹性梁的挠度、弯矩等最终趋于经典弹性梁的静挠度、弯矩,此时,孔隙流体压力为零,梁的固相骨架承担所有的外载荷。     

A stochastic model for the permeability characteristics of saturated cemented porous media undergoing freezing

As the temperature of a saturated porous medium drops, the water in the pores starts to freeze. Since the temperature at which the phase change takes place is dependent on the pore size, the permeability of the medium changes continuously. Simultaneously, due to the expansion of water on freezing, it is forced to migrate through the pore body thus inducing stresses in material matrix. The stresses developed and the consequent frost damage are therefore dependent on the change in the permeability characteristics of the medium on freezing. This paper deals with the numerical prediction of permeability characteristics of porous cemented media saturated with water undergoing progressive freezing.A bond percolation model is used to generate the pore structure according to an assumed poresize distribution. Permeability of the medium at various temperatures is computed by solving the network problem. The computed results are compared with other analytical and experimental results. The proposed model predicts a threshold temperature below which permeability drops to zero. This phenomenon is crucial in developing a deeper understanding of the mechanism of frost damage to cemented porous materials such as bricks, stone, concrete, etc.