نشریه مهندسی مکانیک مدرس
سال نوزدهم شماره 12 (آذر 1398)

Biodegradable polymers have widespread usages in the biomedical field, such as stents, sutures, scaffolds, and implants. Due to the importance of behavior of these materials exposed to environmental effects, whether in nature or the human body, extensive researches have been carried out in the last decade that most of them are experimental results and very few are theoretical results. These researches have mainly been performed for specific loading and temperature conditions and so on. For this purpose, in addition to validating the theoretical and empirical relationships derived through the experimental results, the effects of more complex conditions can be considered using the finite element method and numerical solution. In this paper, an analytical relationship extraction method has been presented, as well as the abilities and weaknesses of biodegradable polymers have been investigated by presenting the experimental results of biodegradable polymers. A numerical and finite element analysis is also provided to analyze the behavior of biodegradable polymers. The theoretical analysis and numerical simulation of biodegradable polymers have been carried out using the neo-Hookean hyperelastic model. First, the relationship of stress, versus the stretch has been derived using the strain energy of neo-Hookean material. Next, by assuming a degradation parameter, changes in the properties of the material exposed to environmental effects, according to the time in Abaqus Umat subroutine have been applied to the model. Finally, the accuracy of the simulation has been studied by a comparison between the experimental results and theoretical analyses with numerical solutions.

In this research, the behaviour of a single droplet of the dielectric field under the effect of the applied external uniform magnetic field has been investigated. Previously, it was thought that no force is applied to the dielectric fluids when exposed to the uniform magnetic field. A stagnant droplet in a quiescent fluid column was considered in order to determination of the net effect of the applied uniform magnetic field. Considering that the droplet behaviour has been investigated in the horizontal plane, the net effect of the gravity on the droplet and the surrounding fluid is also zero. Therefore, any change in the condition of the considered droplet will be due to the effect of the applied magnetic field. Numerical analysis has been used to perform this research. The governing equations of the problem are the continuity, momentum, level set equations for interface simulation, re-initialization and re-construction equations of the level set equations to control the mass dissipation of this method. The governing equations have been discretized and solved by developing code in the Fortran programming environment. The behaviour of the considered droplet in various regimes has been investigated under the different magnitudes of the applied magnetic field. The results of the research in various cases show that stagnant droplet deforms under the effect of the applied magnetic field and starts to vibrate which also induces the motion in the surrounding quiescent fluid.

In this paper, the effect of the Gurney flap parameters such as the angle of attack, Reynolds number, angle and height of the flap and its location from the escape edge on the lift coefficient of a symmetric airfoil is considered with the help of simulation in computational fluid dynamic software of Fluent. The turbulence model k-ε is used for the two-dimensional domain. Also, the value of the lift coefficient is introduced as a function of effective parameters by the design of experiment (DOE) method and using the backward elimination regression model which is a statistical method for selecting the model and estimated error terms. The value of the airfoil lift coefficient can be determined and predicted by the obtained function. The numerical values derived from the function of the lift coefficient resulted from the design of experiment method are in good agreement with other valid papers. The results show that at the constant attack angle by increasing Gurney flap height, lift coefficient increase. On the other hand, at the constant height of the Gurney flap, this coefficient decreases with increasing angle of attack. Moreover, the lift coefficient increased by increasing the distance of the Gurney flap from the airfoil escape edge at a 90-degree angle and 1.5%, as well as increasing the Reynolds number at a constant height of a Gurney flap.

In this paper, the effects of unified temperature loading and Winkler-Pasternak elastic foundation on the vibration of functionally graded curved nanobeam have been studied. The proposed model is based on the modified couple stress theory and the Timoshenko beam model. The continuous distribution of material along the thickness of functionally graded curved nanobeam is achieved by changing the gradient index in the volume fraction. The governing equations and related boundary conditions are obtained using the Hamilton principle. By analyzing the quantitative and qualitative results in the tables and figures, influences of geometrical and thermo-physical parameters such as gradient index, aspect ratio, unified temperature difference, the ratio of thickness to length scale parameter and arc angle of functionally graded curved nanobeam on the natural frequency for different vibration mode have been interpreted. There is an excellent agreement between the present results and the results of the previous works. Applied temperature loading increases the sensitivity of the natural frequency to the changes in the aforementioned parameters and also increases the range of its changes. Also, applying the Pasternak elastic foundation changes the behavior of the natural frequency to the temperature changes.

