多孔介质平板通道强迫对流中热局部非平衡时的热应力

根据微观不可压饱和多孔介质热-力-流相互作用的一般理论,在固相骨架小变形的假定下,考虑固相和流相相互作用的粘性耗散,研究了多孔介质平板通道强迫对流热局部非平衡的热应力问题.建立了问题的热-力数学模型,根据饱和多孔介质的平衡方程,在固相骨架只存在横向位移的假定下,求解了固相骨架的位移和相应的热应力,数值考察了各种物性参数对热应力分布的影响,讨论了热局部平衡模型的适用性.     

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

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

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

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

GURTIN VARIATIONAL PRINCIPLE AND FINITE ELEMENT SIMULATION FOR DYNAMICAL PROBLEMS OF FLUID-SATURATED POROUS MEDIA

Based on the theory of porous media,a general Gurtin variational principle for theinitial boundary value problem of dynamical response of fluid-saturated elastic porous media isdeveloped by assuming infinitesimal deformation and incompressible constituents of the solid andfluid phase.The finite element formulation based on this variational principle is also derived.Asthe functional of the variational principle is a spatial integral of the convolution formulation,thegeneral finite element discretization in space results in symmetrical differential-integral equationsin the time domain.In some situations,the differential-integral equations can be reduced to sym-metrical differential equations and,as a numerical example,it is employed to analyze the reflectionof one-dimensional longitudinal wave in a fluid-saturated porous solid.The numerical results canprovide further understanding of the wave propagation in porous media.     

Modeling and Simulation of Local Thermal Non-equilibrium Effects in Porous Media with Small Thermal Conductivity

The two-equation model in porous media can describe the local thermal non-equilibrium (LTNE) effects between fluid and solid at REV scale, with the temperature differences in a solid particle neglected. A multi-scale model has been proposed in this study. In the model, the temperature differences in a solid particle are considered by the coupling of the fluid energy equation at REV scale with the heat conduction equation of a solid particle at pore scale. The experiments were conducted to verify the model and numerical strategy. The multi-scale model is more suitable than the two-equation model to predict the LTNE effects in porous media with small thermal conductivity. The effects of particle diameter, mass flow rate, and solid material on the LTNE effects have been investigated numerically when cryogenic nitrogen flows through the porous bed with small thermal conductivity. The results indicate that the temperature difference between solid center and fluid has the same trend at different particle diameters and mass flow rates, while the time to reach the local thermal equilibrium is affected by solid diameter dramatically. The results also show that the temperature difference between solid center and surface is much greater than that between solid surface and fluid. The values of \( \rho {\text{c}} \) for different materials have important influence on the time to reach the local thermal equilibrium between solid and fluid.     

Boundary conditions for porous solids saturated with viscous fluid

Boundary conditions are derived to represent the continuity requirements at the boundaries of a porous solid saturated with viscous fluid. They are derived from the physically grounded principles with a mathematical check on the conservation of energy. The poroelastic solid is a dissipative one for the presence of viscosity in the interstitial fluid. The dissipative stresses due to the viscosity of pore-fluid are well represented in the boundary conditions. The unequal particle motions of two constituents of porous aggre~ gate at a boundary between two solids are explained in terms of the drainage of pore-fluid leading to imperfect bonding. A mathematical model is derived for the partial connec- tion of surface pores at the porous-porous interface. At this interface, the loose-contact slipping and partial pore opening/connection may dissipate a part of strain energy. A numerical example shows that, at the interface between water and oil-saturated sandstone, the modified boundary conditions do affect the energies of the waves refracting into the isotropic porous medium.     

The Onset of Convection in an Anisotropic Porous Layer Using a Thermal Non-Equilibrium Model

The stability of a horizontal fluid saturated anisotropic porous layer heated from below and cooled from above is examined analytically when the solid and fluid phases are not in local thermal equilibrium. Darcy model with anisotropic permeability is employed to describe the flow and a two-field model is used for energy equation each representing the solid and fluid phases separately. The linear stability theory is implemented to compute the critical Rayleigh number and the corresponding wavenumber for the onset of convective motion. The effect of thermal non-equilibrium and anisotropy in both mechanical and thermal properties of the porous medium on the onset of convection is discussed. Besides, asymptotic analysis for both very small and large values of the interphase heat transfer coefficient is also presented. An excellent agreement is found between the exact and asymptotic solutions. Some known results, which correspond to thermal equilibrium and isotropic porous medium, are recovered in limiting cases.     

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.     

