Sharif Digital Repository

In this paper, a double-surface plasticity model, based on a combination of a convex yield surface consisting of a failure envelope, such as a Mohr-Coulomb yield surface and, a hardening cap model, is developed for the nonlinear behaviour of powder materials in the concept of a generalized plasticity formulation for the description of cyclic loading. This model reflects the yielding, frictional and densification characteristics of powder along with strain and geometrical hardening which occur during the compaction process. The solution yields details on the powder displacement from which it is possible to establish the stress state in the powder and the densification is derived from... 

A very simple and computationally low cost numerical algorithm is developed to generate a quasi-structured data structure for an unstructured grid. To achieve this purpose, the data structure in the matrices of an unstructured grid is classified to address the element layers and node lines in the computational domain. In this regard, elements and nodes of the unstructured grid are renumbered in a directional ordering-based strategy. The elements and nodes arrangement in each layer and line is accomplished in a unique direction around an interior object of the grid either clockwise or counterclockwise. Furthermore, a new searching scheme is introduced which guarantees a quick search inside... 

This paper presents a numerical algorithm for BIBO stability testing of a certain class of the so-called fractional-delay systems. The characteristic function of the systems under consideration is a multi-valued function of the Laplace variable s which is defined on a Riemann surface with finite number of Riemann sheets where the origin is a branch point. The stability analysis of such systems is not straightforward because there is no universally applicable analytical method to find the roots of the characteristic equation on the right half-plane of the first Riemann sheet. The proposed method is based on the Rouche's theorem which provides the number of the zeros of a given function in a... 

The analytical formulation of the fully coupled and the uncoupled rolling contact mechanics problems for all possible stick/slip regimes are derived in Part I (Alinia et al., 2013). In this part, we focus on the numerical algorithm, the iteration procedure and the numerical results. The coupled and the uncoupled solutions corresponding to each of the assumed stick/slip regimes are provided. The uncoupled solution provides an acceptable approximation to the problem for small coefficient of friction values. However, for high values of the coefficient of friction the problem should be solved in fully coupled form. In addition, the effect of several parameters such as the stiffness ratio, the... 

The objective of this study is to develop a general finite element formulation associated with an incremental adaptive procedure which established for thermoviscoelastic stress analysis of orthotropic cylinders made of polymers. This paper concerned with development of a numerical algorithm for the solution of the quasistatic initial/boundary value problems involving the linear viscoelastic media with thermal and mechanical deformations. The viscoelastic constitutive equations, represented in an integral form and involving relaxation functions, are transformed into an incremental algebraic relation. An incremental relaxation is then developed for the finite element formulation to deal with... 

Problems regarding simulation of field emitter array (FEA) electron guns are discussed. A simple method is proposed to significantly reduce computational requirements such as computation power, system memory, and time of FEA electron gun simulation and modeling. The method can be applied to any numerical solver regardless of its meshing technique. In order to extract field emission parameter from any experimental cathode I-V curve, a partly numerical algorithm, which uses the presented truncation method at the heart of its solver, is proposed. The proposed method and algorithm are applied to a number of examples, including a double-gated FEA problem, and its effectiveness in terms of error... 

In this paper, we consider the bacterial point-to-point communication problem with one transmitter and one receiver by considering the ligand receptor binding process. The most commonly investigated signalling model, referred to as the Level Scenario (LS), uses one type of a molecule with different concentration levels for signaling. An alternative approach is to employ multiple types of molecules with a single concentration level, referred to as the Type Scenario (TS). We investigate the trade-offs between the two scenarios for the ligand receptor from the capacity point of view. For this purpose, we evaluate the capacity using numerical algorithms. Moreover, we derive an upper bound on the... 

Abstract Elastic optical networking has recently been proposed for use in optical transport networks to cope with increasingly heterogeneous and dynamic demand patterns. In this paper, we study the blocking performance of a multi-class elastic optical link for which a demand needs to be allocated a contiguous subset of the entire spectrum. This problem is different than the well-known blocking problem in multi-class multi-server loss systems due to the contiguous allocation constraint. We first propose a non-work-conserving aligned spectrum allocation policy which is shown to outperform the conventional first fit-based work-conserving allocation policy without alignment. Subsequently, for... 

This work focuses on transient thermal behavior of radial fins of rectangular, triangular and hyperbolic profiles with temperature-dependent properties. A hybrid numerical algorithm which combines differential transformation (DTM) and finite difference (FDM) methods is utilized to theoretically study the present problem. DTM and FDM are applied to the time and space domains of the problem, respectively. The accuracy of this method solution is checked against the numerical solution. Then, the effects of some applicable parameters were studied comparatively. Since a broad range of governing parameters are investigated, the results could be useful in a number of industrial and engineering... 

Today, turbomachinery designers are trying to attain zero clearances for optimizing performance and efficiency of machine by designing advanced compliant-segmented seals. This trend involves secondary issues such as rotor-to-stator contact interaction which is to be treated with more technical knowledge in rotordynamic design and condition monitoring programs. The aim of this paper is to gain insight into dynamic interaction between a flexible rotor and a set of compliant arc-shaped segments as stator system. Specifically, the unbalanced-induced forward rubbing response is investigated in tight clearance condition during resonance-passing situations. A Jeffcott rotor modeling is used for... 

