Sharif Digital Repository

In this work, the rising of a single bubble in a quiescent liquid under microgravity condition was simulated. In addition to general studies of microgravity effects, the initiation of hydrodynamic convection, solely due to the variations of interface curvature (surface tension force) and thus the generation of shearing forces at the interfaces, was also studied. Then, the variation of surface tension due to the temperature gradient (Marangoni convection), which can initiate the onset of convection even in the absence of buoyancy, was studied. The related unsteady incompressible full Navier-Stokes equations were solved using a finite difference method with a structured staggered grid. The... 

In this work, a supersonic turbulent flow over a long axisymmetric body was investigated, both experimentally and computationally. The experimental study consisted of a series of wind tunnel tests for the flow over an ogive-cylinder body at a Mach number of 1.6 and at a Reynolds number of 8 × 10 6, at angles of attack between -2 and 6 degrees. It included the surface static pressure and the boundary layer profile measurements. Further, the flow around the model was visualized using a Schlieren technique. All tests were conducted in the trisonic wind tunnel of the Qadr Research Center (QRC). Also, the same flow at zero angle of attack was computationally simulated using a multi-block grid... 

The kinetic-theory-based numerical schemes, such as direct simulation Monte Carlo (DSMC) and information preservation (IP), can be readily used to solve transition flow regimes. However, their high computational cost still promotes the researchers to extend the Navier-Stokes (NS) equations beyond the slip flow and to the transition regime applications. Evidently, a suitable extension would accurately predict both the local velocity profiles and the mass flow rate magnitude as well as the streamwise pressure distribution. The second-order slip velocity model derived from kinetic theory can provide relatively accurate velocity profiles up to a Knudsen (Kn) number of around 0.5; however, its... 

The unsteady flow of blood through stenosed artery, driven by an oscillatory pressure gradient, is studied. An appropriate shape of the time-dependent stenoses which are overlapped in the realm of the formation of arterial narrowing is constructed mathematically. A msathematical model is developed by treating blood as a non-Newtonian fluid characterized by the Oldroyd-B and Cross models. A numerical scheme has been used to solve the unsteady nonlinear Navier-stokes equations in cylindrical coordinate system governing flow, assuming axial symmetry under laminar flow condition so that the problem effectively becomes two-dimensional. Finite difference technique was used to investigate the... 

We use a direct simulation Monte Carlo (DSMC) method to simulate gas heating/cooling and choked subsonic flows in micro/nanoscale channels subject to either constant wall temperature or constant/variable heat flux boundary conditions. We show the effects of applying various boundary conditions on the mass flow rate and the flow parameters. We also show that it is necessary to add a buffer zone at the end of the channel if we wish to simulate more realistic conditions at the channel outlet. We also discuss why applying equilibrium-based Maxwellian distribution on molecules coming from the channel outlet, where the flow is nonequilibrium, will not disturb the DSMC solution. The current... 

Heat generation from Very Large-Scale Integrated (VLSI) circuits increases with the development of high-density integrated circuit technology. One of the efficient techniques is liquid cooling by using a microchannel heat sink. Numerical simulations on the microchannel heat sink in the literature are mainly two dimensional. The purpose of the present study is to develop a three-dimensional procedure to investigate flow and conjugate heat transfer in the microchannel heat sink for electronic packaging applications. A finite volume numerical code with a multigrid technique, based on an additive correction multigrid (AC-MG) scheme, which is a high-performance solver, is developed to solve the... 

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

The accuracy of Reynolds averaged Navier-Stokes (RANS) turbulence models to predict the behavior of two-dimensional (2-D) density current has been examined. In this work, a steady density current is simulated by the κ -ε , κ -ε RNG, two-layer κ -ε and modified v 2̄ - f models. All models are compared with available experimental data. Density current with uniform velocity and concentration enters a channel via a sluice gate into a lighter ambient fluid and moves forward down-slope. The eddy-viscosity concept cannot accurately simulate this flow because of two stress production structures in it. Results show that all isotropic models have a weak outcome for this current, but with improving the... 

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 study seeks to explore the ability of unsteady Reynolds-averaged Navier-Stokes (URANS) simulation approach for resolving two-dimensional (2D) gravity currents on fine computational meshes. A 2D URANS equations closed by a buoyancy-modified k-ε turbulence model are integrated in time with a second-order fractional step approach coupled with a direct implicit method and discretized in space on a staggered mesh using a second-order accurate finite volume approach incorporating a high resolution semi-Lagrangian technique for the convective terms. A series of 2D simulations are carried out for gravity currents from both discontinuous and continuous sources. Comparisons with experimental... 

