JOURNAL OF COMPUTATIONAL METHODS IN ENGINEERING (ESTEGHLAL)
Year:2020 | Volume:38 | Issue:2 |Papers Count:8

Wave-field extrapolation based on solving the wave equation is an important step in seismic modeling and needs a high level of accuracy. It has been implemented through a various numerical methods such as finite difference method as the most popular and conventional one. Moreover, the main drawbacks of the finite difference method are the low level of accuracy and the numerical dispersion for large time intervals (Δ t). On the other hand, the symplectic integrators due to their structure can cope with this problem and act more accurately in comparison to the finite difference method. They reduce the computation cost and do not face numerical dispersion when time interval is increased. Therefore, the aim of the current paper is to present a symplectic integrator for wave-field extrapolation using the Euler method. Then, the extrapolation is implemented for rather large time intervals using a simple geological model. The extrapolation employed for both symplectic Euler and finite difference methods showed a better quality image for the proposed method. Finally the accuracy was compared to the finite difference method.

In this study, we present a numerical solution for the two-phase incompressible flow in the porous media under isothermal condition using a hybrid of the linear lower-order nonconforming finite element and the interior penalty discontinuous Galerkin (DG) method. This hybridization is developed for the first time in the two-phase modeling and considered as the main novelty of this research. The pressure equation and convection dominant saturation equation are discretized using the nonconforming Crouziex-Raviart finite element (CR FEM) and the weighed interior penalty discontinuous Galerkin (SWIP) method, respectively. Utilizing the nonconforming finite element method for solving the flow equation made the pressure and velocity values be consistent with respect to the degrees of freedom arrangement at the midpoint of the neighboring element edges. The boundary condition governing the simulation is the Robin type at entrance boundaries, and the time marching discretization for the governing equations is the sequential solution scheme. An H (div) projection using Raviart-Thomas element is implemented to improve the results’ resolution and preserve the continuity of the normal component of the velocity field. At the end of each time step, the non-physical oscillation is omitted using a slope limiter, namely, modified Chavent-Jaffre limiter, in each element. Also, in this study, the developed algorithm is verified using some benchmark problems and the test cases are considered to demonstrate the efficiency of the developed model in capturing the shock front at the interface of fluid phases and discontinuities.

This paper presents a numerical analysis of granular column collapse phenomenon using a two-dimensional smoothed particle hydrodynamics model and a local constitutive law proposed by Jop et al. This constitutive law, which is based on the viscoplastic behaviour of dense granular material flows, is characterized by an apparent viscosity depending both on the local strain rate and the local pressure. The rheological parameters are directly derived from the experiments. A simple proposed regularization method used in the viscosity relation to reproduce the stopping condition and the free surface of a granular flow where the pressure is disappeared. Pressure oscillation, as the main disadvantage of the weakly compressible SPH method, leads to an inaccurate pressure distribution. In this research, a new algorithm is proposed to remove the nonphysical oscillations by relating the divergence of velocity to the Laplacian of pressure. The simulations based on the proposed SPH algorithm satisfactorily capture the dynamics of gravity-driven granular flows observed in the experiments. The maximum thickness of a granular flowing on a rough inclined plane is obtained based on the local rheology model and compared with the experimental results. The run-out distances and the slopes of the deposits in the simulations showed a good agreement with the values found in the experiments. The results of the simulation proved that the initial column ratio played an important role in spreading the granular mass.

Graphene is one of the nanostructured materials that has recently attracted the attention of many researchers. This is due to the increasing expansion of nanotechnology and the application of this nanostructure in technology and industry owing to its mechanical, electrical and thermal properties. Thermal buckling behavior of single-layered graphene sheets is studied in this paper. Given the failure of classical theories to consider the scale effects and the limitations of the nano-sized experimental investigations of nano-materials, the small-scale effect is taken into account in this study, by employing the modified couple stress theory which has only one scale parameter. On the other hand, the two-variable refined plate theory, which considers the shear deformations in addition to bending deformations, is used to define the displacement field and to formulate the problem. The developed finite strip method formulation is used to evaluate the critical buckling temperature of the nanoplates. The validity of the proposed method is confirmed by comparing the results of this study with the those in the literature. The effects of different boundary conditions, temperature changing patterns, aspect ratio, and the ratio of length parameter to thickness on the critical buckling temperature are considered and the results are presented in the form of Tables and Figures.

With respect to the great application of pipes conveying fluid in civil engineering, presenting a mathematical model for their stability analysis is essential. For this purpose, a concrete pipe, reinforced by iron oxide (Fe2O3) nanoparticles, conveying fluid is considered. The goal of this study is to investigate the structural stability to show the effects of the inside fluid and the nanoparticles. The structure was modeled by a cylindrical shell and using Reddy theory. To obtain the force induced by the inside fluid, the Navier-Stokes equation was used. To assume the effect of the nanoparticles in the pipe, the Mori-Tanaka model was utilized so that the effects of agglomeration of nanoparticles could be considered. Finally, by applying energy method and the Hamilton's principle, the governing equations were derived. For the stability analysis of the structure, differential quadrature method (DQM) was proposed and the effects of different parameters such as volume fraction of the nanoparticles and agglomeration of the nanoparticles inside fluid and geometrical parameters were investigated. The results showed that the existence of the nanoparticles as the reinforcement for the pipe led to the delay in the pipe instability.

As thermalisation loss is the dominant loss process in the quantum dot intermediate band solar cells (QD-IBSCs), it has been investigated and calculated for a QD-IBSC, where IB is created by embedding a stack of InAs(1-x) Nx QDs with a square pyramid shape in the intrinsic layer of the AlPySb(1-y) p-i-n structure. IB, which is an optically coupled but electrically isolated miniband, divides the total band gap of AlPySb(1-y) into two sub-band gaps. To obtain the thermalisation loss of AlPySb(1-y)/InAs(1-x)Nx QD-IBSCs, the position and width of IB in the band gap of AlPySb(1-y) should be calculated. The position of IB, which is equal to the first eigen-energy of a unit cell of QD, is obtained by solving the 3D Schrö dinger equation with a finite-element method and the width of IB is obtained by the absorption characteristics. Then, with the investigation of the effect of nitrogen and phosphorous molar fraction, QDs size and the distance between the QDs on the thermalisation loss, the minimized loss for the optimized structure of AlPySb(1-y)/InAs(1-x)Nx QD-IBSCs is obtained.

Nowadays, outpatient providers are struggling to reduce the current costs and improve the service quality. A part of the outpatient service provider is a hemodialysis department with expensive supplies and equipment. Therefore, in the present paper, the scheduling of hemodialysis patients with their preferences has been studied. The aim of scheduling hemodialysis patients in this study is to minimize the normalized weighted sum of deviations from the patients' preferences and the total completion time. It should be noted that the patient's preferences include beds, treatment combination of days and their turn. To solve the problem, two mathematical models have been presented. Performence of the models in solving the real data of the hemodyalisis department of Imam Khomeini Hospital, in Kermanshah, was investigated. The results showed the efficiency of the proposed models in considering the preferences of patients; however, these preferences in the hospital schedule were considered in few cases, as far as it was possible. So, these preferences has no priority in the hospital schedule. In addition to considering the patients’ preferences, the solution of models reduced the total completion time of the pationts treatment. Also, one of the proposed models in this papercould optimally solve the instances three times larger than the hospital cases.