Sharif Digital Repository

Smoothed Particle Hydrodynamics(SPH) is the oldest Lagrangian method for solving fluid equations. Fluid is approximated with particles that represent a control mass and carry physical properties su as mass, temprature and . . . Simpler formulation relative to other methods like finite volume and finite element, no need for surface traing and ease of adding new physics are benefits of this method. In this method it is necessary to find neighboring particles that are closer than smoothing lenght (h) for ea particle. is is a time consuming operation causing SPH to be a computationaly expensive method. First simulation of Power-law fluids is considered. It has been observed that using of... 

This thesis studies a gas condensate reservoir and considers the effect of different parameters. In this group of gas reservoirs, if the pressure drops below the dew point pressure, liquid phase appears. Appearing and disappearing of phases creates four regions in the reservoir. In the region far from the production well, pressure is still above the dew point pressure and vapor phase is present. By approaching production well, pressure drops below the dew point pressure and liquid phase starts to form but it is immobile. In closer regions to the well the liquid phase saturation grows and this phase becomes mobile. As a result, relative permeability of the vapor phase decreases sharply. In... 

Phase change and deposition of solid particles in liquid flows are undesirable in industrial and natural processes and hazardous in some cases. Considering the intentions to petroleum exploration in deep waters such as Caspian Sea, the Gulf of Mexico and North Sea, precipitation and deposition of wax crystals in petroleum and gas pipelines is a real concern due to the fact that during production, oil is exposed to more intense temperature gradient (the dominant mechanism in wax deposition) and its heavier components change their phase to solid and deposit. Wax deposition in production lines starting at crude oil storage tanks is not accepted due to reduction in production rate and... 

One of the most important issues in the petroleum industry, is how to increase the oil production from hydrocarbon reservoirs. The optimum design of an enhanced oil recovery process needs the knowledge on the physical phenomena of the working flow in oil reservoirs. An important aspect of any oil recovery process is the effectiveness of process fluid in removing oil from the rock pores at microscopic scales, in a way that it determines the success or failure of the recovery process. Meanwhile, one of the most important factors on the fluid mobility and as a result the microscopic mobility efficiency is the pressure difference between different phases in narrow paths of porous media. This... 

Magnetohydrodynamic (MHD) flow is concerned with the interaction between the magnetic field and electrically conducting fluid flow. In this investigation, the MHD flow in the duct of magnetohydrodynamic wave energy conversion (MWEC) generator was simulated by use of OpenFOAM software. In the MWEC generator the reciprocating motion of ocean waves lead to electrically conducting fluid flow in the duct of this generator. By applying magnetic field perpendicular to this flow, electric current generate in the fluid. This generator has a higher efficiency compared with other ocean wave energy conversion generators. Firstly, to evaluate accuracy and performance of mhdFOAM code of OpenFOAM, three... 

Fractured reservoirs show a different behavior from common reservoirs because of the existence of a broad network of fractures. This phenomenon makes it necessary to apply special methods for fractured reservoirs in the procedure of reservoir simulation. Since twenty percent of petroleum content in the world is buried in fractured reservoirs, investigating these reservoirs is of great importance.
The first step in simulation of these kinds of reservoirs is to come up with a geometrical model which can be used to take the fracture network influence into account. In the course for reaching such an objective, various models have been developed which are based on specific assumptions and in... 

In this thesis the oil reservoir simulation is carried out by using a combination of fast Multi-scale finite volume method and Streamline method considering gravity effects. In this simulation, the governing equations for incompressible fluid flow in porous media are solved using the IMPES (implicit pressure- explicit saturation) viewpoint. Multi- scale method is rapid (equivalent on coarse-scale methods) and accurate (equivalent on fine-scale methods) in solving the pressure field. Using the pressure field obtained through this way combined with Darcy equation make us able to extract the conservative velocity field. In this study, analytical and numerical Streamline methods are used to... 

In this work, simulation of 3D three-phase naturally fractured reservoirs is studied. Both Dual Porosity and Dual Porosity-Dual Permeability models which are commonly used in commercial softwares such as ECLIPSE and CMG, are utilized. The governing equations include the mass (in Black-oil model) and linear momentum (in Darcy law form) equations, which are converted into pressure equations and a set of transport equations for the saturations and mole fractions of phases. These equations are solved once for the matrix domain and subsequently for the fracture domain. These two set of equations are coupled through transfer functions. The accuracy of computations significantly depends on how... 

