Sharif Digital Repository

Complicated sedimentary processes control the spatial distribution of geological heterogeneities. This serves to make the nature of the fluid flow in the hydrocarbon reservoirs immensely complex. Proper modeling of these heterogeneities and evaluation of their connectivity are crucial and affects all aspects of fluid flow. Since the natural variability of heterogeneity occurs in a myriad of length scales, accurate modeling of the rock type connectivity requires a very fine scheme, which is computationally very expensive. Hence, this makes other alternative methods such as the percolation approach attractive and necessary. The percolation approach considers the hypothesis that a reservoir can... 

This study describes a 1-D code, which determines the one-dimensional flow and heat transfer distributions in arbitrarily connected flow passages with multiple inlets and outlets. The procedure uses a component solver, which solves the one-dimensional compressible fluid flow and heat transfer equations inside an individual cooling passage. Also, it uses experimental correlations for the heat transfer coefficient and friction factor for each cooling passage. In addition to the component solver, there is a flow network solver subroutine, which applies the mass and energy conservation equations properly at each internal node. In order to validate this 1-D code, both the component and the flow... 

This study describes a 1-D code, which determines the one-dimensional flow and heat transfer distributions in arbitrarily connected flow passages with multiple inlets and outlets. The procedure uses a component solver, which solves the one-dimensional compressible fluid flow and heat transfer equations inside an individual cooling passage. Also, it uses experimental correlations for the heat transfer coefficient and friction factor for each cooling passage. In addition to the component solver, there is a flow network solver subroutine, which applies the mass and energy conservation equations properly at each internal node. In order to validate this 1-D code, both the component and the flow... 

Colloidal Gas Aphron (CGA) are finding increasing application in fields of science and engineering because of their distinctive characteristic. As interest in the application of CGA based fluids grows and in order to select the best procedure for using them in successful petroleum engineering operations, there is a need to gain a better understanding of the factors that affect their properties and behavior. This article discusses the rheological characterization, stability analysis and filtration properties of CGA based fluids for three bio-polymers and two ionic surfactant. The stability and filtration analysis were investigated with the static drain rate technique and API filtration tests,... 

This paper simulates Al2O3-water nanofluid flow and forced convection around a rotating circular cylinder. The governing parameters are Reynolds number (1 ≤ Re ≤ 100), solid volume fraction of nanoparticles (0 ≤ φ ≤ 0.05) and non-dimensional rotation rate (0 ≤ α ≤ 3). The simulations are performed to study the effects of mentioned parameters on the heat transfer rate and fluid flow characteristics. The governing equations including the continuity, momentum, and energy equations are solved with a finite volume method. It is observed that the reduction of heat transfer with increase in rotation rate is in the vicinity of 6.9% and 32% for Re = 5 and 100, respectively at φ = 0.05. Furthermore,... 

The algebraic multigrid (AMG) algorithm as a preconditioner to the Krylov subspace methods has drawn the attention of many researchers in solving fluid flow and heat transfer problems. However, the efficient employment of this solver needs experience, because users have to quantify several important parameters. In this work, we choose a hybrid finite-volume element method and quantify the optimum magnitudes for those parameters. To generalize our results, two sets of fluid flow governing equations, the thermobuoyant flow and confined diffusion flame, are studied and the optimum values are determined. The results indicate that the AMG can be very effective if a proper storage method is chosen... 

This study investigated the lake–aquifer hydraulic interactions in Lake Urmia (LU) as the second largest hypersaline lake in the world. Due to the scarcity of hydrogeological data required for modelling, a method based on Darcy’s Law and lake water budget was used to quantify the lake–aquifer interaction. Long-term ground- and satellite-based hydrological datasets over the time frame 2001–2019 were used. Results indicate that the groundwater flux between LU and the aquifers controls 18.74 ± 1.67% of the lake’s water storage. While 10 out of 14 adjacent aquifers recharge LU at a rate of less than 180 m3/m.month, one phreatic aquifer recharges the LU up to 1400 m3/m.month. Two aquifers are... 

This paper presents an analysis of the combined electro-osmotic and pressure-driven axial flows of viscoelastic fluids in a rectangular microchannel with arbitrary aspect ratios. The rheological behavior of the fluid is described by the complete form of Phan-Thien–Tanner (PTT) model with the Gordon–Schowalter convected derivative which covers the upper convected Maxwell, Johnson–Segalman and FENE-P models. Our numerical simulation is based on the computation of 2D Poisson–Boltzmann, Cauchy momentum and PTT constitutive equations. The solution of these governing nonlinear coupled set of equations is obtained by using the second-order central finite difference method in a non-uniform grid... 

This paper presents an analysis of the combined electro-osmotic and pressure-driven axial flows of viscoelastic fluids in a rectangular microchannel with arbitrary aspect ratios. The rheological behavior of the fluid is described by the complete form of Phan-Thien–Tanner (PTT) model with the Gordon–Schowalter convected derivative which covers the upper convected Maxwell, Johnson–Segalman and FENE-P models. Our numerical simulation is based on the computation of 2D Poisson–Boltzmann, Cauchy momentum and PTT constitutive equations. The solution of these governing nonlinear coupled set of equations is obtained by using the second-order central finite difference method in a non-uniform grid... 

Capacitance-Resistive Model (CRM), as a fast yet efficient proxy model, suffers from some limitations in modeling relatively complex reservoirs. Some current improvements on this proxy made it a more powerful simulator with updating parameters over time. However, the model's intrinsic uncertainty arisen from simplifying fluid-flow modeling by some limited number of constant parameters is not addressed yet. In this study, this structural limitation of CRM has been addressed by introducing a time-correlated model error, including stochastic and non-stochastic parameters, embedded into this proxy's formulation. The error term's non-stochastic parameters have been tuned to be used in forecasting... 

