Effective rheology of bubbles moving in a capillary tube

Steady-State Two-Phase Flow of Compressible and Incompressible Fluids in a Capillary Tube of Varying Radius

We study immiscible two-phase flow of a compressible and an incompressible fluid inside a capillary tube of varying radius under steady-state conditions. The incompressible fluid is Newtonian and the compressible fluid is an inviscid ideal gas. The surface tension associated with the interfaces between the two fluids introduces capillary forces that vary along the tube due to the variation in the tube radius. The interplay between effects due to the capillary forces and the compressibility results in a set of properties that are different from incompressible two-phase flow. As the fluids move towards the outlet, the bubbles of the compressible fluid grows in volume due to the decrease in pressure. The volumetric growth of the compressible bubbles makes the volumetric flow rate at the outlet higher than at the inlet. The growth is not only a function of the pressure drop across the tube, but also of the ambient pressure. Furthermore, the capillary forces create an effective threshold below which there is no flow. Above the threshold, the system shows a weak nonlinearity between the flow rates and the effective pressure drop, where the nonlinearity also depends on the absolute pressures across the tube.

Non-Newtonian Rheology in a Capillary Tube with Varying Radius

Fluid blobs in an immiscible Newtonian fluid flowing in a capillary tube with varying radius show highly nonlinear behavior. We consider here a generalization of previously obtained results to blobs of non-Newtonian fluids. We compute here the yield pressure drop and the mean flow rate in two cases: (i) When a single blob is injected, (ii) When many blobs are randomly injected into the tube. We find that the capillary effects emerge from the non-uniformity of the tube radius and contribute to the threshold pressure for flow to occur. Furthermore, in the presence of many blobs the threshold value depends on the number of blobs and their relative distances which are randomly distributed. For a capillary fiber bundle of identical parallel tubes, we calculate the probability distribution of the threshold pressure and the mean flow rate. We consider two geometries: tubes of sinusoidal shape, for which we derive explicit expressions, and triangular-shaped tubes, for which we find that essential singularities are developed. We perform numerical simulations confirming our analytical results.

The Co-Moving Velocity in Immiscible Two-Phase Flow in Porous Media

We present a continuum (i.e., an effective) description of immiscible two-phase flow in porous media characterized by two fields, the pressure and the saturation. Gradients in these two fields are the driving forces that move the immiscible fluids around. The fluids are characterized by two seepage velocity fields, one for each fluid. Following Hansen et al. (Transport in Porous Media, 125, 565 (2018)), we construct a two-way transformation between the velocity couple consisting of the seepage velocity of each fluid, to a velocity couple consisting of the average seepage velocity of both fluids and a new velocity parameter, the co-moving velocity. The co-moving velocity is related but not equal to velocity difference between the two immiscible fluids. The two-way mapping, the mass conservation equation and the constitutive equations for the average seepage velocity and the co-moving velocity form a closed set of equations that determine the flow. There is growing experimental, computational and theoretical evidence that constitutive equation for the average seepage velocity has the form of a power law in the pressure gradient over a wide range of capillary numbers. Through the transformation between the two velocity couples, this constitutive equation may be taken directly into account in the equations describing the flow of each fluid. This is, e.g., not possible using relative permeability theory. By reverse engineering relative permeability data from the literature, we construct the constitutive equation for the co-moving velocity. We also calculate the co-moving constitutive equation using a dynamic pore network model over a wide range of parameters, from where the flow is viscosity dominated to where the capillary and viscous forces compete. Both the relative permeability data from the literature and the dynamic pore network model give the same very simple functional form for the constitutive equation over the whole range of parameters.

Rheology of Immiscible Two-phase Flow in Mixed Wet Porous Media: Dynamic Pore Network Model and Capillary Fiber Bundle Model Results

Immiscible two-phase flow in porous media with mixed wet conditions was examined using a capillary fiber bundle model, which is analytically solvable, and a dynamic pore network model. The mixed wettability was implemented in the models by allowing each tube or link to have a different wetting angle chosen randomly from a given distribution. Both models showed that mixed wettability can have significant influence on the rheology in terms of the dependence of the global volumetric flow rate on the global pressure drop. In the capillary fiber bundle model, for small pressure drops when only a small fraction of the tubes were open, it was found that the volumetric flow rate depended on the excess pressure drop as a power law with an exponent equal to 3/2 or 2 depending on the minimum pressure drop necessary for flow. When all the tubes were open due to a high pressure drop, the volumetric flow rate depended linearly on the pressure drop, independent of the wettability. In the transition region in between where most of the tubes opened, the volumetric flow depended more sensitively on the wetting angle distribution function and was in general not a simple power law. The dynamic pore network model results also showed a linear dependence of the flow rate on the pressure drop when the pressure drop is large. However, out of this limit the dynamic pore network model demonstrated a more complicated behavior that depended on the mixed wettability condition and the saturation. In particular, the exponent relating volumetric flow rate to the excess pressure drop could take on values anywhere between 1.0 and 1.8. The values of the exponent were highest for saturations approaching 0.5, also, the exponent generally increased when the difference in wettability of the two fluids were larger and when this difference was present for a larger fraction of the porous network.

Effective Rheology of Two-Phase Flow in Three-Dimensional Porous Media: Experiment and Simulation

We present an experimental and numerical study of immiscible two-phase flow of Newtonian fluids in three-dimensional (3D) porous media to find the relationship between the volumetric flow rate (Q) and the total pressure difference ([Formula: see text]) in the steady state. We show that in the regime where capillary forces compete with the viscous forces, the distribution of capillary barriers at the interfaces effectively creates a yield threshold ([Formula: see text]), making the fluids reminiscent of a Bingham viscoplastic fluid in the porous medium. In this regime, Q depends quadratically on an excess pressure drop ([Formula: see text]). While increasing the flow rate, there is a transition, beyond which the overall flow is Newtonian and the relationship is linear. In our experiments, we build a model porous medium using a column of glass beads transporting two fluids, deionized water and air. For the numerical study, reconstructed 3D pore networks from real core samples are considered and the transport of wetting and non-wetting fluids through the network is modeled by tracking the fluid interfaces with time. We find agreement between our numerical and experimental results. Our results match with the mean-field results reported earlier.

Ensemble distribution for immiscible two-phase flow in porous media

We construct an ensemble distribution to describe steady immiscible two-phase flow of two incompressible fluids in a porous medium. The system is found to be ergodic. The distribution is used to compute macroscopic flow parameters. In particular, we find an expression for the overall mobility of the system from the ensemble distribution. The entropy production at the scale of the porous medium is shown to give the expected product of the average flow and its driving force, obtained from a black-box description. We test numerically some of the central theoretical results.

History independence of steady state in simultaneous two-phase flow through two-dimensional porous media

It is well known that the transient behavior during drainage or imbibition in multiphase flow in porous media strongly depends on the history and initial condition of the system. However, when the steady-state regime is reached and both drainage and imbibition take place at the pore level, the influence of the evolution history and initial preparation is an open question. Here, we present an extensive experimental and numerical work investigating the history dependence of simultaneous steady-state two-phase flow through porous media. Our experimental system consists of a Hele-Shaw cell filled with glass beads which we model numerically by a network of disordered pores transporting two immiscible fluids. From measurements of global pressure evolution, histograms of saturation, and cluster-size distributions, we find that when both phases are flowing through the porous medium, the steady state does not depend on the initial preparation of the system or on the way it has been reached.