Sharif Digital Repository

Identification of existing instability modes from experimental pressure measurements of rocket engines is difficult, specially when steep waves are present. Actual pressure waves are often non-linear and include steep shocks followed by gradual expansions. It is generally believed that interaction of these non-linear waves is difficult to analyze. A method of mode identification is introduced. After presumption of constituent modes, they are superposed by using a standard finite difference scheme for solution of the classical wave equation. Waves are numerically produced at each end of the combustion tube with different wavelengths, amplitudes, and phases with respect to each other. Pressure... 

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

SnO2-CuO is one the most promising systems for detection of detrimental H2S gas. Although previous experimental research has suggested a sulfidation reaction to explain selectivity toward H2S, little is known about the origin of change of electrical response of this system by changing the H2S gas concentration. In this study the relation between sensing response of continuous SnO2-CuO bilayer thin film and H2S gas concentration is computed based on changeability of chemical composition of covellite CuxS. For this purpose, chemical activity of sulfur as a function of atomic fraction in covellite copper sulfide is estimated using Gibbs energies of formation and chemical thermodynamics. By... 

The formation damage (FD) caused by the invasion of drilling fluid severely affects reservoir performance during production. Most of the published research studies which address this type of FD have been carried out at the core or field scale. Thus, the main aim of the paper is to investigate the pore-scale mechanisms of FD induced by drilling fluids and their control with silica nanoparticles (NPs) using a microfluidic approach. The proper identification of the mechanisms of FD can lead to the proper selection of NP type and concentration as well as a suitable method to remediate FD. The micromodel was designed in a way to closely simulate the cross-flow at the wellbore surface. A... 

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

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is extended and applied to accurately simulate two-phase liquid-vapor flows with high density ratios. Herein, the He-Shan-Doolen-type lattice Boltzmann multiphase model is used and the spatial derivatives in the resulting equations are discretized by using the fourth-order compact finite-difference scheme and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient two-phase flow solver. A high-order spectral-type low-pass compact nonlinear filter is used to regularize the numerical solution and remove spurious waves generated by flow nonlinearities in smooth regions... 

In this paper, two very simple and accurate equations are introduced that enable geotechnical engineers to determine the undrained ultimate bearing capacity of shallow strip foundations resting on two-layered clays. Results of the analyses show that in the case of weak-overstrong clay, usually general shear failure occurs, and in the case of strong-over-weak clay, the soil will most likely experience punching failure with a large amount of plastic settlement before it reaches its ultimate bearing capacity. This study applied the finite-difference method using a computer program to evaluate the undrained bearing capacity for all cases analyzed. The Mohr-Coulomb failure criterion was used for... 

In this work, the capability and performance of the vorticity confinement (VC) implemented in a high-order accurate flow solver in predicting two-dimensional (2D) compressible mixing layer flows on coarse grids are investigated. Here, the system of governing equations with incorporation of the VC in the formulation is numerically solved by the fourth-order compact finite difference scheme. To stabilize the numerical solution, a low-pass high-order filter is applied, and the nonreflective boundary conditions are used at the farfield and outflow boundaries to minimize the reflections. At first, the numerical results without applying the VC are validated by available direct numerical... 

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

Recent studies in information computation technology (ICT) are focusing on Next-generation networks, SDN (Software-defined networking), 5G, and 6G. Optimal working mode for device-to-device (D2D) communication is aimed at improving the quality of service with the frequency spectrum structure is of research areas in 5G. D2D communication working modes are selected to meet both the predefined system conditions and provide maximum throughput for the network. Due to the complexity of the direct solutions, we formulated the problem as an optimization problem and found the optimal working modes under different parameters of the system through extensive simulations. After determining the links'... 

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

This paper extends the analysis of solute dispersion in electrohydrodynamic flows to the case of band broadening in polyelectrolyte-grafted (soft) capillaries by accounting for the effects of ion partitioning, irreversible catalytic reaction and pulsatile flow actuation. In the Debye-Hückel limit, we present the benchmark solutions of electric potential and velocity distribution pertinent to steady and oscillatory mixed electroosmotic-pressure-driven flows in soft capillaries. Afterwards, the mathematical models of band broadening based on the Taylor-Aris theory and generalized dispersion method are presented to investigate the late-time asymptotic state and all-time evolution of... 

This paper extends the analysis of solute dispersion in electrohydrodynamic flows to the case of band broadening in polyelectrolyte-grafted (soft) capillaries by accounting for the effects of ion partitioning, irreversible catalytic reaction and pulsatile flow actuation. In the Debye-Hückel limit, we present the benchmark solutions of electric potential and velocity distribution pertinent to steady and oscillatory mixed electroosmotic-pressure-driven flows in soft capillaries. Afterwards, the mathematical models of band broadening based on the Taylor-Aris theory and generalized dispersion method are presented to investigate the late-time asymptotic state and all-time evolution of... 

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

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 study, a coupled numerical approach based on Lattice Boltzmann Method (LBM) and Discrete Element Method (DEM) is employed for two-dimensional simulation of fluid flow in deformable particulate media comprising of movable circular particles. The developed LB-DE code is validated against the results of a bi-axial shear test as well as two well-known benchmark problems including settling of a circular particle under gravity force inside a viscous fluid, and motion of a neutrally buoyant particle released in a Poiseuille flow. The verified code is then utilized for simulation of "Sand Production" phenomenon which is of importance for oil producing wells in weakly cemented sandstone... 

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

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