Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 27 February 2023
Sec. Interdisciplinary Physics

Parameterizations of immiscible two-phase flow in porous media

  • PoreLab, Department of Physics, Norwegian University of Science and Technology, Trondheim, Norway

A fundamental variable characterizing immiscible two-phase flow in porous media is the wetting saturation, which is the ratio between the pore volume filled with wetting fluid and the total pore volume. More generally, this variable comes from a specific choice of coordinates on some underlying space, the domain of variables that can be used to express the volumetric flow rate. The underlying mathematical structure allows for the introduction of other variables containing the same information, but which are more convenient from a theoretical point of view. We introduce along these lines polar coordinates on this underlying space, where the angle plays a role similar to the wetting saturation. We derive relations between these new variables based on the Euler homogeneity theorem. We formulate these relations in a coordinate-free fashion using differential forms. Finally, we discuss and interpret the co-moving velocity in terms of this coordinate-free representation.

1 Introduction

Flow of immiscible fluids in porous media [14] is a problem that has been in the hands of engineers for a long time. This has resulted in a schism between the physics at the pore scale and the description of flow at scales where the porous medium may be seen as a continuum, also known as the Darcy scale. This fact has led to the phenomenological relative permeability equations proposed by Wyckoff and Botset in 1936 [5] with the inclusion by Leverett of the concept of capillary pressure in 1940 [6], still being the unchallenged approach to calculate flow in porous media at large scales. The basic ideas of the theory are easy to grasp. Seen from the viewpoint of one of the immiscible fluids, the solid matrix and the other fluid together reduce the pore space in which that fluid can move. Hence, the effective permeability as seen by the fluid is reduced, and the permeability reduction factor of each fluid is their relative permeability. The capillary pressure models the interfacial tension at the interfaces between the immiscible fluids by assuming that there is a pressure difference between the pressure fields in each fluid. The central variables of the theory, relative permeability, and capillary pressure curves are determined routinely in the laboratory under the name special core analysis and then used as input in reservoir simulators [3, 7]. A key assumption in the theory is that the relative permeability and the capillary pressure are functions of the saturation alone. This simplifies the theory tremendously, but it also distances the theory from realism.

Some progress has been made in order to improve on relative permeability theory. Barenblatt et al. [8] recognized that a key assumption in relative permeability theory is that locally, the fluids are in phase equilibrium, even if the flow as a whole is developing. This has as an implication that the central variables of that theory, relative permeability, and the capillary pressure are functions of the saturation alone. They then go on to generalize the theory to flow which is locally out of equilibrium, exploring how the central variables change. Wang et al. [9] went further by introducing dynamic length scales due to the mixing zone variations, over which the spatial averaging is performed.

The relative permeability approach is phenomenological and going beyond relative permeability theory means making a connection between the physics at the pore scale with the Darcy scale description of the flow. The attempts at constructing a connection between these two scales—which may stretch from micrometers to kilometers—have mainly been focused on homogenization, replacing the original porous medium by an equivalent spatially structureless one.

The most famous approach to the scale-up problem along these lines is thermodynamically constrained averaging theory (TCAT) [1014], based on thermodynamically consistent definitions made at the continuum scale based on volume averages of pore-scale thermodynamic quantities, combined with closure relations based on homogenization [15]. All variables in TCAT are defined in terms of pore-scale variables. However, these result in many variables and complicated assumptions are needed to derive useful results.

Another homogenization-approach based on non-equilibrium thermodynamics uses Euler homogeneity to define the up-scaled pressure. From this, Kjelstrup et al. derived constitutive equations for the flow while keeping the number of variables down [1618].

There is also an ongoing effort in constructing a scaled-up theory based on geometrical properties of pore space by using the Hadwiger theorem [1921]. The theorem implies that we can express the properties of the spatial geometry of the three-dimensional porous medium as a linear combination of four Minkowski functionals: volume, surface area, mean curvature, and Gaussian curvature. These are the only four numbers required to characterize the geometric state of the porous medium. The connectivity of the fluids is described by the Euler characteristic, which by the Gauss–Bonnet theorem can be computed from the total curvature.

A different class of theories is based on detailed and specific assumptions concerning the physics involved. An example is local porosity theory [2227]. Another example is the decomposition in prototype flow (DeProf) theory which is a fluid mechanical model combined with non-equilibrium statistical mechanics based on a classification scheme of fluid configurations at the pore level [2830]. A third example is that of Xu and Louge [31] who introduced a simple model based on Ising-like statistical mechanics to mimic the motion of the immiscible fluids at the pore scale, and then scale up by calculating the mean field behavior of the model. In this way, they obtain the wetting fluid retention curve.

The approach we take in this paper is built on the approach to the upscaling problem found in [3236]. The underlying idea is to map the flow of immiscible fluids in porous media, a dissipative and therefore out-of-equilibrium system, onto a system which is in equilibrium. The Jaynes principle of maximum entropy [37] may then be invoked and the scale-up problem is transformed into that of calculating a partition function [35].

In order to sketch this approach, we need to establish the concept of steady-state flow. In the field of porous media, “steady-state flow” has two meanings. The traditional one is to define it as flow where all fluid interfaces remain static [38]. The other one, which we adopt here, is to state that it is flow where the macroscopic variables remain fixed or fluctuate around well-defined averages. The fluid interfaces will move and on larger scales than the pore scale where one will see fluid clusters breaking up and merging. If the porous medium is statistically homogeneous, we will see the same local statistics describing the fluids everywhere in the porous medium [36].

Imagine a porous plug as shown in Figure 1. There is a mixture of two immiscible fluids flowing through it in the direction of the cylinder axis under steady-state conditions.

FIGURE 1
www.frontiersin.org

FIGURE 1. A porous plug in the form of a cylinder. There is a mixture of two immiscible fluids flowing through it in the direction of the cylinder axis under steady-state conditions. Two disks orthogonal to the cylinder axis are further imagined. These are REAs.

A central concept in the following is the representative elementary area (REA) [34, 35, 39]. We pick a point in the porous plug. In a neighborhood of this point, there will a set of streamlines associated with the velocity field v. The overall flow direction is then defined by the tangent vectors to the streamlines. We place an imaginary disk of area A at the point, with orientation such that A is orthogonal to the overall flow direction. We illustrate this as the lower disk in Figure 1. We assume that the disk is small enough for the porous medium to be homogeneous over the size of the plane with respect to porosity and permeability. This disk is an REA.

