Sharif Digital Repository

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

Acidizing is one of the most effective techniques to remove the formation damage and restore (or even increase) the permeability of the near wellbore region. Although this technique is widely used to resolve the skin problem, it could initiate new damages itself that hinder the fluid flow and decrease the well production, consequently. Acid-oil emulsion and sludge formation are known as two major induced formation damage and the main reasons for the oil well acid treatment failures. Despite its critical effects, no comprehensive study has been addressed this specific type of formation damages in details. In this study, the acid-oil emulsion and asphaltic sludge formation were evaluated... 

Capillary imbibition is an important recovery mechanism in naturally fractured reservoirs when water-filled fractures surround water-wet matrix blocks. A large amount of studies of imbibition process is simply total or partial immersion of nonwetting phase saturated rock in aqueous wetting phase. However, water advance in fractures during water flooding or water encroachment from an active aquifer introduces time dependent boundary conditions where invariant exposure of rock surface to water is not representative. In this work, a laboratory simulated matrix-fracture system was used to investigate different aspects of imbibition in the presence of fracture fluid flow (namely dynamic... 

In the last decade, extensive attention is devoted to intelligibly designed materials of macro/micro-structures containing the fluid flow. In this study, intelligent control and vibrational stability of cantilevered fluid conveying macro/micro-tubes utilizing axially functionally graded (AFG) materials are considered. The governing equation of motion of the system is derived based on modified couple stress theory and then is discretized using Galerkin method. A detailed investigation is carried out to elaborate the influence of various parameters such as material properties, axial compressive load, and Pasternak foundation on the dynamical behavior of the system, all of which are influential... 

In this study, an extrinsically enriched element free Galerkin (EFG) method is proposed for the thermo-hydro-mechanical simulation of saturated porous media. By taking advantage of partition of unity property of moving least square shape functions, weak discontinuities such as material interfaces are modeled using the Ridge enrichment function and impermeable strong discontinuities are simulated using the Heaviside function. Some guidelines are proposed for the selection of EFG numerical parameters to ensure the stability and accuracy of the results. Numerical examples are provided to illustrate the capability of the proposed approach for fully coupled THM analysis of discontinuous porous... 

COMSOL Multiphysics is a comprehensive simulation software environment for a wide range of applications. COMSOL has an interactive interface that facilitates the modeling procedure and allows an easy coupling of different physical processes. The Subsurface Flow module extends the COMSOL modeling environment to applications related to fluid flow in saturated and variably saturated porous media. COMSOL is increasingly used in the investigation of geophysical, hydrogeological and environmental phenomena. The main goal of this work is to explore the ability of COMSOL for simulating seawater intrusion (SWI) in fractured coastal aquifers. Numerical modeling of such a problem is of high interest as... 

In this paper, a three-dimensional simulation of fully coupled multiphase fluid flow and heat transfer through deforming porous media is presented in the context of EFG mesh-less method. Spatial discretization of the system of governing equations is performed using EFG and a fully implicit finite difference scheme is employed for temporal discretization. Penalty method is used for imposition of essential boundary conditions. The developed numerical tool is employed to simulate two problems of nuclear waste disposal and CO2 sequestration in deep underground strata. The obtained results demonstrate the capability and robustness of the developed EFG code. © 2018 Elsevier Ltd  

The current study investigates the conjugate heat transfer characteristics for laminar flow in backward facing step channel. All of the channel walls are insulated except the lower thick wall under a constant temperature. The upper wall includes a insulated obstacle perpendicular to flow direction. The effect of obstacle height and location on the fluid flow and heat transfer are numerically explored for the Reynolds number in the range of 10 ≤ Re ≤ 300. Incompressible Navier-Stokes and thermal energy equations are solved simultaneously in fluid region by the upwind compact finite difference scheme based on flux-difference splitting in conjunction with artificial compressibility method. In... 

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 three-dimensional power law fluid flow through rough microchannels has been studied numerically to determine the effects of the topographic structures on the thermal and hydrodynamic characteristics of the system. Rectangular, triangular and sinusoidal element shapes have been considered in order to investigate the effects of roughness height, width, pitch and channel separation on the pressure drop and heat transfer. Uniform wall heat flux boundary condition has been applied for all the peripheral walls. The results indicate that the global heat transfer performance can be improved or reduced by the roughness elements at the expense of pressure head when compared with the smooth... 

Net fluid flow of soap films induced by (ac or dc) electric fields in asymmetric frames is presented. Previous experiments of controllable soap film flow required the simultaneous use of an electrical current passing through the film and an external electric field or the use of nonuniform ac electric fields. Here a single voltage difference generates both the electrical current going through the film and the electric field that actuates on the charge induced on the film. The film is set into global motion due to the broken symmetry that appears by the use of asymmetric frames. If symmetric frames are used, the film flow is not steady but time dependent and irregular. Finally, we study... 

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

Optical and thermal analysis of the most well-known solar concentrator system; parabolic trough collector (PTC) are investigated and analyzed. To evaluate performance of the PTC, four cities of Iran with different weather conditions are chosen as case studies. Effective parameters such as concentration ratio, incident angle correction factor, collector mass flow rate are considered. The main objective of this work is evaluation of the solar energy potential using PTC in under consideration cities with different climates. Numerical modeling of the analysis is done using MATLAB software. Simulation results shows that Shiraz, with an average annual thermal efficiency of 13.91% and annual useful... 

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

In this paper, an extended finite element method is presented for simulation of interaction between hydraulic fracturing and natural fractures in saturated porous media. The well-known u−p formulation is employed in order to obtain the fully coupled set of governing equations. Natural faults are modeled for both opening and closure modes where the fluid inflow and contact conditions are considered at the interface, respectively. The Darcy law is employed in conjunction with an aperture dependent permeability for the fracture channel to describe the interfacial inflow. The contact constraints of both the solid and fluid phases are imposed using the Penalty method. The Heaviside and modified... 

Successful designs for an injection process to improve reservoir recovery seek more accurate values of saturation dependent functions, including relative permeability curves. The ignorance on underlying heterogeneity (e.g. layering) might lead to inaccurate predicting of flow behavior in layered heterogeneous porous media. Analysis of experimental core flood results in composite core systems (the system) along with numerical simulation at the core scale can provide an insight into this problem, especially when the capillary effects are present. In this study unsteady state two phase displacements have been conducted on the system with different internal arrangements. The main contribution is... 

In this work, the electrokinetic properties of asphaltene particles have been investigated. Micro-electrophoresis method by applying DC electric field, was utilized to different mixtures containing asphaltene to determine its electric charge. It was observed that in the case of using n-heptane and its mixture with toluene (heptol), the asphaltene particles were showed to be positively charged however for toluene itself, they expressed no tendency toward the electrodes. While it is expected that larger asphaltene aggregates carry higher electric charge, the results contradictorily showed that they are mainly governed by gravity rather than electro-static force and that “aggregation” reduces... 

The stability of a viscoelastic plate strip, subjected to an axial load with the Kelvin–Voigt fractional order constitutive relationship is studied. Based on the classical plate theory, the structural formulation of the plate is obtained by using the Newton’s second law and the aerodynamic force due to the fluid flow is evaluated by piston theory. The Galerkin method is employed to discretize the equation of motion into a set of ordinary differential equations. To determine the stability margin of plate the obtained set of ordinary differential equations are solved using the Laplace transform method. The effects of variation of the governing parameters such as axial force, retardation time,... 

In this article, three-dimensional (3D) forced-convection heat transfer of magnetic nanofluids in a pipe subject to constant wall heat flux in the presence of single or double permanent magnet(s) or current-carrying wire has been investigated and compared. In this regard, laminar fluid flow and equilibrium magnetization for the ferrofluid were considered. In addition, variations of magnetic field in different media were taken into account and the assumption of having a linear relationship of magnetization with applied magnetic field intensity was also relaxed. Effects of magnetic field intensity, nanoparticle volume fraction, Reynolds number value, and the type of magnetic field source...