Sharif Digital Repository

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

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

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

UF6 gas flow in the total regions of a rotor with the radius of 0.1 m and the length of 1 m in axisymmetric and steady states was simulated using a new hybrid CFD-DSMC method in the OpenFOAM framework. Implicit coupled density-based scheme was performed for CFD method, and Variable Hard Sphere (VHS) and diffuse model were employed in DSMC method. Also, as an initial estimation, the local Knudsen number was applied to determine the interface location between the continuum-rarefied regions (r = 0.0855 m). Then it was modified (r = 0.084 m) to reduce the computational cost. The comparison results of pure CFD and CFD-DSMC methods illustrated that there were large differences between the flow... 

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

Simulation has been shown to be one of the most powerful tools in many areas of the health-care system, a system which attracted the attention of researchers due to its intricacy. The main goal of this study is to evaluate the patients' flow in an emergency department (ED). Discrete event simulation (DES)has been used as the core of this research so that it offers various ways to reduce the length of stay (LOS)of patients in different segments of ED by assessing various scenarios. In particular, a new approach for verification of the simulated model will be discussed. This study reports on the results that extracted from a case study on patient record data in one of the EDs located in Tehran... 

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

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

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

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

In this paper, a unique and applicable method was developed for calculation of efficiency (and mass and heat transfer coefficients) of randomly-packed distillation columns. This method has potential advantages; e.g., unlike all the available methods, it can calculate efficiency without using any empirical mass transfer and hydraulic correlations, and without the need to estimate the operational and hydraulic parameters of an operating column. It, therefore, will be free of errors and limitations of such empirical items and can be used for efficiency calculation of any random packing including new ones. Along with an analysis of the proposed method, the paper also presents a thorough analysis... 

In this study, we use direct simulation Monte Carlo method to simulate subsonic flow in nanochannels and micro/nanoscale backward-facing (BF) step considering a wide range of Knudsen number regimes. The nanochannel flow simulation indicates that the nanoscale flow through the nanochannel resembles unique features such as encountering negative pressure deviation behavior and observing flat velocity profiles at higher Knudsen number regimes. On the other hand, the micro/nano BF step flow simulations demonstrate that the length of separation region considerably decreases as the flow becomes more rarefied and approaches the transition regime. Meanwhile, the variations in the flow properties are... 

Domed roofs have been used in Iran and many other countries to cover large buildings such as mosques, shrines, churches, schools. They have been also employed in other buildings like bazaars or market places in Iran due to their favorable thermal performance. The aim of this research is to study about domed roofs thermal performance in order to determine how they can be helpful in reducing the maximum air temperature of inside buildings during the warm seasons considering all parameters like air flow around them, solar radiation, radiation heat transfer with the sky and the ground as well as some openings on the building. The results of the study show that the thermal performance of the... 

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