Sharif Digital Repository

The aim of this work is to investigate the non-equilibrium effects of phase change in cavitating flows. For this purpose, the concept of phase change thermodynamic probability is used along with homogeneous model to simulate two-phase cavitating flows. For simulation of unsteady behaviors of cavitation, which have practical applications, unsteady PISO algorithm based on the non-conservative approach is utilized. For multi-phase simulation, single-fluid Navier-Stokes equations, along with the volume fraction transport equation, are employed. In this paper, phase change thermodynamics probabilities and cavitation model is briefly summarized. Thus, derivation of the cavitation model, starting... 

In this paper, the effect of different parameters on the thermal operation of a Closed Loop Pulsating Heat Pipe (CLPHP) has been investigated. These parameters include the working fluid, the inclination angle, the filling ratio, and the input heat flux. The effect of nanoparticle mass concentrations has been analyzed as well. It was observed that the CLPHP can decrease thermal resistance up to 11.5 times compared to the same empty copper tube with thermal resistance of 9.4 K/W. Optimum thermal operation for a system with the water-silver nanofluid was achieved at conditions of the 50% filling ratio with thermal resistance of 0.9 K/W, and for the water-titanium oxide system, the optimal... 

Recently, a hybrid Lattice Boltzmann Level Set Method (LBLSM) for two-phase incompressible fluids with large density differences, in cases of negligible or a priori known pressure gradients, has been proposed. In the present work, the mentioned LBLSM method is extended to take into account pressure gradient effects. The lattice Boltzmann method is used for calculating velocities, the interface is captured by the level set function, and the surface tension is replaced by an equivalent body force. The method can be applied to simulate two-phase fluid flows with density ratios up to 1000 and viscosity ratios up to 100. In order to validate the method, the evolution and merging of rising bubbles... 

Detailed geological models of a reservoir may contain many more cells that can be handled by reservoir simulators due to computer hardware limitations. Upscaling is introduced as an effective way to overcome this problem. However, recovery predictions performed on a coarser upscaled mesh are inevitably less accurate than those performed on the initial fine mesh. Dual mesh method is an approach that uses both coarse and fine grid information during simulation. In the reconstruction step of this method, the equations should be solved numerically within each coarse block, which is a time consuming process. Recently, a new coarse-grid generation technique based on the vorticity preservation... 

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is extended and applied to accurately simulate two-phase liquid-vapor flows with high density ratios. Herein, the He-Shan-Doolen-type lattice Boltzmann multiphase model is used and the spatial derivatives in the resulting equations are discretized by using the fourth-order compact finite-difference scheme and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient two-phase flow solver. A high-order spectral-type low-pass compact nonlinear filter is used to regularize the numerical solution and remove spurious waves generated by flow nonlinearities in smooth regions... 

In this study, the unsteady three dimensional nanofluid flow, heat and mass transfer in a rotating system in the presence of an externally applied uniform vertical magnetic field is investigated. This study has different applications in rotating magneto-hydrodynamic (MHD) energy generators for new space systems and also thermal conversion mechanisms for nuclear propulsion space vehicles. The important effects of Brownian motion and thermophoresis have been included in the model of nanofluid. The governing equations are non-dimensionalized using geometrical and physical flow field-dependent parameters. The velocity profiles in radial, tangential and axial directions, pressure gradient,... 

In this paper a locally mass conservative finite volume method is employed to model the one-dimensional, two-phase immiscible flow in a poroelastic media. Since, an appropriate choice of primary variables is critical in simulating multiphase subsurface flow, depending on such a choice, the governing equations can be expressed in different forms. By implementing Picard iteration to a highly nonlinear system of equations, three numerical models including pressure form, mixed form and mixed form with a modified Picard linearization are developed in this study. These models have been evaluated in terms of stability, convergence and mass conservation in various one-dimensional test cases.... 

The process of coarsening the detailed geological model of a reservoir to simulation models is known as upscaling. There are two fundamental steps in the procedure of upscaling, i.e., defining the coarse grid geometry and calculating the average properties for the generated coarse grid. In this paper, the focus will be on investigating the applicability of optimization in the context of coarse grid geometry definition. To do so, different objective function candidates will be defined, and their behavior in terms of predicting the two-phase flow accuracy of coarse grids will be analyzed to determine the proper objective function. A modified objective function employing the idea of analytical... 

In this article, 3D modeling and simulation of bubbling fluidized beds has been conducted using various drag models, and the model predictions were validated against reported experimental data and 2D simulation results. In this regard, different drag models reported in the literature including Gidaspow, Syamlal-O'Brien, Hill-Koch-Ladd, and Wen-Yu were applied. A standard Two-Fluid Model (TFM) closed by the Kinetic Theory of Granular Flows (KTGF) was used to simulate bubbling gas-solid fluidized beds. Excellent agreements between the simulation results and experimental data, concerning bed expansion ratio, gas volume fraction, and time-averaged particles velocity, were found over a wide range... 

In order to safe design and optimize performance of industrial systems which work under two phase flow conditions, it's often needed to categorize flow into different regimes. In present work the experiments of two phase flow were done in a large scale test facility with length of 6m and 5cm diameter. Four main flow regimes were observed in vertical air-water two phase flows at moderate superficial velocities of gas and water: Bubbly, Slug, Churn and Annular. Some image processing techniques were used to extract information from each picture. This information include number of bubbles or objects, area, perimeter, height and width of objects (second phase).Also a texture feature extraction... 

