Sharif Digital Repository

Numerical analysis of the nanofluid flow characteristics of perforated conical rings in a heat exchanger tube has been investigated under constant wall temperature condition. The pitch ratio of the perforated conical rings is 4 and the number of holes is varied from 0 (typical conical ring) to 10. The flow regime is fully turbulent with the Reynolds number is varied from 5000 to 14,000 and Cu-water nanofluid 0<ϕ<1.5% is selected as the working fluid. The main novelty of this paper is to perform a 3D simulation of this problem because some previous studies using similar geometry were restricted to experimental analysis. The Reynolds averaged Navier Stokes (RANS) equations are solved with the... 

Liquid–liquid two-phase flow in microchannels is capable of boosting the heat removal rate in cooling processes. Formation of different two-phase flow patterns which affect the heat transfer rate is numerically investigated here in a T-junction containing water-oil flow. For this purpose, the finite element method (FEM)is applied to solve the unsteady two-phase Navier–Stokes equations along with the level set (LS)equation in order to capture the interface between phases. It is shown that the two-phase flow pattern in microchannels depends on the flow initial condition which causes hysteresis effect in two-phase flow. In this study, the hysteresis is observed in flow pattern and consequently... 

Actuator disc models are frequently used to provide a semi-analytical approach to estimating aerodynamic loads on rotary blades. The basic idea is to distribute the aerodynamic loads on a virtual rotating disc instead of simulating the actual rotating blade. These loads are then imposed to represent the source terms of the Navier-Stokes equations, which can be solved numerically using the computational fluid dynamic methods. The thickness of the actuator disk grid is one important factor considerably affecting calculations of the wind turbine rotor. Past researches generally considered the idea of fixed grid thickness exerting along the blade in their actuator disk modeling. However, this... 

One important functionality of liver cells is ammonia detoxification and urea production. In this study, a numerical model of the urea cycle in hepatocytes was developed. Navier Stokes and convection equations were employed to study the process of ammonia elimination and urea production using a microfluidic channel. The concentration of urea and ammonia throughout the channel was obtained. Furthermore, the urea cycle was modelled with respect to its four main enzymes. This resulted in twelve rate equations that were solved to determine the concentration of each metabolite participating in the urea cycle. Application of results implied common disorders such as hyperammonemia types I and II... 

In this paper, the general analysis of two-dimensional steady-state thermal and mechanical stresses for a hollow thick infinite cylinder made of Functionally Graded Piezoelectric Materials (2D-FGPMs) is performed and developed. The general form of thermal, mechanical, and electrical boundary conditions is considered on the inside and outside surfaces. A direct method is used to solve the heat conduction equation and the non-homogenous system of partial differential Navier equations, using the complex Fourier series and the power law functions method. The material properties are assumed dependent on the radial and circumferential variables and are expressed as power law functions along the... 

One important functionality of liver cells is ammonia detoxification and urea production. In this study, a numerical model of the urea cycle in hepatocytes was developed. Navier Stokes and convection equations were employed to study the process of ammonia elimination and urea production using a microfluidic channel. The concentration of urea and ammonia throughout the channel was obtained. Furthermore, the urea cycle was modelled with respect to its four main enzymes. This resulted in twelve rate equations that were solved to determine the concentration of each metabolite participating in the urea cycle. Application of results implied common disorders such as hyperammonemia types I and II... 

In this paper, the general analysis of two-dimensional steady-state thermal and mechanical stresses for a hollow thick infinite cylinder made of Functionally Graded Piezoelectric Materials (2D-FGPMs) is performed and developed. The general form of thermal, mechanical, and electrical boundary conditions is considered on the inside and outside surfaces. A direct method is used to solve the heat conduction equation and the non-homogenous system of partial differential Navier equations, using the complex Fourier series and the power law functions method. The material properties are assumed dependent on the radial and circumferential variables and are expressed as power law functions along the... 

Fouling in blood flow is very common and may decrease the blood flow in human body and lead to critical health issues. Upon injury in a blood vessel, the body's defensive system triggers a process to create a blood clot called “Thrombus”, which prevents bleeding. Blood clots are formed by a combination of blood cells, platelets, and fibrins. In this study, we investigate a controlled drug delivery using the magnetic nanoparticles in blood vessels under the influence of magnetic fields. For this purpose the Maxwell and the Navier-Stokes equations for the system are solved. In contrary to the previous studies it is assumed that the blood is a non-Newtonian fluid. The number of particles has... 

In this work, a high-order nodal discontinuous Galerkin method (NDGM) is developed and assessed for the simulation of 2D incompressible flows on triangle elements. The governing equations are the 2D incompressible Navier–Stokes equations with the artificial compressibility method. The discretization of the spatial derivatives in the resulting system of equations is made by the NDGM and the time integration is performed by applying the implicit dual-time stepping method. Three numerical fluxes, namely, the local Lax–Friedrich, Roe and AUSM+-up are formulated and applied to assess and compare their accuracy and performance in the simulation of incompressible flows using the NDGM. Several... 

Bioreactor system has been used in bone tissue engineering in order to simulate dynamic nature of bone tissue environments. Perfusion bioreactors have been reported as the most efficient types of shear-loading bioreactor. Also, combination of forces, such as rotation plus perfusion, has been reported to enhance cell growth and osteogenic differentiation. Mathematical modeling using sophisticated infrastructure processes could be helpful and streamline the development of functional grafts by estimating and defining an effective range of bioreactor settings for better augmentation of tissue engineering. This study is aimed to conduct computational modeling for newly designed bioreactors in... 