Today, mathematical models and numerical methods are highly regarded according to their ability to predict and understand the cancer treatment process. In this research, the drug delivery to the solid tumor with considering its normal surrounding tissue has been studied by stimulating the blood flow in the dynamic capillary network and interstitial flow and adding the solute transport equation to the fluid flow equations. In the present study for the first time, drug delivery has been studied by a multi-scale comprehensive model with considering two parent vessels with different branches and different input and output pressures. In this paper, the intravascular flow was simultaneously simulated with the interstitial fluid flow. The distribution of drug concentration has been investigated at different times. The results show the dependence of the drug delivery to the interstitial fluid pressure, the pressure of the parent vessel and in fact, the blood pressure of each patient, and the capillary network structure. In addition, an increase of about 20% of the average drug concentration in the tumor site in the present study compared to the previous study with a parent vessel is evidence of the key role of the capillary network and its dependent parameters.

In this paper, a trajectory tracking control of a nonholonomic wheeled mobile robot is proposed based on terminal sliding mode control, and the proposed method has been implemented on a wheeled mobile robot. A wheeled mobile robot is a nonlinear nonholonomic system, and it has three extended coordinates and a nonholonomic constraint. First, the equation of wheeled mobile robot for the extended chained form is derived by transformation of the nonholonomic system equation to the extended chained form. Then a finite time terminal sliding mode approach for trajectory tracking control of the wheeled mobile robot is presented. Afterward, with a graphical simulation environment which is applicable in the Matlab software, graphical simulations of wheeled mobile robot’s movement are done. The result of the graphical simulation in comparing with sliding mode control show the performance of the proposed method. Finally, the practical results of implementation of the controller for trajectory tracking of the wheeled mobile robot is shown, and the results show good tracking performance of the proposed method.

Maneuvering with the highest speed and low power has always been a challenge to design a satellite and spacecraft control system. In this paper, apart from the complexity of modeling actuators, different control methods were used to control the satellite attitude in the presence of uncertainties and disturbances in satellites, in order to obtain an explicit response to minimize the EULERINT criterion. The EULERINT criterion is the integral of the Euler angles between the body axes and the target around Euler's axis over time and somehow interprets the speed of the satellite maneuver in the three control axes. First, using the proportional-derivative control, the comparison of the EULERINT criterion in the application of different kinematic representations (Euler, quaternion vectors and direction cosine matrix equations) in linear and nonlinear models of the satellite was carried out. Then the comparison of the EULERINT criterion between the different methods was presented using the quaternion kinematic, which has the least amount of EULERINT, through changing the proportional-derivative controller to linear-quadratic regulator controllers, pole placement, adaptive, fuzzy, and adaptive-fuzzy. The comparison was conducted to achieve the best control method in terms of frequency response, the lowest EULERINT and the least control effort to control the attitude of the satellite in the presence of disturbance and uncertainty.

In the present work, the effects of blades number on the performance of two stages axial gas turbine have been investigated numerically. Geometry characteristics of the gas turbine have been chosen based on the F5 model of General Electric Company. First, the blades geometry and fluid passages are initially generated due to the real dimensions of the turbine and the generated geometry is networked. Then, the final model of the turbine is generated by gridding blades which set beside each other. Then, Ansys CFX software is used to solve the 3D Navier-Stokes equations in the generated computational domain. The shear stress transport turbulence model has been employed in order to determine the wall effects on the turbulent flow. Before any change in the main turbine, a numerical study was performed and a comparison was conducted between the numerical results and experimental results measured in the power station which the results show a good level of agreement between them. The number of blades of each row has been changed in order to investigate the effects of blade number on the turbine efficiency. The results show that the power generation of the turbine and its efficiency are increased by 0.83% and 0.81%, respectively by an increase in the number of second-row stator blades from 62 to 71 blades.

Delamination or interlayer cracking is one of the most important imperfections in composite materials. The existence of this defect in a structure reduces the strength and, as a result, disables the structure. To analyze the effective factors in interlayer separation, it is necessary to analyze the effective loading parameters. In this paper, the effect of the change in loading rate on the failure mechanism in I failure mode was analyzed using an acoustic emission for unidirectional samples made of glass fiber/epoxy resin. At first composite, samples were made according to standard and placed at different rates of displacement under loading. Force data, displacement and crack growth rate for different loading rates were used to calculate the exact strain energy release rate. In addition to the extensometer, the Dino camera was used. In this paper, a high-reliability method was proposed to evaluate the separation between the layered composites using acoustic emission method. By comparing mechanical data and acoustic emission signals, the mechanical behavior obtained for each loading rate was determined so that the mechanical behavior of the composite material varied with the change in loading rate. The results show that, with increasing loading rates, the resin lost its elastic properties, and the specimen exhibited a more rigid behavior and is quite rigorous so that the fracture failure process is changed. The failure processes and crack growth rate was validated by use of acoustic emission signals. There was good agreement between the fracture toughness of accretion of acoustic emission signals with the experimental values.