The REA will contain different fluid configurations. By “fluid configuration”, we mean the spatial distributions of the two fluids in the disk and their scalar velocity fields. We also show a second disk in Figure 1, which represents another REA placed further along the average flow. Since the flow is incompressible, we can imagine the z-axis corresponding to a pseudo-time axis. We may state that the fluid configuration in the lower disk of Figure 1 evolves into the fluid configuration in the upper disk under pseudo-time-translation. We now associate entropy to the fluid configurations in the sense of Shannon [40]. Even though the system is producing molecular entropy through dissipation, it is not producing Shannon entropy along the z-axis. This allows us to use the Jaynes maximum entropy principle, and a statistical mechanics based on Shannon entropy ensues [35].

This statistical mechanics scales up the pore-scale physics, represented by the fluid configurations to the Darcy scale, which then is represented by a thermodynamics-like formalism involving averaged velocities.

Rudiments of this thermodynamics-like formalism was first studied by Hansen et al. [32], who used extensivity to derive a set of equations that relate the seepage velocity of each of the two immiscible fluids flowing under steady state conditions through a representative volume element (REV). They introduced a co-moving velocity together with the average seepage velocity v that contained the necessary information to determine the seepage velocities of the more wetting and the less wetting fluids, vw and vn respectively. Stated in a different way, these equations made it possible to make the transformation

v,vmvw,vn.(1)

It is the aim of this paper to present a geometrical interpretation of the transformation in Eq. 1. We wish to provide an intuitive understanding of precisely what the transformation means and to determine the role of the co-moving velocity. Our aim is not to develop the theory presented in [32] further by including new results, but rather consolidate and amend the existing theory with a deeper understanding.

Our geometrical interpretation will rest upon the definition of a space spanned by two central variables: the wetting and non-wetting transversal pore areas. These areas are defined as the area of REAs covered by the wetting and non-wetting fluids, respectively.

Instead of using the two extensive areas directly, one often works with coordinates where one is the (wetting) saturation. If the pore area of the REA is fixed, it is often convenient to let the other variable be this area. A third way—which is new—is to use polar coordinates. This, as we shall see, simplify the theory considerably.

In Section 2, we describe the REA and the relevant variables. We define intensive and extensive variables, tying them to how they scale under scaling of the size of the REV. We also review the central results of Hansen et al. [32] here: the introduction of the auxiliary thermodynamic velocities and the co-moving velocity.

The central concept in our approach to the immiscible two-fluid flow problem is that of the pore areas. In Section 3 we introduce area space, a space that simplifies the analysis considerably. We showcase different coordinate systems on this space, focusing on polar coordinates.

Section 4 expresses the central results of Hansen et al. [32] in polar coordinates. Expressed in this coordinate system, the geometrical structure implied by the extensivity versus intensivity conditions imposes restrictions that simplify the equations. The section goes on to express the total pore velocity, seepage velocities, thermodynamic velocities, and co-moving velocity in terms of polar coordinates. This provides a clear geometrical interpretation of the co-moving velocity.

Section 5 focuses on relations between the different velocities that can be expressed without referring to any coordinate system. We do this by invoking differential forms and exterior algebra [41, 42]. Differential forms constitute the generalization of the infinitesimal line, area, or volume elements used in integrals and exterior calculus is the algebra that makes it possible to handle them. Flanders predicted in the early sixties [41] that within short they would become important tools in engineering. This did not happen, and they have remained primarily within mathematics and theoretical physics. Such objects are also especially prevalent in thermodynamics, and here, we will show that the computational rules of differential forms provide a simple way of expressing the relations between the velocities in the system at hand.

Section 6 focuses the discussion on the co-moving velocity by constructing a coordinate-independent expression that defines it.

We present in Section 7 two concrete systems. 1) A porous medium obeying the Brooks–Corey relative permeability and 2) a capillary fiber bundle model, which we analyze using the methods developed in this paper.

In Section 8, we summarize our main results which may be stated through three Eqs 58, 69, 98. These equations formulate in a coordinate-free way the three relations that exist between the average seepage velocity and the co-moving velocity.

2 Representative elementary area

We will now elaborate on the definition of the REA presented in Section 1.

We define the transversal pore area Ap as previously defined, which defines porosity

ϕ=ApA.(2)

The pore volume contains two immiscible fluids; the more wetting fluid (to be referred to as the wetting fluid) or the less wetting fluid (to be referred to as the non-wetting fluid). The transversal pore area Ap may, therefore, be split into the area of the REA disk cutting through the wetting fluid Aw or cutting through the non-wetting fluid An so that we have

Ap=Aw+An.(3)

We also define the saturations

Sw=AwAp,(4)
Sn=AnAp,(5)

so that

Sw+Sn=1.(6)

We note that the transversal pore areas Ap, Aw, and An are extensive in the REA A, that is ApλAp, AwλAw, and AnλAn when AλA. The porosity ϕ and the two saturations Sw and Sn are intensive in the REA A: ϕλ0ϕ, Swλ0Sw and Snλ0Sn when AλA.

There is a time averaged volumetric flow rate Q through the REA. The volumetric flow rate consists of two components, Qw and Qn, which are the volumetric flow rates of the wetting and the non-wetting fluids. We have

Q=Qw+Qn.(7)

We define the total, wetting, and non-wetting seepage velocities, respectively, as

v=QAp,(8)
vw=QwAw,(9)
vn=QnAn.(10)

Using Eqs 710, we find

v=QAp=AwApQwAw+AnApQnAn=Swvw+Snvn.(11)

The volumetric flow rates Q, Qw, and Qn are extensive in the REA A, and the velocities v, vw, and vn are intensive in the area A.

2.1 Euler homogeneity

As the areas Aw and An, and the volumetric flow rate Q are extensive in the REA A, we have that

QλAw,λAn=λQAw,An.(12)

This scaling was the basis for the theory presented in [32], which we now summarize.

We now apply the Euler homogeneous function theorem to Q. We take the derivative with respect to λ on both sides of Eq. 12 and set λ = 1. This gives

QAw,An=AwQAwAn+AnQAnAw.(13)

We note here that we assume Aw and An are our independent control variables. This makes Ap, Sw, and Sn dependent variables. More precisely, Eq. 13 is the Euler theorem for homogeneous functions applied to a degree-1 homogeneous function Q. By dividing this equation by Ap, we have

v=SwQAwAn+SnQAnAw,(14)

where we have used Eqs 4, 5. We define the two thermodynamic velocities v̂w and v̂n as

