Computational Geomechanics. Manuel Pastor
Чтение книги онлайн.
Читать онлайн книгу Computational Geomechanics - Manuel Pastor страница 33
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,
(2.66)
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
with ps = pw Sw + pa Sa in the case of thermodynamic equilibrium or for incompressible solid grains (2.50).
The sum of (2.67), written for air and water and of (2.68), gives the total stress σ, acting on a unit area of the volume fraction mixture
This is the form of the effective stress (2.51), also called generalized Bishop stress (Nuth and Laloui 2008), employed in the following, as already explained.
Moist air in the pore system is assumed to be a perfect mixture of two ideal gases, dry air and water vapour, with π = ga and π = gw, respectively. The equation of a perfect gas is hence valid
where Mπ is the molar mass of constituent π, R the universal gas constant, and θ the common absolute temperature. Further, Dalton’s law applies and yields the molar mass of moisture
Water is usually present in the pores as a condensed liquid, separated from its vapor by a concave meniscus because of surface tension. The capillary pressure is defined as pc = pg − pw, see Equation (2.43).
The momentum exchange term of the linear momentum balance equation for fluids has the form
(2.72)
where vπs