Computational Geomechanics. Manuel Pastor

Читать онлайн книгу.

Computational Geomechanics - Manuel Pastor


Скачать книгу
mass, linear momentum, and angular momentum are then obtained by systematically applying the averaging procedures to the microscopic balance Equation (2.60) as outlined in Hassanizadeh and Gray (1979a, 1979b, 1980). The balance equations have here been specialized for a deforming porous material, where the flow of water and of moist air (mixture of dryair and vapor) is taking place (see Schrefler 1995).

      The local thermodynamic equilibrium hypothesis is assumed to hold because the time scale of the modeled phenomena is substantially larger than the relaxation time required reaching equilibrium locally. The temperatures of each constituent in a generic point are hence equal. Further, the constituents are assumed to be immiscible and chemically nonreacting. All fluids are assumed to be in contact with the solid phase. As throughout this book, stress is defined as tension positive for the solid phase, while pore pressure is defined as compressively positive for the fluids.

      In the averaging procedure, the volume fractions ηπ appear which are identified as follows: for solid phase ηs = 1 − n, for water ηw = nSw, and for air ηa = nSa.

      The averaged macroscopic mass balance equations are given next. For the solid phase, this equation reads

      where ModifyingAbove bold-italic u With ampersand c period dotab semicolon is the mass averaged solid phase velocity and ρπ is the intrinsic phase averaged density. The intrinsic phase averaged density ρπ is the density of the π phase averaged over the part of the control volume (Representative Elementary Volume, REV) occupied by the π phase. The phase averaged density ρπ, on the contrary, is the density of the π phase averaged over the total control volume. The relationship between the two densities is given by

      (2.62)rho Subscript pi Baseline equals eta Subscript pi Baseline rho Superscript pi

      For water, the averaged macroscopic mass balance equation reads

      For air, this equation reads

      (2.64)StartFraction upper D Overscript a Endscripts left-parenthesis italic n upper S Subscript a Baseline rho Superscript a Baseline right-parenthesis Over italic upper D t EndFraction plus italic n upper S Subscript a Baseline rho Superscript a Baseline v Subscript i slash i Superscript a Baseline equals ModifyingAbove m With ampersand c period dotab semicolon

      where va is the mass averaged air velocity.

      The linear momentum balance equation for the fluid phases is

      where tπ is the partial stress tensor, ρπbπ the external momentum supply due to gravity, ρπaπ the volume density of the inertia force, rho Subscript pi Baseline e Superscript pi Baseline left-parenthesis rho ModifyingAbove bold r With ampersand c period dotab semicolon right-parenthesis the sum of the momentum supply due to averaged mass supply, and the intrinsic momentum supply due to a change of density and referred to the deviation ModifyingAbove bold r With ampersand c period dotab semicolon Superscript pi of the velocity of constituent π from its mass averaged velocity, and ModifyingAbove bold t With ampersand c period circ semicolon Superscript pi accounts for exchange of momentum due to mechanical interaction with other phases. rho Subscript pi Baseline e Superscript pi Baseline left-parenthesis rho ModifyingAbove bold r With ampersand c period dotab semicolon right-parenthesis is assumed to be different from zero only for fluid phases. For the solid phase, the linear momentum balance equation is hence

      (2.66)t Subscript italic j i slash j Superscript s Baseline plus rho Superscript s Baseline left-parenthesis b Subscript i Superscript s Baseline minus a Subscript i Superscript s Baseline right-parenthesis plus rho Superscript s Baseline ModifyingAbove t With ampersand c period circ semicolon Subscript i Superscript s Baseline equals 0

      The average angular momentum balance equation shows that for nonpolar media, the partial stress tensor is symmetric tπ ji = tπ tj at the macroscopic level also and the sum of the coupling vectors of angular momentum between the phases vanishes.

      2.5.4 Constitutive Equations

      Constitutive models are selected here which are based on quantities currently measurable in laboratory or field experiments and which have been extensively validated. Most of them have been obtained from entropy inequality; see Hassanizadeh and Gray (1980, 1990).

      It can be shown that the stress tensor in the fluid is

      where pπ is the fluid pressure, and in the solid phase is