v̂w=QAwAn(15)

and

v̂n=QAnAw(16)

so that we may write Eq. 14 as

v=Swv̂w+Snv̂n.(17)

The thermodynamic velocities v̂w and v̂n are not the physical velocities vw and vn. Rather, the most general relation between vw and v̂w and vn and v̂n, which fulfills both Eqs 11, 17,

v=Swv̂w+Snv̂n=Swvw+Snvn,(18)

can be expressed as

v̂w=vw+Snvm,(19)
v̂n=vnSwvm.(20)

This defines the co-moving velocity vm, which relates the thermodynamic and the physical velocities.

We have up to now used (Aw, An) as our control variables. If we now consider the coordinates (Sw, Ap) instead, we can write

v̂w=QAwAn=QSwApSwAwAn+QApSw=dvdSwSn+v,(21)

where we have used that

SwAwAn=AwAnAwAw+An=SnAp.(22)

Likewise, we find that

v̂n=QAnAw=dvdSwSw+v.(23)

We combine Eqs 21, 23 with Eqs 19, 20 to find

vw=v+SndvdSwvm,(24)
vn=vSwdvdSwvm,(25)

and we see that these two equations give the map (v, vm) → (vw, vn).

We now subtract Eq. 25 from Eq. 24, finding

dvdSw=vwvn+vm,(26)

and we now differentiate Eq. 11 with respect to Sw,

dvdSw=vwvn+SwdvwdSw+SndvndSw.(27)

We compare Eqs 26, 27, finding

vm=SwdvwdSw+SndvndSw.(28)

Eqs 11, 28 constitute the inverse mapping (vw, vn) → (v, vm).

The mapping (v, vm) → (vw, vn) tells us that given the constitutive equations for v and vm, we also have the constitutive equations for vw and vn. It turns out that the constitutive equation for vm is surprisingly simple [34],

vm=a+bdvdSw,(29)

where a and b are coefficients depending on the entropy associated with the fluid configurations [35].

The equations we have used in Section 2 hint at an underlying mathematical structure which may seem complex. As we will now show, it is in fact quite the opposite.

3 Coordinate systems in area space

The transversal pore areas Aw ≥ 0 and An ≥ 0 parametrize the first quadrant of R2. This serves as the “area space” mentioned earlier, which we will continue to denote as such. It is important to realize that this space is not physical space. A value of the transversal pore area Ap corresponds to a point Aw,An in this space. If we consider the entire plane R2 as a whole and simply restrict our attention to the first quadrant, we may treat the area space as a vector space. The point Aw,An may then equivalently be described by a vector,

A=Awew+Anen,(30)

where ew and en form an orthonormal basis set, shown in the leftmost figure in Figure 2. Note that in this picture, the bases are shown as attached to the point Aw,An. We reiterate that both Aw and An are extensive.

FIGURE 2
www.frontiersin.org

FIGURE 2. Three coordinate systems are illustrated and used to parametrize the transversal pore area. (A) The Cartesian coordinate system (Aw, An). A curve of constant Aw and An is indicated as Aw and An, respectively. Both Aw and An are extensive. The orthonormal basis set (ew,en) is shown as attached to the point Aw,An. (B) The saturation coordinate system Ap,Sw. Curves of constant Ap and Sw are denoted by Ap and Sw in the figure. Ap is extensive and Sw is intensive. The non-orthonormal basis (ep,es) is shown. (C) The polar coordinate system (Ar,φ). Curves of constant Ar and φ are denoted by Ar and φ in the figure. Ar is extensive and ϕ is intensive. The basis set (er,eϕ) is orthonormal. In all three figures, the transversal pore area as a vector A is shown.

Every point in the transversal area space corresponds to a given saturation Sw and a transversal pore area Ap and we may view the map (Aw, An) → (Ap, Sw) given by Eqs 3, 4 as a coordinate transformation. This is a more natural coordinate system to work with since Sw is an intensive variable and Ap is in practice kept constant. We name this system the saturation coordinate system. We show the normalized basis vector set in the middle figure in Figure 2. We calculate

up=AApSw=Swew+1Swen,(31)
us=ASwAp=ApewApen.(32)

We normalize the two vectors up and us to find

ep=Swew+1Swen12Sw1Sw1/2,(33)
es=ewen2.(34)

This basis set is not orthonormal, as illustrated in Figure 2. By expressing A in this coordinate system, we find

A=Ap12Sw1Sw1/2ep.(35)

As it will become apparent, the most convenient coordinate system to work with from a theoretical point of view is polar coordinates Ar,ϕ, given by

Ar=Aw2+An2,(36)
ϕ=arctanAnAw,(37)

or vice versa.

Aw=Arcosϕ,(38)
An=Arsinϕ.(39)

We note that ϕ is intensive and Ar is extensive. We show the polar coordinate system in the lower figure in Figure 2.

The basis vector set (er,eϕ) is orthonormal, where

er=cosϕew+sinϕen=ep,(40)
eϕ=sinϕew+cosϕen.(41)

As for the saturation coordinate system, there is one intensive variable, ϕ, and one extensive variable, Ar. However, in contrast to the saturation coordinate system, both variables are varied in practical situations.

We have that

A=Arer.(42)

We see that this is consistent with Eq. 35 since

Ar=Ap12Sw1Sw1/2(43)

and ep=er.

4 Euler homogeneity in polar coordinates

We now focus on Eq. 12 which we repeat here.

QλAw,λAn=λQAw,An.

We may interpret this equation geometrically. If we follow the value of Q along a ray passing through the origin of the two-dimensional space spanned by (Aw, An) keeping the ratio An/Aw constant, it grows linearly with the distance from the Q-axis.

In polar coordinates, this means

QAr,ϕ=v̂rAr,(44)

where

v̂r=QArϕ=v̂rϕ.(45)

The important point here is that v̂r is not a function of Ar.

We may derive Eq. 44 from Eq. 12 as follows. First, we write Aw and An in terms of Ar and ϕ in Eq. 12 so that

1λQλArcosϕ,λArsinϕ=QArcosϕ,Arsinϕ=QAr,ϕ.(46)

Next, we set λ = 1/Ar to find

ArQcosϕ,sinϕ=Arv̂rϕ=QAr,ϕ,(47)

where

v̂rϕQcosϕ,sinϕ,(48)

and we are done. We illustrate in Figure 3 the geometrical meaning of Eq. 12, that is, Euler homogeneity. The graph of the volumetric flow rate Q takes the shape of a “crumpled” cone in the space spanned by (Q, Aw, An).

