The Co-Moving Velocity and Projective Transformations

TL;DR

This paper develops a geometric and thermodynamic framework for two-phase flow in porous media, introducing the co-moving velocity and projective transformations to relate flow properties and saturation.

Contribution

It introduces a novel geometric approach using projective duality and affine geometry to model co-moving velocity in two-phase flow, connecting thermodynamics and geometry.

Findings

Derived an exact constitutive relation using projective invariants.

Expressed co-moving velocity through a simple steady-state advection equation.

Validated the parametrization with pore network model data.

Abstract

In a string of recent papers starting with (Transport in Porous Media, 125, 565 (2018)), a theory of immiscible two-phase flow in porous media based on Euler homogeneity of the total volumetric flow rate has been investigated. The thermodynamic-like theory has an associated statistical mechanics based on a maximum entropy principle. A quantity called the co-moving velocity connects the equations of state of the intensive thermodynamic velocities and the physical seepage velocities of two flowing fluids. The obtained relations have a structure that can be interpreted using affine- and projective geometry. The co-moving velocity can be expressed as a transformation of the saturation using projective duality of points and lines. One obtains an exact constitutive relation depending on a projective invariant, the cross-ratio, which allows the co-moving velocity to be expressed in terms of a simple steady-state advection equation. A kinematic view of the velocity relations is presented, modeled by a well-known non-trivial geometry which turns out to be pseudo-Euclidean. The cross-ratio determines a hyperbolic angle in this space, and can be parametrized in terms of three numbers using a linear fractional transformation. Knowing these parameters, the pore velocity and the derivative of the pore velocity with respect to the saturation, an approximation for the co-moving velocity can be obtained for a range of applied pressures, viscosity ratios and surface tensions. The parametrization is demonstrated using data from a dynamic pore network model and relative permeability data from the literature. This paper only considers the pore areas as extensive variables, however, the geometric principles are general, and the same idea could potentially be used in other systems.