S.m. Naghavi, G.a. Sheikhzadeh,
Volume 39, Issue 1 (8-2020)
Abstract
Lattice Boltzmann method is one of computational fluid dynamic subdivisions. Despite complicated mathematics involved in its background, end simple relations dominate on it; so in comparison to the conventional computational fluid dynamic methods, simpler computer programs are needed. Due to its characteristics for parallel programming, this method is considered efficient for the simulation of complex geometry flows, in which a large amount of computational memories is needed. Because of the curved boundaries in the complex geometries, detecting the proper curved boundary condition is unavoidable for the lattice Boltzmann method. For this purpose, more works have been done, and different curved boundary conditions have been proposed. At the present work, first, some curved boundary conditions have been reviewed; then a simplified curved boundary condition is proposed. A computer program based on the lattice Boltzmann method, in FORTRAN language, has been prepared; in this program, the boundary condition along with some others applied on it is proposed. To verify the accuracy and correctness of the proposed boundary condition, 2D cavity flow has been simulated and compared to the available numerical results. Adaptation of the achieved results with those of previous researchers verifies the prepared program correctness. Also, two fluid flows have been simulated, a flow around a stationary cylinder in a 2D channel and one between two stationary and moving cylinders. The results of simulations with the proposed boundary condition, along with the previous boundary conditions, have been compared to the available results. Comparisons demonstrate that solutions with proper accuracy could be obtained by the proposed boundary condition.
N. Fattahi, M. Reisi-Nafchi, G. Moslehi,
Volume 39, Issue 1 (8-2020)
Abstract
Scheduling in production environments is used as a competitive tool to improve efficiency and respond to customer requests. In this paper, a scheduling problem is investigated in a three-stage flexible flowshop environment with the consideration of blocking and batch processing. This problem has been inspired by the charging and packaging line of a large battery manufacturer. In this environment, the first and third stages involve a single processor machine, and the second one consists of m identical parallel batch processing machines. The objective is to minimize the total weighted tardiness of the received orders.Given the lack of consideration of this problem in the literature, first, a mathematical programming model is presented for the problem. Also, due to the NP-hardness of the problem, a variable neighborhood search algorithm and a memetic algorithm are developed to solve it. The computational results show that the variable neighborhood search algorithm can solve instances up to 1200 orders and 15 machines with an average deviation of about 1.9%, relative to the best solution of the two algorithms, and the memetic algorithm can solve instances up to 1200 orders and 15 machines with an average deviation of about 7.8%, as compared e to the best solution of the two algorithms. In general, computational results show the better performance of the variable neighborhood search algorithm in comparison to the memetic algorithm.
F. Shabani, M. Saghafian, D. Saeidi, F. F. Momennasab ,
Volume 39, Issue 2 (2-2021)
Abstract
Particulate separation has many applications in medicine, biology and industry. In this research, the separation of polystyrene particles with a diameter of 10, 20 and 30 μm in the fluid flow of a microchannel is investigated. The microchannel consists of a spiral region and a straight region under the influence of acoustic waves. In the spiral region, the particles under hydrodynamic effects undergo the initial separation; then the particles enter the straight region of the microchannel, and the final separation of the particles is done by the force generated and exerted through the acoustic waves. The effects of acoustic frequency and the number of spiral region loops on separation are investigated. The results show that for measured dimensions and parameters, at 1 MHz acoustic wave, when the number of loops is 2 for the spiral region, the particles at the end of the path are in a suitable position for separation. In addition, the results show that the separation of particles with this hybrid system is better than that done by the simple methods, and the separation rate can be as high as 100%
Z. Z. Ahangari Sisi, M. Mirzaei, S. Rafatnia, B. Alizadeh,
Volume 39, Issue 2 (2-2021)
Abstract
Active vehicle suspension system is designed to increase the ride comfort and road holding of vehicles. Due to limitations in the external force produced by actuator, the design problem encounters the constraint on the control input. In this paper, a novel nonlinear controller with the input constraint is designed for the active suspension system. In the proposed method, at first, a constrained multi-objective optimization problem is defined. In this problem, a performance index is defined as a weighted combination of the predicted responses of the nonlinear suspension system and control input. Then, this problem is solved by the modified firefly optimization algorithm to find the constrained optimal control input. To evaluate the performance of the proposed method, the results of the unconstrained and constrained controllers are provided and discussed for various road excitations. The results show a remarkable increase in the ride comfort with the limited force, while other suspension outputs including the suspension travel and tire deflection being in the acceptable ranges. In addition, these controllers are compared with Sliding Mode Control (SMC) and Nonlinear Model Predictive Control (NMPC) in the presence of model uncertainty.
M. Abouei Ardakan, S. Talkhabi,
Volume 39, Issue 2 (2-2021)
Abstract
One of the common approaches for improving the reliability of a specific system is to use parallel redundant components in subsystems. This approach, which is known as the redundancy allocation problem (RAP), includes the simultaneous selection of the component type and its level for each subsystem in order to maximize the system reliability.Traditionally, there are two redundancy strategies, namely active and standby, for the redundant components. Recently, a new powerful strategy called mixed strategy has been developed. It has been proved that the mixed strategy has a better performance when compared to both previous strategies. The main issue in utilizing the mixed strategy is its complicated formulation and sophisticated calculations, leading to a time-consuming procedure for solving the problems. Hence, in this paper, a new formulation based on the recursive approach is introduced to ease the calculation of the mixed strategy. In the new formulation, the complex double integral calculations are removed and the calculation times is reduced. The proposed recursive formulation provides a general statement for the mixed strategy formula which is not changed by altering the number of components in each subsystem. This flexibility and stability in the formula can be very important, especially for large scale cases. In order to evaluate the new approach and to compare its performances with the previous formulation, a benchmark problem with 14 subsystems is considered and the results of the two formulation are compared with each other.
Z. Barouei, M. Jabbarzadeh,
Volume 39, Issue 2 (2-2021)
Abstract
In this paper, the nonlinear bending analysis for annular circular nano plates is conducted based on the modified coupled stress and three-dimensional elasticity theories. For this purpose, the equilibrium equations, considering nonlinear strain terms, are calculated using the least energy potential method and solved by the numerical semi-analytical polynomial method. According to the previous works, there have been no studies calculating all boundary conditions numerically based on three-dimensional elasticity. Typically, the research done on three-dimensional elasticity is either finite element or only for a simply-supported boundary condition. In this research, for the first time, the nonlinear analysis of bending is calculated with the help of three-dimensional elasticity for a variety of boundary conditions. Also, with the help of the modified couple stress theory, the results on the nano-scale scale have been studied. In the following, while validating the results, we investigate the changes in the scale parameter for the types of boundary conditions, the effect of changing the parameter of scale in different thicknesses, and the impact of the parameter of scale on the linear and nonlinear results.
R. Keshavarzi, Sh. Hatami, Sh. Hashemi,
Volume 39, Issue 2 (2-2021)
Abstract
Plates made of laminated composite materials with variable stiffness can have wide applications in various branches of engineering due to such advantages as high strength /stiffness to weight ratio. In these composites, curved fibers are used to reinforce each lamina instead of the straight fibers. In this paper, the application of finite strip method for the buckling analysis of moderately thick composite plates with variable stiffness is investigated. For buckling analysis, a semi-analytical finite strip method based on the first-order shear deformation theory is employed. In this method, all displacements are presumed by the appropriate harmonic shape functions in the longitudinal direction and polynomial interpolation functions in the transverse direction. The minimum potential energy method has been used to develop the stability formulations. This analysis examines the effect of using curved fibers instead of straight fibers on the laminate composites. The critical loads obtained from this analysis are compared with those of other researchers and the efficiency and accuracy of the developed finite strip method are confirmed. Comparison of the analysis results of these plates shows that changing the slope of the fibers can lead to a significant change in the buckling response. Also, increasing the number of the terms of shape functions in the longitudinal direction has a significant effect on the convergence to the desired results.
S. Esmizadeh, H. Haftbaradaran, F. Mossaiby,
Volume 39, Issue 2 (2-2021)
Abstract
Experiments have frequently shown that phase separation in lithium-ion battery electrodes could lead to the formation of mechanical defects, hence causing capacity fading. The purpose of the present work has been to examine stress intensity factors for pre-existing surface cracks in spherical electrode particles during electrochemical deintercalation cycling using both analytical and numerical methods. To this end, we make use of a phase field model to examine the time-dependent evolution of the concentration and stress profiles in a phase separating spherical electrode particles. By using a geometrical approximation scheme proposed in the literature, stress intensity factors at the deepest point of the pre-existing surface cracks of semi-elliptical geometry are calculated with the aid of the well-established weight function method of fracture mechanics. By taking advantage of a sharp-interphase core-shell model, an analytical solution for the maximum stress intensity factors arising at the deepest point of the surface cracks during a complete deintercalation half-cycle is also developed. Numerical results for evolution of the concentration profile and the distribution of the hoop stresses in the particle are presented; further, the stress intensity factors found numerically based on the phase field model are compared with those predicted by the analytical core-shell model. The results of the numerical model suggest that the maximum stress intensity factor could significantly vary with changes in the surface flux, increasing potentially by a factor of two within the range of parameters considered here, when the concentration difference between the two phases is decreased.
O. Bateniparvar, N. Noormohammadi, A. M. Salehi,
Volume 39, Issue 2 (2-2021)
Abstract
In this paper, Equilibrated Singular Basis Functions (EqSBFs) are implemented in the framework of the Finite Element Method (FEM), which can approximately satisfy the harmonic PDE in homogeneous and heterogeneous media. EqSBFs are able to automatically reproduce the terms consistent with the singularity order in the vicinity of the singular point. The newly made bases are used as the complimentary enriching part along with the polynomial bases of the FEM to construct a new set of shape functions in the elements adjacent to the singular point. It will be shown that the use of the combined bases leads to the quality improvement of the solution function as well as its derivatives, especially in the vicinity of the singularity.
M. Hashemian, M. Jabbarzadeh,
Volume 40, Issue 1 (9-2021)
Abstract
In this paper, nonlinear bending analysis of functionally graded rectangular and sectorial micro/nano plates is investigated using the modified couple stress theory. For this purpose, a higher-order shear deformation theory and von Kármán geometrically nonlinear theory are employed. The equilibrium equations and the boundary conditions for rectangular and annular sector plates are derived from the principle of minimum total potential energy and solved using the Semi-Analytical Polynomial Method (SAPM). One of the advantages of the implemented shear deformation theory is removing the defects of higher order shear deformation theory, and obtaining the response of the first and the third-order shear deformation theories at the same time. Afterwards, beside investigating the benefits of this theory compared with other ones, the results are verified with those by other researches. At the end, the effects of length scale parameter, boundary conditions, power law index, and geometrical dimensions are investigated
S. Moomivand, H. Davari-Ardakani, H. Mosadegh, M. Abouei Ardakan,
Volume 40, Issue 1 (9-2021)
Abstract
In this paper, a multi-mode resource constrained project selection and scheduling problem is investigated considering the reinvestment strategy in a flexible time horizon. Among a set of available projects, a number of projects are selected and scheduled regarding the constraints on renewable resources and precedence relations. The benefits of project portfolio selection and scheduling are compared in both fixed and flexible time horizons. For this purpose, upper and lower tolerance limits are considered for the predetermined time horizon. If the schedule exceeds the time horizon, a penalty cost will be charged. The objective is to determine the optimal time horizon. A mixed-integer linear programming model is proposed for this problem, and solved by GAMS software/CPLEX solver and also a combination of a proposed heuristic algorithm, Genetic Algorithm, and a local search method. Numerical results show that the proposed approach has an acceptable performance in terms of the quality of the solution and the running time. Also, dealing with the problem in a flexible time horizon is more profitable compared to a fixed time horizon.
P. Rastegar Rajeouni, A. R. Rahmati,
Volume 40, Issue 1 (9-2021)
Abstract
In order to simulate multiphase flow in the presence of dielectric current using the Lattice Boltzmann Method (LBM), three distribution functions are used, two of which for using the He-Chen-Zhang phase field model and one for solving the potential field. Initially, the ability of the code to apply surface tension was tested using the Laplace law and the drop release test. The results show that the present numerical program is capable of modeling well the regulated surface tension force. Then, the Rayleigh–Taylor instability simulation is used to evaluate the code's ability in applying volume forces. The results by the developed numerical program are in good agreement with the numerical results in the references. In this study, for the first time, the effect of electric field on a droplet immersed in another fluid and the presence of droplet in a porous medium is investigated by LBM. For this purpose, first the droplet motion due to the potential difference in the porous and non-porous media is investigated. After modeling the droplet motion due to the potential difference, two electric fields areapplied to the droplet to reverse the droplet deformation. Through various tests, it is shown that at a given potential difference, the droplet breaks down after much deformation and is divided into smaller droplets. The decomposition of droplets in a pre-mixed emulsion is a common technique in the production of monodisperse droplets. The presence of monodisperse droplets in an emulsion improves the physical properties of polymer science experts.
B. Movahedian Attar, M. Sadeghi,
Volume 40, Issue 1 (9-2021)
Abstract
Accurate determination of the response of structures under dynamic loads such as earthquake loads plays an important role in the safe and economical design of structures. The purpose of this paper is to utilize a novel solution method based on the use of exponential basis functions for dynamic analysis of Bernoulli beam subjected to different types of base excitations. This method was firstly introduced for solving scalar wave propagation problems, named as stepwise time-weighted residual method. The proposed method considers the solution as a series of exponential basis functions with unknown constant coefficients; and the problem is solved in time without the need for spatial discretization of the beam and by using an appropriate recursive relation to correct the coefficients of the exponential bases. In order to apply the earthquake excitation, first by using the central finite difference relation, the earthquake acceleration history is converted to displacement history. Moreover, the displacement history is applied to the beam as a time-varying boundary condition. In this study, the capabilities of the proposed method in solving several sample problems of vibration of single and multi-span beams under various stimuli such as earthquake acceleration variations are compared with the results of other existing methods.
M. Azizpooryan, N. Noormohammadi,
Volume 40, Issue 1 (9-2021)
Abstract
In this paper, static analysis of in-plane heterogeneous laminated composite plates is numerically studied. The Mindlin’s theory which considers linear transverse shear deformation has been implemented. The governing partial differential equation is satisfied by a weighted residual integration. Chebyshev polynomials of the first kind are used as basis functions and exponential functions make up the weight functions of the integration. The emerging integrals may be composed of some pre-evaluated 1D normalized ones, which effectively paces up the solution progress. To verify the method, several examples of homogeneous as well as heterogeneous plates with various lamination schemes and boundary conditions have been solved. Results are compared with those from the literature or by commercial codes, which reveal excellent accuracy of the proposed method.
Y. Mollapour, E. Poursaeidi, H. Shayani-Jam, O. Pedram,
Volume 40, Issue 1 (9-2021)
Abstract
Corrosive factors along with mechanical loads on the gas turbine compressor blades, cause phenomena such as pitting corrosion, stress corrosion cracking and corrosion fatigue. Due to erosion of particles in the presence of a corrosive environment, pitting happens on the blade surfaces, which is a source of subsequent cracks. Therefore, it is necessary to get knowledge of its mechanism in order to prevent the phenomena as much as possible. The main purpose of this paper is to investigate the growth of pitting corrosion in CUSTOM 450 stainless steel and to obtain strain values in the growing pits at the maximum bending region. In this regard, a two-point bending specimen was made and subjected to a potentio-static test under the potential of 350 mVSCE in the 3.5 wt% sodium chloride solution. Then the propagated pits were numerically examined. By the digital image correlation method, the local strain was calculated in the pits and a relation was presented to obtain the maximum strain time. Therefore, growth direction of pitting corrosion could be estimated by having maximum strain region. Finally, by simulating the pitting corrosion process of a stress-free sample under the potential of 350 mVSCE in 3.5 wt% sodium chloride solution in COMSOL Multiphysics software, variations in the concentration of ions, electric potential, and corrosion current density were shown in the existing pit. The potential was decreased by moving in-depth and the maximum current density was found at the depth of 18 μm. Thus, without the need of advanced laboratory facilities for surface scanning and analysis, useful information from surface corrosion conditions could be obtained
H. Asadigorji, A. Karami Mohammadi,
Volume 40, Issue 1 (9-2021)
Abstract
Complex nonlinear behaviors such as chaotic motion have devastating effects on dynamic systems. In this study, nonlinear behavior of simply supported rectangular viscoelastic plates was examined during supersonic aerodynamics and compared with the nonlinear elastic plate. Classical plate theory was used to obtain the plate equations, and Von- Kármán strain-displacement relations were used to consider the nonlinear geometric effects. The Kelvin Voigt model was also used to describe the viscoelastic properties and the “first-order piston theory" was used for supersonic aerodynamic flow. The equations of motion of the rectangular plate were extracted using the Lagrangian method and then, discretized by the Rayleigh-Ritz method. Solution of the equations was performed using fourth order Runge Kutta method. To investigate the dynamic behavior of the plates, the eigenvalues of the system, time history curves, phase portraits, Poincaré maps, and bifurcation diagrams were studied and analyzed. The results show that in some aspect ratios, the threshold for the occurrence of the flutter in the viscoelastic plate will be lower than that in the elastic plate. On the other hand, when the control parameter increases, complex nonlinear behavior such as chaos in the elastic plate goes simpler in the viscoelastic plate, such as periodic motion.
A. R. Rahmati , E. Kashai,
Volume 40, Issue 2 (1-2022)
Abstract
A two-phase lattice Boltzmann model considering the interaction forces of nanofluid has been developed in this paper. It is applied to investigate the flow and natural convection heat transfer of Al2O3–H2O nanofluid in an enclosure containing internal heat generation. To understand the heat transfer enhancement mechanism of the nanofluid flow from the particle level, the lattice Boltzmann method is used because of its mesoscopic feature and numerical advantages. By using a two-component lattice Boltzmann model, the heat transfer enhancement of the nanofluid is analyzed through incorporating the different forces acting on the nanoparticles and the base fluid . The effects of interaction forces, nanoparticle volume fractions (0.0-0.05), and internal and external Rayleigh numbers (103-106) on the nanoparticle distributions and heat transfer characteristics are investigated. The average Nusselt number increases with the increase of nanoparticle volume fraction and Rayleigh number. We also compared and analyzed adding internal heat generation on the nanoparticles and the base fluid separately, and it was found that by considering heat generation on the base fluid, it mostly affects the temperature field, and by considering that on nanoparticles, it mostly affects the stream field.
F. Hosseinlou,
Volume 40, Issue 2 (1-2022)
Abstract
Today many complex models, typically finite element models, have been employed in the analysis of jacket offshore structures. However, these comprehensive models are not readily adopted in engineering practice, especially during the preliminary design stage. As the dynamic analysis of jacket platforms is very complicated, it will be very advantageous to make a simplified computational method to assess dynamic performance of such structures. In this work a refined simplified model has been utilized to calculate dynamic responses of jacket platforms. In this regard, the model simplification based on the vibration modal data and Timoshenko’s beam equation has been employed to overcome the uncertainty problem in modeling. According to the curve fitting method on the set of frequency response functions to derive modal parameters, the concept of power spectrum density has been also used to confirm the proposed computational model.In this regard, first the behavior of the physical model in the frequency domainhas been presented and compared with the spectral results obtained from the simplified model based on Timoshenko beam. Because the modal test of the physical model was performed under the force of white noise, the dynamic responses of the simplified model were also extracted under the force of white noise using MATLAB software. In this paper, an applied mathematical model has been produced, and it has been demonstrated that the refined simplified model can reflect the real structural features.
F. F. Heidargheitaghi, M. H. M. H. Ghadiri Rad, M. Kazemi,
Volume 40, Issue 2 (1-2022)
Abstract
Continuously varying cross-section members have found wide applications in engineering for cost and resistance optimization. Since steel structures generally have more slender members compared to concrete structures, buckling analysis of steel members is of more importance. Determining the critical load of functionally varying cross-section columns using the analytical solution is a time-consuming process. In this paper, buckling analysis of non-prismatic steel columns is conducted using the meshless local Petrov-Galerkin (MLPG) method. In meshless methods, the scattered nodes are used rather than the elements to model the problem domain and its boundaries. The change of the inertia moment within the length of a column is characterized by introducing a power function with variable taper ratio and exponent. The radial basis function is used to discretize the differential equation governing the buckling. The penalty method is used for the imposition of the boundary conditions. Numerical examples of the critical buckling load for prismatic and non-prismatic columns using the proposed method are compared with the analytical solution, and the effectiveness of the MLPG method for buckling analysis of non-prismatic columns is validated. Also, buckling analysis of muscle column members subjected to non-uniform axial load is carried out to show the efficiency of the proposed method. The effect of several parameters such as non-uniformity of the load and variation of the cross-section on the buckling load of the column is discussed in details.
M. Nemati, M. Sefid, M. S. Barghi Jahromi, R. Jahangiri,
Volume 40, Issue 2 (1-2022)
Abstract
In the present work, the effect of magnetic field, changes in the angle of inclination of the cavity and the shape of nanoparticles on the flow field and heat transfer of water-alumina with uniform heat generation/absorption is investigated by Lattice Boltzmann method (LBM). The curved wall and the diagonal walls of the cavity are at a constant temperature of hot and cold, respectively. Nanoparticle volume fraction of 0, 0.02 and 0.04, Hartmann number of 0, 15, 30, 45 and 60, heat generation/absorption coefficient of -5, 0 and +5 and inclination angle of 45, 135 and 225 degrees are studied. The high accuracy of the results compared to previous studies confirmed the correctness of the code written in Fortran language. The results shows that in all cases, increasing the Hartmann number leads to a decrease in the maximum value of the streamlines and the average Nusselt number, with the lowest effect at 225 degrees. Also increasing the strength of the magnetic field leads to an average decrease of 28, 23 and 7% of the average Nusselt number for angles of 45, 135 and 225 degrees, respectively. Increasing the heat generation/absorption coefficient is a determining factor in the effectiveness of the magnetic field and adding nanoparticles, and increasing it reduces the amount of heat transfer. On average, heat generation reduces the average Nusselt number by 71, 98, and 145 percent for the angles of 45, 135, and 225 degrees, respectively. In general, the lowest value of the average Nusselt number is related to the angle of 225 degrees, but the effect of adding nanoparticles in increasing the average Nusselt number is the highest at this angle. Generally, an increase in the percentage of nanoparticles leads to an average increase of 12% in the average Nusselt number. The effect of nanoparticle shape is more apparent with increasing their volume fraction. The highest amount of heat transfer is related to the cylindrical nanoparticles, in which the average Nusselt number is on average about 6% higher than the spherical state.
