A Micropolar Theory of Porous Media: Constitutive Modelling

The extension of the classical mixture theory by the concept of volume fractions leads to the theory of porous media. In this article, the theory of porous media is generalised to micropolar constituents. The kinematic relations and the balance equations for a porous medium are developed without restricting the number of constituents. Based on the entropy inequality, the general form of the constitutive equations are derived for a binary medium consisting of a porous elastic skeleton saturated by a viscous pore-fluid. Both constituents are assumed to be compressible. Handling the saturation constraint by a Lagrangian multiplier leads to a compatibility of the proposed model to so-called hybrid and incompressible models.     

A Model for Flow and Deformation in Unsaturated Swelling Porous Media

A thermomechanical theory for multiphase transport in unsaturated swelling porous media is developed on the basis of Hybrid Mixture Theory (saturated systems can also be modeled as a special case of this general theory). The aim is to comprehensively and non-empirically describe the effect of viscoelastic deformation on fluid transport (and vice versa) for swelling porous materials. Three phases are considered in the system: the swelling solid matrix s, liquid l, and air a. The Coleman–Noll procedure is used to obtain the restrictions on the form of the constitutive equations. The form of Darcy’s law for the fluid phase, which takes into account both Fickian and non-Fickian transport, is slightly different from the forms obtained by other researchers though all the terms have been included. When the fluid phases interact with the swelling solid porous matrix, deformation occurs. Viscoelastic large deformation of the solid matrix is investigated. A simple form of differential-integral equation is obtained for the fluid transport under isothermal conditions, which can be coupled with the deformation of the solid matrix to solve for transport in an unsaturated system. The modeling theory thus developed, which involves two-way coupling of the viscoelastic solid deformation and fluid transport, can be applied to study the processing of biopolymers, for example, soaking of foodstuffs and stress-crack predictions. Moreover, extension and modification of this modeling theory can be applied to study a vast variety of problems, such as drying of gels, consolidation of clays, drug delivery, and absorption of liquids in diapers.     

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.     

An elastoplasticity model in porous media theories

In this article, porous media theories are referred to as mixture theories extended by the well-known concept of volume fractions. This approach implies the diverse field functions of both the porous solid matrix and the pore fluid to be represented by average functions of the macroscale.The present investigations are based on a binary model of incompressible constituents, solid skeleton, and pore liquid, where, in the constitutive range, use is made of the second-grade character of general heterogeneous media. Within the framework of geometrically finite theories, the paper offers a set of constitutive equations for the solid matrix, the viscous pore liquid and the different interactions between the constituents. The constitutive model applies to saturated as well as to empty solid materials, taking into account the physical nonlinearities based on elasto-plastic solid deformations. In particular, the constitutive model concentrates on granular materials like soil or concrete, where the elastic deformations are usually small and the plastic range is governed by kinematically hardening properties.     

CONSTITUTIVE RELATION OF UNSATURATED SOIL BY USE OF THE MIXTURE THEORY (Ⅱ)—LINEAR CONSTITUTIVE EQUATIONS AND FIELD EQUATIONS

IntroductionInpart(Ⅰ )ofthework[1],byuseofmixturetheory ,thenonlinearconstitutiveequationsandthefieldequationsofunsaturatedsoilwereconstructed ,andthecompleteequationsforthethermodynamicsystemofunsaturatedsoilwasformed .Inthispart,thelinearconstitutiveequationsandfieldequationsofunsaturatedsoilareobtainedthroughlinearizingnonlinearequations,andthelinearequationsarewrittenintheformssimilartoBiot’sequationsforsaturatedporousmedia .ItisprovedthatDarcy’slawissuitabletodescribethemotionofliquid…     

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.     

A linear dynamic model for a saturated porous medium

A linear isothermal dynamic model for a porous medium saturated by a Newtonian fluid is developed in the paper. In contrast to the mixture theory, the assumption of phase separation is avoided by introducing a single constitutive energy function for the porous medium. An important advantage of the proposed model is it can account for the couplings between the solid skeleton and the pore fluid. The mass and momentum balance equations are obtained according to the generalized mixture theory. Constitutive relations for the stress, the pore pressure are derived from the total free energy accounting for inter-phase interaction. In order to describe the momentum interaction between the fluid and the solid, a frequency independent Biot-type drag force model is introduced. A temporal variable porosity model with relaxation accounting for additional attenuation is introduced for the first time. The details of parameter estimation are discussed in the paper. It is demonstrated that all the material parameters in our model can be estimated from directly measurable phenomenological parameters. In terms of the equations of motion in the frequency domain, the wave velocities and the attenuations for the two P waves and one S wave are calculated. The influences of the porosity relaxation coefficient on the velocities and attenuation coefficients of the three waves of the porous medium are discussed in a numerical example.     

