ONSET CONDITION OF STRAIN LOCALIZATION IN MATRIX OF SATURATED POROUS MEDIA

Introduction Strainlocalizationofgeomaterialsisoneofmostpopularfailuretypesinnature,which canbeshowedaslandslidesandmudflowsinmountainousareasunderincessantorheavy raining,especiallythevegetationisseverelydamagedbywoodsharvest;pipingeffect,a typeoflocalfa…     

Advective fluxes in multiphase porous media under nonisothermal conditions

We consider the case in which more than one fluid phase occupies the void space of a porous medium. The advective flux law is formulated for a fluid phase, under nonisothermal conditions and with the presence of solutes in the fluid phases. The derivation of the flux laws is based on an approximated version of the averaged balance equation for linear momentum. Taking into account momentum transfer through the interface between the fluid phases, leads to coupling between the flow in adjacent phases. Fluxes are also shown to depend on the surface tension at the interface between the adjacent fluid phases. Since the latter depends on temperature and solute concentration in the two phases, the advective flux is shown to depend on both temperature and solute concentration gradients in the two phases. A preliminary order of magnitude analysis gives conditions under which the coupling phenomena are not negligible. The approach is applied to the unsaturated zone, as a typical example of a multiphase porous medium.The main conclusion is that the well known Darcy law for single phase flow, may have to be modified for a multi fluid phase system, especially when temperature and solute concentration are not uniform.     

A thermodynamical gradient theory for deformation and strain localization of porous media

In this work, a thermodynamically consistent gradient formulation for partially saturated cohesive-frictional porous media is proposed. The constitutive model includes a classical or local hardening law and a softening formulation with state parameters of non-local character based on gradient theory. Internal characteristic length in softening regime accounts for the strong shear band width sensitivity of partially saturated porous media regarding both governing stress state and hydraulic conditions. In this way the variation of the transition point (TP) of brittle-ductile failure mode can be realistically described depending on current confinement condition and saturation level. After describing the thermodynamically consistent gradient theory the paper focuses on its extension to the case of partially saturated porous media and, moreover, on the formulation of the gradient-based characteristic length in terms of stress and hydraulic conditions. Then the localization indicator for discontinuous bifurcation is formulated for both drained and undrained conditions.     

Perturbation growth and localization in fluid-saturated inelastic porous media under quasi-static loadings

Instabilities in inelastic saturated porous media are investigated here for general three-dimensional states under quasi-static loadings using a perturbation approach and focussing in particular on the two limiting cases of the onset of growth and of blowing up of perturbations.For associative flow rules for the skeleton, both onset of growth and blowing up of perturbations depend only on the underlying drained properties. Unbounded growth is obtained when the condition of localization for the underlying drained deformation (singularity of the drained acoustic tensor) is approached or just passed. Onset of growth has always a divergence growth character and critical conditions are always associated to the shortwavelength regime leading to the fact that the failure mode is expected to be a localized one.For non-associative behaviour of the skeleton we show in contrast that the onset of growth and unbounded growth may be defined either by the drained or undrained properties. One or the other depends on the details of the constitutive behaviour but also on the type of loadings. In particular, unbounded growth occurs when either the condition of localization under drained or undrained conditions is first passed. Transition from decaying to growing behaviour may have a divergence character or flutter-type character. Here the critical conditions are associated either to the shortwavelength or to the longwavelength regimes and therefore the failure mode may be localized or diffuse.The hierarchy between criticality of drained and undrained properties is analysed for a general class of constitutive equations and the results are fully and explicitly illustrated for saturated porous media with skeleton obeying Drucker-Prager like constitutive model.     

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.     

Surface porosity and permeability of porous media with a periodic microstructure

Various ways of determining the surface porosity, the relation between the porosity and the surface porosity and the representation of the permeability in terms of the characteristics of the microstructure of the porous medium are analyzed with reference to model porous media with a periodic microstructure. It is shown that it is necessary to distinguish between the geometric (scalar) and physical (tensor) suface porosities and that the geometric surface porosity, the physical surface porosity and the porosity are different characteristics of the porous medium.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 79–85, January–February, 1995.     