Numerical approximation based on different forms of the governing partial differential equation can lead to significantly different results for two-phase flow in porous media. Selecting the proper primary variables is a critical step in efficiently modeling the highly nonlinear problem of multiphase subsurface flow. A comparison of various forms of numerical approximations for two-phase flow equations is performed in this work. Three forms of equations including the pressure-based, mixed pressure-saturation and modified pressure-saturation are examined. Each of these three highly nonlinear formulations is approximated using finite difference method and is linearized using both Picard and... 

In this study, a three-dimensional mathematical model, based on a nonuniform droplet concentration distribution, has been developed to simulate gas absorption in a venturi scrubber. The mass transfer process was illustrated by assuming the liquid phase as a combination of droplets and film. The flow, just as the annular two-phase flow, includes a flow of liquid film layer on the walls and a flow of gas and liquid drops in the core. Peclet number was determined using experimental data reported by Viswanathan et al. (1984) for distribution of droplets across the cross section of the scrubber. The mathematical model for gas absorption was justified by comparing the theoretical predictions with... 

In the present study, a theoretical model of a reaction–diffusion within an entrapped-cell photobioreactor packed with gel-granules containing immobilized photosynthetic bacterial cells is presented. The model is based on a system of two coupled nonlinear reaction–diffusion equations under steady-state condition for biochemical reactions occurring in the photobioreactor that describes the substrate and product concentration within the gel-granule. Simple analytical expressions for the concentration of substrate and product have been derived for all values of reaction–diffusion parameters, demonstrating competition between the diffusion and reaction in the gel-granule, using the homotopy... 

In this study, the Lattice Boltzmann Method (LBM) is used to investigate the deformation of two droplets within microfluidic T-junctions (MFTD). In order to increase the accuracy the two immiscible fluids are modeled using the He-Chen-Zhang model. First, this model is applied to ensure that the surface tension effect existing between the droplets and the continuous fluid is properly implemented in the model. Then the collision and merging of the two droplets within the intersection of a T-shaped microchannel is investigated. For generating droplet formation the effects of relevant dimensionless parameters such as the Reynolds, the Weber numbers as well as a collision parameter affecting the... 

This paper deals with investigations of droplet dynamics in rotating flows. In many previous studies droplet dynamics was analyzed in simple unidirectional flows. To fill this gap, the focus of this study is an overview on investigations of droplet dynamics in a complex rotating flow. A Lattice Boltzmann Method with high potential in simulation of two-phase unsteady flows is applied to simulate the physics of the problem in a lid-driven cavity. In spite of its simple geometry, there is a complex rotating flow field containing different vortices and shear regions. The Reynolds number based on the cavity length scale and the upper wall velocity, ReL, is considered to be 1000. We discuss here... 

The hydrodynamics of two immiscible fluids in a rectangular Hele-Shaw cell under the influence of a magnetic field is studied, both theoretically and numerically. A linear stability analysis is conducted to determine the effect of magnetic fields on the formation of viscous fingers. As a result, an analytical solution is found to calculate the growth rate of perturbations. For numerical simulation of the two-phase flow, the interfacial tension is treated as a body force using the continuum surface force model and the interface tracking is performed by the volume of fluid method. The variations of the width and growth rate of fingers in an unstable displacement versus Hartmann number, a... 

A computational tool is developed to capture the interaction of solid object with two-phase flow. The full two-dimensional Navier–Stokes equations are solved on a regular structured grid to resolve the flow field. The level set and the immersed boundary methods are used to capture the free surface of a fluid and a solid object, respectively. A two-step projection method along with Multi-Processing (OpenMP) is employed to solve the flow equations. The computational tool is verified based on numerical and experimental data with three scenarios: a cylinder falling into a rectangular domain due to gravity, transient vertical oscillation of a cylinder by releasing above its equilibrium position,... 

In this article, a standard 2D Two-Fluid Model (TFM) closed by the kinetic theory of granular flow (KTGF) has been applied to simulate the behavior of tapered-in and tapered-out fluidized bed reactors. In this regard, two types of chemical reactions with gas volume reduction and increase were considered to investigate the effects of initial static bed height on the fractional conversion and bed pressure drop. To validate the CFD model predictions, the results of hydrodynamic simulations concerning bed pressure drop and bed expansion ratio were compared against experimental data reported in the literature and excellent agreement was observed. The obtained simulation results clearly indicate... 

In this study, the unsteady three dimensional nanofluid flow, heat and mass transfer in a rotating system in the presence of an externally applied uniform vertical magnetic field is investigated. This study has different applications in rotating magneto-hydrodynamic (MHD) energy generators for new space systems and also thermal conversion mechanisms for nuclear propulsion space vehicles. The important effects of Brownian motion and thermophoresis have been included in the model of nanofluid. The governing equations are non-dimensionalized using geometrical and physical flow field-dependent parameters. The velocity profiles in radial, tangential and axial directions, pressure gradient,... 

In this paper, the effect of several different parameters on the thermal resistance of a Closed Loop Pulsating Heat Pipe (CLPHP) has been investigated. These parameters include the working fluid, the inclination angle, the filling ratio and the heat influx. Also, the impact of using nanofluids with different nano-particle concentrations has been analyzed. It was observed that a CLPHP can increase the heat transfer up to 11.5 times compared to an empty pipe. Optimum performance for a system with the water-silver nanofluid was achieved at conditions of 50% filling ratio and 0.9 K/W of thermal resistance, and for the water-titanium oxide system, these optimal conditions were found to be 40%...