Sharif Digital Repository

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method for simulation of fluid flow. Movable particles are used to interpolate fluid properties instead of fixed Eulerian grids in this method. By essence, the aforesaid particles are movable points which carry the fluid properties such as temperature, enthalpy, density, etc. these variables can also be the stress, strain or the rheological properties of the fluid. Thus, this method needs no grids and the fluid properties are calculated using an interpolation mechanism which classifies this method as meshless. A computational code is designed using SPH for the simulation of multiphase flows. Usage of the Lagrangian formulation has made it... 

Multilateral wells are relatively new alternatives in field development and because they can increase the productivity per well and reduce the cost of field development, they are rapidly replacing old methods. In order to predict the performance of well, three methods are usually used: Analytical methods which use closed-forms models, semi-analytical methods which entails line-source model and numerical methods which is based on techniques such as finite volume or finite difference. In this thesis, line-source model is extended from two dimensional forms to three dimensional forms in order to be used to optimize the path of multilateral wells. Optimization entails direct optimization... 

Building simple and effective models are essential to many applications, such as building performance diagnosis and optimal control. Detailed physical models are time consuming and often not cost-effective. Black box models require large amount of training data and may not always reflect the physical behaviors. In this study, a method is proposed to simplify the building thermal model and to identify the parameters of the simplified model. For instance, predicting the next-24-hour load in a building is essential for the optimal control of heating, ventilating and air-conditioning (HVAC) systems that use thermal or cool storage technology. In the first phase of this study, equations for... 

Physics of the present work is 2D two-phase immiscible flow under gravitational effect through fractured porous media and the most practical formulas for transfer function in dual porosity-dual permeability models have been implemented for numerical simulation of naturally fractured reservoir. The models are: Kazemi et al. (1976), Gilman (1986), Sonier et al. (1988), Quandalle and Sabathier (1989), Lim and Aziz (1995), and Rangel-German and Kovscek (2006), which are compared using a benchmark problem and the effects of gravity, capillary, number of matrix block subdomains and dual continuum media model are studied. The benchmark problem was adopted from Gilman (1986) and involves water... 

In this project we focused on Simulation of two-phase (gas-oil) flow in 2D anisotropic porous medium on unstructured grids by finite volume method.Since the permeability assumes to be full tensor, common discretization methods for pressure and flux equations couldn’t satisfy continuity conditions then a new method of discretization named continuous flux approximation introduced.In this new method, cell center pressure calculated by directly applying flux continuity conditions to control volume boundaries. In this project a ninepoint scheme is used, it means that for making the flux continuity conditions over each control volume in two dimensions all neighboring cells are involved.The reason... 

In this article, we investigate the transient pressure response and production rate of horizontal wells in a dual porosity-dual permeability naturally fractured reservoir.First, the point source solution is derived in Laplace transform space and then pressure distribution and horizontal well response is obtained using the principle of superposition. In our model, despite of dual porosity model, matrix system is directly connected to horizontal wellbore as well as fracture system and has an influence on type curves characteristics. We also showed that type curves characteristics are affected by horizontal well length, well eccentricity, external boundary conditions, permeability and storavity... 

The objective of this thesis is to use proper orthogonal decomposition (POD) method in numerical simulation of a hydrocarbon reservoir to decrease computational time.The algorithm basically consists of two projection steps: (1) projecting the governing equations onto a subspace spanned by a low number of global empirical basis functions extracted from previous time step solutions, and (2) solving the governing equations in this reduced space and projecting the solution back on the original, high dimensional one. After projection of governing equations to new subspace, because of much less dimensions of this subspace, the computational cost of numerical solution of the equations, which mostly... 

Continuous approximation of pressure and flux fields is one of main goals of simulating flow in porous media. To compute the pressure field in heterogeneous and anisotropic media using minimum memory and time, it is necessary to use a conservative flux-continuous scheme. Nowadays, one of the most efficient approaches to minimize the solving time of large-scale problems is parallel computing. In this research, a parallel algorithm for flux-continuous multi-point flux approximation (MPFA) scheme with full pressure support (FPS) is designed and implemented on unstructured triangular grids by cell-centered approach. In this work, 2D single phase flows are investigated. The results are presented... 

Trombosis formation in blood flow is a complex biochemical and very important process for body health which impedes blood loss by creating obstacles in the path of injured wall vessels. Although this phenomenon plays a very significant role in healthy body function, unfortunately it has been observed that it could cause diseases or even stroke. The reason is that thrombosis formation in blood and its detachment from the vessel wall can result in block body capillaries and lead a catastrophic event for the body. One of the factors which motivates blood for forming unnatural thrombosis is the vessel stenosis. It would ange in normal blood flow path and cause thrombosis formation as time is... 

The solution to the hyperbolic conservation laws may involve discontinuities or shocks even if initial condition is smooth. Devising low cost methods which can identify shocks accurately, has always been a challenge. This project focuses on applying high resolution schemes to the solution of two and three phase Buckley-Leverett equation and hyperbolic conservation laws in general.Buckley-Leverett equation is a transport equation which is used to describe two and three phase flows in porous media. Weighted Essentially Non-Oscillatory (WENO) schemes with reconstruction on fluxes in finite difference formulation are used to solve one dimensional systems of equations including Euler equation as... 

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

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 water service processes, tubular heat exchangers are commonly used to provide a heat transfer medium for cooling or heating of which both are subject to fouling of heat transfer surfaces. The deposit would dominantly build up in form of fouling, crystallization, and or particulate fouling depending on the operating and surface conditions. This, in turn, may result in more demand for excessive energy consumption due to higher pumping requirements, as well as for offsetting the impact of fouling. Introducing the subject of the on-line condenser tube cleaning system as a most effective way of encountering with deposits and its components, the difficulties encountered during operation and... 

In this thesis, two types of two-phase flow are studied in one dimension. The first type is oil-gas two-phase flow which is miscible and compressible. The second one is oil-water two-phase flow which is immiscible and incompressible. Oil wet and water wet conditions are studied in the second type of flow. Two approaches, i.e. double porosity and double porosity-double permeability, are used to study naturally fractured reservoirs. In these approaches, the matrix and fracture are modeled as two separate media that exchange mass via a transfer function. Black oil model is applied to simulate the fluid system. The governing equations consist of the conservation of mass equations for each... 

In this thesis an SPH method is used for simulating multiphase two- and three-dimensional Newtonian fluid flows. The weakly compressible algorithm consists of two steps of prediction and correction. In the prediction step, the velocity field is integrated forward in time without enforcing incompressibility. The correction step consists of enforcing incompressibility by solving the pressure Poisson equation which creates a trade-off between the pressure and density variations. The surface tension force is applied in the linear momentum conservation equation by the divergence of the surface stress tensor. Advanced techniques are applied to achieve better convergence and reduce pressure-related... 

Study of two-phase flowin porous media forms a basis for investigating various fluid systemsin the nature and industries. In recent years, more and more researchers use pore-scale mod-els to explore various physical phenomena occurring in complex media such as hydrocarbonreservoirs. After solving the flow equations in pore-scale, this approach can be utilized tocalculate some macroscopic values using particular forms of averaged flow field.In this research, Smoothed Particle Hydrodynamics (SPH) is used to simulate two-phaseflow in pore spaces of two- and three-dimensional porous media. The SPH method is a La-grangian meshless method that is suitable for flows with free surfaces, deformable... 

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