Sharif Digital Repository

Spontaneous imbibition, the capillary-driven process of displacing the nonwetting phase by the wetting phase in porous media, is of great importance in oil/gas recovery from matrix blocks of fractured reservoirs. The question of how properly scaling up the recovery by counter-current spontaneous imbibition has been the subject of extensive research over decades, and numerous scaling equations have been proposed. As a convention, the scaling equations are usually defined analytically by relating the early time squared recovery to squared pore volume. We show this convention does not apply to common scaling practices and, if used, causes nontrivial scatter in the scaling plots. We explain that... 

In this article, the problem of combined forced and free convection in vertical porous and regular channels for both regular fluids and nanofluids has been solved using the CFD technique in the entrance regions of momentum and heat transfer taking into account the influences of viscous heating and inertial force. In this regard, various types of viscous dissipation models reported in the literature such as the Darcy model, the power of the drag force model, and the clear fluid-compatible model were applied. In the case of nanofluid flow, both the Brownian and thermophoresis molecular transfer mechanisms were considered. The dimensionless distributions of velocity, temperature, and the volume... 

Rubble-mound breakwaters are common marine structures that provide a safe area for human coastal activities. The stability of these structures against sea-waves requires their seaward slope to be protected by an armor layer consisting of natural rock or concrete units. To provide a safe breakwater, it is reasonable to establish a relation between the exerted wave loads and the stability of the armor units. However, up to now, the empirical design equations, derived from model tests, relate wave parameters to armor weight, and keeps the effect of wave loads in a black box. In this paper, a new approach, based on numerically-derived wave loads on the armor, is presented to evaluate the... 

The fresh groundwater lenses (FGLs) of small islands can be highly vulnerable to climate change impacts, including sea-level rise (SLR). Many real cases of atoll or sandy islands involve two-layer hydrogeological conceptualizations. In this paper, the influential factors that affect FGLs in two-layer small islands subject to SLR are investigated. An analytical solution describing FGLs in circular islands, composed of two geological layers, is developed for the simplified case of steady-state and sharp-interface conditions. An application of the developed model is demonstrated to estimate the FGL thickness of some real-world islands by comparison with existing FGL thickness data. Furthermore,... 

A precise semi-empirical correlation for the calculation of interfacial tension (IFT) between the carbon dioxide and paraffin group to be used in an enhanced oil recovery process and the chemical industry is introduced. Genetic programming and dimension analysis (GPDA) are combined to create a correlation for the calculation of the equilibrium interfacial tension of the carbon dioxide and paraffin group, based on the explicit functionality of the pressure and temperature. The parameters of the correlation consist of critical temperature, critical pressure, density of paraffin at normal temperature, and diffusion coefficients. The pool of experimental data for developing the correlation... 

Spontaneous water imbibition into the matrix blocks is known as the main mechanism for increased oil recovery from naturally fractured oil reservoirs. The rate of oil recovery and its ultimate value is mostly affected by wettability of the rocks and their pore structure. Oil viscosity also greatly influences the rate of oil recovery. A novel experimental model was utilized to study the imbibition mechanism under different wettability conditions. Matrix blocks made from different grain types and size distributions of glass beads were saturated with two different types of synthetic oil, to mimic the oil-saturated matrixes. The wetting characteristic of the models used in this study were... 

Using hybrid energy storage system is a method for increasing the storage capability of solar thermal energy. If multiple energy storage devices with complementary performance characteristics are used together, the resulting system will be a 'Hybrid Energy Storage System'. In other words, a Hybrid Energy Storage System (HESS) has several media available for storage at any time. In this way, increase in storable energy is obtained without increasing collectors' area. When there are more than one storage mediums, the system should be able to choose the best medium for storing energy according to the conditions. In the previous works, an optimizer program was used to find the proper medium... 

Because of technological restrictions, it is actually impossible to fabricate smooth microchannels. Therefore, exploring the roughness effects on the flow characteristics at microscale is of great importance for scientific communities. The present investigation deals with the effects of corrugated roughness on the fully developed slip flow forced convection in micropipes. The governing equations subject to first-order slip boundary conditions are solved by means of the straightforward perturbation method. Closed-form expressions are obtained for the dimensionless velocity and temperature distributions, for the friction coefficient and pressure drop, and finally for the Nusselt number. The... 

The axisymmetric problem of a concentric set of energetically consistent annular and penny-shaped cracks in an infinite piezoelectric body subjected to uniform far-field electromechanical loading is addressed. With the aid of a robust innovated technique, the pertinent four-part mixed boundary value problem (MBVP) is reduced to a decoupled Fredholm integral equation of the second kind. The results of two limiting cases of a single penny-shaped crack and a single annular crack are recovered. The contour plots of dimensionless intensity factors (IFs) at each crack front provide the stress and electric displacement intensity factors (SIFs and EDIFs, respectively) for all combination of crack... 