Co-firing of biomass and fossil fuels in industrial furnaces is a suitable way to reduce the environmental impact from human activities, with acceptable investment. In this paper, the results of numerical simulation co-firing of sulfide concentrates and three auxiliary fuels including gasoil, kerosene and sawdust biomass are compared in the flash furnace copper smelting. For modeling of turbulent flow and combustion, RNG, k-ε model and probability density function model (pdf) have been used, respectively. This study has been carried out to investigate the furnace temperature and combustion pollutants distribution. The numerical simulation results show that the flame temperature resulting from the combustion of diesel fuel and sawdust as auxiliary fuel is the highest and lowest, respectively. In biomass combustion, despite that the flame temperature is low, but the NOx mass fraction increases because there is nitrogen in the sawdust chemical composition. Also in sawdust combustion that the oxygen content is high, the SO2 and SO3 sulfur pollutants increase in the high temperatures regions of the furnace and the lower temperature of the auxiliary fuel burner, respectively. Because SO2 is formed at high temperatures (> 1273K) and oxygen-rich and SO3 species is produced at relatively low temperatures with excess oxygen. The amount of CO emissions in sawdust combustion is much lower than the amount of combustion of diesel and oil.

Reference trajectory tracking and guarding against system disturbances and uncertainties are the important factors in the realm of lower limb exoskeleton robots control. Sliding mode PID is one of the robust controllers, which has a sliding manifold in the form of the PID controller. Chattering is the substantial predicament of the PIDSMC so that boundary layer around the sliding manifold is applied to eliminate the phenomenon. In this step, not only the chattering phenomenon is not eliminated but the robustness of the controller is also mitigated. In this study, supper twisting PID sliding mode controller (STPIDSMC) was used to eradicate the chattering phenomenon and enhancing controller robustness. The STPIDSMC robustness is protected indigenously and without defining the boundary layer, and the chattering phenomenon is reduced. Furthermore, to meet the external disturbances and uncertainties with unlimited amplitude, adaptive active force control method is combined by STPIDSMC as a modifying input control loop. In the active force control approach, the control input is online modified based on the estimation of moment inertia of the robot links. In order to accomplish maximum performance, control parameters were optimized using harmony search algorithm. In the optimal state, the performance of the proposed controller has been compared with PIDSMC and STPIDSMC that revealed the priority of the proposed controller rather than other controllers. The results indicate that the three error criteria, ITAE, ITASE, and IASE experience significant reduction about 39, 48, and 66 percent respectively compared to STPIDSM.

Acoustic doppler velocimetry (ADV) is a common measurement technique for flow field in open-channel flows. Since ADV is an intrusive measurement method, the presence of ADV probe may causes changes in flow structure and may intensify the turbulence in sampling volume which can affect the experimental results and analysis. To explore these effects, in this study, particle image velocimetry (PIV) was employed to measure the flow field with and without the presence of a side looking ADV probe in a compound channel. The results of this research showed that the intrusion of ADV in the flow field increases the streamwise velocity in the ADV measurement volume by 1.7%. The more enhanced effect is also notified in the secondary currents so that in the present test conditions, the presence of the side looking ADV causes a decrease in lateral velocity by 4.5 times and causes an increase in vertical velocity by 2.7 times. Investigation of the turbulent intensities showed that the presence of the side looking ADV causes an increment in streamwise turbulent intensity, while does not significantly affect the lateral and vertical components. Furthermore, examining vertical Reynolds shear stress measurement data with or without the presence of ADV showed that the measured values differ in the two states and the presence of ADV decreases 30% of Reynolds shear stress in place of the control volume.

