Sharif Digital Repository

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

Absorber is an important component in absorption machines and its characteristics have significant effects on the overall efficiency of absorption machines. This article reports a model of simultaneous heat and mass transfer process in absorption of refrigerant vapor into a lithium bromide solution of water - cooled vertical plate absorber in the Reynolds number range of 5 < Re < 150. The boundary layer assumptions were used for the transport of mass, momentum and energy equations and the fully implicit finite difference method was employed to solve the governing equations in the film flow. Dependence of lithium bromide aqueous properties to the temperature and concentration and film... 

Nonlinear ultrasound modeling is finding an increasing number of applications in both medical and industrial areas where due to high pressure amplitudes the effects of nonlinear propagation are no longer negligible. Taking nonlinear effects into account makes the ultrasound beam analysis more accurate in these applications. One of the most widely used nonlinear models for propagation of 3D diffractive sound beams in dissipative media is the KZK (Khokhlov, Kuznetsov, Zabolotskaya) parabolic nonlinear wave equation. Various numerical algorithms have been developed to solve the KZK equation. Generally, these algorithms fall into one of three main categories: frequency domain, time domain and... 

The characteristic boundary conditions are applied and assessed for the solution of incompressible inviscid flows. The two-dimensional incompressible Euler equations based on the artificial compressibility method are considered and then the characteristic boundary conditions are formulated in the generalized curvilinear coordinates and implemented on both the far-field and wall boundaries. A fourth-order compact finite-difference scheme is used to discretize the resulting system of equations. The solution methodology adopted is more suitable for this assessment because the Euler equations and the high-accurate numerical scheme applied are quite sensitive to the treatment of boundary... 

This study reports on the simulation of temperature distribution of human tooth under a laser beam based on non-Fourier models. The temperature in the tooth depth that directly results from the conduction heat transfer process is caused by the lengthy thermal relaxation time in the tooth layers. A detailed tooth composed of enamel, dentin, and pulp with unstructured shape, uneven boundaries, and realistic thicknesses was considered. A finite difference scheme was separately adopted to solve time-dependent equations in solid layers and soft tissue of the tooth. In this study, a dual-phase-lag non-Fourier heat conduction model was applied to evaluate temperature distribution induced by laser... 

The preconditioned characteristic boundary conditions based on the artificial compressibility (AC) method are implemented at artificial boundaries for the solution of two- and three-dimensional incompressible viscous flows in the generalized curvilinear coordinates. The compatibility equations and the corresponding characteristic variables (or the Riemann invariants) are mathematically derived and then applied as suitable boundary conditions in a high-order accurate incompressible flow solver. The spatial discretization of the resulting system of equations is carried out by the fourth-order compact finite-difference (FD) scheme. In the preconditioning applied here, the value of AC parameter... 

One of the most important problems facing structural engineers is the analysis of dynamic behavior of bridges subjected to moving vehicles. In addition, viscoelastic supports under bridges change their dynamic behavior under passing traffic loads. This paper presents how to model a bridge with viscoelastic supports and how the maximum dynamic stress of bridges changes during the passing of moving vehicles. Furthermore, this paper presents an algorithm to solve the governing equation of the bridge with viscoelastic supports as well as the equation of motion of a real European truck with different speeds, simultaneously. Using viscoelastic supports with appropriate characteristics can make a... 

In this study, the behavior of gavoshan dam was evaluated during construction and the first impounding. A two-dimensional (2D) numerical analysis was conducted based on a finite difference method on the largest cross-section of the dam using the results of instrument measurements and back analysis. These evaluations will be completed in the case that back analysis is carried out in order to control the degree of the accuracy and the level of confidence of the measured behavior since each of the measurements could be controlled by comparing it to the result obtained from the numerical model. Following that, by comparing the results of the numerical analysis with the measured values, it is... 

The present research is devoted to theoretical study of the pull-in performance of double-sided and paddle-type NEMS actuators fabricated from cylindrical nanowire operating in the Casimir regime and in the presence of the centrifugal force. D'Alembert's principle was used to transform the angular velocity into an equivalent static, centrifugal force. Using the couple stress theory, the constitutive equations of the actuators were derived. The equivalent boundary condition technique was applied to obtain the governing equation of the paddle-type actuator. Three distinct approaches, the Duan-Adomian Method (DAM), Finite Difference Method (FDM), and Lumped Parameter Model (LPM), were applied... 

This paper presents an analysis of the combined electro-osmotic and pressure-driven axial flows of viscoelastic fluids in a rectangular microchannel with arbitrary aspect ratios. The rheological behavior of the fluid is described by the complete form of Phan-Thien–Tanner (PTT) model with the Gordon–Schowalter convected derivative which covers the upper convected Maxwell, Johnson–Segalman and FENE-P models. Our numerical simulation is based on the computation of 2D Poisson–Boltzmann, Cauchy momentum and PTT constitutive equations. The solution of these governing nonlinear coupled set of equations is obtained by using the second-order central finite difference method in a non-uniform grid... 

