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

Mixed-convection flow of nanofluids inside a vertical rectangular channel partially filled with open-cell metal foam and subject to a constant wall-heat flux was investigated experimentally and numerically. Al2O3-water nanofluids with different concentrations were prepared and their stability was examined using UV-Vis spectroscopy. Dynamic light scattering method was used to determine particle size distribution of the nanofluid feedstock. The outlet temperature and pressure drop were measured for different nanofluid flow rates (i.e., Reynolds number values). In the numerical section, a two-dimensional volume-averaged form of the governing equations was used. The velocity and temperature... 

Dynamic modeling of a double-pipe heat exchanger is the subject of the current study. The basis of this study is the same velocity of vapor and liquid phases or, in other words, homogeneous phase, in the annulus part of the exchanger. The model can predict the temperature and vapor quality along the axial pipe from the pipe inlet up to a distance where steady state conditions are achieved. The simulation is conducted for two modes of co- and counter-flow in a one dimensional transient system. The physical properties of water are estimated from empirical correlation and a saturated vapor table with cubic spline interpolation. The exchanger model, which is a set of Ordinary Differential... 

In this study, the nonlinear partial differential equations governing two phase flow through porous media are solved using two different methods, namely, finite difference and finite element. The capillary pressure term is considered in the mathematical model. The numerical results on a 2-D test case are then compared with the experimental drainage process and water flooding performed on a glass type micromodel. Based on the obtained results, finite difference technique needs less computational time for solving governing equations of two phase flow, but findings of this method show less agreement with the experimental data. The finite element scheme was found to be more adequate and its... 

Fast fault detection (FD) and reconfiguration is necessary for fault tolerant power electronic converters in safety critical applications to prevent further damage and to make the continuity of service possible. The aim of this study is to minimize the number of the used additional voltage sensors in a fault tolerant three-phase converter. In this paper, first a practical implementation of a very fast FD scheme with reduced sensor number is discussed. Then, an optimization in this scheme is also presented to decrease the detection time. For FD, special time and voltage criterion are applied to observe the error in the estimated phase-to-phase voltages for a specific period of time. The... 

Fault detection (FD) in power electronic converters is necessary in embedded and safety critical applications to prevent further damage. Fast FD is a mandatory step in order to make a suitable response to a fault in one of the semiconductor devices. The aim of this study is to present a fast yet robust method for fault diagnosis in nonisolated dc-dc converters. FD is based on time and current criteria which observe the slope of the inductor current over the time. It is realized by using a hybrid structure via coordinated operation of two FD subsystems that work in parallel. No additional sensors, which increase system cost and reduce reliability, are required for this detection method. For... 

To predict the amount of different phases in gray cast iron by a finite difference model (FDM) on the basis of cooling rate (R), the volume fractions of total γphase, graphite, and cementite were calculated. The results of phase composition were evaluated to find a proper correlation with cooling rate. More trials were carried out to find a good correlation between the hardness and phase composition. New proposed formulas show that the hardness of gray cast iron decreases as the amount of graphite phase increases, and increases as the amount of cementite increases. These formulas are developed to correlate the phase volume fraction to hardness. The results are compared with experimental data... 

In this study, dynamic pull-in instability and snap-through buckling of initially curved microbeams are investigated. The microbeams are actuated by suddenly applied electrostatic force. A finite element model is developed to discretize the governing equations, and Newmark time discretization is employed to solve the discretized equations. The static pull-in behavior is investigated to validate the model. The results of the finite element model are compared with finite difference solutions and their convergence is examined. In addition, the influence of different parameters on dynamic pull-in instability and snap-through buckling is explored  

This paper investigates implementation of upwind compact implicit and explicit high-order finite difference schemes for solution of the level set equation. The upwind compact implicit and explicit high-order finite difference schemes are well-known techniques to descritize spatial derivatives for convection term in hyperbolic equations. Applying of upwind high-order schemes on the level set equation leads to less error and CPU time reduction compared to essential non-oscillatory (ENO), weighted essential non-oscillatory schemes (WENO), and even different particle level set methods. The results indicate the error based on area loss decreases drastically with applying high-order upwind,... 

The aim of this paper is to correlate interfacial heat transfer coefficient (IHTC) to applied external pressure, in which IHTC at the interface between A356 aluminum alloy and metallic mold during the solidification of casting under different pressures were obtained using the inverse heat conduction problem (IHCP) method. The method covers the expedient of comparing theoretical and experimental thermal histories. Temperature profiles obtained from thermocouples were used in a finite difference heat flow program to estimate the transient heat transfer coefficients. The new simple formula was presented for correlation between external pressure and heat transfer coefficient. Acceptable... 

In this study, a coupled numerical approach based on Lattice Boltzmann Method (LBM) and Discrete Element Method (DEM) is employed for 2D simulation of fluid flow in porous media comprising of movable circular particles. The developed model is used for simulation of sand production which is one of the important problems in petroleum industry. The numerical tool has proved to have the capability of investigating the mechanisms involved in sand production problem. The results show that the rate of sand production is strongly affected by flow rate and confining pressure  

The present work reports the outcome of a comprehensive parametric study on mixed electroosmotically and pressure-driven flow in slit microchannels with constant wall heat fluxes. Special attention is given to disclose the applicability ranges of usual assumptions in simplified analyses. The governing equations for fully developed conditions are first made dimensionless and then solved by means of an implicit finite difference method. The results reveal that the assumption of constant thermophysical properties does not leadto significant errors in practical applications. Although the Debye-Huckel linearization may successfully be used to evaluate velocity profiles up to the zeta potentials... 

Normal stress has some role in the deformation analysis of hydroforming processes. In this study, analytical modeling is pursued to evaluate the effect of normal stress on the hydro-mechanical deep drawing (HDD) process. Analyses are carried out for axisymmetric elements of the formed cup-shaped part for increments of the punch travel. The formulations are obtained using mechanical and geometrical relations and the finite difference method, thereby being solved by proper numerical algorithms. Furthermore, in the present work, part thickness is variable, the loading and straining are non-proportional, and bending/unbending effects over the part curvature are considered. The results show that... 

The goal of this study is to evaluate the effects of proportional loading, plane stress, and constant thickness assumptions on hydro-mechanical deep drawing (HDD) by developing analytical models. The main model includes no simplifying assumption, and then each of the mentioned assumptions is considered in a specific model. The interrelationships between geometrical and mechanical variables are obtained in the finite difference form based on the incremental strain theory, thereby being solved by Broyden algorithm. Published experimental and FE results are used for evaluation of the results obtained in the present work. The results of models under proportional loading, plane stress, and... 

Among major components of LiBr-H 2O absorption chillers is the absorber, which has a direct effect on the chillier size and whose characteristics have significant effects on the overall efficiency of absorption machines. In this article, heat and mass transfer process in absorption of refrigerant vapor into a lithium bromide solution of water-cooled incline plate absorber in the Reynolds number range of 5