Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userModelling imbibition processes in heterogeneous porous media
83
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Imbibition is a commonly encountered multiphase problem in various fields, and exact prediction of imbibition processes is a key issue for better understanding capillary flow in heterogeneous porous media. In this work, a numerical framework for describing imbibition processes in porous media with material heterogeneity is proposed to track the moving wetting front with the help of a partially saturated region at the front vicinity. A new interface treatment, named the interface integral method, is developed here, combined with which the proposed numerical model provides a complete framework for imbibition problems. After validation of the current model with existing experimental results of one-dimensional imbibition, simulations on a series of two-dimensional cases are analysed with the presences of multiple porous phases. The simulations presented here not only demonstrate the suitability of the numerical framework on complex domains but also present its feasibility and potential for further engineering applications involving imbibition in heterogeneous media.
rate research
Read More
Imbibition plays a central role in diverse energy, environmental, and industrial processes. In many cases, the medium has multiple parallel strata of different permeabilities; however, how this stratification impacts imbibition is poorly understood. We address this gap in knowledge by directly visualizing forced imbibition in three-dimensional (3D) porous media with two parallel strata. We find that imbibition is spatially heterogeneous: for small capillary number Ca, the wetting fluid preferentially invades the fine stratum, while for Ca above a threshold value, the fluid instead preferentially invades the coarse stratum. This threshold value depends on the medium geometry, the fluid properties, and the presence of residual wetting films in the pore space. These findings are well described by a linear stability analysis that incorporates crossflow between the strata. Thus, our work provides quantitative guidelines for predicting and controlling flow in stratified porous media.
Imbibition, the displacement of a nonwetting fluid by a wetting fluid, plays a central role in diverse energy, environmental, and industrial processes. While this process is typically studied in homogeneous porous media with uniform permeabilities, in many cases, the media have multiple parallel strata of different permeabilities. How such stratification impacts the fluid dynamics of imbibition, as well as the fluid saturation after the wetting fluid breaks through to the end of a given medium, is poorly understood. We address this gap in knowledge by developing an analytical model of imbibition in a porous medium with two parallel strata, combined with a pore network model that explicitly describes fluid crossflow between the strata. By numerically solving these models, we examine the fluid dynamics and fluid saturation left after breakthrough. We find that the breakthrough saturation of nonwetting fluid is minimized when the imposed capillary number Ca is tuned to a value Ca$^*$ that depends on both the structure of the medium and the viscosity ratio between the two fluids. Our results thus provide quantitative guidelines for predicting and controlling flow in stratified porous media, with implications for water remediation, oil/gas recovery, and applications requiring moisture management in diverse materials.
Porous media with hierarchical structures are commonly encountered in both natural and synthetic materials, e.g., fractured rock formations, porous electrodes and fibrous materials, which generally consist of two or more distinguishable levels of pore structure with different characteristic lengths. The multiphase flow behaviours in hierarchical porous media have remained elusive. In this study, we investigate the influences of hierarchical structures in porous media on the dynamics of immiscible fingering during fluid-fluid displacement. By conducting a series of numerical simulations, we found that the immiscible fingering can be suppressed due to the existence of secondary porous structures. To characterise the fingering dynamics in hierarchical porous media, a phase diagram is constructed by introducing a scaling parameter, i.e., the ratio of time scales considering the combined effect of characteristic pore sizes and wettability. The findings present in this work provide a basis for further research on the application of hierarchical porous media for controlling immiscible fingerings.
Electrical conductivity is an inherent property of a hydrophobic porous media (HPM) and has critical applications. This research aims to provide a solution for predicting the electrical conductivity of nanoscale HPM with heterogeneous pore structure. Molecular dynamics (MD) simulations are compared with the modified Poisson-Boltzmann (MPB) model for understanding ionic charge density distributions in nanopores. The effective medium approximation (EMA) participates in calculating the effective conductance and conductivity of the nanoscale HPM. The results show that the surface charge density affects the ionic density profiles in the hydrophobic nanopores. As the pore size increases, the conductance increases. As the molarity of the aqueous electrolyte solution (AES) decreases, the conductance decreases. A phenomenon related to the conductance saturation occurred when the molarity of AES is very low. The effective conductance of an HPM increase as the coordination number increases. Finally, based on the calculated effective conductance and the heterogeneous pore structure parameters, the electrical conductivity of a nanoscale HPM is calculated.
Flow and multicomponent reactive transport in saturated/unsaturated porous media are modeled by ensembles of computational particles moving on regular lattices according to specific random walk rules. The occupation number of the lattice sites is updated with a global random walk (GRW) procedure which simulates the evolution of the ensemble with computational costs comparable to those for a single random walk simulation in sequential procedures. To cope with the nonlinearity and the degeneracy of the Richards equation the GRW flow solver uses linearization techniques similar to the $L$-scheme developed in finite element/volume approaches. Numerical schemes for reactive transport, coupled with the flow solver via numerical solutions for saturation and water flux, are implemented in splitting procedures. Diffusion-advection steps are solved by GRW algorithms using either biased or unbiased random walk probabilities. Since the number of particles in GRW simulations can be as large as the number of molecules involved in chemical reactions, one avoids the cumbersome problem of rescaling particle densities to approximate concentrations. Reaction steps are therefore formulated in terms of concentrations, as in deterministic approaches. The numerical convergence of the new schemes is demonstrated by comparisons with manufactured analytical solutions. Coupled flow and reactive transport problems of contaminant biodegradation described by the Monod model are further solved and the influence of flow nonlinearity/degeneracy and of the spatial heterogeneity of the medium is investigated numerically.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
