No Arabic abstract
In particle-laden flows through porous media, porosity and permeability are significantly affected by the deposition and erosion of particles. Experiments show that the permeability evolution of a porous medium with respect to a particle suspension is not smooth, but rather exhibits significant jumps followed by longer periods of continuous permeability decrease. Their origin seems to be related to internal flow path reorganization by avalanches of deposited material due to erosion inside the porous medium. We apply neutron tomography to resolve the spatio-temporal evolution of the pore space during clogging and unclogging to prove the hypothesis of flow path reorganization behind the permeability jumps. This mechanistic understanding of clogging phenomena is relevant for a number of applications from oil production to filters or suffosion as the mechanisms behind sinkhole formation.
The flow of viscoelastic fluids in porous media is encountered in many practical applications, such as in the enhanced oil recovery process or in the groundwater remediation. Once the flow rate exceeds a critical value in such flows, an elastic instability with fluctuating flow field is observed, which ultimately transits to a more chaotic and turbulence-like flow structure as the flow rate further increases. In a recent study, it has been experimentally shown that this chaotic flow behaviour of viscoelastic fluids can be suppressed by increasing the geometric disorder in a model porous media consisting of a microchannel with several micropillars placed in it. However, the present numerical study demonstrates that this is not always true. We show that it depends on the initial arrangement of the micropillars for mimicking the porous media. In particular, we find that for an initial ordered and aligned configuration of the micropillars, the introduction of geometric order actually increases the chaotic flow dynamics as opposed to that seen for an initial ordered and staggered configuration of the micropillars. We suggest that this chaotic flow behaviour actually depends on the number of the stagnation points revealed to the flow field where maximum stretching of the viscoelastic microstructure happens. Our findings and explanation are perfectly in line with that observed and provided in a more recent experimental study.
Immiscible fluid displacement in porous media is fundamental for many environmental processes, including infiltration of water in soils, groundwater remediation, enhanced recovery of hydrocarbons and carbon geosequestration. Microstructural heterogeneity, in particular of particle sizes, can significantly impact immiscible displacement. For instance, it may lead to unstable flow and preferential displacement patterns. We present a systematic, quantitative pore-scale study of the impact of spatial correlations in particle sizes on the drainage of a partially-wetting fluid. We perform pore-network simulations with varying flow rates and different degrees of spatial correlation, complemented with microfluidic experiments. Simulated and experimental displacement patterns show that spatial correlation leads to more preferential invasion, with reduced trapping of the defending fluid, especially at low flow rates. Numerically, we find that increasing the correlation length reduces the fluid-fluid interfacial area and the trapping of the defending fluid, and increases the invasion pattern asymmetry and selectivity. Our experiments, conducted for low capillary numbers, support these findings. Our results delineate the significant effect of spatial correlations on fluid displacement in porous media, of relevance to a wide range of natural and engineered processes.
We study the transport of inertial particles in water flow in porous media. Our interest lies in understanding the accumulation of particles including the possibility of clogging. We propose that accumulation can be a result of hydrodynamic effects: the tortuous paths of the porous medium generate regions of dominating strain/vorticity, which favour the accumulation/dispersion of the inertial particles. Numerical simulations show that essentially two accumulation regimes are identified: for low and for high flow velocities. When particles accumulate in high-velocity regions, at the entrance of a pore throat, a clog is formed. The formation of a clog significantly modifies the flow, as the partial blockage of the pore causes a local redistribution of pressure. This redistribution can divert the upstream water flow into neighbouring pores. Moreover, we show that accumulation in high velocity regions occurs in heterogeneous media, but not in homogeneous media, where we refer to homogeneity with respect to the distribution of the pore throat diameters.
We investigate the elastoviscoplastic flow through porous media by numerical simulations. We solve the Navier-Stokes equations combined with the elastoviscoplastic model proposed by Saramito for the stress tensor evolution. In this model, the material behaves as a viscoelastic solid when unyielded, and as a viscoelastic Oldroyd-B fluid for stresses higher than the yield stress. The porous media is made of a symmetric array of cylinders, and we solve the flow in one periodic cell. We find that the solution is time-dependent even at low Reynolds numbers as we observe oscillations in time of the unyielded region especially at high Bingham numbers. The volume of the unyielded region slightly decreases with the Reynolds number and strongly increases with the Bingham number; up to 70% of the total volume is unyielded for the highest Bingham numbers considered here. The flow is mainly shear dominated in the yielded region, while shear and elongational flow are equally distributed in the unyielded region. We compute the relation between the pressure drop and the flow rate in the porous medium and present an empirical closure as function of the Bingham and Reynolds numbers. The apparent permeability, normalized with the case of Newtonian fluids, is shown to be greater than 1 at low Bingham numbers, corresponding to lower pressure drops due to the flow elasticity, and smaller than 1 for high Bingham numbers, indicating larger dissipation in the flow owing to the presence of the yielded regions. Finally we investigate the effect of the Weissenberg number on the distribution of the unyielded regions and on the pressure gradient.
We investigate the chemical dissolution of porous media using a network model in which the system is represented as a series of interconnected pipes with the diameter of each segment increasing in proportion to the local reactant consumption. Moreover, the topology of the network is allowed to change dynamically during the simulation: as the diameters of the eroding pores become comparable with the interpore distances, the pores are joined together thus changing the interconnections within the network. With this model, we investigate different growth regimes in an evolving porous medium, identifying the mechanisms responsible for the emergence of specific patterns. We consider both the random and regular network and study the effect of the network geometry on the patterns. Finally, we consider practically important problem of finding an optimum flow rate that gives a maximum increase in permeability for a given amount of reactant.