[page 2323, §1]
[2323.1.1] A macroscopic theory of two phase flow inside a rigid porous medium poses not only challenges to nonequilibrium statistical physics and geometry , but is also crucial for many applications [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. [2323.1.2] Despite its popularity the accepted macroscopic theory of two phase flow seems unable to reproduce the experimentally observed phenomenon of saturation overshoot .
[2323.2.1] Models for twophase flow in porous media can be divided into macroscopic (laboratory or field scale) models popular in engineering, and microscopic (pore scale) models such as network models [12, 13, 14, 15, 16, 17] that are popular in physics. [2323.2.2] As of today no rigorous connection exists between microscopic and macroscopic models [1, 18]. [2323.2.3] In view of the predominantly non-specialist readership with a physics background it is appropriate to remind the reader of the traditional theory, introduced between 1907 and 1941 by Buckingham, Richards, Muskat, Meres, Wyckoff, Botset, Leverett and others [19, 20, 21, 22, 23]a (This is a footnote:) aThe following introductory paragraphs are quoted from Ref. for convenience of the interdisciplinary readership and following an explicit request from the editor.. [2323.2.4] One formulation of the traditional macroscopic theory starts from the fundamental balance laws of continuum mechanics for two fluids (called water and oil ) inside the pore space (called ) of a porous sample with a rigid [page 2324, §0] solid matrix (called ). [2324.0.1] Recall the law of mass balance in differential form
where denote mass density, volume fraction and velocity of phase as functions of position and time . [2324.0.2] Exchange of mass between the two phases is described by mass transfer rates giving the amount of mass by which phase changes per unit time and volume. [2324.0.3] Momentum balance for the two fluids requires in addition
where is the stress tensor in the th phase, is the body force per unit volume acting on the th phase, is the momentum transfer into phase from all the other phases, and denotes the material derivative for phase .
[2324.1.1] Defining the saturations as the volume fraction of pore space filled with phase one has the relation where is the porosity of the sample. [2324.1.2] Expressing volume conservation in terms of saturations yields
[2324.1.3] In order to obtain the traditional theory these balance laws for mass, momentum and volume have to be combined with specific constitutive assumptions for and .
[2324.2.1] Great simplification is afforded by assuming that the porous medium is rigid and macroscopically homogeneous
where are the fluid pressures. [2324.2.4] Realistic subsurface flows have low Reynolds numbers so that the inertial term
can be neglected in the momentum balance equation (2). [2324.2.5] It is further assumed that the body forces
are given by gravity. [2324.2.6] As long as there are no chemical reactions between the fluids the mass transfer rates vanish, so that holds. [2324.2.7] Momentum transfer between the fluids and the rigid matrix is dominated by viscous drag in the form
[page 2325, §0] where are the constant fluid viscosities, is the absolute permeability tensor, and are the nonlinear relative permeabilitiy functions for water and oil b (This is a footnote:) b account for the fact, that the experimentally observed permability of two immiscible fluids deviates from their partial (or mean field) permeabilities obtained from volume averaging of the absolute permeability. .
while the momentum balance eq. (2)
give the generalized Darcy laws for the Darcy velocities [3, p. 155]. [2325.1.2] Equations (9) and (10) together with eq. (3) provide 9 equations for 12 primary unknowns ,. Additional equations are needed.
where is the oil-water interfacial tension and is the mean curvature of the oil-water interface. [2325.2.2] The system of equations is closed with two equations of state relating the phase pressures and densities. [2325.2.3] In petroleum engineering the two fluids are usually assumed to be incompressible
while in hydrology one thinks of water and air setting
where J kgK is the specific gas constant and the temperature is assumed to be constant throughout .
[2325.3.1] When the fluids (water and oil) are incompressible (as in petroleum engineering) eqs. (12a) and (12b) hold. In this case, adding equations (9a) and (9b), using eq. (3) and integrating the result shows
[page 2326, §0] where (with )
are the mobilities , total mobility and fractional flow functions , respectively. [2326.0.1] Multiplying eq. (10a) with , eq. (10b) with and subtracting eq. (10b) from eq. (10a), using eq. (13) and to eliminate gives the result
which can be inserted into eq. (9a) to give
a nonlinear partial differential equation for the saturation field . [2326.0.2] For small or when gravity and capillarity effects can be neglected the last two terms vanish and eq. (17) reduces to the Buckley-Leverett equation 
a quasilinear hyperbolic partial differential equation. [2326.0.3] Equation (17) supplemented with a (quasilinear elliptic) equation obtained from by defining a global pressure in such a way that the total flux obeys a Darcy law with respect to the global pressure provides, for incompressible fluids, an equivalent formulation of eqs. (9)-(12b). [2326.0.4] For compressible fluids the situation is different.
[2326.1.1] When corresponds to water and to air (as for applications in hydrology) eqs. (12c) and (12d) hold. [2326.1.2] The large density difference suggests to consider the case of a very rarified gas or vacuum as a first approximation. [2326.1.3] For eq. (9b) is identically fulfilled, eq. (12d) implies and then eq. (10b) implies . [2326.1.4] In this way the -phase vanishes from the problem and one is left only with the -phase. [2326.1.5] Inserting eq. (10a) into eq. (9a) and using eq. (11) gives the Richards equation 
for saturation or
for pressure after writing with the help of eq. (11). [2326.1.6] To define the nonlinear function the typical sigmoidal shape has been assumed for .
[2326.2.1] The quasilinear elliptic-parabolic Richards equation (19) is the basic equation in hydrology, while the quasilinear hyperbolic Buckley-Leverett equation (18) is fundamental for applications in petroleum engineering. [2326.2.2] Both equations, (18) and (19), differ from the general fractional flow formulation (17) in terms of saturation and global pressure . [2326.2.3] They differ also from the formulation in terms of given by (3),(9),(10) and (11). [2326.2.4] These latter equations appropriately supplemented with initial and boundary conditions and spaces of functions resp. generalized functions for the unknowns constitute the traditional theory of macroscopic capillarity in porous media.
[page 2327, §1] [2327.1.1] The question of domains is important for wellposedness and numerical solution. [2327.1.2] For eq. (18) it is well known that classical solutions, i.e. locally Lipschitz continuous functions, will in general exist only for a finite length of time [28, 29, 30]. [2327.1.3] Hence it is necessary to consider also weak solutions . [2327.1.4] Weak solutions are locally bounded, measurable functions satisfying eq. (18) in the sense of distributions. [2327.1.5] Weak solutions are frequently constructed by the method of vanishing viscosity or the theory of contraction semigroups. [2327.1.6] For the Richards equation (19) a domain of definition in the space of Bochner-square-integrable Sobolev-space-valued functions has been discussed in . [2327.1.7] In many engineering applications formulations such as eqs. (17), (18) or (19) with (11) augmented with appropriate initial and boundary conditions are solved by computer programs [33, 34, 35]. [2327.1.8] This concludes our brief introduction into the traditional theory.