In this study, simulation of two-phase water-oil compressible flow in naturally fractured reservoirs is studied. The studied reservoir has a 1D and horizontal ge-ometry that its properties such as porosity and permeability are homogeneous.Dual Porosity and Dual Porosity-Dual Permeability models are employed for simulation of reservoir. Reservoir descriptive equations are conservation of mass and momentum. Conservation of momentum equation is simplified by Darcy law and two-phase fluid flow is simulated with Black-oil model. Discretization of equations is done by finite volume method. In order to solve the equations of two-phase compressible fluid flow, the so-called IMPES (Implicit Pressure... 

Innovations in laser, microwaves, and similar technologies have significantly advanced thermal treatments for diseases or even injuries concerning skin tissue. For a thorough understanding in the underlying mechanisms of bioheat transfer behavior of skin,a1D unsteady non-dimensional hyperbolic model of heat transfer through this tissue with metabolic heat generation which is subject to specific boundary conditions, is solved numerically using the finite difference method. A thermal shock is generated at the base of the tissue, which moves forward with a finite speed. A Numerical solution for a simple one-layer skin tissue is obtained. Then, the effects of various parameters, time step,... 

In this study, bulk viscosity of suspensions is predicted using numerical simulation.The elaborated numerical method consists of the lattice Boltzmann (LB) method for fluid matrix simulation and the immersed boundary (IB) method for solid boundary implementation. In order to successfully implement the IB method in the frame-work of the LB method, two important challenges, which are body force modeling in the LB method and numerical velocity slip on solid boundaries, are discussed. As a first challenge, body force modeling in the LB method is studied by using a uni-fied framework which has been devised in this study. Using this unified framework,firstly, errors of different body force... 

In this project, a response of a single well in a gas condensate reservoir has been investigated by compositional modeling and implicit numerical method. The most significant goals of this paper are as follows: 1- Examination of the flow produced in a gas condensate reservoir in constant pressure or mass rate. 2- Investigation of variables affecting the productivity of reservoir. 3- The study of effect of injection and gas cycling to the well on the production. A three-zone model has been used in this study as if there is only gas phase in the third and second zones. Condensation occurs since the pressure drops below the dewpoint pressure in the second zone; however, gas phase is only... 

In this project، the immiscible two phase flow of oil-water in the absence of gravity in the two dimensional domain is investigated. The uncertainty quantification of pressure and saturation of fluid in the reservoir are calculated by the methods of statistical moment equation (presented by Tchelepi) and probabilistic collocation method combined with Karhunen Loeve expansion (presented by Heng Lee). Then the results are compared with Monte Carlo simulation. Matlab Reservoir Simulation Toolbox is used for flow simulation. The main problem based on Tchelepi’s work[56[ ، is a horizontal two dimensional problem that there is two injection and production wells at the two end point grids (a... 

In This work, the performance of a finite volume (FV) numerical method is examined for simulating naturally fractured hydrocarbon reservoir. The dual-porosity approach applied for fracture modeling, and the flow transfer formulated based on a generalized black oil model which is based on matrix formulation, leading to a system of equations prone to degeneracy. This FV high-resolution central scheme numerical method is not sensitive to the choice of CFL number like some of other central schemes and can conveniently handle degenerate equations appearing in the reservoir simulation governing equation. In this effort, one dimension multi phase flow simulated in naturally fractured reservoir and... 

Reservoir simulation is a useful tool for efficient reservoir management and protective exploitation of non-renewable hydrocarbon resources as the main supplier of development programs in the country. Due to the nature of underground hydrocarbon reservoirs, the simulation faces several challenges, both in rock and fluid parts. These challenges should be overcome in order to provide reliable results from simulations. One of the major challenges pertaining to rock, is large discontinuous variations accompanied by anisotropy of rock permeability fields which causes numerous problems when governing equations are solved numerically. A viable numerical procedure shall account for complex reservoir... 

The main goal of this study is simulation of the flow-like landslide phenomenon and investigation of the effect of bed material entrainment on its behavior. Debris flows are considered as one of the types of flow-like landslides, and flow of dry granular materials on horizontal and inclined beds as a simple model of these types of flows in nature is numerically simulated. Also, in order to investigate the effect of entrainment on dynamics of flow, a bed containing granular materials is installed on horizontal and inclined beds. The governing equations of the flow are solved in the continuum context using the Smoothed Particle Hydrodynamics (SPH) method and the depth-averaged assumptions are... 

Measuring flow rate in various industries is an operational challenge that requires to be accurate and fast to avoid wasting energy and money. Various flowmeters are employed in different industries using different principles and methods. Sonar flowmeters are one of the acoustic methods of fluid flow measurement, which do this by listening to the sound of the flow via arrays of sensors. Sonar flow measurement can be achieved by finding either the convection velocity of eddies or the speed of sound in the fluid. The first method is suitable for measuring single-phase flows, which is the subject of this thesis. This method needs a detailed study of turbulence. In the present work, this is... 

It has been demonstrated that the key features of fracture networks can be identified with the help of high-resolution aerial images of outcrops (rock formations exposed at the Earth's surface) obtained by drones and photogrammetry. The aim of this study is to utilize these images to improve the geometrical modeling of Naturally Fractured Reservoirs (NFRs). Initially, an algorithm is developed for processing aerial images of outcrops, so that the pattern of fracture network can be obtained. After obtaining the fracture network, the Embedded Discrete Fracture Model (EDFM) implemented in the MATLAB Reservoir Simulation Toolbox (MRST) is employed to study two-phase (oil-water) flows in the... 

When classical Fourier law of heat conduction is not applicable due to finite speed of heat propagation in a medium, non-Fourier heat conduction is used. In non-Fourier heat conduction, heat transfer is studied considering finite speed of heat flux (thermal wave) or temperature gradient or both (dual phase lag). Biological tissues like skin are one of the materials which show non-Fourier behavior during usual heating processes. In this research heat transfer in human skin is modeled by Fourier, thermal wave and dual phase lag (DPL) models. The equations are solved using the finite difference method and the temperature distribution across the tissue is calculated. In thermal wave model,... 

In this research, a novel scheme for flux approximation in finite volume central schemes have been presented, evaluated and used in PARSI reservoir simulator to solve black oil equations with Trangenstein & Bell formulation with gravitational effects. Also, a numerical method, have been proposed to increase the ability of modeling of capillary phenomena to PARSI.Edward's dominant wave-capturing scheme (2005), with Harten-Hymen entropy fix (1983), has been used to solve flow equations with gravity. This scheme relies on the detection of the dominant wave in the system without recourse to characteristic decomposition and upwinding while avoiding the excessive numerical diffusion that is...