饱和多孔介质中的混合有限元法和有限应变下应变局部化分析

对基于Biot理论的饱和多孔介质中动力-渗流耦合分析提出了一个耦合场混合元．固相位移、应变和有效应力以及流相压力、压力梯度和Darcy速度在单元内均处理为独立变量分别插值．基于胡海昌-Washizu三变量广义变分原理给出的饱和多孔介质动力-渗流耦合问题控制方程的单元弱形式，导出了单元公式．进一步导出了考虑压力相关非关联塑性的非线性单元公式和发展了相应的一致性算法．对几何非线性分析，采用了共旋公式途径．数值结果例题显示所发展耦合场混合元模拟大应变下由应变软化引起以应变局部化为特征的渐进破坏现象的性能．     

Elasticoplastic conditions of the filtration of a liquid in porous media

The article discusses questions in the theory of filtration in porous media, taking account of elastic, elasticoplastic, and plastic deformations. Parameters are introduced to evaluate irreversible effects in petroleum- and water-bearing strata, i.e., coefficients of the change in the porosity and the permeability. Equations are derived for filtration under unsteady-state and steady-state working conditions of wells and galleries. Two limiting cases, which allow analytical solutions, are separated out. In the general case, the equations of elasticoplastic filtration conditions are solved on an electronic computer. The numerial calculations show that the predominating effect results from taking account of the irreversible change in the permeability, depending on the change of the pressure in the stratum.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 84–90, September–October, 1973.The author is grateful to V. N. Nikolaevskii for his evaluation of the work.     

On the properties of compressible gas flow in a porous media

The flow of an adiabatic gas through a porous media is treated analytically for steady one- and two-dimensional flows. The effect on a compressible Darcy flow by inertia and Forchheimer terms is studied. Finally, wave solutions are found which exhibit a cut-off frequency and a phase shift between pressure and velocity of the gas, with the velocity lagging behind the pressure.Nomenclature A area of tube for one-dimensional flow - B drag coefficient associated with Forchheimer term - c speed of sound - M Mach number - p ^* gas pressure - p dimensionless gas pressure - s coordinate along the axis of tube - t ^* time variable - t dimensionless time variable - V^* gas velocity in the porous media - V dimensionless gas velocity Greek Letters ratio of specific heat capacities - phase angle between gas pressure and velocity for linear waves - parameter indicating the importance of the inertia term - viscosity - _p natural frequency of the porous media - ^* gas density - dimensionless gas density - parameter indicating the importance of the Forchheimer term - porosity of porous media - velocity potential - stream function     

Flow of gas and water through flooded jointed and porous media under conditions of falling pressure

Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 168–171, May–June, 1990.     

Effect of the microstructure of porous media on relative phase permeabilities

Correlations between the characteristics of a porous medium and the capillary-pressure and phase-permeability curves are discussed.     

Effect of variable porosity on laminar convection in a uniformly heated vertical porous channel

The present investigation is devoted to the study of fully developed mixed convective flow through a vertical porous channel. The lateral variations of porosity and thermal diffusivity in the bed near the wall, are approximated by exponential functions. The correlation between permeability and porosity is brought through Kozney-Carman approximation. The volume averaged one dimensional low speed momentum equation proposed by Vafai is employed for the analysis of the problem. Results are obtained for steady heating of ascending cold fluid and steady cooling of ascending hot fluid. For the above physical situations it is observed that the heat transfer rate, and ratio of friction factor increases with increase in porous parameter, whereas the ratio of mass flow rate decreases with increase in porous parameter. The velocity profiles exhibit hydrodynamic channelling and peak velocity shifts towards the wall for higher values of the porous parameter. For steady heating of ascending could fluid increase in Rayleigh number enhances the heat transfer rate, and mass flow rate, while it reduces the ratio of friction factor. An opposite trend is observed for the case of steady cooling of ascending hot fluid.     

