Acta Polytechnica doi:10.14311/AP.2018.58.0378 Acta Polytechnica 58(6):378–387, 2018 © Czech Technical University in Prague, 2018 available online at http://ojs.cvut.cz/ojs/index.php/ap THE USE OF DEM SIMULATION FOR CONFIRMING THE PROCESS OF PARTICULATE MATERIAL MIXING Peter Peciar∗, Oliver Macho, Roman Fekete, Marián Peciar Slovak University of Technology in Bratislava, Faculty of Mechanical Engineering, Institute of Process Engineering, Námestie Slobody 17, 812 31, Bratislava, Slovakia ∗ corresponding author: peter.peciar@stuba.sk Abstract. At the present time in the research of particulate materials, computer methods that work independently with individual particles are coming to the forefront. One such method is the Discrete Element Method – DEM, which is already widely used. Its usage, however, is complex, mostly due to the input data – what the properties of the researched material are, plus their interaction in processes. And despite the progress, it is still always necessary to validate the experimental equipment and to verify the individual simulations by an experimental measuring or by theoretical knowledge. This study focuses on the verification of the simulation of the mixing of a particulate material with the help of the DEM method, whereby simulations are compared with an experimental measurement and theoretic calculations. The theoretical calculation was carried out by the Novosad model, while the experimental material was granulated polyethylene with strictly defined mechanical-physical properties. Keywords: particulate material, mixing, theoretical knowledge, experimental measurements, DEM method. 1. Introduction Mixing of a particulate material is a very old process, but nonetheless relatively little studied. Research in this area is progressing quite slowly. Studies dealing with the effects of the forces acting during the process of mixing appear only sporadically and the majority of them focuses on the kinetics of mixing and computer simulations without any verification by experiments [1, 2]. Among the pioneers in the field of mixing belong D. Bagster and J. Bridgwater [3–5], who undertook measuring of the forces acting on a flat mixing element. In their first experiments, they used a static blade and the chamber performed a straight-line motion. They later designed a new mixing apparatus with a static glassed chamber in the shape of a prism. The blade moved in a reciprocating motion in a horizontal direction in the material while the station, equipped with measuring elements, recorded forces in both the horizontal and the vertical directions. The results of these researches pointed to a depen- dence of the measured forces on the properties of the particulate material, specifically the bulk density, material particle density, internal friction angle, wall friction angle, material porosity and blade geometry. This dependence was formulated into the relation H ρ(1 − ε)gLB2 = f ( z B , L B , v2 gB , ϕi, ϕ, ε ) , (1) where H is the horizontal force acting on the blade, ρ the particle density of the material, ε the porosity of the material, g the gravitational acceleration, L the blade length, B the blade height, z the immersion of the blade in the mixed material, v the blade speed, ϕi the internal friction angle, ϕ the wall friction angle between particulate material and blade material. Their experiments showed that the size of the hori- zontal force is not dependent on the porosity of the material or blade speed. The authors also found out that the value of the angle of the force from the hori- zontal, by which the blade acts on the mixed material, is approximately the same as the value of the wall friction angle of the particulate material and the blade material. In his work, J. Novosad [6, 7] carried out an in experimental monitoring of a course of mixing of sev- eral kinds of bulk solids with a blade mixer. In his experiments, he used glass cylindrical chambers in a vertical arrangement. His contribution was mainly the usage of mechanical properties of mixed particulate materials in relation to the calculation of mutually strong effects of mixed particulate materials and the mixing element. He found out that the most impor- tant variable characterizing particulate materials was the role of the internal and wall friction and bulk density. K. Kravárik [8] worked on the explanation of the force effects of particulate materials on a flat blade go- ing through the layer of its batch and also on methods of the calculation of the forces acting on a flat blade of rectangular shape moving at a constant speed in a horizontal direction. The investigated materials included granulated fer- tilizer, granulated polyethylene and sand. As the parameters according to which he monitored the de- pendence of the amount of force acting on the blade, he used the internal friction angle, the wall friction 378 http://dx.doi.org/10.14311/AP.2018.58.0378 http://ojs.cvut.cz/ojs/index.php/ap vol. 58 no. 6/2018 The Use of DEM Simulation angle of the blade material, the bulk density, the height and length of the blade, the depth of the blade immersion in the material and the blade speed. For a theoretical clarification of the formation of the shear zone in a particulate material before the blade, he compared various calculation models and, in experi- mental conditions, studied the shape of the shear zone before the advancing blade from the Novosad model and formulated a dimensionless, critical dependence to calculate the force acting on the blade: H ρSgLB2 = k1 ( z B )k2( v2 gB )k3( tanϕi tanϕ )k4 , (2) where H is the horizontal force acting on the blade, ρS the bulk density of the particulate material, g the gravitational acceleration, L the blade length, B the blade height, z the immersion of the blade in the mixed material, v the blade speed, ϕi the internal friction angle, ϕ the wall friction angle between the particulate material and blade material, k1, k2, k3, k4 the coefficients of critical dependencies for the calculation of the dimensionless force. On the basis of the results achieved, the author concluded that the dimensionless force significantly depends on the dimensionless immersion of the blade. The effect of the material’s internal friction angle is relatively significant and the effect of the dimensionless speed of the blade in the given extent has a serious effect on the size of the dimensionless horizontal force. The influence of the wall friction angle was eliminated by creating a static wedge area in front of the blade. At the present time, a simulation of tasks in a vari- ety of programs with the use of several computational methods is gaining prominence. However, for inves- tigating particulate materials and mixing elements, the DEM method remains the most used throughout the world. In 1997, Y. Muguruma, T. Tanaka, S. Kawatake and Y. Tsuji [9] experimentally studied the three-dimensional motion of particles in a rotation mixer with two baffles at different degrees of filling of the mixing chamber. They then simulated the ex- periments with the DEM method. They found out that the mixing curve, which coordinates the level of premixing and time through the DEM method, was in an accordance with the foregoing experiments. The result was the emergence of two distinct stages on the mixing curve. M. S. Siraj, S. Radl, B. J. Glasser and J. G. Khinast in their work [10], point out the influence of turning the blade travelling through the particulate material by the DEM method. The mixing equipment was made up of a prism-shaped static chamber with a moving blade with a variable slope. As the modelling material, glass balls of different diameters were used. The authors performed the simulation for several an- gles of turning the blade and for different particle size ratios. They expected that, to a certain degree, some particles in the bulk layers would show a tendency to rise to the top and others fall to the bottom, which represents segregation. If the larger particles rose to the top layer, they tended to remain there at the be- ginning of the experiment. With further simulations, they registered a significant influence from the mate- rial friction on the chamber walls. A large difference in the mixing of the particulate material was mainly visible with smaller angles of the blade slope (60◦, 90◦), while at an angle 140◦, there was no significant change between the simulation that included wall fric- tion and that which did not. In their simulations, they also tried to identify the normal and the tangential forces acting between the particles. The experiment showed that the extent of the normal and tangential forces was indirectly proportional to the slope angle of the blade, such that the higher forces were for the slope angles of the blade less than 90◦ and the lower forces were for the slope angles above 90◦. In their work, G. R. Chandratilleke, A. B. Yu and J. Bridgwater attempted [11] to verify, with the use of the DEM method, the studies of D. F. Bagster and J. Bridgwater. On the basis of the dimensional analysis of the horizontal force acting on a blade during the mixing of a particulate material, they searched for the influence of dimensionless parameters on the mixing process described in (1). The results of the usage of DEM simulations were fundamental studies on the understanding of the behaviour of round particles during the mixing with a flat blade. The effect of the slope of the head of the blade as well as the blade speed of mixing was studied. They discovered that the slope of the head of the blade has a real influence on the shape and the height of the pile of the material in front of the mixing blade. With a slanted blade head, a high heap in the front of the blade was formed, where a static zone was also created. In the investigation of the turning of the blade head to 135◦, they came to the conclusion that such a rotation significantly reduces the duration and height of the pile in front of the blade and that it is more suitable for mixing round particles. Horizontal forces acting on the blade were strongly affected by the blade head angle and were lessened by the increase of this slope. Due to the simulation, the authors showed that, with a constant blade head slope and with a constant depth of immersion, the horizontal force is not dependent on the blade speed during mixing. In the conclusion of their work, they tried to define the speed fields during the mixing and to apply them to the convective style of mixing, which would significantly reduce the calculation time. The authors S. Radl, D. Brandl, H. Haimburg, B. J. Glasser and J. G. Khinast added this to their research work [12], in which they attempted to com- pare the flow of the fluid of the bulk material on the basis of computer simulations and simple experiments of the movement of the investigated material. For the fluid flow, they used the dimensionless Froud number and applied this number for the flow of the particu- late material, by which they divided the moving bulk material into several areas simply according to the 379 P. Peciar, O. Macho, R. Fekete, M. Peciar Acta Polytechnica Figure 1. The Novosad model [13]. value of the Froud number. The authors came to the conclusion that the behaviour of the grained material can be described by the Froud number in analogy with the fluids. They found out that a recirculation zone is located in front of the blade and that this zone disap- pears when the critical value of the Froud number 0.5 is exceeded. The difference between the area beyond the blade in the mixing of particulate materials and fluids is caused by the rheology of granular materi- als. Furthermore, they came to the conclusion that if the mixing chamber is fuller, it leads to a significant change in the speed profile of the material before and after the moving blade, and that the recirculation zone disappears. The influence of the distance of the blade from the bottom of the vessel on the speed profiles was also significant, while a change in the slope angle of the blade did not have a significant effect on the speed of the particles in the upper part of the batch. 2. Energy aspect of the mixing process The energy aspect of the mixing process can be de- rived from the theory of a stress balance of a par- ticle material and the theory of Mohr’s circle. The Novosad model (Figure 1), which theoretically allows an approximate calculation of the force acting on the blade during mixing, also comes out of these theo- ries. This model takes into consideration the effect of the mechanical-physical properties of the particulate material. The calculation of the horizontal element of the force comes out of relation (3) for the calculation of the normal element of the burdening acting on the front side of the blade during a low blade immersion (z/B ≤ 1), where z is the distance from the surface of the material to the lower edge of the blade and B is the height of the blade) with consideration of the effect of the weight of the particulate material: pmn = FMqme 2ω tan ϕi + σC ( FMe2ω tan ϕi − 1 ) + σCS + 2 3ρSgB ( JMe3ω tan ϕi − SM + TM ( MMe3ω tan ϕi −NM )) , (3) where FM = sinµA2 cosµP 2 cosϕp1 cos(µA1 + ϕp1) sinµP 1 , (4) JM = sin2 µA2 sinµP 2 cosµP 2 cosϕp1 sin(µA1 + µA2) sinµP 1 cos(µA1 + ϕp1) , (5) SM = sin2 µA1 sinµA2 cosϕp1 sin(µA1 + µA2) cos(µA1 + ϕp1) , (6) TM = 1 9 tan2 ϕp1+1 × sin2 µA2 cosϕp1 sin(µA1 + µA2) cos(µA1 + ϕp1) cosϕi , (7) MM = 3 sin(µA1 + ω) tanϕi − cos(µA1 + ω), (8) NM = 3 sinµA1 tanϕi − cosµA1 (9) The calculation of the normal el ment of the burdening acting on the front side of the blade during a medium blade immersion (z/B > 1) with consideration of the effect of the weight of the particulate material comes out of the following relation: pmn = FSq1me 2ω tan ϕi − qp+ + σCS 2 3ρSgB ( RSe 3ω tan ϕi − SS + TS ( MSe 3ω tan ϕi −NS )) , (10) where FS = sinµA2 cos(µP 2 − ϕp2) cosϕp1 cos(µA1 + ϕp1) sinµP 1 , (11) RS = sin2 µA2 sinµP 2 cos(µP 2 − α) cosϕp1 sin(µA1 + µA2) sinµP 1 cos(µA1 + ϕp1) , (12) SS = sin2 µA1 sinµA2 cosϕp1 sin(µA1 + µA2) cos(µA1 + ϕp1) , (13) TS = 1 9 tan2 ϕp1 + 1 × sin2 µA2 cosϕp1 sin(µA1µA2) cos(µA1 + ϕp1) cosϕi , (14) MS = 3 sin(µA1 + ω) tanϕi − cos(µA1 + ω), (15) NS = 3 sinµA1 tanϕi − cosµA1 (16) For non-cohesive particulate materials, relations (3) and (10) also apply, substituting zero cohesion σC = 0Pa and σCS = 0Pa, in order to the pressure in the particulate material to safely rise over the value of the limit load, it is necessary to include in the calculation the safety factor k in the calculation. In the case of a non-cohesive particulate material, k = 1.5. The 380 vol. 58 no. 6/2018 The Use of DEM Simulation Figure 2. Scheme of experimental station. Particulate material Granulated polyethylene Diameter of particle* 3.4mm Mass of particle 2.13 · 10−2 g Volume of particle 2.03 · 10−8 m−3 Density of partilce 932.0 kgm−3 Bulk density (tapped) 607.9 kgm−3 Angle of repose 30.9° Angle of internal friction 31.8° Angle of wall friction (steel) 19.4° Table 1. Mechanical and physical properties of gran- ulated polyethylene [16] (*equivalent diameter of par- ticle by volume) horizontal force component is calculated according to relation HT EOR = kpmn(BL). (17) The values of the horizontal force components acting on the blade found in this way are the average force values acting on the whole blade area. 3. Experimental part Experimental measures were performed on an experi- mental measuring station made up of a static chamber, a frame, the drive, the measuring element and a mix- ing blade. The mixing chamber includes a glass side wall for a visual control of events during the mixing. The mobile measuring element is equipped with two force sensors to record the total horizontal and vertical components of the force. Granulated polyethylene was used for the experi- ment. The properties of granulated polyethylene are shown in Table 1. All the properties were measured by standard procedures [14, 15]. A set of measurements was performed for several blade immersion depths and different motion speeds of homogenizer blade motion. 4. Simulation of DEM method process The DEM is a numerical method for the calculation of the pressure and shifting of the individual elements Figure 3. Detail of measuring element. in a volume that contains a large amount of parti- cles. A granular material is modelled as a set of solid particles and the interactions between the individual particles. The program used for the modelling of this task was EDEM Academic™ from the DEM Solution UK Company. The goal of the experimental simula- tion was to look at the fields of particles in front of a moving mixing element and to establish the average force acting on the blade due to the movement of the particulate material layer. The calculation method was Herz-Mindlin with RVD Rolling Friction. Round particles with the properties of the experimental mate- rial were used for the simulations. For simplification of the simulation, the width of the chamber was reduced, but it was considered as if the interaction between the particles had actually occurred in contact with the walls. The properties of the materials and their interactions are shown in Tables 1 and 2. The simulations were performed for 4 blade immer- sions at a speed of 0.2m s−1. In Figures 4, 5, 6 and 7, a cross-section of the mixing chamber for the inves- tigated blade immersion depths with marking of the individual layers of the particulate material can be seen [17]. Graphic behaviours of the particle movement before and after the mixing element for the investigation of the blade immersion depth are shown in Figures 8, 9, 10 and 11. 381 P. Peciar, O. Macho, R. Fekete, M. Peciar Acta Polytechnica Model material Poisson’s ratio Solids density Shear modulus Steel (mixing chamber) 0.3 7800 kgm−3 8 · 1010 Pa Steel (mixing blade) 0.3 7800 kgm−3 8 · 1010 Pa Granulated polyethylene 0.245 932 kgm−3 1.3 · 108 Pa Interaction Coefficient restitution static rolling friction friction Granulated polyethylene – Granulated polyethylene 0.5 0.62 0.1 Granulated polyethylene – Steel (mixing blade) 0.35 0.35 0.1 Granulated polyethylene – Steel (mixing chamber) 0.35 0.62 0.1 Table 2. Properties of the materials and their interactions entered in the program EDEM Academic™ Figure 4. Cross-section of mixing chamber with dimensionless immersion of blade z/B = 0.5 at blade speed 0.2m s−1 in program EDEM Academic™. Figure 5. Cross-section of mixing chamber with dimensionless immersion of blade z/B = 1.0 at blade speed 0.2m s−1 in program EDEM Academic™. 382 vol. 58 no. 6/2018 The Use of DEM Simulation Figure 6. Cross-section of mixing chamber with dimensionless immersion of blade z/B = 1.5 at blade speed 0.2m s−1 in program EDEM Academic™. Figure 7. Cross-section of mixing chamber with dimensionless immersion of blade z/B = 2.0 at blade speed 0.2m s−1 in program EDEM Academic™. Figure 8. Graphic illustration of the movement of particles during mixing in mixing cross-section with dimensionless immersion of blade z/B = 0.5 at blade speeds 0.2m s−1 in the program EDEM Academic™. 383 P. Peciar, O. Macho, R. Fekete, M. Peciar Acta Polytechnica Figure 9. Graphic illustration of the movement of particles during mixing in mixing cross-section with dimensionless immersion of blade z/B = 1.0 at blade speeds 0.2m s−1 in the program EDEM Academic™. Figure 10. Graphic illustration of the movement of particles during mixing in mixing cross-section with dimensionless immersion of blade z/B = 1.5 at blade speeds 0.2m s−1 in the program EDEM Academic™. Figure 11. Graphic illustration of the movement of particles during mixing in mixing cross-section with dimensionless immersion of blade z/B = 2.0 at blade speeds 0.2m s−1 in the program EDEM Academic™. 384 vol. 58 no. 6/2018 The Use of DEM Simulation Figure 12. Dependency of horizontal force on dimen- sionless blade immersion. 5. Results 5.1. Influence of blade speed on horizontal force From the measured data, it follows that the force by which the particulate material acts on the moving blade at a given speed depends only slightly on its value for all blade immersions. This assumption was confirmed experimentally by Š. Gužela in the work [18] and K. Kravárik in the work [8] where they stated that this influence is also more expressive for speeds higher than the speed of the blade. The dependency of the horizontal force on the blade immersion for the investigated speeds is shown in Figure 12. 5.2. Influence on depth of blade immersion on horizontal force From the measured data, it follows that the force by which the particulate material acts on the mov- ing blade at a given speed depends, to a significant degree, on the immersion of the blade in the batch of the particulate material. A graphic illustration of the dependency of the horizontal force acting on the blade during the mixing on the speed for the given immersion in a batch is shown in Figure 13. From the experiments, it is clear that the depth of the blade immersion is a significant parameter affecting the size of the horizontal force acting on the blade. With the increasing depth, this force component increases. 5.3. Comparison of resulting horizontal force acting on a blade The comparison of the resulting horizontal force act- ing on a blade for individual blade immersion depths during mixing at blade speed 0.2m s−1 is shown in Figure 14. The calculated value of the horizontal force according to the Novosad model was experimen- tally measured on a HBM U9B sensor. The reworked measured values were also compared on ceramic mem- brane pressure sensors, with this procedure carried out according to [19] and stated in [13]. Figure 13. Dependency of horizontal force on blade speed during mixing in investigated dimensionless blade immersions in batches of particulate material. Figure 14. Dependency of horizontal force on blade speed during mixing in investigated dimensionless blade immersions in batches of particulate material. 6. Conclusion The experiments performed on the mixing of partic- ulate materials confirmed that the behaviour of a material during the mixing is in conformity with the theoretical knowledge and simulations by the DEM method in the EDEM Academic™ program. The horizontal force acting on a homogenizer blade was used for the comparison of theoretical calculations, experimental measurements and simulation. The mea- sured values of the horizontal force were in excellent conformity with the theory and these experiments can be considered successful, as Figure 14 shows. The DEM simulations performed for the support of the theoretical knowledge and experimental mea- sures demonstrated the contemporary trend of using calculation models, while the DEM method for mon- itoring the movement of individual particles showed to be suitable. Through the use of this method, it was possible to observe the movement and behaviour of the particles in front of and behind the mixing element over time, as it is presented in Figures 8, 9, 10 and 11. For confirming the accuracy of the move- 385 P. Peciar, O. Macho, R. Fekete, M. Peciar Acta Polytechnica Figure 15. Photographic illustration of particle movement during mixing in a mixing chamber with dimensionless blade immersion z/B = 1.0 at blade speed 0.2m s−1 in program EDEM Academic™ (legend: areas delineated by red lines show the individual slip areas in front of a moving blade constructed according to the Novosad model, and the thick black line illustrates the limit slip line constructed according to the Novosad model). ment of the individual particles in simulations, images were also created by photographic methods that com- pared the movement of the particles in front of and behind the mixing element. A comparison of particle movement in the DEM method and in experiments at z/B = 1.0 and blade movement speed 0.2m s−1 is shown in Figure 15. List of symbols µA1 Deflection angle of slide lines in the active state of limit equilibrium from the load plane [°] µA2 Deflection angle of slide lines in the active state of limit equilibrium from the load plane [°] µP 1 Deflection angle of slide lines in the passive state of limit equilibrium from the load plane [°] µP 2 Deflection angle of slide lines in the passive state of limit equilibrium from the load plane [°] ρS Bulk density [kg m−3] σC Normal cohesion of experimental material [Pa] σCS Normal cohesion contact between experimental ma- terial and construction material [Pa] ϕ Angle of wall friction between experimental particulate material and steel [°] ϕi Angle of internal friction [°] ϕp1 Virtual angle of wall friction particulate material contained in the active state of limit equilibrium [°] ω Rotation angle of slide lines in the transition area [°] g Acceleration due to gravity [m s−2] k Safety factor [1] pmn Normal component of the limit uniform load [Pa] rw Radius of logarithm spiral [m] qp Uniform surface load [Pa] q1m Virtual surrounding load acting on the passive areain particulate material [Pa] z Immersion of the blade [m] B Blade height [m] FM Dimensionless coefficient for the shallow immersion of the blade [1] FS Dimensionless coefficient for the medium immersion of the blade [1] HT EOR Horizontal force acting on the blade calculated by theory [N] JM Dimensionless coefficient for the shallow immersion of the blade [1] L Blade length [m] MM Dimensionless coefficient for the shallow immersion of the blade [1] MS Dimensionless coefficient for the medium immersion of the blade [1] NM Dimensionless coefficient for the shallow immersion of the blade [1] NS Dimensionless coefficient for the medium immersion of the blade [1] RS Dimensionless coefficient for the medium immersion of the blade [1] SM Dimensionless coefficient for the shallow immersion of the blade [1] SS Dimensionless coefficient for the medium immersion of the blade [1] TM Dimensionless coefficient for the shallow immersion of the blade [1] TS Dimensionless coefficient for the medium immersion of the blade [1] Acknowledgements This article was created with the support of the Ministry of Education, Science, Research and Sport of the Slovak Republic within the Research and Development Opera- tional Programme for the project "University Science Park of STU Bratislava", ITMS 26240220084, co-funded by the European Regional Development Fund. The authors wish to acknowledge the Scientific Grant Agency of the Ministry of Education, Science, Research and Sport of the Slovak Republic and the Slovak Academy of Sciences for the financial support of this research by grant VEGA 1/0276/17. References [1] J. Nečas, J. Hlosta, D. Žurovec, et al. Feeder type optimisation for the plain flow discharge process of an underground hopper by discrete element modelling. Advances in Science and Technology 11(3):246–252, 2017. doi:10.12913/22998624/76100. 386 http://dx.doi.org/10.12913/22998624/76100 vol. 58 no. 6/2018 The Use of DEM Simulation [2] J. Rozbroj, J. Nečas, D. Žurovec, et al. Simulation of material flow through a sample divider. Advances in Science and Technology 12(1):194–199, 2018. doi:10.12913/22998624/85663. [3] D. F. Bagster. The prediction of the force needed to move blades through a bed of cohesionless granules. Powder Technology 3:153–162, 1969/70. doi:10.1016/0032-5910(69)80067-7. [4] D. F. Bagster, J. Bridgwater. The flow of granular material over a moving blade. Powder Technology 3:323–338, 1969/70. doi:10.1016/0032-5910(69)80104-X. [5] D. F. Bagster, J. Bridgwater. The measurement of the force needed to move blades through a bed of cohesionless granules. Powder Technology 1:189–198, 1969. doi:10.1016/0032-5910(67)80036-6. [6] J. Novosad. Mechanika sypkých hmot. Praha: Institut pro výchovu vedoucích pracovniků ministerstva průmyslu, 1983. (In Czech). [7] J. Novosad. Studies on Granular Materials. I. Kinematics of granular materials mixed by mechanical impeller. Praha: Československá akadémia vied, 1963. [8] K. Kravárik. Energetické aspekty miešania partikulárnych látok. Dissertations thesis. STU, Faculty of Mechanical Engineering, Bratislava, 1999. (In Slovak). [9] Y. Muguruma, T. Tanaka, S. Kawatake, Y. Tsuji. Discrete particle simulation of a rotary vessel mixer with baffles. Powder Technology 93:261–266, 1997. doi:10.1016/S0032-5910(97)03281-6. [10] M. S. Siraj, S. Radl, B. J. Glasser, J. G. Khinast. Effect of blade angle and particle size on powder mixing performance in a rectangular box. Powder Technology 211:100–113, 2011. doi:10.1016/j.powtec.2011.04.004. [11] G. R. Chandratilleke, A. B. Yu, J. Bridgwater. A dem study of the mixing of particles induced by a flat blade. Chemical Engineering Science 79:54–74, 2012. doi:10.1016/j.ces.2012.05.010. [12] S. Radl, D. Brandl, H. Heimburg, et al. Flow and mixing of granular material over a single blade. Powder Technology 226:199–212, 2012. doi:10.1016/j.powtec.2012.04.042. [13] P. Peciar, M. Peciar, R. Fekete. Analysis of mixing of particulate material by a single blade. Solid State Phenomena 244:40–47, 2016. ISSN 1012-0394. [14] D. Geldart, E. Abdullah, A. Verlinden. Characterisation of dry powders. Powder Technology 190:70–74, 2009. doi:10.1016/j.powtec.2008.04.089. [15] D. McGlinchey. Characterisation of Bulk Solids. Wiley-Blackwell, 2009. ISBN: 978-1-405-14363-9. [16] P. Peciar. Pressure Distribution on the Particulate Material Homogenizer Blade. Dissertations thesis. STU, Faculty of Mechanical Engineering, Bratislava, 2014. ISBN 978-80-227-4146-0 (In Slovak). [17] P. Peciar. Research of the selected processes of particulate materials processing. Habilitation thesis. STU, Faculty of Mechanical Engineering, Bratislava, 2018. ISBN 978-80-227-4801-8 (In Slovak). [18] Š. Gužela, M. Peciar. Analýza napätia pôsobiaceho na lopatku postupujúcu vrstvu partikulárnej látky. STU, Faculty of Mechanical Engineering, Bratislava, 2006. ISBN: 80-227-2443-2, (In Slovak). [19] P. Peciar, M. Peciar, R. Fekete. Mobile universal measuring system of pressure in particulate material and calibration method of pressure sensors. Utility model no. 7283, Application no. 5068-2014. Industrial Property Office of the Slovak Republic, Banská Bystrica, Slovakia, 2015. 387 http://dx.doi.org/10.12913/22998624/85663 http://dx.doi.org/10.1016/0032-5910(69)80067-7 http://dx.doi.org/10.1016/0032-5910(69)80104-X http://dx.doi.org/10.1016/0032-5910(67)80036-6 http://dx.doi.org/10.1016/S0032-5910(97)03281-6 http://dx.doi.org/10.1016/j.powtec.2011.04.004 http://dx.doi.org/10.1016/j.ces.2012.05.010 http://dx.doi.org/10.1016/j.powtec.2012.04.042 http://dx.doi.org/10.1016/j.powtec.2008.04.089 Acta Polytechnica 58(6):378–387, 2018 1 Introduction 2 Energy aspect of the mixing process 3 Experimental part 4 Simulation of DEM method process 5 Results 5.1 Influence of blade speed on horizontal force 5.2 Influence on depth of blade immersion on horizontal force 5.3 Comparison of resulting horizontal force acting on a blade 6 Conclusion List of symbols Acknowledgements References