FIGURE 3
www.frontiersin.org

FIGURE 3. Geometrical meaning of the function Q being degree-1 Euler homogeneous in the extensive variables, meaning it fulfils equation (12), using polar coordinates. We here use the same symbol Q as both a variable and a function. Moving along a ray in the space spanned by the variables Q, Aw and An, the volumetric flow rate increases linearly with the distance from the Q-axis. For a constant value Ar of the radial coordinate Ar (A) or a constant value Q of the volumetric flow rate (B), the graph of Q appears as a “crumpled” cone originating at the origin, with the level set (shown in red) imagined as an “edge”. The level sets are different in the two cases. In particular, the volumetric flow rate will in (A) have different limiting values Q1 and Q2 in the limit of single-phase flow of either the wetting or non-wetting fluid.

We now express the different fluid velocities, namely the pore velocity, seepage velocities, thermodynamic velocities, and the co-moving velocity, in terms of polar coordinates.

4.1 Seepage velocities

The area space spanned by Aw,An is simply (a subset of) the real plane. We can therefore identify it with the tangent space of velocities, which is the same Cartesian plane as that of the areas Aw,An. We can therefore regard the seepage velocity as a vector in the area space,

v=vwew+vnen.(49)

The volumetric flow rate Q is then given by the scalar product

Q=vA=vwAw+vnAn=v̂wAw+v̂nAn.(50)

In polar coordinates, the seepage velocity is given by

v=vrer+vϕeϕ.(51)

Hence, the volumetric flow rate in polar coordinates is given by

Q=vA=vrer+vϕeϕArer=vrAr.(52)

Comparing with Eq. 44, we find

vr=v̂r=v̂rϕ.(53)

By using Eqs 49, 51 combined with Eqs 40, 41, we find

vr=vwcosϕ+vnsinϕ=Swvw+1Swvn12Sw1Sw1/2,(54)
vϕ=vwsinϕ+vncosϕ=1Swvw+Swvn12Sw1Sw1/2.(55)

We have here expressed cos ϕ and sin ϕ in terms of Sw.

We now compare Eq. 44 with Eq. 8, rewritten as

Q=vAp.(56)

Using Eq. 43, we find that their equality demands

v=v̂r12Sw1Sw1/2.(57)

Hence, v is not the norm of v.

4.2 Thermodynamic velocities

We may write Eq. 13 as

Q=v̂wAw+v̂nAn=v̂A,(58)

where we have defined, in the same manner as for Eq. 49,

v̂=v̂wew+v̂nen.(59)

We use Eqs 40, 41 to express this equation in polar coordinates,

v̂=v̂rer+v̂ϕeϕ,(60)

where

v̂r=v̂wcosϕ+v̂nsinϕ,,(61)
v̂ϕ=v̂wsinϕ+v̂ncosϕ.(62)

4.3 Co-moving velocity

We defined the co-moving velocity vm as the velocity function that relates the physical seepage velocities and the thermodynamic velocities shown in Eqs 19, 20. We use these two equations together with Eqs 49, 59 to form the vector

v̂v=v̂wvwew+v̂nvnen=vmApAnewAwen=vmApArsinϕewArcosϕen,(63)

where we have used Sw = Aw/Ap and Eqs 38, 39. We rewrite this equation in terms of eϕ, Eq. 41, finding

v̂v=vmcosϕ+sinϕeϕ.(64)

We define

v̂m=vmcosϕ+sinϕ(65)

so that Eq. 63 may be written as

v=v̂+v̂m,(66)

where

v̂m=v̂meϕ.(67)

In polar coordinates, the relation in Eq. 66 is isolated to the ϕ-coordinate. Using Eqs 51, 62, we can express the relation between the vector components v̂ϕ and vϕ as

vϕ=v̂ϕ+v̂m.(68)

We illustrate this relation in Figure 4. This figure also demonstrates why Q=v̂A=vA; it is due to (vv̂)=v̂mA. Hence, we have that

v̂mA=0.(69)

FIGURE 4
www.frontiersin.org

FIGURE 4. Vectorial representation for eqs. 52, 58, where Ar=|A| is illustrated. The difference in the vectors v̂ and v, v̂m, is orthogonal to A, and cannot be detected by the properties of A alone.

5 Coordinate-free representation

Since the difference between the thermodynamic- and seepage velocities in the previous section was shown to sit in the tangent vector components, we can benefit from examining the problem from a position where these components are easier to work with. In this section, we therefore introduce differential forms and exterior algebra [41, 42] to make this treatment more economic. In addition, this makes it easier to formulate the equations met so far in a way that does not depend on the coordinate system used.

Eq. 44, Q=v̂rAr, may be differentiated to give

dQ=v̂rdAr+dv̂rAr=v̂rdAr+v̂rArdϕ,(70)

where v̂r=dv̂r/dϕ. Here, dQ, dAr, and are one-forms. Furthermore, we have that

dQ=QArϕdAr+QϕArdϕ=v̂rdAr+v̂ϕArdϕ.(71)

Comparing Eqs 70, 71 gives

v̂ϕ=v̂r.(72)

We do the same using the coordinate system (Aw, An). From Eq. 50, we find

dQ=v̂wdAw+v̂ndAn+Awdv̂w+Andv̂n=v̂wdAw+v̂ndAn.(73)

In the second equality, we assumed that dQ is an exact differential, and in the first equality, we applied the differential to Q=v̂wAw+v̂nAn. This implies the relation

Awdv̂w+Andv̂n=0.(74)

This equation corresponds to the Gibbs–Duhem equation in thermodynamics, which is a statement of the dependency amongst the intensive variables.

We express v̂w and v̂n in terms of v̂r and v̂ϕ by inverting Eqs 61, 62 to find

v̂w=v̂rcosϕv̂ϕsinϕ,,(75)
v̂n=v̂rsinϕ+v̂ϕcosϕ.(76)

By applying the exterior derivative, we get the relations

dv̂w=cosϕdv̂rsinϕdv̂ϕv̂ndϕ,(77)
dv̂n=sinϕdv̂r+cosϕdv̂ϕ+v̂wdϕ.(78)

Combining Eqs 7578 with Gibbs–Duhem Eq. 74 gives

dv̂rdϕ=v̂ϕ,(79)

which is identical to Eq. 72. Hence, this is the Gibbs–Duhem equation in polar coordinates.