In some diseases there is a focal pattern of velocity in regions of bifurcation, and thus the dynamics of bifurcation has been investigated in this work. A computational model of blood flow through branching geometries has been used to investigate the influence of bifurcation on blood flow distribution. The flow analysis applies the time-dependent, three-dimensional, incompressible Navier-Stokes equations for Newtonian fluids. The governing equations of mass and momentum conservation were solved to calculate the pressure and velocity fields. Movement of blood flow from an arteriole to a venule via a capillary has been simulated using the volume of fluid (VOF) method. The proposed simulation... 

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

In this paper, the truly Meshless Local Petrov-Galerkin (MLPG) method is extended for computation of steady incompressible flows, governed by the Navier-Stokes equations (NSE), in vorticity-stream function formulation. The present method is a truly meshless method based on only a number of randomly located nodes. The formulation is based on two equations including stream function Poisson equation and vorticity advection-dispersion-reaction equation (ADRE). The meshless method is based on a local weighted residual method with the Heaviside step function and quartic spline as the test functions respectively over a local subdomain. Radial basis functions (RBF) interpolation is employed in shape... 

An appropriate combination of the thin-layer Navier-Stokes (TLNS) and parabolized Navier-Stokes (PNS) solvers is used to accurately and efficiently compute hypersonic flowfields of equilibrium air around blunt-body configurations. The TLNS equations are solved in the nose region to provide the initial data plane needed for the solution of the PNS equations. Then the PNS equations are employed to efficiently compute the flowfield for the afterbody region by using a space marching procedure. Both the TLNS and PNS equations are numerically solved by using the efficient implicit non-iterative finite-difference algorithm of Beam and Warming. A shock fitting technique is used in both the TLNS and... 

In this study, the permeability for interdendritic liquid flow through equiaxial mushy zone has been modeled and experimentally measured. In the present work, by applying a virtual dendrite in a micro-domain and solving Navier-Stokes equation, flow pattern around the dendrite has been obtained and then by applying Darcy's law to this 3D domain the permeability has been determined. In this micro-model the influence of solid fraction and geometry of dendrites have been assessed. Numerically determined values of permeabilities have been analyzed by the use of SPSS statistical software. Then an experimental method is used to measure the permeability for flow through equiaxial mushy zone of Pb-Sn... 

An incompressible-smoothed particle hydrodynamics (I-SPH) formulation is presented to simulate impulsive waves generated by landslides. The governing equations, Navier-Stokes equations, are solved in a Lagrangian form using a two-step fractional method. Landslides in this paper are simulated by a submerged mass sliding along an inclined plane. During sliding, both rigid and deformable landslides mass are considered. The present numerical method is examined for a rigid wedge sliding into water along an inclined plane. In addition solitary wave generated by a heavy box falling inside water, known as Scott Russell wave generator, which is an example for simulating falling rock avalanche into... 

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

In this research, the centrifugal compressor of a turbocharger is investigated experimentally and numerically. Performance characteristics of the compressor were obtained experimentally by measurements of rotor speed and flow parameters at the inlet and outlet of the compressor. Three dimensional flow field in the impeller and diffuser was analyzed numerically using a full Navier-Stokes program with SST turbulence model. The performance characteristics of the compressor were obtained numerically, which were then compared with the experimental results. The comparison shows good agreement. Furthermore, the effect of area ratio and tip clearance on the performance parameters and flow field was... 

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

A new finite-volume-based finite-element method using the quadratic elements is developed in the present study, to analyze the flow in micro and nano sizes with higher-order slip boundary conditions. The method is applied to gaseous flow in micro and nanoscale-channels. The developed method is carried out over a wide range of Knudsen numbers, which cover not only the continuum slip flow regime with 0≤Kn≤0.1 but also it entire the range of transient regime with 0.1