Sharif Digital Repository

Most of numerical approaches describing the behavior of miscible displacement in porous media through the solution of advection-dispersion equation suffer from numerical dispersion. Random Walk (RW) method is a good candidate to overcome this problem especially in heterogeneous media. In addition, how treating dispersion coefficient as a tensor might affect the accuracy of RW simulation results is not well understood. Here, a RW model was developed and used for simulating miscible displacement experiments performed on heterogeneous micromodels including single fracture/flow barrier. Dispersion coefficient was treated as a tensor and a hybrid scheme was used for velocity interpolation. The... 

The numerical errors associated with explicit upstream finite difference solutions of two-dimensional advection - Dispersion equation with linear sorption are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation. The numerical truncation errors are defined using Peclet and Courant numbers in the X and Y direction, a sink/source dimensionless number and new Peclet and Courant numbers in the XY plane. The effects of these truncation errors on the explicit solution of a two-dimensional advection-dispersion equation with a first-order reaction or degradation are demonstrated by comparison with an analytical solution in... 

The main objective of this paper is to provide comparative quantitative examinations on the capabilities of two-dimensional horizontal and pseudo-three-dimensional (3D) modelling approaches for simulating spatial and temporal variability of the flow and salinity in Lake Urmia, Iran. The water quality in the lake has been an environmentally important subject partly because this shallow hypersaline aquatic ecosystem is considered to be one of the largest natural habitats of a unique multicellular organism, Artemia urmiana. This brine shrimp is the major food source for many of the protected and rare shorebirds that visit the lake. A.urmiana can grow and survive in certain ranges of salinity,... 

The efficiency of low salinity waterflooding, particularly during tertiary mode injection, is highly controlled by in situ mixing between the stagnant regions holding high salinity water (HSW) and the flowing regions containing low salinity water (LSW) because it impacts directly the electrokinetics of wettability alteration and the time scale of the low salinity effect. This study aims to address the effects of oil polarity and charged rock surfaces on the time scale of mixing and transport under two-phase flow conditions. A systematic series of micromodel experiments were performed. The micromodels were first saturated with high salinity formation brine and oil (both model and crude oil);... 

Gravity drainage is known as the controlling mechanism of oil recovery in naturally fractured reservoirs. The efficiency of this mechanism is controlled by block-to-block interactions through capillary continuity and/or reinfiltration processes. In this study, at first, several free-fall gravity drainage experiments were conducted on a well-designed three-block apparatus and the role of tilt angle, spacers’ permeability, wettability and effective contact area (representing a different status of the block-to-block interactions between matrix blocks) on the recovery efficiency were investigated. Then, an experimental-based numerical model of free-fall gravity drainage process was developed,... 

The present study proposes an exact solution for the release kinetic of a solute from inside a hollow cylindrical polymeric matrix into an infinite medium when the initial concentration of the solute (A) is greater than the solubility limit (Cs). A combination of analytical and numerical methods was used to calculate the solute concentration profile and the release rate. The model was developed for two different conditions including: (1) The release medium was flowing through the hollow cylinder in which the boundary layer may be neglected, and (2) The release medium inside the hollow cylinder was stagnant where the boundary layer needed to be considered. The results indicated that the... 

In this study a new iterative algorithm is developed to evaluate dispersivity in fracture and matrix, distinctly. The novelty of proposed algorithm is using mathematical model of solute transport in fractured porous media coupled with experimental data iteratively. A fractured glass micromodel has been designed to visualize the interaction between fracture and matrix during displacement of n-Decane by n-Octane at constant rate. The similarity between numerical and experimental model has been enhanced by reducing the assumptions which were applied in previous related studies. The iteration is performed on velocity components of solute transport and longitudinal as well as transversal... 

In this study a new iterative algorithm is developed to evaluate dispersivity in fracture and matrix, distinctly. The novelty of proposed algorithm is using mathematical model of solute transport in fractured porous media coupled with experimental data iteratively. A fractured glass micromodel has been designed to visualize the interaction between fracture and matrix during displacement of n-Decane by n-Octane at constant rate. The similarity between numerical and experimental model has been enhanced by reducing the assumptions which were applied in previous related studies. The iteration is performed on velocity components of solute transport and longitudinal as well as transversal... 

Interaction between intracellular dynamics and extracellular matrix (ECM) generally occurred into very thin fragment of moving cell, namely lamellipodia, enables all movable cells to crawl on ECM. In fast-moving cells such as fish Keratocytes, Lamellipodia including most cell area finds a fan-like shape during migration, with a variety of aspect ratio function of fish type. In this work, our purpose is to present a novel and more complete two-dimensional continuum mathematical model of actomyosin-cytosolic two-phase flow of a self-deforming Keratocyte with circular spreaded to steady fan-like shape. In the new approach, in addition to the two-phase flow of the F-actin and cytosol, the... 

Existing closed-form solutions of contaminant transport problems are limited by the mathematically convenient assumption of uniform flow. These solutions cannot be used to investigate contaminant transport in coastal aquifers where seawater intrusion induces a variable velocity field. An adaptation of the Fourier-Galerkin method is introduced to obtain semi-analytical solutions for contaminant transport in a confined coastal aquifer in which the saltwater wedge is in equilibrium with a freshwater discharge flow. Two scenarios dealing with contaminant leakage from the aquifer top surface and contaminant migration from a source at the landward boundary are considered. Robust implementation of... 