The thermodynamic velocity, Eq. 59, is, therefore, given by

v̂=v̂rer+v̂reϕ.(80)

We now have one of the main results of Hansen et al. [32] in a very compact form.

v=v̂rer+v̂r+v̂meϕ.(81)

We note that since v̂r only depends on ϕ, the same must be true for v̂m. Hence, v̂m=v̂m(ϕ). Here, we see that the relation between the thermodynamic- and seepage velocities are not captured in the extensive structure tied to the radial coordinate, but instead in the intensive quantities which only has a dependency on ϕ.

We now differentiate Eq. 56, Q=vA written in the saturation coordinate system, to find

dQ=vdAp+dvdSwApdSw.(82)

We note that

dAp=dAw+dAn,(83)

from differentiating Eq. 3. Likewise, we note that

dSw=SwAwAndAw+SwAnAwdAn=1Ap2AndAwAwdAn,(84)

using Eq. 4. We combine Eqs 8284 and find

dQ=v+1SwdvdSwdAw+vSwdvdSwdAn.(85)

By using Eq. 13, definitions of Eqs 15, 16, and transformations Eqs 19, 20, we may write

dQ=v̂wdAw+v̂ndAn=vw1SwvmdAw+vn+SwvmdAn.(86)

Equating Eq. 85 and Eq. 86 gives Eqs 24, 25.

Let us now rewrite the differential dQ as follows:

dQ=dApv=dApSwv̂w+1Swv̂n=Apv̂wv̂ndSw+vdAp,(87)

where we have used Eq. 74. We now combine this equation with Eq. 82 to get

dvdSw=v̂wv̂n.(88)

By using Eqs 19, 20, we may rewrite this equation in terms of the seepage velocities, resulting in Eq. 26.

We differentiate Q a second time, which is zero according to the Poincaré lemma [41, 42]. We see that this is indeed so by using Eq. 70,

d2Q=dv̂rdAr+v̂rArdϕ=v̂rdϕdAr+v̂rdArdϕ=v̂rv̂rdϕdAr=0,(89)

where the wedge ∧ signifies the antisymmetric exterior product.

The same equation expressed in the coordinate system (Aw, An) gives

d2Q=dv̂wdAw+v̂ndAn=v̂nAwAnv̂wAnAndAwdAn=0.(90)

We rewrite the coefficient in terms of the saturation coordinate system,

v̂nAwAnv̂wAnAn=dv̂ndSwSwAwAndv̂wdSwSwAnAw=Swdv̂wdSw+1Swdv̂ndSw=0,(91)

which is nothing but the Gibbs–Duhem equation yet again. Writing this equation in terms of the seepage velocities gives us Eq. 28.

We see that the different equations relating the velocities can all be found from the differential geometric structure of the space spanned by the areas, with addition to Euler homogeneity, see Figure 3. All relations are either a consequence of writing dQ using different coordinate systems or a consequence of the Poincaré lemma, expressed as d2Q = 0.

6 The co-moving velocity in coordinate-free representation

The co-moving velocity appears in several equations in the previous sections. Common to all of them is that they all are written out in terms of a given coordinate system. We may write the representation of the co-moving velocity in the saturation coordinate system, Eq. 28, as

Apvm=AdvdSw.(92)

We may write Eq. 65 in polar coordinates as

Arv̂m=Advdϕ.(93)

Expressed in terms of the seepage velocities, we have that

dQ=vwdAw+vndAn+Awdvw+Andvn.(94)

In order to obtain an expression for the co-moving velocity which is independent of the coordinate system, we construct

Awdvw+Andvn=vmAndAwAwdAnAw+An=vmApdSw=v̂mArdϕ,(95)

where the second line only reflects Eqs 92, 93.

We write dQ as

dQ=d[Av̂]=v̂idAi,(96)

using Eq. 74. Here, Einstein summation convention has been applied, with the index i running over w,n in the basis Aw,An. In terms of the seepage velocities, this becomes

dQ=d[Av]=vidAi+Aidvi.(97)

Combining Eqs 96, 97, 66 gives us

v̂midAi+Aidvi=0,(98)

where v̂mi denotes the components of v̂m. This holds in any coordinate system, since it is just obtained by differentiating a function, namely Av. We have, therefore, obtained an expression for the co-moving velocity which is independent of the chosen coordinate system.

We can introduce vector-valued differential forms [42] to make a connection between the interpretation in terms of (tangent) vectors used up until now and the language of differential forms. Let ei be an arbitrary basis for vectors on the area space as we have discussed until now, be it area-, saturation-, or polar coordinates. We now take A as an example. We write

A=βiei,(99)

where we have applied the Einstein summation convention. A is now regarded as a vector-valued differential form (thus seeing ei as components written out in the βi basis, i.e., a reversal of viewpoint), the βi are 0-forms, just functions, and the index i runs over the coordinate indices. We can still treat A as a vector, but the difference in viewpoints is apparent when we apply the exterior derivative d to A, which gives

dA=ieidβi=dβiei,(100)

where we have suppressed the tensor product of the forms and basis vectors since we are working over a vector space, as is standard notation in for example [42]. Strictly speaking, the exterior derivative d should be replaced by the exterior covariant derivative in this setting, a generalization of d which is defined on both tangent vectors and forms (there is currently no term in Eq. 100 that reflect the changes in the basis ei). The terms Aidvi are then possible to interpret in terms of connection forms. However, this is outside the scope of the current discussion, and we leave this topic for future work.

If we write out Eq. 100 in the coordinate system Aw,An, we get

dA=dAwew+dAnen.(101)

With this formalism in mind, we now return to polar coordinates. Eq. 98 may be written in these coordinates by using Eqs 40, 41 to define the vector valued 0-forms,

ωr=cosϕew+sinϕen,,(102)
ωϕ=sinϕew+cosϕen,(103)

which are equivalent to the corresponding basis vectors in Eqs 40, 41, but interpreted as 0-forms. We then compute

dωr=eϕdϕ,,(104)
dωϕ=ωrdϕ(105)

so that

dA=dArωr=dArωr+Ardωr=dArer+Areϕdϕ,(106)

where in the second equality, we have used ωr which is just another way of writing the basis vector er that emphasizes its role as a form.

We can similarly define vector valued 1-forms dv̂ and dv from Eqs 59, 49 respectively. We then straightforwardly obtain

dv̂=dv̂wew+dv̂nen,,(107)
dv=dvwew+dvnen.(108)

We express dv in polar coordinates as

