Sharif Digital Repository

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is proposed and applied to accurately compute steady and unsteady incompressible flows. Herein, the spatial derivatives in the lattice Boltzmann equation are discretized by using the fourth-order compact FD scheme, and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient incompressible flow solver. A high-order spectral-type low-pass compact filter is used to stabilize the numerical solution. An iterative initialization procedure is presented and applied to generate consistent initial conditions for the simulation of unsteady flows. A sensitivity study is also... 

We report the existence of interfacial instability in the two-dimensional channel flow of a sediment suspension whose particles diffuse in the carrier fluid due to shear-induced collisions. We derive partial differential equations that govern the deformations of the interface between the sediment suspension and the clear fluid, and devise a perturbation method that preserves the positivity of the particle volume fraction. We solve perturbed momentum, particle transport and deforming interface equations to show that a Kelvin-Helmholtz-type unstable wave develops at the interface for wavelengths longer than a critical value. Short-wavelength oscillations of the interface are damped due to... 

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

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

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 efficiency of proton exchange membrane (PEM) fuel cell is straightly correlated to the bipolar plate design and fluid channel arrangements. Higher produced energy can be attained by optimal design of type, size, or patterns of the channels. Previous researches showed that the bipolar plate channel design has a considerable effect on reactant distribution uniformity as well as humidity control in PEM fuel cells. This paper concentrates on enhancements in the fuel cell performance by optimization of bipolar plate design and channels configurations. A numerical model of flow distribution based on Navier-Stokes equations using individual computer code is presented. The results gained from... 

The present study has been conducted to investigate the effect of inlet flow regimes on vertical velocity and concentration profiles of the turbidity current. Experiments carried out on small scale channel and Kaolin with a density of 2.65 kg/m3 and a mean particle diameter of 4.5 μm, was used to generate the currents. Velocities were obtained at a rate of 10 Hz using an acoustic Doppler velocimeter. Using the siphon sampling approach, the acoustic backscattering intensity was calibrated and concentration was determined using a backscattering analysis. Results showed that the ratio of maximum velocity height to the height of the current was around 0.4-0.5 in sub- and super-critical regimes.... 

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

A comprehensive set of experiments was performed on weirs placed obliquely in a rectangular open channel. The main objectives were to investigate the effect of the weir angle relative to the approach flow direction on the behaviour and the hydraulic flow characteristics and to examine methods such as using upstream guide vanes and inclined aprons to increase the effective oblique weir length. Results indicate that by increasing the oblique angle, the effective length of the oblique weir increases significantly. Relationships are provided for the threshold submergence and the effective length of the oblique weir. Design recommendations are also presented for the use of guide vanes and... 

High bed-shear stress resulting from secondary flows and velocity increase at a channel junction cause local erosion and deposition in natural rivers and earth-lined channels. Herein, the scour patterns at the junction of two loose bed channels were studied under clear-water conditions. The main dimensionless variables are the angle between the two approach flow branches, the discharge and width ratios of the tributary to the downstream channel branches, and the ratio of the mean downstream velocity to the threshold velocity. The results indicate that the position of the maximum scour depth temporally moves to the outer wall and upstream to the main channel, as affected by the dimensionless... 

We present a novel unstructured coarse grid generation technique based on vorticity for upscaling two-phase flow in permeable media. In the technique, the fineness of the gridblocks throughout the domain is determined by vorticity distribution such that where the larger is the vorticity at a region, the finer are the gridblocks at that region. Vorticity is obtained from single-phase flow on original fine grid, and is utilized to generate a background grid which stores spacing parameter, and is used to steer generation of triangular and finally Voronoi grids. This technique is applied to two channelized and heterogeneous models and two-phase flow simulations are performed on the generated... 

The effect of an intermediate standing baffle on the flow structure in a rectangular open channel has been investigated by a three-dimensional acoustic Doppler velocimeter. Investigation of time-averaged velocity profiles at different streamwise positions reveals that the approach flow is fully developed upstream of the baffle. By analysing the space-averaged power spectra of streamwise velocity, a peak structure was observed in the upstream baffle region. Downstream of the baffle this peak structure has been alleviated by the baffle. The same analysis for the vertical component indicates the existence of a peak structure both up- and downstream of the baffle. Consequently, a baffle affects... 

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

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

Side weirs are hydraulic control structures widely used in irrigation, drainage networks and waste water treatment plants. These structures can be used for adjusting and diverting of flow with minimum energy loss. In spite of many studies were carried out on rectangular side weirs, the studies on oblique and labyrinth side weirs are scarce. In this study, based on the experimental data from more than 210 laboratory tests and through using the multivariable nonlinear partial least square (PLS) method, two nonlinear equations are presented for discharge coefficient CM of triangular labyrinth side weirs with one and two cycles. The obtained empirical equations relating CM with the relevant... 

The capability of backward extrusion (BE) method was assessed to achieve modified structures in AZ80 magnesium alloy. At first, 3D-Deform was employed to simulate the deformation flow through the deformed cup which gives an evidence from the flow behavior of the material. The material was processed via BE method at various temperatures of 250, 350, and 450 °C. Metallographic investigations were conducted in three different regions of the BE-processed cup (wall, bottom, and flow channel). The main feature observed at the wall of the BE cup was the presence of mechanical twins, the frequency of which was reduced by raising the process temperature. The flow localization in the form of shear... 

In this work, the Chebyshev collocation spectral lattice Boltzmann method is implemented in the generalized curvilinear coordinates to provide an accurate and efficient low-speed LB-based flow solver to be capable of handling curved geometries with non-uniform grids. The low-speed form of the D2Q9 and D3Q19 lattice Boltzmann equations is transformed into the generalized curvilinear coordinates and then the spatial derivatives in the resulting equations are discretized by using the Chebyshev collocation spectral method and the temporal term is discretized with the fourth-order Runge–Kutta scheme to provide an accurate and efficient low-speed flow solver. All boundary conditions are... 

In this work, a high-order weighted essentially nonoscillatory (WENO) finite-difference lattice Boltzmann method (WENOLBM) is developed and assessed for an accurate simulation of incompressible flows. To handle curved geometries with nonuniform grids, the incompressible form of the discrete Boltzmann equation with the Bhatnagar-Gross-Krook (BGK) approximation is transformed into the generalized curvilinear coordinates and the spatial derivatives of the resulting lattice Boltzmann equation in the computational plane are solved using the fifth-order WENO scheme. The first-order implicit-explicit Runge-Kutta scheme and also the fourth-order Runge-Kutta explicit time integrating scheme are... 

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