This paper presented a numerical study that predicts critical mass flow rate, pressure, vapor quality, and void fraction along a very long tube with small diameter or capillary tub under critical condition by the drift flux model. Capillary tubes are simple expansion devices and are necessary to design and optimization of refrigeration systems. Using dimensional analysis by Buckingham's π theory, some generalized correlations are proposed for prediction of flow parameters as functions of flow properties and tube sizes under various critical conditions. This study is performed under the inlet pressure in the range of 0.8 ≤ pin ≤ 1.5Mpa, subcooling temperature between 0 ≤ ΔTsub ≤10 °C. The... 

The drainage of oil under gravity forces is an efficient method because it can farther reduce the remaining oil saturation to below that obtained after water flooding. This paper describes a series of visual experiments under forced gas invasion with special attention to the effects of wettability. From oil production history and image analysis, we examine a threshold criterion for displacement stability that is consistent with the results of gradient percolation theory. The effect of the destabilized front velocity on oil recovery and residual saturation is investigated for both wettability conditions. Different recovery rates occur with different fluid morphologies, which depend on the... 

In order to relate the wetting properties at the pore scale to the macroscale prevailing forces, a series of experiments was performed in vertical porous media under forced gas invasion at various wettability conditions with partially spreading oil. To describe the dynamics of oil recovery in a three-phase flow condition, the downward gas flood experiments were continued by water injection from the bottom. Experimental results obtained in situations where the magnitudes of viscous, capillary, and gravity forces are comparable. We study the transition from flow configurations where the interface is stable with respect to viscous instability to flow configurations where viscous fingering... 

In this work, an experimental and theoretical study has been conducted to investigate the permeability reduction due to CaSO4 scale deposition in packed column porous media. Permeability reduction by calcium sulphate deposition follows a systematic trend considering various important parameters that are affected in this complex process. Hence, a novel dimensionless model has been proposed for the prediction of permeability reduction rate with high accuracy. The developed model is based on the data obtained from glass bead and carbonate grain packed column at low pressure. The proposed model was validated with Berea sandstone cores data at high pressure (100-20,678 kPa), various flow rates... 

Including gravity and wettability effects, a full analytical solution for the frontal flow period for 1D counter-current spontaneous imbibition of a wetting phase into a porous medium saturated initially with non-wetting phase at initial wetting phase saturation is presented. The analytical solution applicable for liquid-liquid and liquid-gas systems is essentially valid for the cases when the gravity forces are relatively large and before the wetting phase front hits the no-flow boundary in the capillary-dominated regime. The new analytical solution free of any arbitrary parameters can also be utilized for predicting non-wetting phase recovery by spontaneous imbibition. In addition, a new... 

In this study, the Lattice Boltzmann Method (LBM) is used to investigate the deformation of two droplets within microfluidic T-junctions (MFTD). In order to increase the accuracy the two immiscible fluids are modeled using the He-Chen-Zhang model. First, this model is applied to ensure that the surface tension effect existing between the droplets and the continuous fluid is properly implemented in the model. Then the collision and merging of the two droplets within the intersection of a T-shaped microchannel is investigated. For generating droplet formation the effects of relevant dimensionless parameters such as the Reynolds, the Weber numbers as well as a collision parameter affecting the... 

The morphology of symmetric and asymmetric membranes obtained by precipitation can be rationalized by thermodynamic and kinetic parameters. Such parameters which the former relates to thermodynamic effect and the latter controls the time of separation of the cast solution film immersing in the non-solvent bath from the glass plate were calculated as dimensionless parameters for further analysis. Phase diagram was used to consider the role of thermodynamic. The kinetic properties of the membranes could be investigated by the rate of mass transfer between solvent and non-solvent. By adding additives to the polymer solution the desire for becoming two phase and also viscosity of the solution... 

In this research, effective length of one-dimensional combustion in a dilute monopropellant spray, constant area and fixed volume chamber is analytically predicted. A new evaporation rate in the form of d k+1 relation is introduced. In the case of controlling vaporization by radiative heat transfer, k is equal to zero, and when molecular processes control the vaporization, k will be equal to one and in some cases vaporization data need the value of k greater than one to fit properly to related equation. Development of this approach can be used in the design of combustion chambers with optimum length and with using vaporization rate of R = R0〈r〉 0 k/〈r〉k. Spray equation and distribution... 

In this research, we have investigated the three-dimensional elastic collision of two balls, based on friction in the tangential plane. Our aim is to offer analytical closed form relations for post collision parameters such as linear and angular velocities, collision time and tangential and normal impulse in three dimensions. To simplify the problem, stick regime is not considered. In other words, balls have a low tangential coefficient of restitution. Sliding, sliding then rolling, and rolling at the beginning of contact are three cases that can occur during impact which have been considered in our research. The normal interaction force is described by the Hertz contact force and... 

This paper presents a new procedure for determining the fracture toughness of rock-like specimens using the diametric compression test with the center-cracked horseshoe disk (CCHD) method. Using finite element analysis, a dimensionless stress intensity factor was obtained and a formula was rendered for determining mode I fracture toughness. To evaluate the accuracy of the measurement results produced by the CCHD method, fracture toughness experiments were conducted on the same rock-like material using the notched Brazilian disk (NBD) method. The CCHD tests were simulated using a two-dimensional particle flow code for validation of the experimental results, and a great agreement between the...