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

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

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

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

Karun Bridge is part of a national highway project underway in southern Iran. The bridge is a 336 m steel arch structure that weighs more than 2500 tons and crosses the Karun River 270 meters above the river valley. The bridge will be the largest suspended bridge in the Middle East. Cantilever construction is to be used to construct the bridge from both sides simultaneously. The weight of the bridge is to be carried by two abutments and four piers that are anchored to exposed weathered rock mass. Loads on the abutments and piers include significant cantilever loads and moments during bridge erection and wind loads. Daily rapid temperature changes also will impart significant thermal loads on... 

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 obtain the wave and boundary layer characteristics such as bed shear stress. A mixing-length theory is used for turbulence modeling. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 by ASME  

In this paper, the numerical errors associated with the finite difference solutions of two-dimensional advection-dispersion equation with linear sorption are obtained from a Taylor analysis and are removed from numerical solution. The error expressions are based on a general form of the corresponding difference equation. The variation of these numerical truncation errors is presented as a function of Peclet and Courant numbers in X and Y direction, a Sink/Source dimensionless number and new form of Peclet and Courant numbers in X-Y plane. It is shown that the Crank-Nicolson method is the most accurate scheme based on the truncation error analysis. The effects of these truncation errors on... 

The numerical errors associated with explicit upstream finite difference solutions of two-dimensional advection - Dispersion equation with linear sorption are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation. The numerical truncation errors are defined using Peclet and Courant numbers in the X and Y direction, a sink/source dimensionless number and new Peclet and Courant numbers in the XY plane. The effects of these truncation errors on the explicit solution of a two-dimensional advection-dispersion equation with a first-order reaction or degradation are demonstrated by comparison with an analytical solution in... 

The transient heat transfer analysis of a layer has been studied much less than the steady state. However, transient temperature distribution resulted from including radiation and conduction simultaneously, is significantly different from those obtained by considering conduction alone. In order to include the effect of radiation heat transfer, we must insert the gradient of radiative flux in the energy equation. For this purpose, a variety of multi-flux methods have been suggested. A simplified procedure is the two-flux method, which is the one used in the present paper. This paper is focused on one-dimensional transient heat transfer of a layer using Finite Difference Method. To this end, a... 

In this paper, the development of a neutron noise simulator for hexagonal-structured reactor cores using both the forward and the adjoint methods is reported. The spatial discretisation of both 2-D 2-group static and dynamic equations is based on a developed box-scheme finite difference method for hexagonal mesh boxes. Using the power iteration method for the static calculations, the 2-group neutron flux and its adjoint with the corresponding eigenvalues are obtained by the developed static simulator. The results are then benchmarked against the well-known CITATION computer code. The dynamic calculations are performed in the frequency domain which leads to discarding of the time... 

Mesh point positions in Finite Difference Method (FDM) of discretization for the neutron diffusion equation can remarkably affect the averaged neutron fluxes as well as the effective multiplication factor. In this study, by aid of improving the mesh point positions, an enhanced finite difference scheme for the neutron diffusion equation is proposed based on the neutron importance function. In order to determine the neutron importance function, the adjoint (backward) neutron diffusion calculations are performed in the same procedure as for the forward calculations. Considering the neutron importance function, the mesh points can be improved through the entire reactor core. Accordingly, in... 

Flow movement in unsaturated soil can be expressed by a partial differential equation, named Richards equation. The objective of this study is the finding of an appropriate implicit numerical solution for head based Richards equation. Some of the well known finite difference schemes (fully implicit, Crank Nicolson and Runge-Kutta) have been utilized in this study. In addition, the effects of different approximations of moisture capacity function, convergence criteria and time stepping methods were evaluated. Two different infiltration problems were solved to investigate the performance of different schemes. These problems include of vertical water flow in a wet and very dry soils. The... 

When cells exposed to an electric field, localized changes in the distribution of the electric field will be induced and these changes in turn lead to electrical stresses on cell surface. The electrical stresses play a key role in the cell membrane structural changes which leads to important phenomena like hydrophilic pores formation on the cell membrane resulting in the cell permeability. In this work, protoplast cell interaction with direct current (DC) electric field is investigated. The electrical stresses acted on the cell membrane in the presence of electric field are investigated numerically by a modified finite difference method, fast Immersed Interface Method (IIM). Exact solution... 

The numerical solution to 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 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 simultaneously...