[126.96.36.199] The approach presented here is based on the traditional theory presented in Section 1 combined with the distinction between percolating and nonpercolating phases introduced in Refs. [14, 16, 15] and discussed in Section 2. [188.8.131.52] The volume fractions of the subsets and are denoted as . [184.108.40.206] Let denote the porosity (volume fraction of ). [220.127.116.11] Volume conservation implies the relations
where are volume fractions, and are saturations. [18.104.22.168] The water saturation is defined as and the oil saturation as .
[22.214.171.124] As in the traditional theory the porous medium is again assumed to be macroscopically homogeneous so that eq. (4) holds. [126.96.36.199] The flows are slow, and hence also eq. (7) continues to hold without change. [188.8.131.52] For incompressible fluids one has now
analogous to eq. (5a).
[184.108.40.206] The momentum transfer into phase from all the other phases is assumed to be a simple viscous drag,
where the resistance coefficient quantifies the viscous coupling between phase and . [220.127.116.11] For the rigid rock matrix [page 4, §0] . [18.104.22.168] Hence is the momentum transfer from the wall into phase . [22.214.171.124] Then
where and was used because there is no common interface and hence no direct viscous interaction between these phase pairs. [126.96.36.199] Each is a -matrix.
[188.8.131.52] In the present approach the body forces are given by gravity plus capillary forces
contrary to eqs. (8). [184.108.40.206] The capillary body forces are responsible for keeping the trapped fluids inside the medium. [220.127.116.11] They are obtained as gradients of capillary potentials
where the capillary potentials are defined as
with constants and exponents .
[18.104.22.168] Next the stress tensor for percolating phases can be specified in analogy with eq. (6) as
where and are the fluid pressures. [22.214.171.124] The stress tensor for the nonpercolating phases cannot be specified in the same way because the forces cannot propagate in nonpercolating phases. [126.96.36.199] Here it is assumed that these stresses are given by the pressure in the surrounding percolating phase modified by the energy density stored in the common interface with the surrounding percolating phases. [188.8.131.52] This suggests an Ansatz 
where is the oil-water interfacial tension, and the unknowns are the interfacial areas per unit volume of porous medium between phases and , resp. and . [184.108.40.206] (The factor arises from the definition of per unit volume of porous medium.) [220.127.116.11] To avoid equations of motion for the unknowns and it is assumed that geometrical relations of the form
hold, where are prefactors assumed to be constant. [18.104.22.168] Thus for the nonpercolating phases
will be used below. [22.214.171.124] Here
[126.96.36.199] The mass transfer rates must depend on rates of saturation change. [188.8.131.52] They are here assumed to be
where are constants. [184.108.40.206] The parameters , are defined by
[page 5, §0] where are limiting saturations for and are equilibration time scales for reaching capillary equilibrium. [220.127.116.11] For simplicity
will be assumed below.