Lattice Boltzmann Pore Scale Simulation of Natural Convection in a Differentially Heated Enclosure Filled with a Detached or Attached Bidisperse Porous Medium

In this research, pore scale simulation of natural convection in a differentially heated enclosure filled with a conducting bidisperse porous medium is investigated using the thermal lattice Boltzmann method. For the first time, the effect of connection of the bidisperse porous medium to the enclosure walls is studied by considering the attached geometry in addition to the detached one. Effect of most relevant parameters on the streamlines and isotherms as well as hot wall average Nusselt number is studied for two of the bidisperse porous medium configurations. It is observed that effect of geometrical and thermo-physical parameters of the bidisperse porous medium on the heat transfer characteristics is more complicated for the attached configuration. To assess the validity of the local thermal equilibrium condition in the micro-porous media, the pore scale results are used to compute the percentage of the local thermal non-equilibrium for two of the bidisperse porous medium configurations. It is concluded that for the detached configuration, the local thermal equilibrium condition is confirmed in the entire micro-porous media for the ranges of the parameters studied here. However, for the attached geometry, it is shown that departure from the local thermal equilibrium condition is observed for the higher values of the Rayleigh number, micro-porous porosity, solid–fluid thermal conductivity ratio, and the smaller values of the macro-pores volume fraction.     

A Local Thermal Non-Equilibrium Analysis of Fully Developed Forced Convective Flow in a Tube Filled with a Porous Medium

A local thermal non-equilibrium model has been considered for the case of thermally fully developed flow within a constant heat flux tube filled with a porous medium. Exact temperature profiles for the fluid and solid phases are found after combining the two individual energy equations and then transforming them into a single ordinary differential equation with respect to the temperature difference between the solid phase and the wall subject to constant heat flux. The exact solutions for the case of metal-foam and air combination reveal that the local thermal equilibrium assumption may fail for the case of constant heat flux wall. The Nusselt number is presented as a function of the Peclet number, which shows a significant increase due to both high stagnant thermal conductivity and thermal dispersion resulting from the presence of the metal-foam.     

Fixed domain methods for free and moving boundary flows in porous media

A survey is presented concerning fixed domain methods used to solve mathematical models of free and moving boundary flow problems in porous media. These include the following: variational inequality or quasi-variational inequality formulations; general inequality formulations which have been set and solved in fixed domains; and the residual flow procedure. Finally, some parallel computing methods and mesh adaptation methods are discussed to demonstrate how these fixed domain formulations can be solved with current technology.The fixed domain methods that are referenced herein can be classified into two groups: the variational inequality method and the extended pressure head method. Baiocchi was the first to apply the variational inequality method to free boundary problems of flows through porous media. This method in general also uses an extension of the pressure head but adds an application of an integral transformation (a Baiocchi transformation) to the problem. The method possesses a beautiful mathematical structure for its theory and yields simple numerical solution algorithms. However, application of the method is difficult if not impossible in some cases depending upon the regularity of the seepage domain.The extended pressure head method is based on the concept that the pressure is extended smoothly across the free or moving boundary into the unsaturated region from the flow domain. The extension of the pressure head to the entire porous medium yields an extended coefficient of permeability of the medium which is equal to the saturated coefficient in the seepage region and is equal to zero or some small value (for computational purposes) in the unsaturated region.     

流体饱和多孔介质的动力学Gurtin型变分原理和有限元模拟

基于多孔介质理论,在两相不可压和小变形的假设下,建立了流体饱和弹性多孔介质的动力学Gurtin型变分原理,并导出了以此变分原理为基础的有限元离散公式.由于Gurtin型变分原理是卷积型的空间积分泛函,空间的有限元离散导致一个关于时间的对称微分-积分方程组.在一般条件下,该积分-微分方程组可转化为对称的微分方程组,这组方程有别于标准Galerkin有限元的非对称离散方程组.作为数值例子,分析了流体饱和弹性多孔介质中一维纵向波的传播和反射,其结果进一步揭示了饱和多孔介质中波的传播特性.     

Thermo-Osmosis Effect in Saturated Porous Medium

