Sharif Digital Repository

Taghizadeh Manzari, M. (In this work, a modified Smoothed Particle Hydrodynamics (SPH) method, with a new moving solid boundary treatment approach, is utilized to simulate the particulateflow problems. The renormalized first and second derivative schemes which lead tothe consistency of the method, are also used along with a modification to the continuityequation which prevents the spurious pressure oscillations. The proposed methodis validated by solving benchmark problems of solid body motion in channel flows.There is a good agreement between the obtained results and those reported in theliterature. The convergence of solutions for different domain discretizations is alsoassessed. In order... 

The goal of this research is to introduce and improve the new multiscale finite volume method for simulation of multiphase flow in porous media. This method has recently presented as an efficient approach in computational fluids dynamics field. This thesis contains the followings: In the first step a brief introduction is presented about petroleum reservoir engineering and different methods of oil production. Then the importance of new simulation methods to anticipate the behavior of these reservoirs with detailed geographical information is discussed. In the second section black-oil formulation in petroleum reservoirs is descritized with the aid of finite volume method. Using general... 

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

In this thesis, an Incompressible SPH (ISPH) method is used to simulate suspensions. A suspension is assumed to be a mixture of solid particles in a Newtonian fluid. First, renormalization tensors are applied to the discretized equations of SPH method. This leads to a higher accuracy in this method. It should be noted that the projection method is used for ISPH. After that, the Immersed Boundary Method (IBM) is used to model and simulate the presence of solid particles in fluid flow. Some modifications were necessary to use IBM with Lagrangian approach. Then, some benchmark problems are solved to validate the presented method. The main objective of the present thesis is simulation of... 

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

If the reservoir temperature lies between the critical temperature and cricondentherm of the reservoir fluid, the reservoir is classified as a retrograde gas-condensate reservoir. When the pressure is decreased on these mixtures, instead of expanding or vaporizing as might be expected, they vaporize instead of condensing. If the reservoir pressure is above the upper dew-point pressure, the hydrocarbon system exists as a single phase in the reservoir (region 3). As the reservoir pressure declines isothermally during production from the initial pressure to the upper dew-point pressure, the attraction between the molecules of the light and heavy components move further apart. As this occurs,... 

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

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

In recent years, oil and gas reservoirs are one of the most valuable natural resources and are the main economical basis of the country. Therefore any research which can help to optimize the efficiency and enhance their recovery is of great importance. In the past 30 years, reservoir simulation has evolved from a research field to one of the most flexible tools in reservoir engineering. Simulation is usually more quick, cost effective and reliable than other methods in predicting reservoir performance. Because of complexity of the great amount of computations, research in oil reservoirs field is usually done by mathematical/computer programs, named as simulators. For this reason, various... 

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

Since coronary artery diseases are increasing every day, simulation of blood flow in blood vessels and their stenosis is one of the most important topics nowadays. Many efforts have been made to achieve numerical simulations using different methods such as Finite-Difference. In this thesis, an SPH method is used to simulate pulsatile blood flow in arteries. 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... 

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

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

Multi-Scale finite volume (MSFV) method have been developed and applied for multiphase flow in porous media. The most important advantage of MSFV method is its computational efficiency. In this thesis we apply and implement this method to porous media simulation. Then, computational parameters such as coarsening ratio, effect of threshold on adaptive computation and different boundary conditions for computing basis functions are investigated. The two difficulties of sMSFV method are problems with high permeability contrast and problems with high course grid aspect ratio which are considered. In standard MSFV (sMSFV) method reduced boundary condition is used to determine the basis and...