Compressible turbulent flow through the abrupt enlargement in pipes is studied numerically by means of Advection Upstream Splitting Method (AUSM+-up). In low Mach numbers, a pressure correction equation of elliptic type is derived. This equation is compatible with the nature of governing equations and retrieves hyperbolic characteristic at higher Mach numbers. It is shown that the proposed numerical algorithm is computationally more efficient than the preconditioned density-based methods. The flow parameters such as reattachment length, pressure loss coefficient and wall shear stress are predicted. It is found that the loss coefficient of the compressible flow rises drastically with... 

In this study, a numerical model of unsteady potential flow around submerged marine propellers has been developed. The boundary element approach in combination with time stepping method to model free wake dynamics has been implemented. An important feature of this method in the simulation of pressure-dominant problems is a proper balance between time and accuracy in the numerical process. Another advantage of time stepping method is that there is no need to define wake geometry before modeling. Due to inherent instability of boundary integral equations, a smoothing function to damp the effect of singularities is imposed to the solution. The main innovative idea of this work is that the... 

In this paper, we propose an analytical approach to compute the probability of connectivity for one-dimensional ad hoc networks. The proposed analysis gives the exact probability of connectivity for an arbitrary distribution of nodes, provided that nodes are independently and identically distributed. We conduct separate analyses for two cases; in the first case, the number of nodes varies by time under a stationary distribution and in the second case, there is a fixed (known) number of nodes in the network. Using the approaches presented in this work, we are able to derive closed-form formulas for the probability of connectivity for some spatial distributions, while for more complicated... 

The present study focuses on the numerical simulation of cavitation around the NACA 0015. The unsteady behaviors of cavitation which have worthwhile applications are investigated. The cavitation patterns, velocity fields and frequency of the cavitating flow around hydrofoil is obtained. For multi phase simulation, single-fluid Navier-Stokes equations, along with a volume fraction transport equation, are employed. The bubble dynamics model is utilized to simulate phase change. SIMPLE algorithm is used for velocity and pressure computations. For discretization of equations the finite-volume approach written in body fitted curvilinear coordinates, on collocated grid, is used. In this study,... 

In this paper, an extended finite element method is employed to simulate the presence of discontinuities caused by the contact surface. In X-FEM, the need for mesh adaption to interface is neglected and the process is accomplished by partitioning the domain with some triangular sub-elements whose Gauss points are used for integration of the elements. The modified level set technique and the Heaviside enrichment function are employed to approximate the discontinuous displacement field of elements located on the contact surface. The penalty method is used to impose the contact constraints and establish the non-penetration condition. An efficient numerical algorithm is employed to model the... 

Size distribution of aerosol particles is prevalently obtained through electrical characterization techniques and study of charged particles' dynamics in the presence of electric field. In this work, a wire-cylinder corona charger is presented, redesigned, and aero dynamically optimized. An initial 2D axisymmetric geometry of the charger was employed for the simulations using the Computational Fluid Dynamics (CFD) commercial code FLUENT 6.3.26. Through successive attempts, a new geometry was obtained by streamlining the walls to eliminate the undesired vortices produced in the flow field of the previous ones. The process optimized the charger by minimizing losses and dilutions of the... 

In this manuscript, an advanced damage-plasticity model is utilized to simulate response of concrete structures under high-velocity impact. Due to presence of large deformations, it is necessary to incorporate the damage-plasticity model into the finite deformation framework. In high-velocity impact, severe numerical problems could be encountered while updating stress values due to the sudden increase in their levels, especially for the complicated material model used in this study. To overcome these obstacles, an enhancement is made in the nonlinear system of equations of stress- updating procedure. In addition, an adaptive multi-step numerical algorithm is introduced which improves the... 

This article aims at developing a semi-analytic approach for studying the free and forced vibrations of a pseudoelastically behaving shape memory alloy beam. Based on the Euler–Bernoulli beam theory, equations of motion were derived through Hamilton principle, and the obtained partial differential equations were decomposed by applying the Galerkin approach and were solved using Newmark integration method. A three-dimensional phenomenological model of shape memory alloy, which is capable of identifying the main properties of the shape memory alloy, was employed to model the behavior of the shape memory alloy beam. A closed-form numerical algorithm was introduced to simulate the governing... 

This research proposes a large deformation, time-dependent viscoplastic constitutive relationship to enhance the prediction of the compaction degree of asphalt concrete materials under laboratory and field conditions. A large-deformation thermodynamic-based framework is presented. The Helmholtz free energy and rate of energy dissipation functions were assumed to derive rate-dependent constitutive relationships to relate multi-axial state of stresses to the recoverable and non-recoverable deformation response of asphalt concrete during compaction. A straightforward method that allows the calibration of the proposed model against laboratory compaction data (e.g., data from Superpave Gyratory...