Sharif Digital Repository

This paper presents a model to analyze pull-in phenomenon and dynamics of multi layer microplates using coupled finite element and finite difference methods. Firstorder shear deformation theory is used to model dynamical system using finite element method, while Finite difference method is applied to solve the nonlinear Reynolds equation of squeeze film damping. Using this model, Pull-in analysis of single layer and multi layer microplates are studied. The results of pull-in analysis are in good agreement with literature. Validating our model by pull-in results, an algorithm is presented to study dynamics of microplates. These simulations have many applications in designing multi layer... 

In this paper, two very simple and accurate equations are introduced that enable geotechnical engineers to determine the undrained ultimate bearing capacity of shallow strip foundations resting on two-layered clays. Results of the analyses show that in the case of weak-overstrong clay, usually general shear failure occurs, and in the case of strong-over-weak clay, the soil will most likely experience punching failure with a large amount of plastic settlement before it reaches its ultimate bearing capacity. This study applied the finite-difference method using a computer program to evaluate the undrained bearing capacity for all cases analyzed. The Mohr-Coulomb failure criterion was used for... 

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is extended and applied to accurately simulate two-phase liquid-vapor flows with high density ratios. Herein, the He-Shan-Doolen-type lattice Boltzmann multiphase model is used and the spatial derivatives in the resulting equations are discretized by using the fourth-order compact finite-difference scheme and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient two-phase flow solver. A high-order spectral-type low-pass compact nonlinear filter is used to regularize the numerical solution and remove spurious waves generated by flow nonlinearities in smooth regions... 

Large-stencil schemes which their spectral properties are acceptable in the vicinity of ω = π are analyzed for the first time. A machine independent model for evaluating the efficiency of generalized time-marching finite-difference algorithms over periodic domains is developed. This model which is based on operation count reveals that for small values of Total Computational Cost(TCC), the previous low-order small-stencil schemes are more efficient while for moderate TCC, the efficiency of optimized large-stencil schemes abruptly increases. This important result is the motivation for developing optimized large-stencil schemes. The current schemes are successfully implemented in a full... 

SnO2-CuO is one the most promising systems for detection of detrimental H2S gas. Although previous experimental research has suggested a sulfidation reaction to explain selectivity toward H2S, little is known about the origin of change of electrical response of this system by changing the H2S gas concentration. In this study the relation between sensing response of continuous SnO2-CuO bilayer thin film and H2S gas concentration is computed based on changeability of chemical composition of covellite CuxS. For this purpose, chemical activity of sulfur as a function of atomic fraction in covellite copper sulfide is estimated using Gibbs energies of formation and chemical thermodynamics. By... 

The formation damage (FD) caused by the invasion of drilling fluid severely affects reservoir performance during production. Most of the published research studies which address this type of FD have been carried out at the core or field scale. Thus, the main aim of the paper is to investigate the pore-scale mechanisms of FD induced by drilling fluids and their control with silica nanoparticles (NPs) using a microfluidic approach. The proper identification of the mechanisms of FD can lead to the proper selection of NP type and concentration as well as a suitable method to remediate FD. The micromodel was designed in a way to closely simulate the cross-flow at the wellbore surface. A... 

In this study, the numerical simulation of the fluid flow and acoustic field of a supersonic jet is performed by using high-order discretization and the vorticity confinement (VC) method on coarse grids. The three-dimensional Navier-Stokes equations are considered in the generalized curvilinear coordinate system and the high-order compact finite-difference scheme is applied for the space discretization, and the time integration is performed by the fourth-order Runge-Kutta scheme. A low-pass high-order filter is applied to stabilize the numerical solution. The non-reflecting boundary conditions are adopted for all the free boundaries, and the Kirchhoff surface integration is utilized to... 

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

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