Sharif Digital Repository

In theories of multiphase flow through porous media, capillary pressure-saturation and relative permeability-saturation curves are assumed to be intrinsic properties of the medium. Moreover, relative permeability is assumed to be a scalar property. However, numerous theoretical and experimental works have shown that these basic assumptions may not be valid. For example, relative permeability is known to be affected by the flow velocity (or pressure gradient) at which the measurements are carried out. In this article, it is suggested that the nonuniqueness of capillary pressure-relative permeability-saturation relationships is due to the presence of microheterogeneities within a laboratory... 

The objective of this paper is to provide optimization of falling film Li/Br solution on a horizontal single tube based on minimization of entropy generation. Flow regime is considered to be laminar, the effect of boiling has been ignored and wall temperature is constant. Velocity, temperature and concentration distributions are numerically determined and dimensionless correlations are obtained for predicting the average heat transfer coefficient and average evaporation factor on the horizontal tube. Thermodynamic imperfection due to passing lithium bromide solution is attributed to non-isothermal heat transfer; fluid flow friction and mass transfer irreversibility. Scale analysis shows that... 

Computational simulations of two-phase flow in porous media are used to investigate the feasibility of replacing a porous medium containing heterogeneities with an equivalent homogeneous medium. Simulations are performed for the case of infiltration of a dense nonaqueous phase liquid (DNAPL) in a water-saturated, heterogeneous porous medium. For two specific porous media, with periodic and rather simple heterogeneity patterns, the existence of a representative elementary volume (REV) is studied. Upscaled intrinsic permeabilities and upscaled nonlinear constitutive relationships for two-phase flow systems are numerically calculated and the effects of heterogeneities are evaluated. Upscaled... 

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... 

Water can be pumped in nanochannels by limiting it between the surfaces with different hydrophobicities and exerting a spinning electric field. The asymmetrical hydrophobicity combined with the spinning electric field and the fact that the water molecules have a dipole moment create a situation in which the angular momentum of water molecules is transformed into a linear momentum and the water is pumped into the nanochannel. The hydrophobicity of the surfaces can be manipulated by using nanostructures to reduce friction. In this study, two types of nanostructures have been used which are hierarchical nanostructures and polymer nanostructures made of Poly(N-isopropylacrylamide). The walls of... 

Minimization of drop contact time is so important and critical for applications such as self-cleaning and anti-corrosion. In recent years, surface acoustic waves are presented as a powerful method for the manipulation of the drops. In this manuscript, a numerical study of drop shedding under the effect of acoustic waves is presented, which may have potential use in the anti-icing systems. Therefore, the effects of different parameters, such as acoustic wave frequency, amplitude, and direction of the wave on the water removal, are investigated. For this purpose, a color gradient lattice Boltzmann method (LBM) is developed and used in these simulations. The acoustic actuator effect is added as... 

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,... 

In this paper, a computational technique is presented for the isothermal and non-isothermal water injection into naturally fractured oil reservoirs. A remarkable number of naturally fractured reservoirs contain relatively heavy oils that could not be extracted economically; hence, the thermal recovery methods are extensively used for such reservoirs. In this study, the effectiveness of hot water injection over cold (isothermal) water injection in oil production is quantified. The influence of long and short fractures and their alignments on oil recovery are discussed. To this end, a 2D model for two-phase fluid flow and heat transfer is presented. The medium is assumed to be partially... 

In this paper, the dissolution and mobilization of non-aqueous phase liquid (NAPL) blobs in the Surfactant-Enhanced Aquifer Remediation (SEAR) process were upscaled using dynamic pore network modeling (PNM) of three-dimensional and unstructured networks. We considered corner flow and micro-flow mechanisms including snap-off and piston-like movement for two-phase flow. Moreover, NAPL entrapment and remobilization were evaluated using force analysis to develop the capillary desaturation curve (CDC) and predict the onset of remobilization. The corner diffusion mechanism was also applied in the modeling of interphase mass transfer to represent NAPL dissolution as the dominant mass transfer... 

Accurate determination of relative permeability curves and their hysteresis is vital for reliable prediction of the performance of oil and gas reservoirs under enhanced recovery processes. Two out of the three available approaches to simulate three-phase relative permeability hysteresis are based on two-phase hysteresis. A few options (e.g., Carlson, Killough and Jargon models) are available in commercial reservoir simulators to account for hysteresis in relative permeability curves under two-phase flow. These models are based on the assumptions of water-wet state of the rocks, while most of the reservoir rocks are mixed-wet. As a result the aim of the present work is to evaluate the... 

Steam/water stratified flow would occur in transient condition (e.g. LOCA) in light water Nuclear Power Plants (NPPs). Due to high gradient of flow characteristics at the interface of steam/water flow, the prediction of flow characteristics (e.g. temperature, pressure, velocity, and Turbulent Kinetic Energy (TKE)) requires further attention and special interfacial models. Also, accurate simulation of these mentioned characteristics needs fine spatial mesh and very small time steps based on Computational Fluid Dynamics (CFD) standard criteria. In order to reduce the computational cost, the combination of thermal–hydraulic modelling and soft computing is considered as a new strategy in this... 

The dynamics of droplet detachment from a surface is a fundamental topic studied in coating engineering, fluid mechanics, and subsurface engineering applications. This topic has direct relevance to wettability alteration using the modified ionic composition of water in contact with oil droplet, low salinity waterflooding (LSWF). Previous experimental studies of LSWF have shown a very long timescale in wettability alteration which cannot be explained using bulk diffusion coefficient. In the present study, we address both the time scale of detachment, as well as the impact of buoyancy and interfacial forces (referred to as Bond number) on droplet detachment by proposing an advanced...