On mechanics of saturated granular materials

A continuum theory of saturated granular materials is formulated. The basic balance laws for the solid phase as well as for the fluid phase are presented. The constitutive equations are derived and the basic equations of motion of the solid and fluid continua are obtained. Several cases of interest, such as incompressible granules saturated with liquids are discussed. It is shown that the theory contains, as its special cases, the Mohr-Coulomb criterion for a granular material as well as Darcy's law of flow through porous media.     

An analytical homogenization method for heterogeneous porous materials

The objective of this work is to develop an analytical homogenization method to estimate the effective mechanical properties of fluid-filled porous media with periodic microstructure. The method is based on the equivalent inclusion concept of homogenization applied earlier for solid–solid mixture. It is assumed that porous media are described by the poroelastic constitutive law developed by Biot where porosity is a material parameter. By solving the governing equations of poroelasticity in Fourier transformed domain, the relation between periodic strain and eigenstrain in porous media is established. This relation is subsequently used in an average consistency condition involving both solid and fluid phase stresses and strains. The geometry of the porous microstructure is captured in the g-integral. This homogenization method can also be applied to estimate the equivalent properties of solid–fluid mixture where a pure solid and fluid can be modeled by assuming very low and high porosity, respectively. Several examples are considered to establish this new method by comparing with other existing analytical and numerical methods of homogenization. As an application, poroelastic properties of cortical bone fibril are estimated and compared with previously computed values.     

两相混合流体的非线性本构理论

本文从连续介质力学的基本原理出发,建立了微极流体与经典流体两相流动的非线性扩散理论。给出了混合流体本构方程的一般形式。对单相流体、单相微极流体及稀悬浮体三种特殊情形,得到了具体形式的二阶非线性本构方程,并同已有的理论进行了比较。     

流体饱和两相多孔介质拟静态问题的有限元解法

给出基于混合物理论的流体饱和两相多孔介质模型,该模型由一可变形固体 一流体相组成。采用Ｇａｌｅｒｋｉｎ加权残值法导出求解拟静态问题的有限元公式,并编制了二维有限元程序。用程序分析了一维和二维问题,得到合理的结果。     

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

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

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

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

INSTABILITY AND DISPERSIVITY OF WAVE PROPAGATION ININELASTIC SATURATED/UNSATURATED POROUS MEDIA

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｒａｎｓｉｅｎｔｐｈｅｎｏｍｅｎａｏｆｄｙｎａｍｉｃｐｒｏｂｌｅｍｓａｎｄｔｈｅｍａｎｎｅｒｓｔｏｔｒｅａｔｔｈｅｍａｒｅｇｅｎｅｒａｌｌｙｒｅｌａｔｅｄｔｏｔｈｅｓｐｅｃｔｒａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｅｘｃｉｔａｔｉｏｎ .Ｏｎｅｍａｙｏｍｉｔｐｒｏｐａｇａｔｉｎｇｗａｖｅｓｉｆｒｅｌａｔｉｖｅｌｙｌｏｗｅｒｆｒｅｑｕｅｎｃｉｅｓｇｏｖｅｒｎｔｈｅｒｅｓｐｏｎｓｅｓｉｎｔｈｅｐｒｏｂｌｅｍｓｓｕｃｈａｓｓｅｉｓｍｉｃｒｅｓｐｏｎｓｅｓａｎｄｒｅｓｐｏｎ…     

Mixture modeling of jointed fluidsaturated porous media

