Cluster evolution in steady-state two-phase flow in porous media

Morphogenesis of Biofilms in Porous Media and Control on Hydrodynamics

The functioning of natural and engineered porous media, like soils and filters, depends in many cases on the interplay between biochemical processes and hydrodynamics. In such complex environments, microorganisms often form surface-attached communities known as biofilms. Biofilms can take the shape of clusters, which alter the distribution of fluid flow velocities within the porous medium, subsequently influencing biofilm growth. Despite numerous experimental and numerical efforts, the control of the biofilm clustering process and the resulting heterogeneity in biofilm permeability is not well understood, limiting our predictive abilities for biofilm-porous medium systems. Here, we use a quasi-2D experimental model of a porous medium to characterize biofilm growth dynamics for different pore sizes and flow rates. We present a method to obtain the time-resolved biofilm permeability field from experimental images and use the obtained permeability field to compute the flow field through a numerical model. We observe a biofilm cluster size distribution characterized by a spectrum slope evolving in time between -2 and -1, a fundamental measure that can be used to create spatio-temporal distributions of biofilm clusters for upscaled models. We find a previously undescribed biofilm permeability distribution, which can be used to stochastically generate permeability fields within biofilms. An increase in velocity variance for a decrease in physical heterogeneity shows that the bioclogged porous medium behaves differently than expected from studies on heterogeneity in abiotic porous media.

Pore Network Modeling of the Effects of Viscosity Ratio and Pressure Gradient on Steady-State Incompressible Two-Phase Flow in Porous Media

We perform steady-state simulations with a dynamic pore network model, corresponding to a large span in viscosity ratios and capillary numbers. From these simulations, dimensionless steady-state time-averaged quantities such as relative permeabilities, residual saturations, mobility ratios and fractional flows are computed. These quantities are found to depend on three dimensionless variables, the wetting fluid saturation, the viscosity ratio and a dimensionless pressure gradient. Relative permeabilities and residual saturations show many of the same qualitative features observed in other experimental and modeling studies. The relative permeabilities do not approach straight lines at high capillary numbers for viscosity ratios different from 1. Our conclusion is that this is because the fluids are not in the highly miscible near-critical region. Instead they have a viscosity disparity and intermix rather than forming decoupled, similar flow channels. Ratios of average mobility to their high capillary number limit values are also considered. Roughly, these vary between 0 and 1, although values larger than 1 are also observed. For a given saturation, the mobilities are not always monotonically increasing with the pressure gradient. While increasing the pressure gradient mobilizes more fluid and activates more flow paths, when the mobilized fluid is more viscous, a reduction in average mobility may occur.

Velocity distributions in trapped and mobilized non-wetting phase ganglia in porous media

Understanding the mobilisation of trapped globules of non-wetting phase during two-phase flow has been the aim of numerous studies. However, the driving forces for the mobilisation of the trapped phases are still not well understood. Also, there is little information about what happens within a globule before, at the onset and during mobilization. In this work, we used micro-particle tracking velocimetry in a micro-fluidic model in order to visualise the velocity distributions inside the trapped phase globules prior and during mobilisation. Therefore, time-averaged and instantaneous velocity vectors have been determined using fluorescent microscopy. As a porous medium, we used a polydimethylsiloxane (PDMS) micro-model with a well-defined pore structure, where drainage and imbibition experiments were conducted. Three different geometries of trapped non-wetting globules, namely droplets, blobs and ganglia were investigated. We observed internal circulations inside the trapped phase globules, leading to the formation of vortices. The direction of circulating flow within a globule is dictated by the drag force exerted on it by the flowing wetting phase. This is illustrated by calculating and analyzing the drag force (per unit area) along fluid-fluid interfaces. In the case of droplets and blobs, only one vortex is formed. The flow field within a ganglion is much more complex and more vortices can be formed. The circulation velocities are largest at the fluid-fluid interfaces, along which the wetting phase flows and decreases towards the middle of the globule. The circulation velocities increased proportionally with the increase of wetting phase average velocity (or capillary number). The vortices remain stable as long as the globules are trapped, start to change at the onset of mobilization and disappear during the movement of globules. They reappear when the globules get stranded. Droplets are less prone to mobilization; blobs get mobilised in whole; while ganglia may get ruptured and get mobilised only partially.

Influence of phase connectivity on the relationship among capillary pressure, fluid saturation, and interfacial area in two-fluid-phase porous medium systems

Multiphase flows in porous medium systems are typically modeled at the macroscale by applying the principles of continuum mechanics to develop models that describe the behavior of averaged quantities, such as fluid pressure and saturation. These models require closure relations to produce solvable forms. One of these required closure relations is an expression relating the capillary pressure to fluid saturation and, in some cases, other topological invariants such as interfacial area and the Euler characteristic (or average Gaussian curvature). The forms that are used in traditional models, which typically consider only the relationship between capillary pressure and saturation, are hysteretic. An unresolved question is whether the inclusion of additional morphological and topological measures can lead to a nonhysteretic closure relation. Relying on the lattice Boltzmann (LB) method, we develop an approach to investigate equilibrium states for a two-fluid-phase porous medium system, which includes disconnected nonwetting phase features. A set of simulations are performed within a random close pack of 1964 spheres to produce a total of 42 908 distinct equilibrium configurations. This information is evaluated using generalized additive models to quantitatively assess the degree to which functional relationships can explain the behavior of the equilibrium data. The variance of various model estimates is computed, and we conclude that, except for the limiting behavior close to a single fluid regime, capillary pressure can be expressed as a deterministic and nonhysteretic function of fluid saturation, interfacial area between the fluid phases, and the Euler characteristic. To our knowledge, this work is unique in the methods employed, the size of the data set, the resolution in space and time, the true equilibrium nature of the data, the parametrizations investigated, and the broad set of functions examined. The conclusion of essentially nonhysteretic behavior provides support for an evolving class of two-fluid-phase flow in porous medium systems models.

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.

