Sharif Digital Repository

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

The prediction of flow behavior in porous media can provide useful insights into the mechanisms involved in CO2 sequestration, petroleum engineering and hydrology. The multi-phase flow is usually simulated by solving the governing equations over an efficient model. The geostatistical (or fine grid) models are rarely used for simulation purposes because they have too many cells. A common approach is to coarsen a fine gird realization by an upscaling method. Although upscaling can speed up the flow simulation, it neglects the fine scale heterogeneity. The heterogeneity loss reduces the accuracy of simulation results. In this paper, the relation between heterogeneity loss during upscaling and... 

In this study, the Chebyshev collocation spectral lattice Boltzmann method (CCSLBM) is developed and assessed for the computation of low-speed flows. Both steady and unsteady flows are considered here. The discrete Boltzmann equation with the Bhatnagar-Gross-Krook approximation based on the pressure distribution function is considered and the space discretization is performed by the Chebyshev collocation spectral method to achieve a highly accurate flow solver. To provide accurate unsteady solutions, the time integration of the temporal term in the lattice Boltzmann equation is made by the fourth-order Runge-Kutta scheme. To achieve numerical stability and accuracy, physical boundary... 

von Willebrand factor (VWF) is a naturally collapsed protein that participates in primary haemostasis and coagulation events. The clotting process is triggered by the adsorption and conformational changes of the plasma VWFs localized to the collagen fibres found near the site of injury. We develop coarse-grained models to simulate the adsorption dynamics of VWF flowing near the adhesive collagen fibres at different shear rates and investigate the effect of factors such as interaction and cooperativity of VWFs on the success of adsorption events. The adsorption probability of a flowing VWF confined to the receptor field is enhanced when it encounters an adhered VWF in proximity to the... 

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

This Paper presents the simulation of electroosmotic flow in nanochannels using the dissipative particle dynamics (DPD) method. Most of the past electroosmotic phenomenon studies have been carried out using the continuum flow assumptions. However, there are many electroosmotic applications in nanoscales NEMS and microscales MEMS, which need to be treated using non-continuum flow assumptions. We simulate the electroosmotic flow within the mesoscopic scale using the DPD method. Contrary to the ordinary molecular dynamics method, the DPD method provides less computational costs. We will show that the current DPD results are in very good agreement with other available non-DPD results. To expand... 

Contraction flow is one of important geometries in fluid flow both in Newtonian and non-Newtonian fluids. In this study, flow of a viscoelastic fluid through a planar 4:1 contraction with rounded corners was investigated. Six different rounding ratios (RR =0, 0.125, 0.25, 0.375, 0.438, 0.475, 0.488) was examined using the linear PTT constitutive equation at creeping flow and isothermal condition. Then the resulting PDE set including continuity, momentum, and PTT constitutive equations were implemented to the OpenFOAM software. The results clearly show vortex deterioration with increasing rounding diameter, so that when rounding corner exceeds a critical value, the vortex disappears... 

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

Forced convective of a nanofluid that consists of water and Al2O3 in horizontal tubes has been studied numerically. Computed results were validated with existing well established correlation. Two-phase Eulerian model has been implemented for the first time to study such a flow field. A single-phase model and two-phase mixture model formulations were also used for comparison. The comparison of calculated results with experimental values shows that the mixture model is more precise. It is illustrated that the single-phase model and the two-phase Eulerian model underestimates the Nusselt number. Effects of nanoparticles concentration on the thermal parameters are also discussed  

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

This paper investigates two-dimensional, time-independent elecroosmotic pressuredriven flow generated by a direct current electric potential with asymmetrical and symmetrical zeta potential distributions along the microchannel walls. Fluid flow through the horizontal microchannel is simulated using a numerical method. Two different cases are proposed to study the effect of electric potential on the flow field. First, negative electric potential is applied on the microchannel walls. In this case, large segments with negative electric potential are initially placed on the first half of the microchannel walls with two different arrangements. Afterward, smaller segments with negative electric... 

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

Presence of fracture roughness and occurrence of nonlinear flow complicate fluid flow through rock fractures. This paper presents a qualitative and quantitative study on the effects of fracture wall surface roughness on flow behavior using direct flow simulation on artificial fractures. Previous studies have highlighted the importance of roughness on linear and nonlinear flow through rock fractures. Therefore, considering fracture roughness to propose models for the linear and nonlinear flow parameters seems to be necessary. In the current report, lattice Boltzmann method is used to numerically simulate fluid flow through different fracture realizations. Flow simulations are conducted over a... 

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

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

Applying upscaling techniques is an undeniable demand in reservoir simulation, considering the difference between level of details in a geological model and level of details that can be handled by reservoir simulators. Upscaling reservoir model involves first constructing a coarse grid by employing gridding algorithms and then computing average properties for coarse grid blocks. Although various techniques have been proposed for each of these steps, one has to be aware of strengths and weaknesses of each technique before attempting to apply them. In this paper, we focus on different gridding methods and evaluate their performances. Three main grid generation techniques are considered:...