Sharif Digital Repository

Production logging is used to evaluate wells production performance. Interpretation of production log data provides velocity profile and contribution of each zone on total production. In multi-phase flow conditions, production log interpretation can be challenging since producing fluids do not have similar densities and travel with different speed depending on fluids properties and wellbore deviation. Production log interpretation in multi-phase producing wells requires identifying downhole flow regimes and determining velocity profile for each phase. There are different flow regimes and velocity models available, which are being used in production log interpretation to determine wells flow... 

The present study focuses on the numerical simulation of cavitation around the NACA 0015. The unsteady behaviors of cavitation which have worthwhile applications are investigated. The cavitation patterns, velocity fields and frequency of the cavitating flow around hydrofoil is obtained. For multi phase simulation, single-fluid Navier-Stokes equations, along with a volume fraction transport equation, are employed. The bubble dynamics model is utilized to simulate phase change. SIMPLE algorithm is used for velocity and pressure computations. For discretization of equations the finite-volume approach written in body fitted curvilinear coordinates, on collocated grid, is used. In this study,... 

Reservoir flow simulation involves subdivision of the physical domain into a number of gridblocks. This is best accomplished with optimized gridpoint density and a minimized number of gridblocks, especially for coarse-grid generation from a fine-grid geological model. In any coarse-grid generation, proper distribution of gridpoints, which form the basis of numerical gridblocks, is a challenging task. We show that this can be achieved effectively by a novel grid-generation approach based on a background grid that stores gridpoint spacing parameters. Spacing parameter (L) can be described by Poisson's equation (▽2L = G), where the local density of gridpoints is controlled by a variable source... 

Three-dimensional simulation of incompressible flow in rotating tubes for both laminar and turbulent flows has been performed using a finite-volume method for elliptic flows. The influence of Reynolds number on the velocity field and the effects of temperature gradient on temperature profiles have been presented by numerical simulations. Also the effects of velocity field, flow regime, and temperature distribution along the tube have been studied from different points of view. The results have been calculated for rotational Reynolds numbers ranging from 1000 to 320,000. The comparisons between numerically calculated velocity field and the Nusselt number have shown satisfactory agreement with... 

In the present paper, conceptual duct shape design for kinetic energy extraction with hydrokinetic turbines is discussed. The goal is to find a single-passage axisymmetric geometry that holds stable flow with maximum kinetic energy flux at duct throat. For finding the optimum duct shape, the fluid flow was numerically simulated in a wedge shaped space with Flow-Simulation Software. In a multi-stage conceptual design, tabulated configurations were employed to study each geometrical characteristic separately. These include curvature of profile camber, trailing edge shape, profile tip shape, and duct exit cross sectional area. The revolved profile of each duct consists of a well constrained... 

In this article, geometry optimization of a jet stirred reactor (JSR) combustor has been carried out for minimum NO emissions in methane oxidation using a combined numerical algorithm based on computational fluid dynamics (CFD) and differential evolution (DE) optimization. The optimization algorithm is also used to find a fairly accurate reduced mechanism. The combustion kinetics is based on a five-step mechanism with 17 unknowns which is obtained using an optimization DE algorithm for a PSR-PFR reactor based on GRI-3.0 full mechanism. The optimization design variables are the unknowns of the five-step mechanism and the cost function is the concentration difference of pollutants obtained... 

In this work, the Chebyshev collocation spectral lattice Boltzmann method is implemented in the generalized curvilinear coordinates to provide an accurate and efficient low-speed LB-based flow solver to be capable of handling curved geometries with non-uniform grids. The low-speed form of the D2Q9 and D3Q19 lattice Boltzmann equations is transformed into the generalized curvilinear coordinates and then the spatial derivatives in the resulting equations are discretized by using the Chebyshev collocation spectral method and the temporal term is discretized with the fourth-order Runge–Kutta scheme to provide an accurate and efficient low-speed flow solver. All boundary conditions are... 

The Eulerian methods are susceptible to generate the nonphysical spurious currents in the multiphase flow simulations near the interfaces. This paper presents a new Eulerian method to accurately simulate the velocity fields, especially near the multiphase flow interfaces and prevents the numerical results from generating the nonphysical currents. A Eulerian central difference finite-volume scheme equipped with the suitable numerical dissipation terms is used to simulate incompressible multiphase flows. The interface is captured by Flux Corrected Transport-Volume of Fluid method (FCT-VOF). Increasing the accuracy near the sharp gradients, such as interface, the conservative form of... 

In this article, a v2f model is employed to conduct a series of computations of incompressible flow in a periodic array of square cylinders simulating a porous media. A Galerkin/least-squares finite element formulation employing equal order velocity-pressure elements is used to discretize the governing equations. The Reynolds number is varied from 1000 to 84,000 and different values of porosities are considered in the calculations. Results are compared to the available data in the literature. The v2f model exhibits superior accuracy with respect to κ - ε results and is closer to LES calculations. The macroscopic pressure gradients for all porosities studied showed a good agreement with... 

This paper presents a novel coarse grid generation technique based on using vorticity as a measure of the impact of heterogeneity on flow. Vorticity is a maximum when the total flow is high and perpendicular to a large permeability gradient. Maps of vorticity were generated from single-phase flow simulations and used to generate coarse simulation grids from finely gridded geological models. The resulting grid was more refined in areas of high vorticity and coarser in areas of low vorticity. The method is first demonstrated on a simple five-layered model before being applied to three, 2D models of geologically realistic heterogeneity. The homogenized model generated from vorticity maps shows... 