dv=vrerdϕ+vreϕdϕ+vϕeϕdϕvϕerdϕ=vrvϕer+vr+vϕeϕdϕ,(109)

giving

Arv̂mdϕ=Arvϕvrdϕ,(110)

or

v̂m=vϕvr.(111)

By combining Eqs 54, 68, 72, we obtain the same equation. Eq. 98 written in the saturation coordinate system yields Eq. 28.

7 Applications

We give in this section a couple of examples of practical use of the ideas that have been presented. We consider first the Brooks–Corey relative permeability model [43] and then a capillary fiber bundle model [44, 45].

7.1 The brooks–corey model

We assume an irreducible wetting fluid saturation Sw,i and residual non-wetting fluid saturation Sn,r so that Sw,iSw ≤ 1 − Sn,r and Sn,rSn ≤ 1 − Sw,i. We renormalize the saturations to be

Sw*=SnwSn,r1Sw,iSn,r,(112)
Sn*=SwSw,i1Sw,iSn,r.(113)

The Brooks–Corey relative permeability is [43]

krw=krw0Sw*nw,(114)
krn=krn0Sn*nn,(115)

The seepage velocities are

vw=Kkrw0φμwSw*nw1PMwSw*nw1,(116)
vn=Kkrn0φμnSn*nn1PMnSn*nn1.(117)

Here, K is the absolute permeability, φ is the porosity, and μw and μn are the viscosities. We assume no gradients in the saturation so that no capillary pressure enters, only the pressure P.

We find the average seepage velocity to be

v=MwSw*nw+MnSn*nn(118)

and the co-moving velocity to be

vm=nw1MwSw*nw1nn1MnSn*nw1.(119)

We also note that

dvdSw*=nwMwSw*nw1nnMnSn*nw1.(120)

Hence, if nw = nn = n, we get

vm=n1ndvdSw*,(121)

a result which is consistent with the observations in Roy et al. [34], see Eq. 29. In the following, we assume that nw = nn = n. This is in accordance with the findings in [46, 47], where in addition was suggested that nw = nn = 2. We will simplify the notation for the saturation Sw*Sw and Sn*Sn.

We find from Eqs 54, 55 that the seepage velocities are in polar coordinates.

vr=Mwcosϕcosϕ+sinϕncosϕ+Mnsinϕcosϕ+sinϕnsinϕ,(122)
vϕ=Mwcosϕcosϕ+sinϕnsinϕ+Mnsinϕcosϕ+sinϕncosϕ.(123)

We find the thermodynamic velocities in polar coordinates using Eqs 53, 79,

v̂r=vr=Mwcosϕcosϕ+sinϕncosϕ+Mnsinϕcosϕ+sinϕnsinϕ(124)

and

v̂ϕ=dv̂rdϕ=ncscϕsecϕcosϕ+sinϕvϕ.(125)

Finally, from Eq. 81, we find the co-moving velocity in polar coordinates,

v̂m=1ncscϕsecϕcosϕ+sinϕvϕ.(126)

7.2 Capillary fiber bundle model

We now turn to the capillary fiber bundle model [44, 45] which serves as a rudimentary model of a porous medium. We will revisit a version of the model that has already been discussed by Hansen et al. [32]. Consider a bundle of N capillary tubes, all of length L. Some fraction of the tubes has a cross-sectional area as and the rest has a cross-sectional area al. We assume al > as. The total pore areas of the small and large capillary tubes are denoted As and Al, respectively. These two areas sum up to the total pore area Ap,

As+Al=Ap.(127)

We assume that the thinnest capillaries are so narrow that the non-wetting fluid may not penetrate them. Hence, the saturation due to the wetting fluid in these capillaries is irreducible, and we have

As=Sw,iAp.(128)

Furthermore, we assume that each capillary is either fully filled with the wetting fluid or the non-wetting fluid, but never a mixture.

In the following, we sketch the analysis of Hansen et al. [32]. The wetting area Aw is then

Aw=Sw,iAp+SwSw,iAp=As+Alw,(129)

where we define Alw ≡ (SwSw,i)Ap, the total wetting area in the larger capillaries. The seepage velocity of the wetting fluid in the smaller and larger capillaries is denoted vsw and vlw respectively. The non-wetting fluid only flows in the larger capillaries. Its seepage velocity is denoted by vn, and the saturation of the non-wetting fluid is given by Sn. The total seepage velocity v is then

v=Sw,ivsw+SwSw,ivlw+Snvn,(130)

where we have 0 ≤ Sn ≤ 1 − Sw,i. We multiply both sides of the equation by Ap to obtain the associated volumetric flow rate Q,

Q=Apv=Asvsw+Alwvlw+Anvn.(131)

To find the co-moving velocity vm, one can explicitly compute the thermodynamic velocities in Eq. 15 or Eq. 16, and use the relations between the thermodynamic- and seepage velocities, Eq. 19 and Eq. 20. The resulting co-moving velocity for this system can then be shown [32] to be

vm=Sw,iSwvlwvsw.(132)

We note that these expressions are fairly complicated. This is due to the choice of origin of the coordinate system, namely Sw = 0 and Ap = 0. It is convenient to shift the origin from (Sw, Ap) = (0, 0) to (Sw,i, As), as conducted in Section 7.1. Hence, we define Sw*=SwSw,i, Ap*=ApAs, and Aw*=AwAs. We immediately drop the star superscripts in these expressions.

We find the thermodynamic velocities easily,

v̂w=QAwAn=vlw,(133)
v̂n=QAnAw=vn.(134)

The wetting volumetric flow rate is

Qw=ApSw,ivsw+Awvlw.(135)

We subtract the volumetric flow rate due to the irreducible wetting fluid.

Qw*=QwApSw,ivsw=AsvlwQw.(136)

The seepage velocity vector then becomes

v=vlwew+vnen.(137)

The thermodynamic velocity vector is

v̂=vlwew+vnen,(138)

and from Eq. 66, we have

v̂m=0.(139)

We note the large difference between Eqs 132, 139; one is a complicated expression, the other one is zero—but both are correct. The difference between the two is where the origin of parameter space is placed.

8 Conclusion and discussion

The aim of this work has been to formulate the immiscible two-phase flow in porous media problem as a geometrical problem in area space. We did this by

• defining the two pore area variables Aw and An, and considering their span,

• endowing this space with different coordinate systems, (Aw, An), (Ap, Sw), and (Ar, ϕ), and pointing out expressions using the polar coordinate system which simplifies the discussion considerably,

