Sharif Digital Repository

Fluid flow modeling in porous media has many applications in waste treatment, hydrology and petroleum engineering. In any geological model, flow behavior is controlled by multiple properties. These properties must be known in advance of common flow simulations. When uncertainties are present, deterministic modeling often produces poor results. Percolation and Random Walk (RW) methods have recently been used in flow modeling. Their stochastic basis is useful in dealing with uncertainty problems. They are also useful in finding the relationship between porous media descriptions and flow behavior. This paper employs a simple methodology based on random walk and percolation techniques. The... 

During miscible displacements in fractured porous media, one of the most important factors that plays a significant role in oil production is the matrix-fracture interaction. In this work, a series of hydrocarbon injection experiments have been performed on a fractured glass micromodel that was designed specifically to study matrix-fracture interaction. A high quality image analysis method was used to determine the fluid flow behavior, solvent front movement, and viscous fingering associated with solvent movement in matrix and fractures. Observations showed that in the case of unit viscosity ratio, the injection rate increased the slope of recovery curve and consequently improved the final... 

During the past two decades, the lattice Boltzmann (LB) method has been introduced as a class of computational fluid dynamic methods for fluid flow simulations. In this method, instead of solving the Navier Stocks equation, the Boltzmann equation is solved to simulate the flow of a fluid. This method was originally developed based on uniform grids. However, in order to model complex geometries such as porous media, it can be very slow in comparison with other techniques such as finite differences and finite elements. To eliminate this limitation, a number of studies have aimed to formulate the lattice Boltzmann on the unstructured grids. This paper deals with simulating fluid flow through a... 

Finger initiation/development at fluid-fluid interface during miscible floods can cause poor displacement efficiency, which is undesirable in enhanced oil recovery processes. In this work, a series of hydrocarbon injection experiments performed on 5-spot glass micromodels that were initially saturated with the heavy crude oil. The fractured micromodels with different fracture geometrical characteristics were used in the tests. High quality image analysis was applied to determine the fluid flow behavior, solvent front movement, and viscous fingering associated with solvent movement in matrix and fractures. Observations showed that higher solvent dispersion in the fractures rather than matrix... 

A brass vortextube is used to carry out a series of experiments. The main objective of the present research is to investigate the separation performances of a vortex tube (VT) for a hydrocarbon mixture. Examination is also applied to study the effects of nozzle intakes number and cold fraction on the gas species separation at specific inlet pressure 236.37 kPa in a VT with two gas mixtures (LPG as a hydrocarbon mixture and LPG-N2). A two-dimensional computational fluid dynamic (CFD) model simulation of a VT is presented. CFD code after validation is also applied to investigate the role of cold fraction and nozzle intakes number on the gas species separation. The highly rotating flow field... 

Body force modeling is studied in the Generalized Newtonian (GN) fluid flow simulation using a single relaxation time lattice Boltzmann (LB) method. First, in a shear thickening Poiseuille flow, the necessity for studying body force modeling in the LB method is explained. Then, a parametric unified framework is constructed for the first time which is composed of a parametric LB model and its associated macroscopic dual equations in both steady state and transient simulations. This unified framework is used to compare the macroscopic behavior of different forcing models. Besides, using this unified framework, a new forcing model for steady state simulations is devised. Finally, by solving a... 

Aerodynamically alleviated marine vehicle (AAMV) is a high speed craft equipped with aerodynamic surfaces that operating in ground effect zone provides this craft with the ability to achieve much higher cruising speeds. Reducing the take-off mode of an AAMV is highly desirable. Additionally, it is seen where there is a considerable reserve thrust take-off can occur in the lower get-away speeds that shorten the take-off run and, therefore, is favourable. Accordingly, in this study an attempt has been made to develop a nonlinear mathematical model for an AAMV to simulate accelerations in take-off and landing phases, using semi-empirical equations mainly proposed for mono-hull high-speed craft,... 

Respiratory drug delivery has been under the research spotlight for the past few decades, mainly due to the high incidence of pulmonary diseases and the fact that this type of delivery offers the highest efficiency for treatment. Despite its invaluable benefits, there are some major drawbacks to respiratory drug delivery, the most important of which being poor delivery efficiency and relatively high drug deposition in undesirable regions, such as the mouth cavity. One way to improve the efficiency of respiratory drug delivery with metered-dose inhalers is placing a respiratory spacer between the inhaler exit and the mouth. It is argued that high drug deposition in the immediate airways of... 

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

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is proposed and applied to accurately compute steady and unsteady incompressible flows. Herein, the spatial derivatives in the lattice Boltzmann equation are discretized by using the fourth-order compact FD scheme, and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient incompressible flow solver. A high-order spectral-type low-pass compact filter is used to stabilize the numerical solution. An iterative initialization procedure is presented and applied to generate consistent initial conditions for the simulation of unsteady flows. A sensitivity study is also... 