Capacitance-Resistive Model (CRM), as a fast yet efficient proxy model, suffers from some limitations in modeling relatively complex reservoirs. Some current improvements on this proxy made it a more powerful simulator with updating parameters over time. However, the model's intrinsic uncertainty arisen from simplifying fluid-flow modeling by some limited number of constant parameters is not addressed yet. In this study, this structural limitation of CRM has been addressed by introducing a time-correlated model error, including stochastic and non-stochastic parameters, embedded into this proxy's formulation. The error term's non-stochastic parameters have been tuned to be used in forecasting... 

A mathematical model is developed to assess the solidification behaviour of the weld pools. To do so, during gas tungsten arc welding of commercial pure aluminium, equations of conversation of mass, energy and momentum are numerically solved considering three-dimensional steady state heat transfer and fluid flow conditions. The weld pool geometry, weld thermal cycles and various solidification parameters are calculated using temperature and velocity fields acquiring from the utilised model. The solidification behaviour of the weld pool at the weld centreline and the fusion line is then studied using the solidification parameters including temperature gradient G, solidification rate R and the... 

Microfluidic devices have been widely used for biological and cellular studies. Microbioreactors for three-dimensional (3D) multicellular spheroid culture are now considered as the next generation in in vitro diagnostic tools. The feasibility of using 3D cell aggregates to form multicellular spheroids in a microbioreactor with U-shaped barriers has been demonstrated experimentally. A barrier array is an alternative to commonly used microwell traps. The present study investigates oxygen and glucose concentration distributions as key parameters in a U-shaped array microbioreactor using finite element simulation. The effect of spheroid diameter, inlet concentration and flow rate of the medium... 

In this study, the thermal cycles and the grain structure in the weld heat-affected zone (HAZ) are predicted. At the first stage, a combined heat transfer and fluid flow model is employed to assess the temperature fields during and after welding of 304 stainless steel and then, the evolution of grain structure is conducted using the predicted temperature distribution and an analytical model of grain growth. The grain sizes of the CGHAZ (coarse grain heat affected zone) achieved from the model are basically in agreement with those obtained from experimental measurement under different heat inputs in the range of 0.33-1.07 MJ/m. Both the experimental data and the calculated results show that... 

The current study aims at introducing a 2D and fast-running code for the issues pertinent to design, analysis and safety in modular high temperature reactors. While the porous media approach is only applied to pebble bed type, the analysis in this paper covers both pebble bed and prismatic reactor. A time-dependent mass equation along with energy conservation equation for the cooling gas and a time-dependent energy conservation equation for the solid was solved. Appropriate series of constitutive equations (e.g. heat transfer coefficient, effective heat conductivity of solid, heat transfer coefficient, pressure drop etc.) has been recruited as well. In addition a finite-volume method is... 

A comprehensive study was performed to analyze the unsteady laminar flow characteristics around a porously covered, a fully porous, and a solid squared section cylinder located in the middle of a plane channel. In order to simulate fluid flow inside porous media and porous–fluid interface accurately (minimizing modeling error), the porous region was analyzed in pore scale, using LBM. Additionally, to minimize the LBM-related compressibility error through the porous region, a multi-block multiple relaxation time lattice Boltzmann method (MRT-LBM) was used. Also, to decrease CPU time, a Navier–Stokes flow solver, based on finite volume method and SIMPLE algorithm, was coupled with MRT-LBM to... 

Since the mechanical properties and the integrity of the weld metal depend on the solidification behaviour and the resulting microstructural characteristics, understanding weld pool solidification is of importance to engineers and scientists. Thermal and fluid flow conditions affect the weld pool geometry and solidification parameters. During solidification of the weld pool, a columnar grain structure develops in the weld metal. Prediction of the formation of the microstructure during welding may be an important and supporting factor for technology optimization. Nowadays, increasing computing power allows direct simulations of the dendritic and cell morphology of columnar grains in the... 

Valveless piezoelectric micropumps are in wide practical use due to their ability to conduct particles with absence of interior moving mechanical parts. The objective of this paper is to obtain the fluid flow response to actuation frequency of a passive diffuser valve under harmonic pressures. In this regards a 2D model of a micropump valves and chambers is analyzed. The analysis is performed for 10Kpa back pressure on micropump chamber and actuation frequencies within the range of 1Hz to 10 KHz. Results show the highest velocity in the direction of diffuser axis occurs at the narrow diffuser neck while flow direction reverses every half period. For low frequencies, a parabolic velocity... 

The casing treatment and flow injection upstream of the rotor tip are two effective approaches in suppressing instabilities or recovering from a fully developed stall. This paper presents numerical simulations for a high-speed transonic compressor rotor, NASA Rotor 37, applying a state-of-the-art design for the blade tip injection. This is characterized by introducing a jet flow directly into the casing treatment machined into the shroud. The casing treatment is positioned over the blade tip region and exceeds the impeller axially by ∼30 per cent of the tip chord both in the upstream and in the downstream directions. To numerically solve the governing equations, the three-dimensional finite... 

Improvements in sealing mechanism between the rotating and the stationary parts of a rurbomachine can extensively reduce the endwall leakage flow. In this regard, abradable seals are incorporated into compressor and turbine blade-tip region. In a gas turbine, equipped with abradable seals, tip of the rotor blade is designed to cut into the material coating of the casing and to form a close fitted circumferential groove for the movement of the blade tip. As a result, the resistance to the leakage flow in the tip gap region increases due to smaller tip clearances (available without any rub-induced damages). Minimizing the tip clearance size can lead to an increase in performance and stability....