Flow movement in unsaturated soil can be expressed by a partial differential equation, named Richards equation. The objective of this study is the finding of an appropriate implicit numerical solution for head based Richards equation. Some of the well known finite difference schemes (fully implicit, Crank Nicolson and Runge-Kutta) have been utilized in this study. In addition, the effects of different approximations of moisture capacity function, convergence criteria and time stepping methods were evaluated. Two different infiltration problems were solved to investigate the performance of different schemes. These problems include of vertical water flow in a wet and very dry soils. The... 

This paper extends the analysis of solute dispersion in electrohydrodynamic flows to the case of band broadening in polyelectrolyte-grafted (soft) capillaries by accounting for the effects of ion partitioning, irreversible catalytic reaction and pulsatile flow actuation. In the Debye-Hückel limit, we present the benchmark solutions of electric potential and velocity distribution pertinent to steady and oscillatory mixed electroosmotic-pressure-driven flows in soft capillaries. Afterwards, the mathematical models of band broadening based on the Taylor-Aris theory and generalized dispersion method are presented to investigate the late-time asymptotic state and all-time evolution of... 

This paper extends the analysis of solute dispersion in electrohydrodynamic flows to the case of band broadening in polyelectrolyte-grafted (soft) capillaries by accounting for the effects of ion partitioning, irreversible catalytic reaction and pulsatile flow actuation. In the Debye-Hückel limit, we present the benchmark solutions of electric potential and velocity distribution pertinent to steady and oscillatory mixed electroosmotic-pressure-driven flows in soft capillaries. Afterwards, the mathematical models of band broadening based on the Taylor-Aris theory and generalized dispersion method are presented to investigate the late-time asymptotic state and all-time evolution of... 

Joining of dissimilar light metals by friction stir welding (FSW) is of interest to reduce weight and fuel consumption in the transport sector. Such coupled metals may need protective surface treatments, e.g. against wear or corrosion, for some applications. In this work, the formation of plasma electrolytic oxidation (PEO) coatings in a silicate-based electrolyte for corrosion protection of FSW AZ31B magnesium alloy-titanium joints has been studied. The joints, if unprotected, may be susceptible to severe galvanic corrosion in chloride-containing environments. The coatings were characterized by scanning electron microscopy, energy-dispersive spectroscopy and X-ray diffraction. Mg2SiO4 and... 

In this paper, a numerical model is developed based on the X-FEM technique to simulate the transport of dense solute in a single fluid phase through the fractured porous media. The governing equation is based on the mass conservation law which is applied to the fluid phase and the solute in both matrix and fracture domain. The integral governing equations of the mass exchange between the fracture and the surrounding matrix is derived. The extended finite element method (X-FEM) is applied by employing appropriate enrichment functions to model the fractured porous domain. The superiority of the X-FEM is that the FE mesh is not necessary to be conformed to the fracture geometry, so the regular... 

DeePore2 is a deep learning workflow for rapid estimation of a wide range of porous material properties based on the binarized micro–tomography images. By combining naturally occurring porous textures we generated 17,700 semi–real 3–D micro–structures of porous geo–materials with size of 2563 voxels and 30 physical properties of each sample are calculated using physical simulations on the corresponding pore network models. Next, a designed feed–forward convolutional neural network (CNN) is trained based on the dataset to estimate several morphological, hydraulic, electrical, and mechanical characteristics of the porous material in a fraction of a second. In order to fine–tune the CNN design,... 

The dispersion of a neutral solute band by electrokinetic flow in polyelectrolyte layer (PEL)-grafted rectangular/slit microchannels is theoretically studied. The flow is assumed to be both steady and fully developed and a first-order irreversible reaction is considered at the wall to account for probable surface adsorption of solutes. Considering low electric potentials, analytical solutions are obtained for electric potential, fluid velocity and solute concentration. Special solutions are also obtained for the case without wall adsorption. To track the dispersion properties of the solute band, the generalized dispersion model is adopted by considering the exchange, the convection and the... 

The majority of literature regarding optimized Latin hypercube sampling (OLHS) is devoted to increasing the efficiency of these sampling strategies through the development of new algorithms based on the combination of innovative space-filling criteria and specialized optimization schemes. However, little attention has been given to the impact of the initial design that is fed into the optimization algorithm, on the efficiency of OLHS strategies. Previous studies, as well as codes developed for OLHS, have relied on one of the following two approaches for the selection of the initial design in OLHS: (1) the use of random points in the hypercube intervals (random LHS), and (2) the use of... 

Practical designs of non-aqueous phase liquids (NAPLs) remediation strategies require reliable modeling of interphase mass transfer to predict the retraction of NAPL during processes such as dissolution. In this work, the dissolution process of NAPL during two-phase flow in heterogeneous porous media is studied using pore-network modeling and micromodel experiments. A new physical-experimental approach is proposed to enhance the prediction of the dissolution process during modeling of interphase mass transfer. In this regard, the normalized average resident solute concentration is evaluated for describing the dissolution process at pore-level. To incorporate the effect of medium... 

Interaction between intracellular dynamics and extracellular matrix (ECM) generally occurred into very thin fragment of moving cell, namely lamellipodia, enables all movable cells to crawl on ECM. In fast-moving cells such as fish Keratocytes, Lamellipodia including most cell area finds a fan-like shape during migration, with a variety of aspect ratio function of fish type. In this work, our purpose is to present a novel and more complete two-dimensional continuum mathematical model of actomyosin-cytosolic two-phase flow of a self-deforming Keratocyte with circular spreaded to steady fan-like shape. In the new approach, in addition to the two-phase flow of the F-actin and cytosol, the...