Heat generation from very large-scale integrated (VLSI) circuits increases with the advent of high-density integrated circuit technology. One of the promising techniques is liquid cooling by using microchannel heat sink. Numerical works on the microchannel heat sink in the literature are mostly two dimensional. The purpose of the present study is to develop a three-dimensional analysis procedure to investigate flow and conjugate heat transfer in the microchannel-based heat sink for electronic packaging applications. The micro-heat sink model consists of a 10 mm long silicon substrate, with rectangular microchannels, 57 μm wide and 180 μm deep, fabricated along the entire length. A finite... 

In this investigation, the permeability of interdendritic liquid flow through the equiaxial mushy zone has been modelled for Pb-Sn alloys based on experimental measurements. In the present work by solving Navier-Stokes equation, the flow pattern around the equiaxed dendrite has been obtained and then permeability has been determined by applying Darcy's law. Numerical determined values of permeabilities have been analysed by the use of Statistical Package for the Social Sciences (SPSS) statistical software. Then an experimental method has been used to measure the permeability for flow through equiaxial mushy zone of Pb-Sn alloys. Results show that increasing the solid fraction and... 

This work is conducted with evaluation of different turbulence models capabilities in predicting three dimensional jet-into-crossflow (JICF) interactions. For this purpose, first of all, comprehensive discussions on the near wall flow complexities due to discharge of a jet into a crossflow are presented. In this regards, large scale coherent structures such as: counter rotating vortex pairs (CRVP's), near wall secondary motions, horseshoe vortices, and wall jets like are discussed. Secondly, the abilities of different turbulence models in predicting such flows (JICF) are evaluated. Our evaluation is based on three points of view including: 1) JICF characteristics, 2) computed location, and... 

Catalytic partial oxidation of methane in short residence time rhodium coated monolithic reactors offers an attractive route for syngas production. The plug flow and boundary layer flow approximations are considered as computationally efficient substitutes for the full Navier-Stokes model of the reactor while including detailed heterogeneous and homogeneous chemistry. The one dimensional plug flow model has trivial computational demands but only a limited range of application. The boundary layer model provides an excellent, computationally manageable substitute for the full Navier-Stokes model over a wide range of operating conditions. Using the 38-step elementary surface reaction mechanism... 

The numerical solution to the parabolized Navier-Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth-order compact finite-difference method. The PNS and IPNS equations in the general curvilinear coordinates are solved by using the implicit finite-difference algorithm of Beam and Warming type with a high-order compact accuracy. A shock-fitting procedure is utilized in both compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are simultaneously... 

Oil trapping behavior during the pre-flush stage is critically important to evaluate the effectiveness of matrix acidizing for the oil well stimulation. In this study, the visco-capillary behavior of the two-phase flow in the pore-scale is analyzed to investigate the influence of wetting properties for a natural rock sample. A two-dimensional model, based on Cahn-Hilliard phase-field and Navier-Stokes equations, was established and solved using the finite element method. A stability phase diagram for log capillary number (Ca)-log viscosity ratio (M) was constructed and then compared with the reported experimental works. The maximum and minimum ranges of capillary number and viscosity ratio... 

Using our recently developed model, for the first time in the literature, the effect of fluid/fluid and rock/fluid interactions on the performance of Low Salinity Waterflooding (LSWF, as an Enhanced Oil Recovery process) at pore-doublet scale is investigated. The model is incorporated into OpenFOAM and both the Navier-Stokes equation for oil/water two-phase flow and the advection-diffusion equation for ion transport (at both fluid/fluid and rock/fluid interface) are solved via direct numerical simulation (DNS). The model is validated against imbibition and drainage pore-doublet experiments reported in the literature, and then applied to investigate the sole effect of wettability alteration... 

Low salinity waterflooding (LSWF) is a process in which by lowering the ionic strength and/or manipulation of the composition of the injection water, the long term equilibrium in oil/brine/rock system is disturbed to reach a new state of equilibrium through which the oil production will be enhanced due to fluid/fluid and/or rock/fluid interactions. In spite of recent advances in the simulation of the LSWF at core scale and beyond, there are very few works that have modelled and simulated this process at the pore scale specially using direct numerical simulation (DNS). As a result the effects of wettability alteration and/or Interfacial Tension (IFT) change on the distribution of the phases... 

Using our recently developed model, for the first time in the literature, the effect of fluid/fluid and rock/fluid interactions on the performance of Low Salinity Waterflooding (LSWF, as an Enhanced Oil Recovery process) at pore-doublet scale is investigated. The model is incorporated into OpenFOAM and both the Navier-Stokes equation for oil/water two-phase flow and the advection-diffusion equation for ion transport (at both fluid/fluid and rock/fluid interface) are solved via direct numerical simulation (DNS). The model is validated against imbibition and drainage pore-doublet experiments reported in the literature, and then applied to investigate the sole effect of wettability alteration... 

Low salinity waterflooding (LSWF) is a process in which by lowering the ionic strength and/or manipulation of the composition of the injection water, the long term equilibrium in oil/brine/rock system is disturbed to reach a new state of equilibrium through which the oil production will be enhanced due to fluid/fluid and/or rock/fluid interactions. In spite of recent advances in the simulation of the LSWF at core scale and beyond, there are very few works that have modelled and simulated this process at the pore scale specially using direct numerical simulation (DNS). As a result the effects of wettability alteration and/or Interfacial Tension (IFT) change on the distribution of the phases...