• recognizing the meaning of the volumetric flow rate being a degree-1 homogeneous function when expressing it in polar coordinates, and deriving a number of properties related to the different velocities in the problem; the seepage velocities, thermodynamic velocities, and the co-moving velocity,

• using differential forms to derive relations between the different velocities, and lastly,

• formulate an expression for the co-moving velocity which is independent of the coordinates on the underlying space.

We remind the reader of the following equations. The difference between the thermodynamic and seepage velocities may be expressed by combining Eqs 96, 97,

dQ=v̂idAi=Aidvi+vidAi.(140)

The co-moving velocity is given by Eq. 66.

v=v̂+v̂m.

These equations lead us to the three central equations summarizing the central results of this paper. They are Eq. 58,

vA=Q,

Eq. 69,

v̂mA=0,

and Eq. 98,

v̂midAi+Aidvi=0.

It is clear from this discussion that the co-moving velocity vm, which together with the average seepage velocity v, determines vw and vn through the transformations [24, 25], cannot be found by any measurement of the volumetric flow rate Q or any expression obtained from it and it alone such as dQ. The fundamental question which remains open is the following: is it possible to measure the co-moving velocity vm without explicitly measuring vw and vn?

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions

AH wrote the first draft based on the initial work by AH and provided the applications in Section 7. HP wrote the second draft and added geometrical considerations and figures.

Funding

This work was partly supported by the Research Council of Norway through its Centre of Excellence funding scheme, project number 262644.

Acknowledgments

The authors thank Dick Bedeaux, Carl Fredrik Berg, Eirik Grude Flekkøy, Magnus Aa. Gjennestad, Signe Kjelstrup, Marcel Moura, Knut Jørgen Måløy, Santanu Sinha, Per Arne Slotte, and Ole Torsæter for interesting discussions.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Bear J. Dynamics of fluids in porous media. Mineola: Dover (1988). doi:10.1097/00010694-197508000-00022

CrossRef Full Text | Google Scholar

2. Sahimi M. Flow and transport in porous media and fractured rock: From classical methods to modern approaches. New York: Wiley (2011). doi:10.1002/9783527636693

CrossRef Full Text | Google Scholar

3. Blunt MJ. Multiphase flow in permeable media. Cambridge: Cambridge Univ. Press (2017). doi:10.1017/9781316145098

CrossRef Full Text | Google Scholar

4. Feder J, Flekkøy EG, Hansen A. Physics of flow in porous media. Cambridge: Cambridge Univ. Press (2022). doi:10.1017/9781009100717

CrossRef Full Text | Google Scholar

5. Wyckoff RD, Botset HG. The flow of gas-liquid mixtures through unconsolidated sands. Physics (1936) 7:325–45. doi:10.1063/1.1745402

CrossRef Full Text | Google Scholar

6. Leverett MC. Capillary behavior in porous sands. Trans AIMME (1940) 12:152.

Google Scholar

7. Dandekar Ay. Petroleum reservoir rock and fluid properties. 2. ed. Boca Raton: CRC Press (2013).

Google Scholar

8. Barenblatt GI, Patzek TW, Silin BB. The mathematical model of non-equilibrium effects in water-oil displacement. SPE-75169-MS (2002). doi:10.2118/75169-MS

CrossRef Full Text | Google Scholar

9. Wang Y, Aryana SA, Allen MB. An extension of Darcy?s law incorporating dynamic length scales. Adv Water Res (2019) 129:70–9. doi:10.1016/j.advwatres.2019.05.010

CrossRef Full Text | Google Scholar

10. Hassanizadeh SM, Gray WG. Mechanics and thermodynamics of multiphase flow in porous media including interphase boundaries. Adv Water Res (1990) 13:169–86. doi:10.1016/0309-1708(90)90040-B

CrossRef Full Text | Google Scholar

11. Hassanizadeh SM, Gray WG. Toward an improved description of the physics of two-phase flow. Adv Water Res (1993) 16:53–67. doi:10.1016/0309-1708(93)90029-F

CrossRef Full Text | Google Scholar

12. Hassanizadeh SM, Gray WG. Thermodynamic basis of capillary pressure in porous media. Water Resour Res (1993) 29:3389–405. doi:10.1029/93WR01495

CrossRef Full Text | Google Scholar

13. Niessner J, Berg S, Hassanizadeh SM. Comparison of two-phase Darcy’s law with a thermodynamically consistent approach. Transp Por Med (2011) 88:133–48. doi:10.1007/s11242-011-9730-0

CrossRef Full Text | Google Scholar

14. Gray WG, Miller CT. Introduction to the thermodynamically constrained averaging theory for porous medium systems. Berlin: Springer-Verlag (2014). doi:10.1007/978-3-319-04010-3

CrossRef Full Text | Google Scholar

15. Whitaker S. Flow in porous media II: The governing equations for immiscible, two-phase flow. Transp Por Med (1986) 1:105–25. doi:10.1007/BF00714688

CrossRef Full Text | Google Scholar

16. Kjelstrup S, Bedeaux D, Hansen A, Hafskjold B, Galteland O. Non-isothermal transport of multi-phase fluids in porous media. the entropy production. Front Phys (2018) 6:126. doi:10.3389/fphy.2018.00126

CrossRef Full Text | Google Scholar

17. Kjelstrup S, Bedeaux D, Hansen A, Hafskjold B, Galteland O. Non-isothermal transport of multi-phase fluids in porous media. Constitutive equations. Front Phys (2019) 6:150. doi:10.3389/fphy.2018.00150

CrossRef Full Text | Google Scholar

18. Bedeaux D, Kjelstrup S. Fluctuation-dissipation theorems for multiphase flow in porous media. Entropy (2022) 24:46. doi:10.3390/e24010046

CrossRef Full Text | Google Scholar

19. McClure JE, Armstrong RT, Berrill MA, Schlüter S, Berg S, Gray WG, et al. Geometric state function for two-fluid flow in porous media. Phys Rev Fluids (2018) 3:084306. doi:10.1103/physrevfluids.3.084306

CrossRef Full Text | Google Scholar

20. Armstrong RT, McClure JE, Robins V, Liu Z, Arns CH, Schlüter S, et al. Porous media characterization using Minkowski functionals: theories, applications and future directions. Transp Por Med (2019) 130:305–335. doi:10.1007/s11242-018-1201-4

CrossRef Full Text | Google Scholar

