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

The numerical solution of 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 the 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... 

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

This study is concerned with numerical modeling of viscous surface wave motion using boundary element method (BEM). The equations of motion for thin boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM-finite difference technique is used to obtain the viscosity-related quantities such as wave damping rate, shear stress, and velocity distribution inside the boundary layer. The technique is presented for standing surface wave motion. An excellent agreement is obtained between the numerical predictions and the previous results. The extension to other free surface problems is straightforward. © 2005 Elsevier Ltd. All rights reserved  

A numerical model for simulation of liquid/gas phase flow during mold filling is presented. The incompressible Navier-Stokes equations are discretized on a fixed Cartesian mesh with finite difference method. The fractional-step scheme is employed for enforcing incompressibility constraint. The free surface effects are calculated using the volume of fluid method based on the piecewise-linear interface reconstruction and split Lagrangian advection of volume fraction field. Adding limited compressibility to the gas phase led to improvement in convergence rate of Poisson equation solver (about 2-fold). This new concept permits simulation of two-phase incompressible free surface flow during mold... 

Polymer films possess excellent optical properties, such as high transparency, and thermal characteristics, like low heat conductivity, as well as further polymer specific advantages. Consequently, polymer films have an outstanding potential for many solar applications. They are already used for encapsulation of photovoltaic (PV) cells, as convection barrier in solar collectors and as substrate or adhesive layers for glazing. In translucent polymers, energy can be transferred internally by radiation in addition to conduction. Since radiant propagation is very rapid, it can provide energy within the layer more quickly than diffusion by heat conduction. Thus, the transient thermal response of... 

The Adomian Decomposition Method (ADM) for solving the highly non-linear vorticity-stream function formulation of 2D incompressible Navier-Stokes equations has been implemented. The analysis is accompanied by numerical boundary conditions. Also, numerical simulation, using finite difference method (FDM), is performed for comparison purposes. The obtained results only for few terms of the expansion are presented. Because present software such as Mathematica/Maple can not calculate many terms (for example: up to 10 terms) of solution and then ADM approach of this problem is an open problem case  

A multiscale FD-KMC model has been developed to simulate the cyclic voltammetry test of a copper electrode in simple copper sulfate bath. In this coupled model, the FD code provides the cupric ion concentration on OHP for KMC code, while the KMC code provides the electrochemical properties of the copper electrode (surface activity and rate constants of redox reactions) as an input data for FD code. The changes in the electrode properties due to the atomistic phenomena (deposition dissolution and surface diffusion) have been studied for the present potentiodynamic system. The results showed that the CV process consists of some distinct stages, so that the electrode exhibits a specific... 

An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM-finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing-length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd  

A novel method is used to link a 2D kinetic Monte Carlo code to a 1D finite difference code to construct a more realistic and efficient tool for simulating various electrochemical processes. This multiscale model is able to simulate the long-scale mass transfer of electroactive species in bath along with electrode surface phenomena at atomic scale simultaneously. An embedded atom method (EAM) has been used to evaluate the barrier energies of diffusion and redox reactions on electrode surface. The FD code provides the ion concentration on OHP for KMC code, while the KMC code provides the surface activity and rate constants of redox reactions as an input data for FD code. The electrochemical... 

The present study employs Discrete Element Method (DEM) to establish a rheological model that relates the apparent viscosity of a granular material to shear rate, normal stress, and water saturation. In addition, a theoretical model was developed to determine water distribution and water-induced forces between particles for different saturations. The resulting forces were embedded in a 3D shear cell as a numerical rheometer, and a wet specimen was sheared between two walls. A power law rheological model was then obtained as a function of inertia number and saturation. It was found that up to a critical saturation, the apparent viscosity increased with saturation that was higher than that of... 

A finite difference-based numerical model simulating the cone penetration process in unsaturated sands is presented. Mohr–Coulomb model (MCM) with simple modifications and Sun model (SM) were implemented to capture the unsaturated sand behaviour. It was shown that the cone tip resistance values resulting from the two models were fairly comparable. Predicted cone tip resistance values in dry, saturated and unsaturated sands using MCM were validated by the results of field and calibration chamber tests. Sensitivity analyses were performed, and the influence of parameters including relative density, mean effective stress and apparent cohesion due to suction on the tip resistance was... 

We introduce two-dimensional space plus time (2D+1) structure and numerically investigate it using a developed multilayer simulation framework, for the first time. The new structure is consisting of crossed grating with time-varying permittivity which is inspired from1D+1. In this regard, we extend FourierModal Method (FMM) in a general approach for spacetime multilayer states. Our proposed framework is fast, robust, and powerful compared to various finite difference methods. We use the scattering matrix technique to develop the proposed spacetime simulation method for multilayer structures using a non-uniform stack of layers. The method is perfectly suitable to investigate the... 

Utilization of heat storage units in solar energy systems can resolve the challenge of fluctuation and uncertainty of the solar energy. Phase change materials (PCMs) are used as the storage media for solar energy storage systems. In this research, a system including of a solar collector and a PCM-based cascaded energy storage unit was numerically investigated. Air was used as the heat transfer fluid (HTF) and three paraffin-based materials (RT50, RT65, and RT80) were used as PCM for the energy storage unit. The investigated system mainly operates between 15 °C and 90 °C and considering different PCMs, the selected PCMs were appropriate. Paraffin-based PCMs also present acceptable thermal... 

A model of heat and mass transfer of two dimensional MHD micropolar fluid over a cone is constructed. Similarity transformation is adopted for the conversion of partial differential equations into ordinary differential equations have been linearized by employing the Newton's linearization technique and then new sets of equations are discretized using the finite difference method. The impact of non-dimensional parameters is further analyzed and the numerical results for profiles of velocity, temperature and concentration are expressed graphically and the results are discussed in detail. For the higher values of Dufour number, temperature field is enhanced graphically but show the opposite... 

The Computational fluid dynamics (CFD)–discrete element method (DEM) numerical simulation may be applied to predict the hydrodynamic behavior of dense particle–fluid flows. The main drawback of this simulation is the long computational time required owing to the large number of particles and the minute time-step required to maintain a stable solution. In this work, a new method to improve the efficiency and accuracy of CFD–DEM simulations is presented. The particle stiffness coefficient is used as a flexible parameter to improve the accuracy and efficiency of the model. The particle concentration distribution results are compared with experimental one’s to derive the optimum effective... 

A model of heat and mass transfer of two dimensional MHD micropolar fluid over a cone is constructed. Similarity transformation is adopted for the conversion of partial differential equations into ordinary differential equations have been linearized by employing the Newton's linearization technique and then new sets of equations are discretized using the finite difference method. The impact of non-dimensional parameters is further analyzed and the numerical results for profiles of velocity, temperature and concentration are expressed graphically and the results are discussed in detail. For the higher values of Dufour number, temperature field is enhanced graphically but show the opposite...