We introduce a theoretical method for simultaneous measurement of refractive index and thickness of multilayer systems using Fourier domain optical coherence tomography (FD-OCT) without any auxiliary arrangement. The input data to the formalism are the FD-OCT measured optical path lengths (OPLs) and properly selected spectral components of FD-OCT interference spectrum. The outputs of the formalism can be affected significantly by uncertainty in measuring the OPLs. An optimization method is introduced to deal with the relatively large amount of uncertainty in measured OPLs and enhance the final results. Simulation result shows that by using the optimization method, indices can be extracted... 

We introduce a theoretical framework for simultaneous refractive index and thickness measurement of multilayer systems using the Fourier domain optical coherence tomography (FD-OCT) system without any previous information about the item under investigation. The input data to the new formalism are the FD-OCT measured optical path lengths and properly selected spectral components of the FD-OCT interference spectrum. No additional arrangement, reference reflector, or mechanical scanning is needed in this approach. Simulation results show that the accuracy of the extracted parameters depends on the index contrast of the sample while it is insensitive to the sample's thickness profile. For... 

In this paper, we aim to account for the partitioning of finite sized ions and electric double layer (EDL) overlapping effects on the electrostatics and hydrodynamics of soft nanofluidics by stablishing a modified Poisson-Boltzmann (MPB) equation enjoying mean field approach. The application of the present MPB equation enables us to describe the interaction between the steric effect and electrostatic repulsion of EDL ions due to permittivity difference of polyelectrolyte layer (PEL) and electrolyte solution. Utilizing the Debye-Hückel approximation pertinent to low surface potentials, we analytically derive the solutions of electric potential and velocity profiles of mixed electroosmotic and... 

A high-order compact finite-difference total Lagrangian method (CFDTLM) is developed and applied to nonlinear structural dynamic analysis. The two-dimensional simulation of thermo-elastodynamics is numerically performed in generalized curvilinear coordinates by taking into account the geometric and material nonlinearities. The spatial discretization is carried out by a fourth-order compact finite-difference scheme and an implicit second-order accurate dual time-stepping method is applied for the time integration. The accuracy and capability of the proposed solution methodology for the nonlinear structural analysis is investigated through simulating different static and dynamic benchmark... 

The characteristic boundary conditions are applied and assessed for the solution of incompressible inviscid flows. The two-dimensional incompressible Euler equations based on the artificial compressibility method are considered and then the characteristic boundary conditions are formulated in the generalized curvilinear coordinates and implemented on both the far-field and wall boundaries. A fourth-order compact finite-difference scheme is used to discretize the resulting system of equations. The solution methodology adopted is more suitable for this assessment because the Euler equations and the high-accurate numerical scheme applied are quite sensitive to the treatment of boundary... 

In this work, a preconditioned high-order weighted essentially non-oscillatory (WENO) finite-difference lattice Boltzmann method (WENO-LBM) is applied to deal with the incompressible turbulent flows. Two different turbulence models namely, the Spalart–Allmaras (SA) and k−ωSST models are used and applied in the solution method for this aim. The spatial derivatives of the two-dimensional (2D) preconditioned LB equation in the generalized curvilinear coordinates are discretized by using the fifth-order WENO finite-difference scheme and an implicit–explicit Runge–Kutta scheme is adopted for the time discretization. For the convective and diffusive terms of the turbulence transport equations, the... 

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

In this paper, a three-dimensional simulation of fully coupled multiphase fluid flow and heat transfer through deforming porous media is presented in the context of EFG mesh-less method. Spatial discretization of the system of governing equations is performed using EFG and a fully implicit finite difference scheme is employed for temporal discretization. Penalty method is used for imposition of essential boundary conditions. The developed numerical tool is employed to simulate two problems of nuclear waste disposal and CO2 sequestration in deep underground strata. The obtained results demonstrate the capability and robustness of the developed EFG code. © 2018 Elsevier Ltd  

In this paper, a fast and practical method is proposed for open-circuit (OC) fault diagnosis (FD) in a three-phase quasi-Z-source inverter (q-ZSI). Compared to the existing fast OC FD techniques in three-phase voltage-source inverters (VSIs), this method is more cost-effective since no ultra-fast processor or high-speed measurement is required. Additionally, the method is independent of the load condition. The proposed method is only applicable to Z-source family inverters and is based on observing the effect of shoot-through (SH) intervals on the system variables during switching periods. The proposed algorithm includes two consecutive stages: OC detection and fault location identification.... 

In this paper, we propose a new delayed data offloading scheme by exploiting full-duplex (FD) device-to-device (D2D) communications. To this end, we consider a scenario in which a common file is requested by a subset of users in different times and with different maximum tolerable delays. Moreover, a maximum specific bandwidth is dedicated for data transmission. In order to send the file in the specific bandwidth, we use multicasting alongside FD D2D communications. Thanks to the fact that users have different deadlines, we divide them into non-overlapping groups, i.e., coalitions, to receive the file in distinct intervals. In order to guarantee the tolerable delay of all users, a coalition...