In this paper, the thermo-poroelasticity theory is used to investigate the quasi-static response of temperatures, pore pressure, stress, displacement, and fluid flux around a cylindrical borehole subjected to impact thermal and mechanical loadings in an infinite saturated poroelastic medium. It has been reported in literatures that coupled flow known as thermo-osmosis by which flux is driven by temperature gradient, can significantly change the fluid flux in clay, argillaceous and many other porous materials whose permeability coefficients are very small. This study presents a mathematical model to investigate the coupled effect of thermo-osmosis in saturated porous medium. The energy balance equations presented here fulfill local thermal non-equilibrium condition (LTNE) which is different from the local thermal equilibrium transfer theory, accounting for that temperatures of solid and fluid phases are not the same and governed by different heat transfer equations. Analytical solutions of temperatures, pore pressure, stress, displacement, and fluid flux are obtained in Laplace transform space. Numerical results for a typical clay are used to investigate the effect of thermo-osmosis. The effects of LTNE on temperatures, pore pressure, and stress are also studied in this paper.     

Self-similar problem of decomposition of gas hydrates in a porous medium upon depression and heating

We consider the specifics of decomposition of gas hydrates under thermal and depressive action on a porous medium completely filled with a solid hydrate in the initial condition. The existence of volumetric-expansion zones, in which the hydrate coexists in equilibrium with water and gas, is shown to be possible in high-permeable porous media. The self-similar problems of hydrate decomposition upon depression and heating are studied. Ii is shown that there are solutions according to which hydrate decomposition can occur both on the surface of phase transitions and in the volumetric region. We note that, in the first case, decomposition is possible without heat supply to a medium and even with heat removal. Institute of Mechanics of Multiphase Systems, Siberian Division, Russian Academy of Sciences, Tyumen' 625000. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 3, pp. 111–118, May–June, 1998.     

A Local Model for the Energy Transfer in a Saturated Static Mixture

In the present work the transient energy transfer in a nonsaturated porous medium is studied, using a mixture theory viewpoint. The porous matrix is assumed homogeneous, rigid and isotropic, while the fluid is a Newtonian incompressible one and both are assumed static. Since the homogeneous matrix is not saturated, gradients of concentration are present. The porous medium and the fluid (a liquid) will be regarded as continuous constituents of a mixture that will have also a third constituent, an inert gas, assumed with zero mass density and thermal conductivity. The problem is described by a set of two partial differential equations which represent the energy balances for the fluid and the solid constituents. Isovalues for these two constituents are plotted, considering representative time instants and selected values for the energy equations coefficients and for the saturation.     

Three Pressures in Porous Media

In a thermodynamic setting for a single phase (usually fluid), the thermodynamically defined pressure, involving the change in energy with respect to volume, is often assumed to be equal to the physically measurable pressure, related to the trace of the stress tensor. This assumption holds under certain conditions such as a small rate of deformation tensor for a fluid. For a two-phase porous medium, an additional thermodynamic pressure has been previously defined for each phase, relating the change in energy with respect to volume fraction. Within the framework of Hybrid Mixture Theory and hence the Coleman and Noll technique of exploiting the entropy inequality, we show how these three macroscopic pressures (the two thermodynamically defined pressures and the pressure relating to the trace of the stress tensor) are related and discuss the physical interpretation of each of them. In the process, we show how one can convert directly between different combinations of independent variables without re-exploiting the entropy inequality. The physical interpretation of these three pressures is investigated by examining four media: a single solid phase, a porous solid saturated with a fluid which has negligible physico-chemical interaction with the solid phase, a swelling porous medium with a non-interacting solid phase, such as well-layered clay, and a swelling porous medium with an interacting solid phase such as swelling polymers.     

CONSOLIDATION HEORY OF UNSATURATED SOIL BASED ON THE THEORY OF MIXTURE(Ⅰ)

Unsaturated soil is a three-phase media and is composed of soil grain,water andgas.In this paper,the consolidation problem of unsaturated soil is investigated basedon the theory of mixture.A theoretical formula of effective stress on anisotropicporous media and unsaturated soil is derived.The principle of effective stress and theprinciple of Curie symmetry are taken as two fundamental constitutive principles ofunsaturated soil.A mathematical model of consolidation of unsaturated soil isproposed,which consists of25 partial differenfial equations with25 unknowns.Withthe help of increament linearizing method,the model is reduced to5 governingequations with5 unknowns,i.e.,the three displacement components of solid phase,thepore water pressure and the pore gas pressure.7 material parameters are involved inthe model and all of them can be measured using soil tests.It is convenient to use themodel to engineering practice.The well known Biot’s theory is a special case of themodel.