The quasi-static equations of motion are studied for bi-laminated fluid-saturated porous media within the framework of non-phenomenological mixture theories. The flow-deformation coupled behavior of the media is governed by Biot's theory for which all constituents are considered compressible. The asymptotic analysis for a periodic microstructure with multiple scales, developed by Hegemier and Murakami, is adopted to obtain the equations of equilibrium and mass conservation in a binary saturated porous medium. The multiscale analysis appears to be advantageous for dealing with consolidation phenomena because it is capable of transforming a coupled, transient problem into two decoupled, steady-state ones. Various models with different degrees of approximation are generated, and among them a theory for saturated rocks with a single joint system is described. Mixture properties are expressed explicitly in terms of characteristics of intact and joint material. The most distinctive feature of this model comes from the fact that some cross terms, that have not been included in previous models, appear in the constitutive equations for fluid mass change and fluid flux. These cross terms are physically understood because they simply take into account effects occurring on the local level: the deformation-flow coupled phenomenon, the stress continuity and displacement compatibility conditions. These novel results may have far-reaching consequences for future theoretical modeling and experimental programs in two-phase fluid-filled porous media.     

Fully-developed conjugate heat transfer in porous media with uniform heating

We propose a computational method for approximating the heat transfer coefficient of fully-developed flow in porous media. For a representative elementary volume of the porous medium we develop a transport model subject to periodic boundary conditions that describes incompressible fluid flow through a uniformly heated porous solid. The transport model uses a pair of pore-scale energy equations to describe conjugate heat transfer. With this approach, the effect of solid and fluid material properties, such as volumetric heat capacity and thermal conductivity, on the overall heat transfer coefficient can be investigated. To cope with geometrically complex domains we develop a numerical method for solving the transport equations on a Cartesian grid. The computational method provides a means for approximating the heat transfer coefficient of porous media where the heat generated in the solid varies “slowly” with respect to the space and time scales of the developing fluid. We validate the proposed method by computing the Nusselt number for fully developed laminar flow in tubes of rectangular cross section with uniform wall heat flux. Detailed results on the variation of the Nusselt number with system parameters are presented for two structured models of porous media: an inline and a staggered arrangement of square rods. For these configurations a comparison is made with literature on fully-developed flows with isothermal walls.     

Wave propagation in saturated ground using the volume fractions

The purpose of this paper is to study the dynamic behavior of soft ground including a porous layer by considering the porosity change. In order to take the porosity change into account, the concept of the volume fraction, which has been proposed in continuum mechanics, is introduced. The constitutive equations presented by Bowen are applied to the analysis of the porous media. According to Bowen's theory, the porosity is considered as a variable called the volume fraction and has its own constitutive equation. The constitutive equation of the volume fraction has thermoelastic equation coefficients and is determined by the strains of the solid and the fluid. This means that the compressibilities of the solid and the fluid are considered. When the special condition is assumed, Bowen's theory can contain Biots's theory, which has been applied in earthquake engineering. The wave propagation in the ground including a porous layer, modeled by Bowen's theory, is studied and compared with that of Biot's theory. One-dimensional attenuation and surface amplitude are calculated. The effect of the volume fraction is discussed with respect to the compressibilities of the solid and the fluid.     

Saturated Elastic Porous Solids: Incompressible,Compressible and Hybrid Binary Models

Constitutive models for a general binary elastic-porous media are investigated by two complementary approaches. These models include both constituents treated as compressible/incompressible, a compressible solid phase with an incompressible fluid phase (hybrid model of first type), and an incompressible solid phase with a compressible fluid phase (hybrid model of second type). The macroscopic continuum mechanical approach uses evaluation of entropy inequality with the saturation condition always considered as a constraint. This constraint leads to an interface pressure acting in both constituents. Two constitutive equations for the interface pressure, one for each phase, are identified, thus closing the set of field equations. The micromechanical approach shows that the results of Didwania and de Boer can be easily extended to general binary porous media.     

CONSTITUTIVE RELATION OF UNSATURATED SOIL BY USE OF THE MIXTURE THEORY (Ⅱ)-LINEAR CONSTITUTIVE EQUATIONS AND FIELD EQUATIONS

The linear constitutive equations and field equations of unsaturated soils were obtained through linearizing the nonlinear equations given in the first part of this work. The linear equations were expressed in the forms similar to Biot’s equations for saturated porous media. The Darcy’s laws of unsaturated soil were proved. It is shown that Biot’s equations of saturated porous media are the simplification of the theory. All these illustrate that constructing constitutive relation of unsaturated soil on the base of mixture theory is rational.