The production routing problem (PRP) merges the lot-sizing problem and the vehicle routing problem, two classical problems that have been the focus of comprehensive studies for over half a century. Solving the PRP is an effort to optimize decisions about the production, inventory, distribution, and routing in an integrated manner. In the literature of the recent decade, due to economic changes and regulatory issues, reverse logistics has become a focal point. Subsequently, the vehicle routing problem with simultaneous pickup and delivery (VRPSPD) has drawn more and more attention for its considerable effect on the reverse logistics problem. In addition, one of the major arguments in supply... 

In this study, silicon microchannel heat sink (MCHS) performance using nanofluids as coolants was analyzed. The nanofluid was a mixture of nanoscale Cu particles and pure water with various volume fractions. Based on theoretical models and experimental correlations, the heat transfer and friction coefficients required in the analysis were used. In the theoretical model, nanofluid was treated as a single-phase fluid. In the experimental correlation, thermal dispersion due to particle random motion was included. The microchannel heat sink performances for a specific geometries with Wch = W fin = 100 μm and Lch =300 μm is examined. In this study, flow in laminar and turbulent regimes using the... 

We derive a new non-linear two dimensional model for melt flows and interface instability in aluminum reduction cells. This model is based on non-linear de St. Venant shallow water equations and contains the main features of an aluminum reduction cell. In this model we consider linear friction terms but in a new way that has not been considered in previous works. Our results are in good agreement with the results of simulation of viscous flow. This model is applicable both in determination of melt flows in molten aluminum and cryolite layers and also in finding the extreme limit for stability of interfacial waves in an aluminum reduction cell  

Multiphase flow is a challenging area of computational fluid dynamics (CFD) due to their potential large topological change and close coupling between the interface and fluid flow solvers. As such, Lagrangian meshless methods are very well suited for solving such problems. In this paper, we present a new fully explicit incompressible Smoothed Particle Hydrodynamics approach (EISPH) for solving multiphase flow problems. Assuming that the change in pressure between consecutive time-steps is small, due to small time steps in explicit solvers, an approximation of the pressure for following time-steps is derived. To verify the proposed method, several test cases including both single-phase and... 

The Computational fluid dynamics (CFD)–discrete element method (DEM) numerical simulation may be applied to predict the hydrodynamic behavior of dense particle–fluid flows. The main drawback of this simulation is the long computational time required owing to the large number of particles and the minute time-step required to maintain a stable solution. In this work, a new method to improve the efficiency and accuracy of CFD–DEM simulations is presented. The particle stiffness coefficient is used as a flexible parameter to improve the accuracy and efficiency of the model. The particle concentration distribution results are compared with experimental one’s to derive the optimum effective... 

Decision on the location of new wells through infill drilling projects is a complex problem that depends on the reservoir rock and fluid properties, well and surface facilities specifications, and economic measures. Conventional approach to address this is a direct optimization that uses the numerical flow simulation. However, this is computationally very extensive. In this study the authors use a hybrid genetic algorithm (HGA) optimization technique based on the genetic algorithm (GA) with helper functions based on the polytope algorithm and the neural network. This hybridization introduces hill-climbing into the stochastic search and makes use of proxies created and calibrated iteratively... 

Decision on the location of new wells through infill drilling projects is a complex problem that depends on the reservoir rock and fluid properties, well and surface facilities specifications, and economic measures. Conventional approach to address this is a direct optimization that uses the numerical flow simulation. However, this is computationally very extensive. In this study the authors use a hybrid genetic algorithm (HGA) optimization technique based on the genetic algorithm (GA) with helper functions based on the polytope algorithm and the neural network. This hybridization introduces hill-climbing into the stochastic search and makes use of proxies created and calibrated iteratively... 

Simulation of the in-situ combustion process is one of the most complex simulations amongst other reservoir flow simulations. Accurate simulation of the process is critical to obtain a successful implementation of the in-situ combustion process. Several factors impact performance of the simulation of this process. First are all the numerical models used for different sub-processes, such as reactions, fluid phase behavior, heat loss to surrounding formations and fluid physical properties. In the previous numerical models of the in-situ combustion process, very simplified models were used for the phase behavior of the fluid. Recent studies show that the fluid phase behavior model has a great... 

We investigate heat transfer between parallel plates separated by liquid argon using two-dimensional molecular dynamics (MD) simulations incorporating with 6-12 Lennard-Jones potential between molecule pairs. In molecular dynamics simulation of nanoscale flows through nanochannels, it is customary to fix the wall molecules. However, this approach cannot suitably model the heat transfer between the fluid molecules and wall molecules. Alternatively, we use thermal walls constructed from the oscillating molecules, which are connected to their original positions using linear spring forces. This approach is much more effective than the one which uses a fixed lattice wall modeling to simulate the... 

The accuracy of Reynolds Averaged Navier-Stokes (RANS) turbulence models to predict the behavior of 2-D density currents has been examined. In this work, a steady density current is simulated by the k - ε, k - ε RNG, two-layer k - ε and modified v̄2 - f model, all of which are compared with the experimental data. Density currents, with a uniform velocity and concentration, enter a channel via a sluice gate into a lighter ambient fluid and move forward down-slope. The eddy-viscosity concept cannot accurately simulate this flow because of two stress production structures found within it. Results show that all isotropic models have a weak outcome on this current, but by improving the ability of...