No Arabic abstract
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.
Immiscible fluid-fluid displacement in porous media is of great importance in many engineering applications, such as enhanced oil recovery, agricultural irrigation, and geologic CO2 storage. Fingering phenomena, induced by the interface instability, are commonly encountered during displacement processes and somehow detrimental since such hydrodynamic instabilities can significantly reduce displacement efficiency. In this study, we report a possible adjustment in pore geometry which aims to suppress the capillary fingering in porous media with hierarchical structures. Through pore-scale simulations and theoretical analysis, we demonstrate and quantify combined effects of wettability and hierarchical geometry on displacement patterns, showing a transition from fingering to compact mode. Our results suggest that with a higher porosity of the 2nd-order porous structure, the displacement can keep compact across a wider range of wettability conditions. Combined with our previous work on viscous fingering in such media, we can provide a complete insight into the fluid-fluid displacement control in hierarchical porous media, across a wide range of flow conditions from capillary- to viscous-dominated modes. The conclusions of this work can benefit the design of microfluidic devices, as well as tailoring porous media for better fluid displacement efficiency at the field scale.
Recent experiments demonstrate how a soluble body placed in a fluid spontaneously forms a dissolution pinnacle -- a slender, upward pointing shape that resembles naturally occurring karst pinnacles found in stone forests. This unique shape results from the interplay between interface motion and the natural convective flows driven by the descent of relatively heavy solute. Previous investigations suggest these structures to be associated with shock-formation in the underlying evolution equations, with the regularizing Gibbs-Thomson effect required for finite tip curvature. Here, we find a class of exact solutions that act as attractors for the shape dynamics in two and three dimensions. Intriguingly, the solutions exhibit large but finite tip curvature without any regularization, and they agree remarkably well with experimental measurements. The relationship between the dimensions of the initial shape and the final state of dissolution may offer a principle for estimating the age and environmental conditions of geological structures.
A reactive fluid dissolving the surface of a uniform fracture will trigger an instability in the dissolution front, leading to spontaneous formation of pronounced well-spaced channels in the surrounding rock matrix. Although the underlying mechanism is similar to the wormhole instability in porous rocks there are significant differences in the physics, due to the absence of a steadily propagating reaction front. In previous work we have described the geophysical implications of this instability in regard to the formation of long conduits in soluble rocks. Here we describe a more general linear stability analysis, including axial diffusion, transport limited dissolution, non-linear kinetics, and a finite length system.
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.
The dissolution of porous materials in a flow field shapes the morphologies of many geologic landscapes. Identifying the dissolution front, the interface between the reactive and the unreactive regions in a dissolving medium, is a prerequisite for studying dissolution kinetics. Despite its fundamental importance, the dynamics of a dissolution front in an evolving natural microstructure has never been reported. Here we show an unexpected spontaneous migration of the dissolution front against the pressure gradient of a flow field. This retraction stems from the infiltration instability induced surface generation, which can lead to a reactive surface dramatically greater than the ex situ geometric surface. The results are supported by a very good agreement between observations made with real time X-ray imaging and simulations based on static images of a rock determined by nanoCT. They both show that the in situ specific surface area of natural porous media is dependent on the flow field and reflects a balancing between surface generation and destruction. The reported dynamics challenge many long-held understanding of water-rock interactions and shed light on reconciling the discrepancies between field and laboratory measurements of reaction kinetics.