Analysis and computation of gravity-induced migration in porous media

Motivated by the problem of gravity segregation in an inclined porous layer, we present a theoretical analysis of interface evolution between two immiscible fluids of unequal density and mobility, both in two and three dimensions. Applying perturbation theory to the appropriately scaled problem, we derive the governing equations for the pressure and interface height to leading order, obtained in the limit of a thin gravity tongue and a slightly dipping bed. According to the zeroth-order approximation, the pressure profile perpendicular to the bed is in equilibrium, a widely accepted assumption for this class of problems. We show that for the inclined bed two-dimensional problem, in the reference frame moving with the mean gravity-induced advection velocity, the interface motion is dictated by a degenerate parabolic equation, different from those previously published. In this case, the late-time behaviour of the gravity tongue can be derived analytically through a formal expansion of both the solution and its two moving boundaries. In three dimensions, using a moving coordinate along the dip direction, we obtain an elliptic–parabolic system of partial differential equations where the fluid pressure and interface height are the two dependent variables. Although analytical results are not available for this case, the evolution of the gravity tongue can be investigated by numerical computations in only two spatial dimensions. The solution features are identified for different combinations of dimensionless parameters, showing their respective influence on the shape and motion of the interface.