Effect of the microinhomogeneity of the medium on the flow law in porous media

A lattice percolation model of an inhomogeneous medium [6] introducing f(r), the probability density function of the microcapillaries with respect to their radius r, is employed. The calculation scheme described makes it possible to determine the dependence of the permeability K on the pressure gradient G for media with arbitrary f(r). It is shown that for inhomogeneous media the behavior of K(G) is mainly determined by the form of f(r). The question of the effect of the state of stress on the permeability of the medium is considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 84–94, March–April, 1989.     

Numerical experiments on convective heat transfer in water-saturated porous media at near-critical conditions

Fluid and heat flow at temperatures approaching or exceeding that at the critical point (374 °C for pure water, higher for saline fluids) may be encountered in deep zones of geothermal systems and above cooling intrusives. In the vicinity of the critical point the density and internal energy of fluids show very strong variations for small temperature and pressure changes. This suggests that convective heat transfer from thermal buoyancy flow would be strongly enhanced at near-critical conditions. This has been confirmed in laboratory experiments. We have developed special numerical techniques for modeling porous flow at near-critical conditions, which can handle the extreme nonlinearities in water properties near the critical point. Our numerical simulations show strong enhancements of convective heat transfer at near-critical conditions; however, the heat transfer rates obtained in the simulations are considerably smaller than data reported from laboratory experiments by Dunn and Hardee. We discuss possible reasons for this discrepancy and develop suggestions for additional laboratory experiments.     

Micro-mechanical modelling of strain induced porosity under generally compressive stress states

This paper is concerned with the nucleation and growth of voids in a titanium alloy undergoing high temperature deformations under generally compressive stress states typical of forging processes. A micro-mechanical model for void nucleation has been developed based on a debonding process between primary alpha particles and the beta matrix. The finite element model developed has been used to examine in detail the stress state sensitivity of void nucleation within the particle-matrix system. The results obtained are compared with other phenomenological approaches showing good agreement for most stress states, but giving different results for a range of compressive stress states. A continuum-level representation of the micro-mechanical results has been obtained and implemented into a finite element model. Cylindrical specimen compression tests have been carried out over the strain rate range 0.005-5.0s⁻¹ and temperature range 925–975°C under conditions of high specimen-die friction.Regions of stress triaxiality that are tensile in nature were therefore generated, and the specimens tested to an overall strain of 0.5 were sectioned, polished and etched. The resulting distributions of voids were quantified, and compared with those predicted using the finite element model discussed above. Good quantitative agreement was obtained both in terms of the magnitude of the area fractions of voids and their distributions. The model also captures reasonably well the strain rate and temperature dependence of the voiding. However, the model assumptions of uniform distributions of alpha particles which are all perfectly spherical and with identical interfacial bond strengths are overly simple, and need to be improved.     

Effect of dilatancy on the strain localization of water-saturated elasto-viscoplastic soil

It is well known that geomaterials such as soils exhibit an increase in volume during shearing deformation, referred to as dilatancy. Dilatancy is a typical property of such granular materials as soils and is closely related to changes in the microstructure. Normally consolidated clay exhibits negative dilatancy or contractancy, namely, a decrease in volume during shearing. On the other hand, overconsolidated clay shows positive dilatancy, namely, an increase in volume during shearing. The aim of the present paper is to study the effects of the microstructure, such as dilatancy and permeability, on the strain localization of water-saturated clay using an elasto-viscoplastic constitutive model. Based on the non-linear kinematic hardening theory and a Chaboche type of viscoplasticity model, an elasto-viscoplastic model for both normally consolidated and overconsolidated clays is proposed; the model can address both negative and positive dilatancies. Firstly, the instability of the model under undrained creep conditions is analyzed in terms of the accelerating creep failure. The analysis shows that clay with positive dilatancy is more unstable than clay with negative dilatancy. Secondly, a finite element analysis of the deformation of water-saturated clay is presented with focus on the numerical results under plane strain conditions. From the present numerical analysis, it is found that both dilatancy and permeability prominently affect shear strain localization behavior.