Sharif Digital Repository

The deformation and breakup dynamics of a compound ferrofluid droplet under shear flow and uniform magnetic field are numerically studied in this paper. Utilizing magnetic field provides the possibility to obtain better control over the compound droplet morphology and breakup in a simple shear flow. To solve the governing equations for interfaces motion and hydrodynamics, the conservative phase field lattice Boltzmann model is employed, and a finite difference approach is applied for calculating the magnetic field. To verify the accuracy of present simulations, the results are validated with those of four relevant benchmarks including liquid lens between two stratified fluids, three-phase... 

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

A high-order compact finite-difference scheme is applied and assessed for the numerical simulation of structural dynamics. The two-dimensional elastic stress-strain equations are considered in the generalized curvilinear coordinates and the spatial derivatives in the resulting equations are discretized by a fourth-order compact finite-difference scheme. For the time integration, an implicit second-order dual time-stepping method is utilized in which a fourth-order Runge-Kutta scheme is used to integrate in the pseudo-time level. The accuracy and robustness of the solution procedure proposed are investigated through simulating different two-dimensional benchmark test cases in structural... 

In the present study, shock-disturbances interaction in high-speed compressible inviscid flow is simulated utilizing the weighted essentially non-oscillatory (WENO) scheme by employing the shock-capturing technique. For this aim, the two-dimensional Euler equations in strong conservative form are discretized by using the explicit third-order TVD Runge–Kutta scheme in time and the fifth-order WENO finite difference scheme in space. The main advantage of using the WENO scheme is its capability for properly solving the discontinuities in the domain without needing any artificial viscosity, limiter function or filter. Hence, this scheme is stable, and thus, it is suitable for simulating very... 

In this study, a coupling of a hydrologic and hydraulic model was utilized to assess the impacts of river geometry data resolution on the flood inundation characteristics in a data-scarce environment. Hydrological modeling incorporates soil conservation service curve-number (SCS-CN) and the geomorphologic based instantaneous unit hydrograph model (GIUH) to compute the direct runoff hydrograph in Kashkan river basin located in western Iran. 1D HEC-Geo-RAS model was used and performed to simulate inundation extent of 100-yr floods (~1800 m3/s) along 40 km reach of Kashkan river with a ground survey of river cross-section (2000 cross-sections, each including 500 data-points). The effect of... 

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

Freeze desalination (FD) works upon the separation of impurities from pure water during ice crystals formation. The required cold source could be supplied by the cold energy of liquefied natural gas (LNG). In the current study, freeze desalination of seawater is explored by directly exploiting the cold energy of LNG within an appropriate range of temperature after producing work in a power generation cycle. A detailed discussion has been given on the inlet temperature of LNG to the FD unit for the first time. The direct utilization has the privilege of eliminating the addition of a secondary refrigerant and its refrigeration cycle to the FD process. A multi-objective optimization is... 

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

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

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

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

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

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

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