Sharif Digital Repository

The sloshing phenomena in a partially filled tank can affect its stability. Modifications of tank instability due to the movement of the tank carrier, are key design points for the stability of a carrier. Even though the sloshing phenomenon has already been investigated using the BEM-FDM technique, the research in this paper covers this phenomenon in a porous media, which is new in 2-D coordinates. For this purpose, a Laplace equation has been used for potential flow, and kinematic and dynamic boundary conditions have been applied to the free surface. Also, a formulation has been developed for a free surface in porous media. BEM has been used for solving the governing equation and FDM... 

In non-orthogonal multiple access (NOMA) systems, serving multiple users in shared resource blocks can allow untrusted users to overhear the messages of other users. In this context, we study a network consisting of a base station (BS), a near user and a far user, where the latter attempts to overhear the message of the former. The near user is a full-duplex (FD) node that can also act as a relay. Two operating scenarios are considered: 1) friendly jammer (FJ), where the FD node broadcasts noise for degrading the channel between the BS and the far user, while receiving data from the BS; and 2) friendly jammer relay (FJR), where, in addition to degrading the channel between the BS and the far... 

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

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

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

Shape memory polymers (SMPs) are a class of smart materials which can recover their shape even after many shape changes in application of an external stimulus. In this paper, flexural behavior of a composite beam, constructed of a corrugated part filled with SMPs, is studied. This composite beam is applicable in sensor and actuator applications. Since the corrugated profiles display higher stiffness-to-mass ratio in the transverse to the corrugation direction, the beams with a corrugated part along the transverse direction are stiffer than ones with a corrugated part along the length. Employing a developed constitutive model for SMPs and the Euler–Bernoulli beam theory, the behavior of the... 

The particle shape is an important factor playing critical role in evaluation of the interactions between particles in high-concentration particle-fluid flows. In this paper, the well-known multisphere (MS) approximation approach and the novel rolling resistance approach are utilized to examine their performance in order to simplify the generalized shaped particle’s interactions within the framework of discrete element method (DEM) and computational fluid dynamics (CFD). The performance of two approaches are compared with the perfect particle’s shape geometry, for the limited cases of cubic-shaped and disk-shaped particle flows in a horizontal well drilling process as a reference scenario.... 

The particle shape is an important factor playing critical role in evaluation of the interactions between particles in high-concentration particle-fluid flows. In this paper, the well-known multisphere (MS) approximation approach and the novel rolling resistance approach are utilized to examine their performance in order to simplify the generalized shaped particle’s interactions within the framework of discrete element method (DEM) and computational fluid dynamics (CFD). The performance of two approaches are compared with the perfect particle’s shape geometry, for the limited cases of cubic-shaped and disk-shaped particle flows in a horizontal well drilling process as a reference scenario.... 

We investigate the effect of particle shape on the transportation mechanism in well-drilling using a three-dimensional model that couples computational fluid dynamics (CFD) with the discrete element method (DEM). This numerical method allows us to incorporate the fluid-particle interactions (drag force, contact force, Saffman lift force, Magnus lift force, buoyancy force) using momentum exchange and the non-Newtonian behavior of the fluid. The interactions of particle-particle, particle-wall, and particle-drill pipe are taken into account with the Hertz-Mindlin model. We compare the transport of spheres with non-spherical particles (non-smooth sphere, disc, and cubic) constructed via the... 

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

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

Proppants transport is an advanced technique to improve the hydraulic fracture phenomenon, in order to promote the versatility of gas/oil reservoirs. A numerical simulation of proppants transport at both hydraulic fracture (HF) and natural fracture (NF) intersection is performed to provide a better understanding of key factors which cause, or contribute to proppants transport in HF–NF intersection. Computational fluid dynamics (CFD) in association with discrete element method (DEM) is used to model the complex interactions between proppant particles, host fluid medium and fractured walls. The effect of non-spherical geometry of particles is considered in this model, using the multi-sphere... 

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

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

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

Numerical approximation based on different forms of the governing partial differential equation can lead to significantly different results for two-phase flow in porous media. Selecting the proper primary variables is a critical step in efficiently modeling the highly nonlinear problem of multiphase subsurface flow. A comparison of various forms of numerical approximations for two-phase flow equations is performed in this work. Three forms of equations including the pressure-based, mixed pressure-saturation and modified pressure-saturation are examined. Each of these three highly nonlinear formulations is approximated using finite difference method and is linearized using both Picard and... 

A two-dimensional Boussinesq-type model is presented accurate to O(μ)6 , μ = h0/l0, in dispersion and all consequential order for non-linearity with arbitrary bottom boundary, where h0 is the water depth and l0 is the characteristic wave length. The mathematical formulation is an extension of (4,4) the Padé approximant to include varying bottom boundary in two horizontal dimensions. A higher order perturbation method is used for mathematical derivation of the presented model. A two horizontal dimension numerical model is developed based on derived equations using the Finite Difference Method in higher-order scheme for time and space. The numerical wave model is verified successfully in... 

A two-dimensional depth-integrated numerical model is developed using a fourth-order Boussinesq approximation for an arbitrary time-variable bottom boundary and is applied for submarine-landslide-generated waves. The mathematical formulation of model is an extension of (4,4) Padé approximant for moving bottom boundary. The mathematical formulations are derived based on a higher-order perturbation analysis using the expanded form of velocity components. A sixth-order multi-step finite difference method is applied for spatial discretization and a sixth-order Runge-Kutta method is applied for temporal discretization of the higher-order depth-integrated governing equations and boundary... 

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