Vibration analysis of the plate is an important topic in high-speed train body design. Because of the dynamic loads on plates which are used in the wagon body of the train, vibration analysis and determination of the amount of deflection and bending of the structure is important in wagon design. A plate which is used in the high-speed train is composite plate. Composite plates are considered because of many advantages relative to the other plates, such as low weight, high strength and cost-effective. In this paper, the nonlinear free vibration analysis of the used plate in the wagon body of high-speed trains has been presented. First, a three layers sandwich plate used for car body of high-speed trains has been transformed into a single layer equivalent orthotropic plate. Von-Karman theory and the Galerkin method have been employed to solving the equations of motion of the equivalent orthotropic plate. The nonlinear natural frequencies of the first four modes of the system have been determined using the numerical and variational iteration methods (VIM). Then the effect of different parameters on the value of nonlinear frequencies of the first four modes has been studied. The Difference lower than 0.1% is observed between the determined natural frequencies by VIM, with initial condition limited to zero, and natural frequencies determined by linear vibration. The results show that natural frequency is increased by increasing elasticity modulus of the face, the thickness of the core and the thickness of the face of the sandwich plate. In addition, because of nonlinearity of plate vibration equations, natural frequencies of composite plate are increased by increasing initial condition.

In this article, technical-economical investigation of heating and cooling loads in faculty of engineering of Arak university, Iran has been investigated. The tradiational and modern conditioner systems such as CAV, VAV, fan coil systems (all three systems are tested with direct fired absorption chiller and screw chiller), VRF system, split system and evaporative cooling system (in total 9 different systems) are designed for this faculty and are compared technically-economically. In this article, the costs of water, electricity and gas, purchase of systems and equipment, maintenance and repair of equipment of the examined systems are calculated and are compared with each other and then this work is done in different cooling loads, so that the results are not limited to the technical faculty of Arak university and can be used for all existing systems throughout the country. In the end, it is shown, that the evaporative cooling system has the lowest current and initial cost among the other systems. however, since this systems do not have the ability to determination of the amount of humidity and ideal temperature for the desired space can not be considered as an ideal and standard system. Also it is shown, that the current and initial costs of compression chiller are less than of absorption chiller and VAV systems have the better performance than CAV and fan coil systems, and the VRF system, after the evaporative cooling system, has the lowest current and initial cost among the examined systems, especially the split system. Due to the intelligent and optimal control of this system, it can be selected as an ideal system.

Quenching heat treatment is one of the most used processes in the industry, which has a great influence on the properties of materials. Accurate understanding of the effects of this process on the behavior of materials can be effective in better use of this process. In this research, the effect of quenching media on the mechanical behavior of wide used steel of AISI 1045 has been investigated and the residual stress created in the structure has been studied using the contour method. In this regard, three cooling environments of water, oil, and molten salt have been utilized, and after examining the strength and contour of hardness resulted by each cooling environment, the residual stresses have been investigated by the contour method. Also, the uncertainty of residual stresses in the environment with the most influencing factor has been evaluated. Investigation of the results shows that quenching in water can create higher hardness and strength, and more excessive compressive residual stress with greater penetration depth than the other environments. But cooling media of water creates more heterogeneous of the structure between the surface and the center of the piece, while quenching in a molten salt environment, with maintaining a structural homogeneity close to the annealing state, can increase the hardness and strength, and generate compressive residual stresses with a penetration depth of about 1.3 mm. Investigation of uncertainty for quenching in the water environment showed that the greatest error in the residual stresses was about 9%, and the error resulting from data smoothing had the most effect on the measurement of residual stresses by the cantor method.

In this paper, an atomic force microscope is modeled based on non-classical nonlocal theory and nonlinear vibration of the system is analyzed and controlled. In this modeling, the Hamilton principle is used to derive the governing equation of Euler-Bernoulli nanocantilever based on the Eringen nonlocal elasticity theory considering Von-Karman geometric non-linearity. In the next step, using the Galerkin method, the governing dynamics differential equation of the atomic force microscope is obtained in the presence of attractive and repulsive van der Waals forces. The governing nonlinear equation is solved by employing multiple time scales method, and primary and secondary resonance of the atomic force microscope is studied. In this regard, the frequency response and excitation amplitude curves of primary, superharmonic and subharmonic resonances are plotted for different values ​​of the nonlocal parameter. Accordingly, it is shown that primary, superharmonic and subharmonic resonances of atomic force microscope are significantly affected by the nonlocal parameter. The results show that the use of nonlocal theory is a fundamental necessity for analyzing nonlinear vibrations of the atomic force microscope. Then, in addition to dynamic analysis, the chaotic vibrations are completely controlled and removed in the nonlocal model of the atomic force microscope by designing and implementing the robust adaptive fuzzy controller. For this task, the robust adaptive fuzzy controller which is considered as a powerful method of chaos controlling is used in the nonlocal model of atomic force microscope. The obtained results are used in the design and control process of the atomic force microscope.