21. McClure JE, Armstrong RT, Berg S. Geometric evolution as a source of discontinuous behavior in soft condensed matter. arXiv:1906.04073. doi:10.48550/arXiv.1906.04073

CrossRef Full Text

22. Hilfer R, Besserer H. Macroscopic two-phase flow in porous media. Phys B (2000) 279:125–9. doi:10.1016/S0921-4526(99)00694-8

CrossRef Full Text | Google Scholar

23. Hilfer R. Capillary pressure, hysteresis and residual saturation in porous media. Phys A (2006) 359:119–28. doi:10.1016/j.physa.2005.05.086

CrossRef Full Text | Google Scholar

24. Hilfer R. Macroscopic capillarity and hysteresis for flow in porous media. Phys Rev E (2006) 73:016307. doi:10.1103/PhysRevE.73.016307

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Hilfer R. Macroscopic capillarity without a constitutive capillary pressure function. Physica A (2006) 371:209–25. doi:10.1016/j.physa.2006.04.051

CrossRef Full Text | Google Scholar

26. Hilfer R, Döster F. Percolation as a basic concept for macroscopic capillarity. Transp Por Med (2010) 82:507–19. doi:10.1007/s11242-009-9395-0

CrossRef Full Text | Google Scholar

27. Döster F, Hönig O, Hilfer R. Horizontal flow and capillarity-driven redistribution in porous media. Phys Rev E (2012) 86:016317. doi:10.1103/PhysRevE.86.016317

CrossRef Full Text | Google Scholar

28. Valavanides MS, Constantinides GN, Payatakes AC. Mechanistic Model of steady-state two-phase flow in porous media based on ganglion dynamics. Transp Por Med (1998) 30:267–99. doi:10.1023/A:1006558121674

CrossRef Full Text | Google Scholar

29. Valavanides MS. Steady-state two-phase flow in porous media: review of progress in the development of the DeProF theory bridging pore-to statistical thermodynamics-scales. Oil Gas Sci Technol (2012) 67:787–804. doi:10.2516/ogst/2012056

CrossRef Full Text | Google Scholar

30. Valavanides MS. Review of steady-state two-phase flow in porous media: independent variables, universal energy efficiency map, critical flow conditions, effective characterization of flow and pore network. Transp Porous Media (2018) 123:45–99. doi:10.1007/s11242-018-1026-1

CrossRef Full Text | Google Scholar

31. Xu J, Louge MY. Statistical mechanics of unsaturated porous media. Phys Rev E (2015) 92:062405. doi:10.1103/PhysRevE.92.062405

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Hansen A, Sinha S, Bedeaux D, Kjelstrup S, Gjennestad MA, Vassvik M. Relations between seepage velocities in immiscible, incompressible two-phase flow in porous Media. Transp Porous Media (2018) 125:565–87. doi:10.1007/s11242-018-1139-6

CrossRef Full Text | Google Scholar

33. Roy S, Sinha S, Hansen A. Flow-area relations in immiscible two-phase flow in porous media. Front Phys (2020) 8:4. doi:10.3389/fphy.2020.00004

CrossRef Full Text | Google Scholar

34. Roy S, Pedersen H, Sinha S, Hansen A. The Co-moving Velocity in immiscible two-phase flow in porous media. Transp Porous Media (2022) 143:69–102. doi:10.1007/s11242-022-01783-7

CrossRef Full Text | Google Scholar

35. Hansen A, Flekkøy EG, Sinha S, Slotte PA. A statistical mechanics framework for immiscible and incompressible two-phase flow in porous media. Adv Water Res (2022) 171:104336. doi:10.1016/j.advwatres.2022.104336

CrossRef Full Text | Google Scholar

36. Fyhn H, Sinha S, Hansen A. Local statistics of immiscible and incompressible two-phase flow in porous media. arXiv:2209.00030. doi:10.48550/arXiv.2209.00030

CrossRef Full Text

37. Jaynes ET. Information theory and statistical mechanics. Phys Rev (1957) 106:620–30. doi:10.1103/PhysRev.106.620);

CrossRef Full Text | Google Scholar

38. Ekrann S, Aasen JO. Steady-state upscaling. Transp Porous Media (2000) 41:245. doi:10.1023/A:10067654

CrossRef Full Text | Google Scholar

39. Bear J, Bachmat Y. Introduction to modeling of transport phenomena in porous media. Berlin: Springer (2012). doi:10.1007/978-94-009-1926-6

CrossRef Full Text | Google Scholar

40. Shannon CE. A Mathematical theory of communication. Bell Syst Tech J (1948) 27:379–423. doi:10.1002/j.1538-7305.1948.tb01338.x

CrossRef Full Text | Google Scholar

41. Flanders H. Differential forms. New York: Academic Press (1963).

Google Scholar

42. Misner C: W, Thorne KS, Wheeler JA. Gravitation. Princeton: Princeton University Press (2017).

Google Scholar

43. Lake LW. Enhanced oil recovery. Old Tappan: Prentice-Hall (1989).

Google Scholar

44. Scheidegger AE. Theoretical models of porous matter. Prod Monthly (1953) 17:17.

Google Scholar

45. Scheidegger Ae. The physics of flow through porous media. Toronto: University of Toronto Press (1974).

Google Scholar

46. Picchi D, Battiato I. The impact of pore-scale flow regimes on upscaling of immiscible two-phase flow in porous media. Water Resour Res (2018) 54:6683–707. doi:10.1029/2018WR023172

CrossRef Full Text | Google Scholar

47. Picchi D, Battiato I. Relative permeability scaling from pore-scale flow regimes. Water Resour Res (2019) 55:3215–33. doi:10.1029/2018WR024251

CrossRef Full Text | Google Scholar

Keywords: immiscible two-phase flow, porous media, co-moving velocity, seepage velocity, thermodynamics, differential geometry, differential forms, Euler homogeneity

Citation: Pedersen H and Hansen A (2023) Parameterizations of immiscible two-phase flow in porous media. Front. Phys. 11:1127345. doi: 10.3389/fphy.2023.1127345

Received: 19 December 2022; Accepted: 06 February 2023;
Published: 27 February 2023.

Edited by:

Eric Josef Ribeiro Parteli, University of Duisburg-Essen, Germany

Reviewed by:

Humberto Carmona, Federal University of Ceara, Brazil
Muhammad Sahimi, University of Southern California, United States

Copyright © 2023 Pedersen and Hansen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Håkon Pedersen, hakon.pedersen@ntnu.no

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.