Sharif Digital Repository

Labyrinth side weirs provide a potentially effective way to shorten the length of a channel opening with the desired water height, but the discharge coefficient needs to be found experimentally. Experiments were carried out for triangular labyrinth side weirs with one and two cycles in order to obtain the discharge coefficient. The test variables included channel opening, head angle, number of cycles, weir height, upstream water depth and discharge in a subcritical situation. The resultsshow that the discharge coefficient ratio for a triangular labyrinth to a conventional side weir could rise to more than 2 with the same flow and geometric conditions. Relationships for the De Marchi... 

In this study, the fluid dynamics and convective heat transfer for turbulent flows through a 3-sub channel of a rod bundle, which is representative of those used in VVER-1000, are examined. The rod bundle is constructed from parallel rods in a hexagonal array. The rods are on constant pitch by spacer grids spaced axially along the rod bundle. The geometry details of the bundle and heat flux from the fuel rod are similar to that of the Iranian nuclear reactor under construction. A numerical study using Computational Fluid Dynamics (CFD) was carried out to estimate the flow field, pressure loss and heat transfer coefficients in spacer grids and rod bundles. Turbulence has been modeled using... 

Dense underflows are continuous currents which move down-slope due to the fact that their density is heavier than that ambient water. In this work, 2-D and 3-D density current in a channel were investigated by a set of experimental studies and the data were used to simulate the density current. The velocity components were measured using Acoustic Doppler Velocimetry (ADV). The height of density current (current's depth) was also measured. In this study, the density current with a uniform velocity and concentration enters the channel via a sluice gate into a lighter ambient fluid and moves forward down-slope. A low-Reynolds number turbulent model (Launder and Sharma, 1974) has been applied to... 

Dense underflows are continuous currents that move downslope due to their density being heavier than that of the ambient water. In this work, a steady density current with a uniform velocity and concentration from a narrow sluice gate enters into a wide channel of lighter ambient fluid and moves forward downslope. Experiments varying inlet velocity and concentration and hence inlet Richardson numbers were conducted. Numerical simulations were also performed with a low-Reynolds number k-ε model. The results of numerical simulation agree well with the experimental data. © 2009 ASCE  

An analytical solution is presented to study the heat transfer characteristics of the combined pressure - electroosmotically driven flow in planar microchannels. The physical model includes the Joule heating effect to predict the convective heat transfer coefficient in two dimensional microchannels. The velocity field, which is a function of external electrical field, electroosmotic mobility, fluid viscosity and the pressure gradient, is obtained by solving the hydrodynamically fully-developed laminar Navier-Stokes equations considering the electrokinetic body force for low wall zeta potentials. Then, assuming a thermally fully-developed flow, the temperature distribution and the Nusselt... 

The physical processes responsible for bubble formation in microchannels are not well understood and lack fundamental understanding. Experimental results are not exactly in agreement with each other and there are no definite theories to explain the possible effects of different parameters. Among different parameters the microchannel hydraulic diameter can affect the bubble formation mechanism in microchannels strongly. In this paper the effect of the time dependent pressure difference between inside and outside of the bubble on bubble dynamics in microchannels have been investigated. The source of this time dependency can be the emergence of bubble embryos which produces a density... 

Aero-thermoelastic analysis of a simply supported functionally graded truncated conical shell subjected to supersonic air flow is performed to predict the flutter boundaries. The temperature-dependent properties of the FG shell are assumed to be graded through the thickness according to a simple rule of mixture and power-law function of volume fractions of material constituents. Through the thickness steady-state heat conduction is considered for thermal analysis. To perform the stability analysis, the general nonlinear equations of motion are first derived using the classical Love's shell theory and the von Karman-Donnell-type of kinematic nonlinearity together with the linearized... 

This study investigates the effect of stress magnitude and stress history on porosity and permeability values of anhydride and carbonate rocks. Porosity and permeability properties are measured for twelve anhydride and carbonate core samples under stress loading and unloading conditions. The results of permeability measurements show that tighter core samples are more stress dependent while the anhydride samples are generally more sensitive to the stress. The gap between stress loading and unloading (hysteresis) is more considerable at lower effective stress values. The results also indicate that the hysteresis is more noticeable in the anhydride core samples. The gas slippage factor is also... 

We numerically solve the Navier-Stokes equations to study the rarefied gas flow in short micro-and nanoscale channels. The inlet boundary conditions play a critical role in the structure of flow in short channels. Contrary to the classical inlet boundary conditions, which apply uniform velocity and temperature profiles right at the real channel inlet, we apply the same inlet boundary conditions, but at a fictitious position far upstream of the real channel inlet. A constant wall temperature incorporated with suitable temperature jump is applied at the channel walls. Our solutions for both the classical and extended inlet boundary conditions are compared with the results of other available... 

In the present study, a high-order compact finite-difference lattice Boltzmann method is applied for accurately computing 3-D incompressible flows in the generalized curvilinear coordinates to handle practical and realistic geometries with curved boundaries and nonuniform grids. The incompressible form of the 3-D nineteen discrete velocity lattice Boltzmann method is transformed into the generalized curvilinear coordinates. Herein, a fourth-order compact finite-difference scheme and a fourth-order Runge-Kutta scheme are used for the discretization of the spatial derivatives and the temporal term, respectively, in the resulting 3-D nineteen discrete velocity lattice Boltzmann equation to... 

In the present study, a high-order compact finite-difference lattice Boltzmann method is applied for accurately computing 3-D incompressible flows in the generalized curvilinear coordinates to handle practical and realistic geometries with curved boundaries and nonuniform grids. The incompressible form of the 3-D nineteen discrete velocity lattice Boltzmann method is transformed into the generalized curvilinear coordinates. Herein, a fourth-order compact finite-difference scheme and a fourth-order Runge-Kutta scheme are used for the discretization of the spatial derivatives and the temporal term, respectively, in the resulting 3-D nineteen discrete velocity lattice Boltzmann equation to... 

Euler conservation equations, ideal gas state equations and simplified chemical kinetics models were used to simulate two-dimensional straight and baffled shock tubes. In a straight channel, detonation waves were initiated by a strong shock wave and allowed to travel down the channel to reach a CJ wave condition. It has been shown that a two-step reaction, kinetics model with an induction time delay, resulted in a physically plausible transient solution. The one-step kinetics model solution is only valid at the limit of a steady state CJ wave condition and should not be used for transient problems. The two-step kinetics model was then used to simulate a detonation initiation in a baffled... 

Proton exchange membrane (PEM) fuel cell performance is directly related to the bipolar plate design and their channels pattern. Power enhancements can be achieved by optimal design of the type, size, or patterns of the channels. It has been realized that the bipolar plate design has significant role on reactant transport as well as water management in a PEM Fuel cell. Present work concentrates on improvements in the fuel cell performance by optimization of flow-field design and channels configurations. A three-dimensional, multi-component numerical model of flow distribution based on Navier-Stokes equations using individual computer code is presented. The simulation results showed excellent... 

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

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

In this paper an analytical model is proposed to study the nonlinear dynamic behavior of rolling element bearing systems including surface defects. Various surface defects due to local imperfections on raceways and rolling elements are introduced to the proposed model. The contact force of each rolling element described according to nonlinear Hertzian contact deformation and the effect of internal radial clearance has been taken into account. Mathematical expressions were derived for inner race, outer race and rolling element local defects. To overcome the strong nonlinearity of the governing equations of motion, a modified Newmark time integration technique was used to solve the equations... 

The use of multi-element airfoils has been known as a major approach to boost up the lift of wing without dramatic increase in its drag. In fact, the configuration helps to reduce the chance of flow separation over the airfoil. However, the use of a complicated geometry such as multi-element airfoil would normally cause complexity in flow behavior. The experience has shown that the flow field complexities cannot be properly modeled using standard two-equation k-epsilon turbulence model. Therefore, it is important to improve the accuracy of general turbulence models in specific applications and complex computational domains. In this work, we extend a suitable objective function based on... 

We investigate the stability of the pressure-driven, low-Reynolds-number flow of Brownian suspensions with spherical particles in microchannels. We find two general families of stable/unstable modes: (i) degenerate modes with symmetric and antisymmetric patterns; (ii) single modes that are either symmetric or antisymmetric. The concentration profiles of degenerate modes have strong peaks near the channel walls, while single modes diminish there. Once excited, both families would be detectable through high-speed imaging. We find that unstable modes occur in concentrated suspensions whose velocity profiles are sufficiently flattened near the channel centreline. The patterns of growing unstable...