Strong influence of geometrical heterogeneity on drainage in porous media

We present an experimental study of drainage in two-dimensional porous media exhibiting bimodal pore size distributions. The role of the pore size heterogeneity is investigated by measuring separately the desaturation curves of the two pore populations. The displaced wetting fluid remains trapped in small pores at low capillary numbers and is swept only above a critical capillary number proportional to the permeability of the big pores network. Based on this observation, we derive a simple criterion for phase trapping based on the balance of viscous to capillary forces. Numerical implementation of this theory in a pore network model quantitatively fits our experimental results. This combination of approaches demonstrates quantitatively the influence of geometrical heterogeneities on drainage in porous media.

Local wettability reversal during steady-state two-phase flow in porous media

We study the effect of local wettability reversal on remobilizing immobile fluid clusters in steady-state two-phase flow in porous media. We consider a two-dimensional network model for a porous medium and introduce a wettability alteration mechanism. A qualitative change in the steady-state flow patterns, destabilizing the percolating and trapped clusters, is observed as the system wettability is varied. When capillary forces are strong, a finite wettability alteration is necessary to move the system from a single-phase to a two-phase flow regime. When both phases are mobile, we find a linear relationship between fractional flow and wettability alteration.

Oscillation-induced displacement patterns in a two-dimensional porous medium: a lattice Boltzmann study

We present a numerical study of the statistical behavior of a two-phase flow in a two-dimensional porous medium subjected to an oscillatory acceleration transverse to the overall direction of flow. A viscous nonwetting fluid is injected into a porous medium filled with a more viscous wetting fluid. During the whole process sinusoidal oscillations of constant amplitude and frequency accelerates the porous medium sideways, perpendicular to the overall direction of flow. The invasion process displays a transient behavior where the saturation of the defending fluid decreases, before it enters a state of irreducible wetting fluid saturation, where there is no net transport of defending fluid toward the outlet of the system. In this state the distribution of sizes of the remaining clusters are observed to obey a power law with an exponential cutoff. The cutoff cluster size is found to be determined by the flow and oscillatory stimulation parameters. This cutoff size is also shown to be directly related to the extracted amount of defending fluid. Specifically, the results show that the oscillatory acceleration of the system leads to potentially a large increase in extracted wetting fluid.

Steady-state, simultaneous two-phase flow in porous media: an experimental study

We report on experimental studies of steady-state two-phase flow in a quasi-two-dimensional porous medium. The wetting and the nonwetting phases are injected simultaneously from alternating inlet points into a Hele-Shaw cell containing one layer of randomly distributed glass beads, initially saturated with wetting fluid. The high viscous wetting phase and the low viscous nonwetting phase give a low viscosity ratio M=10(-4). Transient behavior of this system is observed in time and space. However, we find that at a certain distance behind the initial front a "local" steady-state develops, sharing the same properties as the later "global" steady state. In this state the nonwetting phase is fragmented into clusters, whose size distribution is shown to obey a scaling law, and the cutoff cluster size is found to be inversely proportional to the capillary number. The steady state is dominated by bubble dynamics, and we measure a power-law relationship between the pressure gradient and the capillary number. In fact, we demonstrate that there is a characteristic length scale in the system, depending on the capillary number through the pressure gradient that controls the steady-state dynamics.

Flux-dependent percolation transition in immiscible two-phase flows in porous media

Using numerical simulations, we study immiscible two-phase flow in a pore network reconstructed from Berea sandstone under flow conditions that are statistically invariant under translation. Under such conditions, the flow is a state function which is not dependent on initial conditions. We find a second-order phase transition resembling the phase inversion transition found in emulsions. The flow regimes under consideration are those of low surface tension-hence high capillary numbers Ca-where viscous forces dominate. Nevertheless, capillary forces are imminent, we observe a critical stage in saturation where the transition takes place. We determine polydispersity critical exponent tau=2.27+/-0.08 and find that the critical saturation depends on how fast the fluids flow.

Steady-state two-phase flow in porous media: statistics and transport properties

We study experimentally the case of steady-state simultaneous two-phase flow in a quasi-two-dimensional porous media. The dynamics is dominated by the interplay between a viscous pressure field from the wetting fluid and bubble transport of a less viscous, nonwetting phase. In contrast with more studied displacement front systems, steady-state flow is in equilibrium, statistically speaking. The corresponding theoretical simplicity allows us to explain a data collapse in the cluster size distribution as well as the relation |nablaP| proportional, sqrt[Ca] between the pressure gradient in the system and the capillary number.

Capillary-driven instability of immiscible fluid interfaces flowing in parallel in porous media

When immiscible wetting and nonwetting fluids move in parallel in a porous medium, an instability may occur at sufficiently high capillary numbers so that interfaces between the fluids initially held in place by the porous medium are mobilized. A boundary zone containing bubbles of both fluids evolves, which has a well-defined thickness. This zone moves at constant average speed toward the nonwetting fluid. A diffusive current of bubbles of nonwetting fluid into the wetting fluid is set up.