We outline a comprehensive numerical procedure for modeling of species transport and surface reaction kinetics in electrokinetically actuated microfluidic devices of rectangular cross section. Our results confirm the findings of previous simplified approaches that a concentration wave is created for sufficiently long microreactors. An analytical solution, developed for the wave propagation speed, shows that, when normalizing with the fluid mean velocity, it becomes a function of three parameters comprising the channel aspect ratio, the relative adsorption capacity, and the kinetic equilibrium constant. Our studies also reveal that the reactor geometry idealized as a slit, instead of a... 

In this paper, the Immersed Moving Boundary-Lattice Boltzmann (IMB-LB) method is compared with the single relaxation time and multiple-relaxation-time versions of the Immersed Boundary-Lattice Boltzmann (IB-LB) method in terms of the amount of numerical velocity slip produced on solid boundaries. The comparisons are performed for both straight and curved boundaries based on the effects of thickness of virtual domain used in the IB method for the first time, and relaxation time parameter(s) of the LB method. For the straight boundaries, a shear flow problem is studied while for the curved boundaries, a falling circular cylinder in an infinite channel is investigated. First, sensitivities of... 

Modeling the behavior of suspended particles in gaseous phase is important for diverse reasons; e.g. aerosol is usually the main subject of CFD simulations in clean rooms. Additionally, to determine the rate and sites of deposition of particles suspended in inhaled air, the motion of the particles should be predicted in lung airways. Meanwhile there are two basically different approaches to simulate the behavior of particles suspension, Lagrangian and Eulerian approaches. This study compares the results of these two approaches on simulating the same problem. An in-house particle tracking code was developed to simulate the motion of particles with Lagrangian approach. In order to simulate the... 

Mixed-convection flow of nanofluids inside a vertical rectangular channel partially filled with open-cell metal foam and subject to a constant wall-heat flux was investigated experimentally and numerically. Al2O3-water nanofluids with different concentrations were prepared and their stability was examined using UV-Vis spectroscopy. Dynamic light scattering method was used to determine particle size distribution of the nanofluid feedstock. The outlet temperature and pressure drop were measured for different nanofluid flow rates (i.e., Reynolds number values). In the numerical section, a two-dimensional volume-averaged form of the governing equations was used. The velocity and temperature... 

A comprehensive mathematical model has been developed for the simulation of simultaneous chemical absorption of carbon dioxide and hydrogen sulfide by means of Monoethanolamine (MEA) aqueous solution in hollow fiber membrane reactors is described. In this regard, a perfect model considering the entrance regions of momentum, energy, and mass transfers was developed. Computational Fluid Dynamics (CFD) techniques were applied to solve governing equations, and the model predictions were validated against experimental data reported in the literature and excellent agreement was found. Effects of different disturbances on the dynamic behavior of the reactor were investigated. Moreover, effects of... 

The inflow pattern of liquid into a droplet is studied experimentally using a surface active dye and compared with results of CFD simulations. The results show visual agreement between experiments and simulations. The CFD simulations show also good agreement with the surface tension measured by drop profile analysis tensiometry (PAT). The inflow of the surfactant induces a Marangoni instability caused by the local arrival of the surfactant at the drop surface. The onset of this Marangoni instability observed experimentally has a delay of about 10. s when compared with the simulation results. Different scenarios are discussed, including a boundary layer barrier, a kinetic-controlled... 

Performance of polymer electrolyte membrane fuel cells (PEMFC) at high current densities is limited to transport reactants and products. Furthermore, large amounts of water are generated and may be condensed due to the low temperature of the PEMFC. Development of a two-phase flow model is necessary in order to predict water flooding and its effects on the PEMFC performance. In this paper, a multiphase mixture model (M2) is used, accurately, to model two-phase transport in porous media of a PEMFC. The cathode side, which includes channel, gas diffusion layer (GDL), microporous layer (MPL), and catalyst layer (CL), is considered as the computational domain. A multidomain approach has been used... 

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

In this paper, we study the frost formation and growth at the walls of a duct with uniform wall temperature variation. The simulation is performed for laminar flow regime considering suitable semi-empirical models incorporated with computational fluid dynamics (CFD) method. The frost growth is considered to be normal to the duct surface. Since the duct aspect ratio is high, we perform our simulations in two-dimensional zones. To simulate the frost layer properly, we solve both the energy and mass balance equations implementing some semi-empirical correlations on the frost side. At this stage, we suitably predict the required heat flux value at the solid boundary and the heat transfer... 

In this paper, we present our results about simulation of 2D-EOF in Nano/Micro scales porous media using lattice Boltzmann method (LBM) in micro-channel for EOF. The high efficient numerical code use strongly high nonlinear Poisson Boltzmann equation to predicate behavior of EOF in complex geometry. The results are developed with precisely investigation of several effective parameters on permeability of EOF, such as geometry (channel height and number and location of charge), external electric field, thickness of Debye length (ionic concentration), and zeta potential. Our results are in excellent agreement with available analytical results. Our results show that for certain external electric...