Sharif Digital Repository

In this paper, a numerical model is developed for the assessment of carbon dioxide transport through naturally fractured cap-rocks during CO2 sequestration in underground aquifers. The cap-rock contains two types of fracture with different length scales: micro-cracks (fissures) and macro-cracks (faults). The effect of micro-cracks is incorporated implicitly by modifying the intrinsic permeability tensor of porous matrix, while the macro-cracks are modeled explicitly using the extended finite element method (X-FEM). The fractured porous medium is decomposed into the porous matrix and fracture domain, which are occupied with two immiscible fluid phases, water and CO2. The flow inside the... 

Fluid flow in porous media is affected by surface characteristics such as roughness and topography. In this work, to simulate the surface of natural porous structures in transparent interconnected media like micromodels, various degrees of roughness have been artificially created on flat glass substrates via different methods of laser ablation, cream etching, combination of laser ablation and cream etching, and hydrofluoric acid (HF) etching. The obtained surfaces by each method were characterized in detail via field emission scanning electron microscopy (FESEM), atomic force microscopy (AFM), energy-dispersive X-ray spectroscopy (EDX/EDS), and surface profilometry. The impact of high... 

In this study, in order to investigate the effect of the underlying pore-scale processes on continuum scale simulations of porous media dissolution, we improve the standard Lattice Boltzmann method using Quadtree grid refinement approach to simulate fluid flow and reactive transport through large domain sizes. Our results have shown considerable computational improvements up to 80% in simulation time together with increased numerical accuracy. The results and the added value of the new approach are discussed using comparison of our model with the conventional LBM. Moreover, we have applied a systematic analysis by increasing complexity levels and starting from fluid flow and continuing with... 

This paper presents a numerical multiscale formulation for analysis of the transient heat and fluid flow in deformable heterogeneous porous media. Due to the heterogeneity of the media, the direct numerical simulation of the micro-structures leads to high computational costs. Hence, the multi-scale method can provide an efficient computational procedure. To this end, the first-order computational homogenization is adopted for two-scale simulation of THM problems. The governing equations of the problem contain a stress equilibrium equation, a mass continuity equation and an advection–diffusion equation in a fully coupled manner. Accordingly, the proper virtual power relations are defined as a... 

The present article aimed to examine the pore-scale examination of a fluid capability to block pores and fractures deduced by return permeability and the blocking ability of a novel lightweight Colloidal Gas Aphron Nano-Fluid (CGANF) in heterogeneous fractured/un-fractured porous medium. Silica (S) and fumed silica (FS) nanoparticles, Sodium Dodecyl Sulfate (SDS) anionic surfactant, and Xanthan Gum (XG) bio-polymer were employed for following up CGANF bubbles' movements in incongruous micro-models’ inclusive crack. The transparent nature of the glass micro-model provides an intuition into the lenticular behavior of CGANF bubbles in the heterogeneous porous media. Differential pressure... 

In this study, the heat transfer and fluid flow characteristics of a three-dimensional microchannel that is partially filled with a layer of porous medium at its bottom solid wall is investigated. The microchannel is consisted of a clear fluid flow region, solid walls and a porous layer that is attached to its solid bottom wall. A constant heat flux is applied to the bottom wall of the microchannel. Darcy-Brinkman-Forchheimer model is used to simulate the fluid flow inside the porous medium. The novelty of this work is to investigate thoroughly and precisely the effect of using of porous layer configuration in MCHSs on hydraulic and thermal performances. The effect of porous layer thickness,... 

Induced Polarization (IP) is a non-intrusive geophysical method to monitor Dense Non-Aqueous Phase Liquid (DNAPL) contamination and remediation processes underground. In this study, an advanced numerical code simulating DNAPL flow and complex electrical resistivity is presented. The model was validated against existing IP results and image measurements that were carried out previously in a series of 2D tank experiment. Multiphase flow modeling in porous media is coupled with electrical current modeling to simulate the process of DNAPL migration and the associated IP response. This brings a broader view of the contamination in space and time compared to surface and borehole measurements,... 

It is believed that it is going to be a sizeable mismatch between supply and demand when it comes to renewable resources. Lately, researchers are on course to compensate for the unpredictabilityof such resources by the employment of phase change materials (PCMs). Having multiple advantages, PCMs generally suffer from inadequate thermal conductivity which causes prolonged transition procedures. To tackle this issue, this study is fixated on two parameterswhich are linked to fins addition and porous media incorporation in a melting process within a triple concentric tube heat exchanger (TCTHX). The results provided by multiple cases underlined the significance of natural convection in the bare... 

The process of fluid flow displacement in porous media has recently gained great prominence owing to its widespread usage in a variety of industries, especially in the case of pore scale investigations. Although, many studies have been conducted to address pore-scale investigations in both modeling and experimental approaches, the role of interfacial tension and contact angle on pore-scale phenomena is less focused. In this work, direct pore-scale modeling was used to precisely examine the effect of interfacial tension and contact angle on the fluid flow at the microscale. Also, several pore-scale mechanisms, including Haines jump and dynamic breakup mechanisms, were observed. Therefore, the... 

One of the most popular models that has been applied to predict the fluid velocity inside the fracture with impermeable walls is the cubic law. It highlights that the mean flux along the fracture is proportional to the cubic of fracture aperture. However, for a fractured porous medium, the normal and tangential interface conditions between the fracture and porous matrix can change the velocity profile inside the fracture. In this paper, a correction factor is introduced for flow equation along the fracture by imposing the continuity of normal and tangential components of velocity at the interface between the fracture and porous matrix. As a result, the mean velocity inside the fracture... 