Air pollution is one of the consequences of industrial development that its severity is increasing day to day, due to the population growth and expanding urbanization, development of transport and fuel consumption increases. Awareness of air quality and trends of pollutants changes in different locations in a city can play an important and effective role in urban health management and macro policymaking. The first practical step in reducing the level of pollutants is the adequate knowledge of the pollutants details, including the type, amount and frequency of events throughout the year to determine the type of pollution and its source. In this paper, the results of hourly, daily, monthly and annual analysis of the various pollutants have been studied and scrutinized using the experimental data obtained from air quality measurement stations in the city of Saveh. In the end, In addition to the general solutions and suggestions for reducing the air pollution level in Saveh city,  two precise solutions include the construction of an underpass and construction of a high way have been presented which the geometry and detailed features of each one has been mentioned in the article.

Regarding the water and energy crisis, improving the efficiency of thermal systems and heat recovery, along with the use of desalination process, has attracted the attention of many researchers in recent years. For this purpose, thermal desalination process and solar collectors were used in steam power plants. In this study, an integrated structure for simultaneous generation of fresh water and power has been developed using a combination of solar collectors, steam power plant for power generation, ORC cycle, and thermal multi-effect desalination cycle. The integrated structure has the capacity of producing 762.6 kg / s of fresh water, 104.1 MW of power in the rankine cycle and 306.7 MW of power in a steam power plant. In this integrated structure, the efficiency of the steam power plant is 37.24% and the total exergy efficiency is 78.54%. Exergy analysis of the integrated structure shows that the highest destruction of exergy in solar collectors and heat exchangers are equal to 45.2% and 37.27%, respectively. The economic analysis of the developed integrated structure shows that the period of return is 3.838 years, and the prime cost of the product is 0.0325 $/kWh. Moreover, the impact of various parameters on the performance of the integrated structure was investigated using sensitivity analysis.

Today, the use of ultrasonic waves is expanding to the separation of particles or cells. One of the effective factors in the separating is cell deformation caused by ultrasonic waves. The most popular models used for deformation are the elastic and viscoelastic models. In this research, the cell has been modeled in a fluid environment under the influence of ultrasonic waves and deformations has been obtained. For this purpose, the Helmholtz equation that is a combine of the disturbance equations of sound waves and Navier-Stokes equation is solved and acoustic pressure is obtained. This pressure is then applied to the cell as deformation agent and the deformation is obtained using fluid-solid interactions modeling. Initially, deformation of the cell with elastic properties has been presented and validation has been conducted using comparison with the previous experimental researches. Finally, the deformation for the viscoelastic cell, which has so far not been used for deformation modeling in the acoustic field, has been obtained and presented. The results show that the viscoelastic model has the most compatibility with the experiment results. Also, the effect of frequency on the aspect ratio has been investigated. As the frequency ranges increased from 2 to 8 MHz, the aspect ratio is increased to 0.3.

Casing collapse is major problem of the oil fields which causes increase of costs to oil companies. This problem can be seen not only at drilling time in some formations but also after the completion and production may cause problems. exact prediction of collapse pressure is a very important factor in the casing design. Casing Collapse generaliy is a function of the geomechanical properties of the formation and the properties of the solid mechanics of the casing. One of the properties of solid mechanics that affects the testing of the collapse can be the ovality of the casing and the difference in the thickness of the casing and the existence of residual stress during the construction of the casing. In this paper, using numerical methods, the effect of each of the above-mentioned solid mechanics parameters and formation creep on the collapse of the casing  has been investigated. The results of this study indicate that pipe defects, such as casing ovality, eccentricity and the presence of residual stress, reduce the strength of the casing.
This reduces the resistance to the extent that the casing at the time of installation due to high plastic strain will collapse and also it was found that the pipe imperfections is more effective than rock salt creep in casing collapse.

In this research, micro-nano bubbles formation inside the venturi tube has been investigated using image processing technique and high-speed photography. For the purpose, two models of venturi tubes with different dimensions were made of light and transparent plexiglass and then they were tested with various water flow discharges. After the injection of air into the venturi tube, a high-speed camera has been used to capture images of two-phase flows passing through the venturi tube. The captured images were processed by MATLAB software. After investigating the results and obtaining the average diameter of the bubbles, the number of Micro-Nano bubbles and the velocity in the center of venturi tube have been calculated and the related graphs have been analyzed. The results show that reducing bottleneck length and increasing flow discharge inside the venturi tube lead to the formation of smaller bubbles and the number of Micro-Nano bubble increases.

