In this paper, we study the influence of sample geometry on the measurement of pressuresaturation relationships, by analyzing the drainage of a two-phase flow from a quasi-2-D random porous medium. The medium is transparent, which allows for the direct visualization of the invasion pattern during flow, and is initially saturated with a viscous liquid (a dyed glycerol-water mix). As the pressure in the liquid is gradually reduced, air penetrates from an open inlet, displacing the liquid which leaves the system from an outlet on the opposite side. Pressure measurements and images of the flow are recorded and the pressure-saturation relationship is computed. We show that this relationship depends on the system size and aspect ratio. The effects of the system's boundaries on this relationship are measured experimentally and compared with simulations produced using an invasion percolation algorithm. The pressure build up at the beginning and end of the invasion process are particularly affected by the boundaries of the system whereas at the central part of the model (when the air front progresses far from these boundaries), the invasion happens at a statistically constant capillary pressure. These observations have led us to propose a much simplified pressure-saturation relationship, valid for systems that are large enough such that the invasion is not influenced by boundary effects. The properties of this relationship depend on the capillary pressure thresholds distribution, sample dimensions, and average pore connectivity and its applications may be of particular interest for simulations of two-phase flow in large porous media.

The intermittent burst dynamics during the slow drainage of a porous medium is studied experimentally. We have shown that this system satisfies a set of conditions known to be true for critical systems, such as intermittent activity with bursts extending over several time and length scales, self-similar macroscopic fractal structure and 1/f α power spectrum. Additionally, we have verified a theoretically predicted scaling for the burst size distribution, previously assessed via numerical simulations. The observation of 1/f α power spectra is new for porous media flows and, for specific boundary conditions, we notice the occurrence of a transition from 1/f to 1/f 2 scaling. An analytically integrable mathematical framework was employed to explain this behavior.

We study the effects of connectivity enhancement due to film flow phenomena on the drainage of an artificial porous medium and the relative influence of gravity for such effects. The medium is initially fully saturated with a liquid (wetting phase), which is displaced by air (nonwetting phase). Our setup allows us to directly visualize the dynamics of the flow and, in particular, to pinpoint which pore invasion events are due to film flow phenomena. We have observed the formation of an active zone behind the liquid-air interface, inside which film flow drainage events are more likely to occur. Understanding the basic mechanisms of film flow in artificial porous media is of relevance for the analysis of real three-dimensional porous systems, in which gravity cannot be neglected and film flow is bound to be present.

In this Letter we give experimental grounding for the remarkable observation made by Furuberg et al. [Phys. Rev. Lett. 61, 2117 (1988)PRLTAO0031-900710.1103/PhysRevLett.61.2117] of an unusual dynamic scaling for the pair correlation function N(r,t) during the slow drainage of a porous medium. Those authors use an invasion percolation algorithm to show numerically that the probability of invasion of a pore at a distance r away and after a time t from the invasion of another pore scales as N(r,t)∝r^{-1}f(r^{D}/t), where D is the fractal dimension of the invading cluster and the function f(u)∝u^{1.4}, for u≪1 and f(u)∝u^{-0.6}, for u≫1. Our experimental setup allows us to have full access to the spatiotemporal evolution of the invasion, which is used to directly verify this scaling. Additionally, we connect two important theoretical contributions from the literature to explain the functional dependency of N(r,t) and the scaling exponent for the short-time regime (t≪r^{D}). A new theoretical argument is developed to explain the long-time regime exponent (t≫r^{D}).

The intermittent dynamics of slow drainage flows in a porous medium is studied experimentally. This kind of two-phase flow is characterized by a rich burst activity and our setup allows us to characterize those bursts directly via images of the flow and pressure measurements. Two different boundary conditions were analyzed: controlled withdrawal rate (CWR) and controlled imposed pressure (CIP). We have characterized geometrical and statistical properties of the bursts from images and pressure measurements. We have shown that in spite of leading to similar final invasion patterns, some dynamical features of the invasion differ considerably between the CWR and CIP boundary conditions. In particular, their pressure signatures are very distinct, which then translates into very distinct features on the power spectrum density of the pressure signals. A fully integrable analytical framework is presented which successfully describes the scaling features of the power spectrum for the CIP case.

We present a theoretical and experimental investigation of drainage in porous media. The study is limited to stabilized fluid fronts at moderate injection rates, but it takes into account capillary, viscous, and gravitational forces. In the theoretical framework presented, the work applied on the system, the energy dissipation, the final saturation and the width of the stabilized fluid front can all be calculated if we know the dimensionless fluctuation number, the wetting properties, the surface tension between the fluids, the fractal dimensions of the invading structure and its boundary, and the exponent describing the divergence of the correlation length in percolation. Furthermore, our theoretical description explains how the Haines jumps’ local activity and dissipation relate to dissipation on larger scales.

The motion of a pair of counter-rotating point vortices placed in a uniform flow around a circular cylinder forms a rich nonlinear system that is often used to model vortex shedding. The phase portrait of the Hamiltonian governing the dynamics of a vortex pair that moves symmetrically with respect to the centerline---a case that can be realized experimentally by placing a splitter plate in the center plane---is presented. The analysis provides new insights and reveals novel dynamical features of the system, such as a nilpotent saddle point at infinity whose homoclinic orbits define the region of nonlinear stability of the so-called F\"oppl equilibrium. It is pointed out that a vortex pair properly placed downstream can overcome the cylinder and move off to infinity upstream. In addition, the nonlinear dynamics resulting from antisymmetric perturbations of the F\"oppl equilibrium is studied and its relevance to vortex shedding discussed.Comment: 21 pages, 6 figure

The flow of multiple fluids through a rock is a dynamic process involving fluid-fluid-rock interactions over multiple spatial and temporal scales. For instance, to understand the trapping of CO 2 underground, research spans from the sub mm-scale of high-resolution X-ray imaging experiments to observe pore scale phenomena