In this article, a computational model is presented for the analysis of coupled thermo-hydro-mechanical process with phase change (evaporation/condensation) in fractured porous media in order to model multiphase fluid flows, heat transfer, and discontinuous deformation by employing the extended finite element method. The ideal gas law and Dalton's law are employed to consider vapor and dry air as miscible gases. To take into account the phase change, latent heat and specific vapor enthalpy are incorporated into the physical model. The set of governing equations consists of linear momentum for the solid-phase, energy balance equation and mass conservation equations of water species (liquid... 

In this paper, a computational modeling tool is developed for fully coupled multiphase flow in deformable heterogeneous porous medium that consists of complex and non-uniform micro-structures using the dual continuum scales based on the computational homogenization approach. The first-order homogenization technique is employed to perform the multi-scale analysis. The governing equations of two-phase flow of immiscible fluids, including an equilibrium equation and two mass continuity equations, are considered based on the appropriate main variables. According to the well-known Hill–Mandel principle of macro-homogeneity, the proper energy types are defined instead of conventional stress power... 

Oil trapping behavior during the pre-flush stage is critically important to evaluate the effectiveness of matrix acidizing for the oil well stimulation. In this study, the visco-capillary behavior of the two-phase flow in the pore-scale is analyzed to investigate the influence of wetting properties for a natural rock sample. A two-dimensional model, based on Cahn-Hilliard phase-field and Navier-Stokes equations, was established and solved using the finite element method. A stability phase diagram for log capillary number (Ca)-log viscosity ratio (M) was constructed and then compared with the reported experimental works. The maximum and minimum ranges of capillary number and viscosity ratio... 

A coupled statistical and conditional generative adversarial neural network is used for 3D reconstruction of both homogeneous and heterogeneous porous media from a single two-dimensional image. A statistical approach feeds the deep network with conditional data, and then the reconstruction is trained on a deep generative network. The conditional nature of the generative model helps in network stability and convergence which has been optimized through a gradient-descent-based optimization method. Moreover, this coupled approach allows the reconstruction of heterogeneous samples, a critical and serious challenge in conventional reconstruction methods. The main contribution of this work is to... 

In this study, two classes of advanced finite volume schemes, including Multi-Point Flux Approximation (MPFA) and Dual Discrete Finite Volume (DDFV) method, have been employed in conjunction with the finite element (FE) geomechanics simulator to model the coupled fluid-solid system. Fully saturated porous media with poroelastic behavior, random field permeability and elastic modulus are considered as parameters. The performance of the proposed hydro-mechanical models, including MPFA O-FEM and DDFV-FEM, is examined through different test cases. First, the models are validated and compared with the closed-form solutions in the homogeneous domain. Second, the methods' stability and convergence... 

The formation of liquid bridges can maintain capillary continuity between matrix blocks during gas/oil gravity drainage in fractured reservoirs. A travelling oil drop draining into a fracture either forms a liquid bridge or breaks into detached drops. However, the different characteristics of a travelling drop during its elongation and required conditions for transformation into a liquid bridge are not well described in the published literature. In this work, a one-dimensional model based on slender-drop theory is employed that holds gravity, viscosity, and surface tension forces but ignores inertia. This model, together with Young-Laplace equation, gives the fracture capillary pressure.... 

The formation of liquid bridge is pertinent to many fields including seepage into underground fractured rocks as an environmental issue and capillary continuity between matrix blocks which controls oil recovery in fractured reservoirs. Evaporation from the surface of liquid bridge into the surrounding gas could affect the stability of liquid bridge and fracture capillary pressure, which is not well discussed in the available literatures. In this research, by the aid of analogy between the diffusive flux and electrostatic potential, a new model for predicting evaporation rate from a liquid bridge inside a horizontal fracture is presented. The proposed model is then coupled with Young-Laplace... 

Low salinity waterflooding (LSWF) is a process in which by lowering the ionic strength and/or manipulation of the composition of the injection water, the long term equilibrium in oil/brine/rock system is disturbed to reach a new state of equilibrium through which the oil production will be enhanced due to fluid/fluid and/or rock/fluid interactions. In spite of recent advances in the simulation of the LSWF at core scale and beyond, there are very few works that have modelled and simulated this process at the pore scale specially using direct numerical simulation (DNS). As a result the effects of wettability alteration and/or Interfacial Tension (IFT) change on the distribution of the phases... 

Static and dynamic properties of porous media are highly dependent on its internal geometry. CT scan images are generally used to characterize porous media geometry. Direct simulation of fluid flow on CT scan images is possible but considerably time-consuming. In this study, a new method was developed for extracting a simplified representation known as “pore network model” by utilizing a rigorous algebraic-analytical method. By using a moving frame in the 3D matrix of the CT scan image and stepwise identifying-removing of image components, running time for a 4003 voxels sample in a typical computer system decreased to less than 350 s. The identification of throats was based on a new... 

Natural convection in a porous enclosure in the presence of thermal dispersion is investigated. The Fourier–Galerkin (FG) spectral element method is adapted to solve the coupled equations of Darcy’s flow and heat transfer with a full velocity-dependent dispersion tensor, employing the stream function formulation. A sound implementation of the FG method is developed to obtain accurate solutions within affordable computational costs. In the spectral space, the stream function is expressed analytically in terms of temperature, and the spectral system is solved using temperature as the primary unknown. The FG method is compared to finite element solutions obtained using an in-house code...