The plug-in hybrid electric powertrain is a new technology and a suitable option to reduce the volume of pollutants in the city. The batteries in these systems can be plugged into an external source in addition to recharging by the combustion engine and regenerative braking. These vehicles have larger full electric range because of their relatively large batteries. As a result, the fuel consumption of these vehicles is low. The aim of this research is the development of a smart sizing and simulation software for plug-in hybrid electric vehicles. The size of the combustion engine, electric motor and battery are calculated in the software according to the input based on the performance of components. Then fuel consumption and emissions of pollutants are estimated in a standard drive cycle using a predictive controller. To verify the result from the software, a produced plug-in hybrid electric vehicle has been used. The software outputs correspond to the determined values for the vehicle. The proposed software is a useful tool for the early phase of the plug-in hybrid electric vehicle development stage.

In this research, nanostructured TiAlN coatings were applied on HSS substrate using cathodic arc evaporation method (CAE) in the different duty cycle values. Then the effect of duty cycle on the coating surface properties including surface morphology and structure, coating thickness and mechanical behavior of nanostructured coatings were investigated. X-ray diffraction (XRD) and scanning electron microscopy (SEM) were used to characterize the surface coatings. Also, micro indentation and adhesion test were utilized to evaluate the mechanical behavior. The results show that by changing the duty cycle, the macro-particles size and amount change which is effective on the roughness and morphology of the coatings. It is attributed to the electrical charge of macro-particles that are produced in the process which can be influenced by the structure. Also, the changes in grain size depend on the changes of duty cycle value. Furthermore, the mechanical properties of the coatings are affected by altering the duty cycle related to the deposition mechanism. The hardness value of TiAlN coatings increases from 3168 HV to 3817 HV when the duty cycle increases from 25% to 50%. But whit an increase in duty cycle from 50% to 75%, hardness reduced to 3582 HV. Consequently, it can be possible to find an optimum duty cycle value to achieve the best mechanical properties. Also, the minimum friction coefficient (0.44) and the minimum wear rate were determined for the TiAlN coating with the duty cycle of 75%, which it can be attributed to better smoothness and higher density of the coating.

Gasification technology is an important part of the clean coal technology. For further development of this technology, the processes inside the gasifier and the interactions of injected gas and solid fuel particles need to be understood. In the current study, a numerical simulation of an entrained flow coal gasifier is investigated using experimental operating conditions. The reactions and kinetic parameters of the gasification process are fitted using coal gasification data obtained from similar published papers. Comparison of the simulation results, experimental data and two other similar studies confirm the accuracy of the developed model. The focus of this study is on the accuracy of the models presented for the devolatilization process and the effect of the oxidizer change from oxygen to air on the gasifier performance. Four devolatilization models are investigated including chemical percolation devolatilization, single rate, Kobayashi and constant rate models. The predicted trends of species changes are similar in different devolatilization models but the amount of produced syngas is somewhat different depending on the accuracy of each model. The Kobayashi and constant rate models predicted the devolatilization rate lower than the other two models. The results obtained from the chemical percolation devolatilization model are more consistent with experimental data compared to other models, but require higher computational times. The use of air oxidizing agents rather than oxygen reduces the species mole concentration of produced syngas and hence reduces the gasifier efficiency.

Particle pollutants in indoor environment are a serious threat for human health. Therefore, it seems necessary to recognize and control the distribution of these particles. In the present study, the effect of air inlet angle of swirling diffusers in UFAD systems has been investigated on micron particles pattern distribution by considering the thermal comfort condition. For evaluating the flow and particle distribution, OpenFoam developed solver by the Eulerian-Lagrangian method and for two node model are used for predicting the flow and thermal conditions. Inlet angles are set in three cases: 30, 45 and 60. Based on the results, in all three cases, the TSENS index is in the comfort zone. But, by changing the inlet angle from 30 to 60, the vertical temperature difference can be reduced about 1 . The results showed that at inlet angle 30 and 60, the percentage of particles exited with 2.5 micrometers diameter were 32% and 55% of the total particles. In other words, increasing the inlet angle can leads to exit more amount of any size of particles from the room. In addition, by increasing particles size, larger particles removed faster from the breathing zone, and smaller particles will remain longer time in the air. In other words, smaller particles have a greater impact on the indoor air quality.