ACCEPTED MANUSCRIPT This is an early electronic version of an as-received manuscript that has been accepted for publication in the Journal of the Serbian Chemical Society but has not yet been subjected to the editing process and publishing procedure applied by the JSCS Editorial Office. Please cite this article as D. Jaćimovski, K. Šućurović, M. Đuriš, Z. Arsenijević, S. Krstić and N. Bošković-Vragolović, J. Serb. Chem. Soc. (2023) https://doi.org/10.2298/JSC230116016J This “raw” version of the manuscript is being provided to the authors and readers for their technical service. It must be stressed that the manuscript still has to be subjected to copyediting, typesetting, English grammar and syntax correc- tions, professional editing and authors’ review of the galley proof before it is published in its final form. Please note that during these publishing processes, many errors may emerge which could affect the final content of the manuscript and all legal disclaimers applied according to the policies of the Journal. https://doi.org/10.2298/JSC230116016J J. Serb. Chem. Soc.00(0)1-16 (2023) Original scientific paper JSCS–12233 Published DD MM, 2023 1 Mass transfer in inverse fluidized beds DARKO JAĆIMOVSKI1, KATARINA ŠUĆUROVIĆ1, MIHAL ĐURIŠ1*, ZORANA ARSENIJEVIĆ1, SANJA KRSTIĆ2 AND NEVENKA BOŠKOVIĆ-VRAGOLOVIĆ3 1Institute of Chemistry, Technology and Metallurgy-National Institute of the Republic of Serbia, University of Belgrade, Belgrade, Serbia, 2Vinča institute of Nuclear Sciences- National Institute of the Republic of Serbia, University of Belgrade, Belgrade, Serbia and 3Faculty of Technology and Metallurgy, University of Belgrade, Belgrade, Serbia (Received 16 January; Revised 27 February; Accepted 22 March 2023) Abstract: In this work, the coefficient of fluid-wall mass transfer in an inverse fluidized bed was determined using the adsorption method. The experiments were carried out in a column with a diameter of 45 mm with spherical and non- spherical particles of polypropylene and polyethylene with a diameter of 3.3-4.9 mm and a density of about 930 kg m-3. A diluted solution of methylene blue was used as a fluidization medium, which was adsorbed on part of the surface of the column on silica gel. The obtained results showed that the presence of particles during inverse fluidization does not contribute significantly to mass transfer compared to the influence of particles on transfer in conventional fluidized beds. Therefore, the pseudofluid concept was introduced into the analysis and an empirical correlation was performed to determine the mass transfer coefficient. The obtained results were compared with literature correlations for inverse and conventional fluidized beds. Keywords: inverse fluidization; fluid-wall mass transfer; pseudofluid. INTRODUCTION Fluidized bed contactors, because of their efficiency and transfer intensity, are often used in systems where contact between liquids and solid particles is required. When the solid phase has a lower density than the liquid phase, fluidization can be achieved by the liquid flowing into the column from the top and forming the bed in the opposite way compared to a conventional bed. Such beds are inverse fluidized beds. Inverse fluidization is most commonly used in wastewater treatment in bioreactors or in some adsorption processes. The beds formed in this way have been shown to be effective for formation and maintenance of biofilm. For the *Corresponding author E-mail: mihal.djuris@ihtm.bg.ac.rs; Tel.:+381 11 33 70 408 https://doi.org/10.2298/JSC230116016J A cc ep te d m an us cr ip t mailto:mihal.djuris@ihtm.bg.ac.rs https://doi.org/10.2298/JSC230116016J 2 JAĆIMOVSKI et al. practical application of inverse fluidized beds in various contactors, it is necessary to know the fluid dynamics of these systems and the mass transfer achieved. Wang et al.1 studied the removal of oil from an emulsion in water by inverse fluidization with a hydrophobic air gel. The authors used nanogel sizes: 0.5 mm-0.85 mm; 0.7 mm-1.2 mm; 1.7 mm-2.35 mm, the density was 64 kg m-3. The experiments were performed in columns with a diameter of 7.6 cm and lengths of 1.47 m and 0.77 m. The oil concentration is monitored using the COD (chemical oxygen demand). The authors have proposed a model that is consistent with the experimental data. It has been shown that the most important parameters affecting oil removal and thus mass transfer are granule size, bed height and fluid velocity. Inverse fluidization is applied in bioreactors where the necessary oxygen is supplied to the microorganisms in a three-phase fluidized bed. The application of inverse fluidization in FBBR reactors (Fluidized bed biofilm reactors) is shown in the work of Begum and Radha2 in which the aerobic biodegradation of phenol was studied using the microorganism Pseudomonas fluorescens. The authors performed the experimental tests in a column with a height of 105 cm and a diameter of 10 cm. Polystyrene particles with a diameter of 3.5 mm and a density of 863 kg m-3 were used as biofilm carriers. Mass transfer is monitored by the parameter COD. It was shown that the value of COD removal increases with the increase of gas velocity and the ratio of the volume of settled bed to the working volume (Vb/ VR). In the best case, COD reduction is 98.5%. Similarly, biological treatment of wastewater in an inverse fluidized bed system was performed by Sokol et al.3. These authors performed the removal of phenol, cresol, isopropylphenol, dimethylphenol, benzene and toluene from wastewater in a column with an internal diameter of 20 cm and a height of 6 m with polypropylene particles of density 910 kg m-3 as biofilm carriers in an inverse fluidized bed. Mass transfer was monitored using the COD parameter. It is shown that the operating conditions that give the best results are Vb/ VR = 0.55 and a gas velocity of 0.024 m s-1. The inverse fluidized bed was used for the treatment of wastewater containing starch4. The experiments were performed in a column with a diameter of 9.2 mm and a height of 1.6 m. Irregularly shaped polypropylene particles with a density of 870 kg m-3 were used as carrier particles. Mass transfer was monitored using the COD parameter. The reduction of the COD parameter increased with increasing air flow rate and with the duration of the experiment. Wastewater treatment in an inverse fluidized system was investigated by Karmanev and Nikolov5. Spherical polystyrene particles with an average diameter of 2.5 mm and density of 200 kg m-3 were used as carriers. It can be seen that biochemical oxygen demand (BOD) decreases A cc ep te d m an us cr ip t MASS TRANSFER IN INVERSE FLUIDIZED BEDS 3 with time. It has been confirmed that inverse fluidization systems are an effective system for treating wastewater with low pollutant concentrations. Knowledge of mass transfer is necessary for the application of inverse fluidized beds. One of the few papers dealing with the study of mass transfer in an inverse fluidized bed is the work of Nikov and Karamanev6 who showed that by increasing the density difference between fluid and particles, better mass transfer is achieved. The experimental data were compared with literature correlations for mass transfer and the following correlation was proposed: 0.33 Sh 0.28( Sc)GaMv= (1) The experiments were carried out with polystyrene and polyethylene particles with a diameter of 2.2-7.1 mm and a density of 80-930 kg m-3. The fluidization medium was water and an aqueous solution of polyethylene glycol with a concentration of 1.9% by mass. In their work, Kumar et al.7 monitored mass transfer in an inverse fluidized system using electrochemical methods. Mass transfer of lead ions was monitored by an electrochemical method in an inverse fluidized bed, where the fluidization medium was an aqueous solution of ferrous salts and ferricyanide in the presence of sodium hydroxide. Cylindrical particles with a diameter of 2.151 mm and a height of 5 mm and a density of 877.6 kg m-3 were used. Inverse fluidization was performed in a column with a diameter of 25.4 mm. As a result of these investigations, a correlation was established for the determination of the mass transfer factor: 0.30 Re 0.18 1 p D j   −   =   −  (2) The aim of this work is to investigate wall-fluid mass transfer in the presence of inverse fluidized particles using the adsorption method. The experimental results are presented as a function of particle and bed parameters and compared with mass transfers in conventional fluidized beds. EXPERIMENTAL The scheme of the experimental apparatus is shown in Figure 1. A cylindrical column with an inner diameter of 45 mm was used. At the top of the column there is a ring, on the inside of which there is a thin aluminum foil coated with silica gel. The column's inner diameter was not affected by the presence of the foil. The fluidization medium was a dilute aqueous solution of methylene blue with a concentration of 2·10-3 g L-1. The experiments were performed at a temperature of 20˚C. Spherical and non-spherical polypropylene and polyethylene particles were used for the experiments, and their properties are given in Table I. A cc ep te d m an us cr ip t 4 JAĆIMOVSKI et al. Figure 1. Scheme of the experimental apparatus (1- cylindrical column for inverse fluidization; 2-distributor; 3-grid; 4-fluidized bed particles; 5-inlet of the solution into the column; 6-vent; 7-piezometers; 8-exit of the solution from the column; 9, 11 -valves; 10-flow meter; 12-reservoir; 13-pump,’A’-ring, a-silica gel foil, b-particles and methylene blue solution) During the fluid flow from the top of the column, an inverse fluidized bed was formed in which methylene blue adsorbed on silica gel during the experiment, i.e. fluid-wall mass transfer occurred in the presence of inert particles. One of the conditions for the application of the method is a short adsorption time to avoid saturation of the adsorption surface, therefore the experiments lasted about 5 minutes8-10. TABLE I. Properties of particles Material Diameter mm Sphericity Density kg m-3 Polypropylene 3.3 1 935 Polypropylene 3.7 0.86 935 Polypropylene 4.4 0.85 937 Polyethylene 3.9 0.87 926 A cc ep te d m an us cr ip t MASS TRANSFER IN INVERSE FLUIDIZED BEDS 5 At the end of each experiment, the coloration of the surface of the silica gel foil, which is equivalent to the adsorbed amount of methylene blue, was analyzed using SigmaScan Pro software10-12. Based on the surface concentration, the mass transfer coefficient was determined according to the following equation: 0 p c k c t = (3) RESULTS AND DISCUSION Figure 2 shows the dependence of the mass transfer coefficient (equation 3) on the fluid velocity in the packed and in the inverse fluidized bed. The dependence shows a significant increase in the mass transfer coefficient at lower velocities, i.e. in a packed bed of particles. U / m s -1 0.00 0.01 0.02 0.03 0.04 0.05 0.06 k / m s -1 0 5e-6 1e-5 2e-5 2e-5 3e-5 3e-5 3.3 mm (PP) 3.7 mm (PP) 3.9 mm (PE) 4.4 mm (PP) single-phase flow Figure 2. Mass transfer coefficient in inverse fluidized bed and single-phase flow depending on fluid velocity After reaching the fluidized state, the mass transfer coefficient increases slightly or is approximately constant. There are no significant differences in mass transfer in the bed with different particles because the differences in diameters and sphericity are relatively small. For comparison, the mass transfer coefficient is also shown in the same diagram for a single-phase flow, i.e., when only liquid (without A cc ep te d m an us cr ip t 6 JAĆIMOVSKI et al. particles) flows through the column. It can be seen that the presence of inert fluidized particles has no significant effect on mass transfer, probably due to the fact that the particle density is close to the liquid density. A somewhat clearer dependence of the mass transfer coefficient on the bed parameters can be observed from the dependence of the Sherwood number Sh on the Reynolds number Re (Figure 3). From Figure 3, it can be seen that the Sherwood number Sh in the fixed bed increases sharply with the increase of the Reynolds number Re (i.e., with the increase of the velocity U). In the fluidized bed, the Sherwood number Sh also increases with the increase in Reynolds number Re, but somewhat less than in the packed bed. It is also evident that there is no effect of particle size and shape on the mass transfer coefficient in the inverse fluidized bed, since all the particles studied show a very similar dependence. Figure 3. Dependence of Sherwood on Reynolds number Figure 4 shows the dependence of the mass transfer coefficient on the bed porosity ε. Since the bed porosity is a function of the fluid velocity U, the dependence of k=f(ε) practically follows the dependence of k=f(U) (Figure 2). In a packed bed, where the porosity is constant, the mass transfer coefficient increases. During the transition to the fluidized bed state up to a bed porosity of about 0.7, the mass transfer coefficient continues to increase, while for higher values of bed porosity it becomes almost constant. Figure 4 shows the appearance A cc ep te d m an us cr ip t MASS TRANSFER IN INVERSE FLUIDIZED BEDS 7 of a slight maximum on the curve, which also occurs with conventional fluidization10. It is also evident that there is no visible difference in the intensity of mass transfer in the inverse fluidized bed with different particle types and sizes. 0.4 0.5 0.6 0.7 0.8 0.9 k / m s- 1 0.0 2.0e-6 4.0e-6 6.0e-6 8.0e-6 1.0e-5 1.2e-5 3.3 mm (PP) 3.7 mm (PP) 3.9 mm (PE) 4.4 mm (PP) Figure 4. Dependence of mass transfer coefficient on bed porosity for inverse fluidized bed (packed and fluidized bed) The adsorption method used in this study to determine the mass transfer coefficient is also useful for flow visualization. Figure 5 shows flow visualization images for different conditions (flows) of 4.4 mm diameter polypropylene particles in the bed. The chromatograms are shown for the same experimental conditions (time duration, temperature, methylene blue concentration, etc.) but for different liquid flow rates. Figures 5(a)-5(c) are characteristic of the liquid flow where the particles are stable, i.e., the packed bed. The flow patterns of the individual particles in the bed can be clearly seen. As the flow increases, an increase in the average color intensity of the silica gel foil surfaces can be seen. At the minimum fluidization velocity Umf (Figure 5(d)), the coloration of the silica gel foil surface is approximately uniform as the particles begin to oscillate slightly and mix the fluid. Figure 5(e) and Figure 5(f) are chromatograms taken in the fluidized bed and are characterized by the uniform coloration of the silica gel foil surface due to the intense mixing of the fluid by the particle fluidization. It can also be seen that the A cc ep te d m an us cr ip t 8 JAĆIMOVSKI et al. coloration of the surface is approximately the same at a higher flow rate in the fluidized bed. Figure 5. Flow chromatogram in packed bed (a-c), at minimum fluidization (d) and in fluidized bed (e-f) A comparison of the data obtained in inverse fluidized beds with those of a conventional fluidized bed as a function of bed porosity is shown in Figure 6. The results presented for the conventional fluidized bed were carried out in a bed of spherical glass particles with a diameter of 3 mm and a density of 2,500 kg m-3 with water as the fluidizing medium. These experiments were carried out under the same experimental conditions and in the same column. From the dependencies shown in Figure 6, it can be seen that the influence of the particles on mass transfer in conventional fluidization is significantly greater than the influence of the particles in inverse fluidization. This can be explained by the fact that significantly higher fluid velocities were used in conventional fluidization compared to inverse fluidization. Figure 6 shows the occurrence of a maximum at a porosity of about 0.7, which is consistent with our previous studies 9,13. A cc ep te d m an us cr ip t MASS TRANSFER IN INVERSE FLUIDIZED BEDS 9 0.4 0.5 0.6 0.7 0.8 0.9 k / m s- 1 0 1e-5 2e-5 3e-5 4e-5 5e-5 6e-5 3.3 mm (PP) 3.7 mm (PP) 3.9 mm (PE) 4.4 mm( PP) conventional fluidized bed (3mm) Figure 6. Comparison of the mass transfer coefficient in inverse and conventional fluidized beds as a function of porosity Based on the experimentally determined results for the mass transfer coefficient, the mass transfer factor jD was calculated. Figure 7 shows the dependence of the mass transfer factor jD on the Reynolds number Re (Figure 7(a)) and the porosity of the bed (Figure 7(b)). Both dependences show a decrease in the mass transfer factor jD with an increase in the Reynolds number Re and the porosity of the bed which is consistent with our previous study10. A cc ep te d m an us cr ip t 10 JAĆIMOVSKI et al. a) b) Figure 7. Dependence of the mass transfer factor on Reynolds number (a) and bed porosity (b) The diagram of the dependence of the mass transfer factor on the Reynolds number also shows the data obtained for a single-phase fluid flow (Figure 7(a)). It can be observed that the influence of the inverse fluidized particles on the mass transfer intensification is very small. This can be explained by the fact that in inverse fluidization the fluid velocity is low and therefore the collisions between the particles in the system are lower, so the influence on mass transfer is smaller. The drag force in conventional fluidized beds is defined as: ( ) 3 6 p d b g p f d F F F g   = − = − (4) in inverse fluidized bed is defined as: ( ) 3 , 6 p d inv g b f p d F F F g   = − = − (5) The ratio of the drag force in these two systems is: 3 3 , , p p fd d inv p inv f p dF F d     − =  − (6) Comparing conventional fluidization, in which particles with a diameter of 3 mm and a density of 2,500 kg m-3 were used, with inverse fluidization, in which particles with a diameter of 3.7 mm and a density of 935 kg m-3 were used (the fluid density in both cases is 1000 kg m-3), then: A cc ep te d m an us cr ip t MASS TRANSFER IN INVERSE FLUIDIZED BEDS 11 , 3.63d d inv F F = (7) It is found that the friction in the inverse fluidized bed is 3.6 times lower than in the conventional fluidized bed. Due to the small effect of particle-fluid friction in the inverse fluidized bed, the entire bed can be treated as a pseudofluid. All this indicates that although the conventional and inverse fluidized beds are described by the same equations, the frictional force is lower in the case of the inverse fluidized bed, so it is mainly used in biofilm reactors where the frictional forces do not disturb the formed biofilm, which is important for mass transfer in such reactors. The system behaves like a pseudofluid, and the frictional forces do not damage the biofilm. On the other hand, the particles are still present and provide gentle mixing, which is ideal for biosystems where microorganisms are involved in the reaction14,15. Re p pf pf pf pf d U  = (8) where pseudofluid density is defined as: ( )1pf f p   = + − (9) and pseudofluid dynamic viscosity as: 5(1 ) exp 3 pf f     −  =     (10) superficial fluid velocity can be calculated: f p pf f p G G U A A  = + (11) Since in a particulate fluidized bed the total particle motion in the column is zero (Gp=0), the velocity of the pseudofluid is equal to the superficial fluid velocity: f pf f G U U A = = (12) Based on the experimental results, the equation for the mass transfer factor was established: 0.48 Re , 13 Re 980 D pf pf j − =   (13) A cc ep te d m an us cr ip t 12 JAĆIMOVSKI et al. Figure 8 shows the dependence of the mass transfer factor jD on the Reynolds number of a pseudofluid Repf. The obtained equation agrees well with the experimental data. The mean absolute deviation of the experimental data from the data calculated by the established equation is 14.1%, while the relative deviation is -4%. Figure 8. Dependence of the mass transfer factor on the Reynolds number of a pseudofluid Table II shows some selected correlations for the determination of mass transfer factors in conventional and inverse fluidized beds and compares them with our experimentally obtained data. The comparison of the experimental data with the data calculated according to the above correlations is shown in Fig. 9. TABLE II. Comparison of literature correlations with experimental data and the obtained model Autor Fluidization Model σR, % σA, % Yutani et al.16 (Eqn. 14) conventional 0.40.4 Re D p j  − = +71.5 73.5 Marooka et al.17 (Eqn. 15) conventional 0.5 Re0.6 1 p D j   −   =   −  -6.93 14.9 A cc ep te d m an us cr ip t MASS TRANSFER IN INVERSE FLUIDIZED BEDS 13 Kunii & Levenspil18 (Eqn. 16) conventional 0.52 0.51 (1 ) Re Re Re D p p p j  −  = + −  -0.89 14.6 Nikov & Karamanev6 (Eqn. 17) inverse ( ) 0.330.28 Re D j GaMv= -82.4 81.5 Our corellation (Eqn. 13) inverse 0.48 Re D pf j − = -4.0 14.1 Figure 9. Comparison of the experimental data with the obtained model and the correlations in the literature It is interesting to note that the correlation of Yutani et al.16 (equation 14), was originally derived for conventional fluidized beds, and the correlation of Nikov and Karamanev6 (equation 17) for inverse fluidized beds, show significant deviation from the experimental data. Very good agreement with the experimental data is shown by the correlations for conventional fluidization of Marooka et al.17 (equation 15) and Kunii and Levenspiel18 (equation 16), as can be seen in Тable II. CONCLUSION In the paper, the fluid-wall mass transfer coefficient was experimentally determined in an inverse fluidized bed. The results show a significant increase in the mass transfer coefficient with increasing flow velocity in the packed bed, while a slight increase is observed in the fluidized bed. Compared to the mass transfers in a fluid flow without particles, it was found that the presence of particles does not contribute significantly to the intensification of the transfer, which is due to the low frictional force between particles and fluid. Based on the obtained results, the inverse fluidized bed was classified as a pseudofluid and a new correlation was presented, which represents the dependence of the mass transfer factor on the Reynolds number of a pseudofluid: 0.48 Re D pf j − = 13 < Repf < 980 (equation 13) The obtained experimental results deviate from the proposed correlation by less than 14.1%, with the mean relative deviation being -4%. When comparing the experimentally obtained data with the correlations from the literature, it was found that the correlations of Marooka et al.17 (equation 15) and Kunii and Levenspiel18 (equation 16) for conventional fluidization show very good agreement with the experimental data. NOMENCLATURE BOD biochemical oxygen demand COD chemical oxygen demand A cc ep te d m an us cr ip t 14 JAĆIMOVSKI et al. A cross section area, m2 cp surface concentration of methylene blue on adsorbent layer, kg m-2 c0 bulk concentration of methylene blue, kg m-2 dp,inv particle diameter in inverse fluidized bed, m dp particle diameter, m Dc column diametar, m FD drag force, N Fb buoyancy force, N Fg gravity force, N g gravitational acceleration, m s-2 Gf mass flow fluids, kg s-1 Gp mass flow particle, kg s-1 jD mass transfer factor k coefficient mass transfer, m s-1 t time, s U superficial fluid velocity, m s-1 Upf superficial pseudofluid velocity, m s-1 Vb/VR ratio of the volume of the settled bed to the working volume ε bed porosity μf viscosity fluid, Pa∙s μpf viscosity pseudofluids, Pa∙s ρf fluid density, kg m-3 ρp particle density, kg m-3 ρpf pseudofluid density, kg m-3 Ga ( 3 2 2 p f f d g /  ) , Galileo number Mv ( p f f ( ) / −  ), relative density Re ( f c f f U D /  ), Reynolds number Rep ( f p f f U d /  ), Reynolds number for particle Repf Reynolds number for pseudofluid Sh ( c AB kD / D ), Sherwoods number Sc ( f f AB / D  ), Schmidts number σA exp exp N cal X X N X  −        1 1 , absolute deviations σR exp exp N cal X X N X  −        1 1 , relative devations Acknowledgements: This work was financially supported by the Ministry of Science, Technological Development and Innovation of the Republic of Serbia (Grant No. 451-03- 47/2023-01/200026 and 451-03-47/2023- 01/200135]. A cc ep te d m an us cr ip t MASS TRANSFER IN INVERSE FLUIDIZED BEDS 15 ИЗВОД ПРЕНОС МАСЕ У ИНВЕРЗНО ФЛУИДИЗОВАНОМ СЛОЈУ ДАРКО ЈАЋИМОВСКИ1, КАТАРИНА ШУЋУРОВИЋ1, МИХАЛ ЂУРИШ1, ЗОРАНА АРСЕНИЈЕВИЋ1, САЊА КРСТИЋ2 И НЕВЕНКА БОШКОВИЋ-ВРАГОЛОВИЋ3 Институт за хемију, технологију и металургију -Институт од националног значаја за Републику Србију, Универзитет у Београду, Београд, Србија, 2Винча Институт за нуклеарне науке- Институт од националног значаја за Републику Србију, Универзитет у Београду, Београд, Србија и 3Технолошко- металуршки факултет, Универзитет у Београду, Београд, Србија У овом раду је одређиван коефицијент преноса масе флуид-зид у инверзно- флуидизованом слоју применом адсорпционе методе. Експерименти су вршени у колони пречника 45 мм са сферичним и несферичним честицама полипропилена и полиетилена пречника 3.3-4.9 мм и густине око 930 кг м-3. Као флуидизациони медијум коришћен је разблажени раствор метиленски плавог који је адсорбован на делу површине колоне на силикагелу. Добијени резултати показали су да присуство честица при инверзној флуидизацији не доприноси значајно преносу масе у поређењу са утицајем честица на пренос масе у конвенционално флуидизованим слојевима. Због тога је у анализу уведен концепт псеудофлуида и изведена је емпиријска корелација за одређивање коефицијента преноса масе. Извршено је поређење добијених резултата са литературним корелацијама за инверзну и конвенционалну флуидизације. (Примљено 16. јануара, ревидирано 27. фебруара, прихваћено 22. марта 2023.) REFERENCES 1. D. Wang, T. Silbaugh, R. Pfeffer Y.S. Lin, Powder Technol. 203 (2010) 298. (https://doi.org/10.1016/j.powtec.2010.05.021) 2. S.S.Begum, K.V. Radha, Korean J. Chem. Eng. 31 (2014) 436 (https://doi.org/10.1007/s11814-013-0260-z) 3. W. Sokol, A. Ambaw, B. Woldeyes, Chem. Eng. J. 150 (2009) 63 (https://doi.org/10.1016/j.cej.2008.12.021) 4. M. Rajasimman, C.Karthikeyan, Int. J. Environ Re. 3 (2009) 569 (https://doi.org/10.22059/IJER.2010.72) 5. D.G. Karamanev, L.N. Nikolov, Environ. Prog. 15 (1996) 3. (https://doi.org/10.1002/ep.670150319) 6. I.Nikov, D. Karamanev, ALChe J. 37 (1991) 781 (https://doi.org/10.1002/aic.690370515) 7. K.A. Kumar, G.V.S. Sarma, M. Vijay, K.V. Ramesh, Test Eng. Manage. 83 (2020) 14318 (http://www.testmagzine.biz/index.php/testmagzine/article/view/9657/7397) 8. S. Končar-Durđević, Nature 172, 878 (1953) 858 (https://doi.org/10.1134/S0036024409090246) 9. D. Jaćimovski, Diskontinualni granični sloj i analogije prenosa u pakovanim, fluidizovanim i transportnim sistemima tečnost-čestice, Doctoral dissertation, Faculty of Technology and Metallurgy, Belgrade, 2017 (http://phaidrabg.bg.ac.rs/o:17299) A cc ep te d m an us cr ip t https://doi.org/10.1016/j.powtec.2010.05.021 http://dx.doi.org/10.1007/s11814-013-0260-z https://doi.org/10.1016/j.cej.2008.12.021 https://doi.org/10.22059/IJER.2010.72 https://doi.org/10.1002/ep.670150319 https://doi.org/10.1002/aic.690370515 http://www.testmagzine.biz/index.php/testmagzine/article/view/9657/7397 https://doi.org/10.1134/S0036024409090246 http://phaidrabg.bg.ac.rs/o:17299 16 JAĆIMOVSKI et al. 10. N. Bošković-Vragolović, R. Garić-Grulović, Ž.Grbavčić, R.Pjanović, Russ. J. Phys. Chem. 83 (2009) 1550 (https://doi.org/10.1134/S0036024409090246) 11. SigmaScan Software, Jandel Scientific, USA ,1999. 12. M. Đuriš, T. Kaluđerović Radoičić, R. Garić-Grulović, Z. Arsenijević, Ž. Grbavčić, Powder Technol. 246 (2013) 98 (https://doi.org/10.1016/j.powtec.2013.05.009) 13. D.Jaćimovski, R.Garić-Grulović, N.Vučetić, R. Pjanović, N. Bošković- Vragolović, Powder Technol. 303 (2016) 68 (https://doi.org/10.1016/J.powtec.2016.09.025) 14. Ž.Grbavčić, Z. Arsenijević, R.Garić-Grulović, Powder Technol. 190 (2009) 283 (https://doi.org/10.1016/j.powtec.2008.08.005) 15. R. Garić-Grulović, Ž. Grbavčić, Z. Arsenijević, J. Serb. Chem. Soc. 70 (2005) 775 (http://dx.doi.org/10.2298/JSC0505775G) 16. N. Yutani, N. Ototake , L.T. Fan, Ind. Eng. Chem. Res. 26 (1987) 343 (https://doi.org/10.1021/ie00062a028) 17. S. Marooka, K. Kusakabe, Y. Kato, Int. Chem. Eng. 20 (1980) 433 18. D. Kunii, O. Levenspiel, Fluidisation Engineering, Wiley, New York, 1969, 195. A cc ep te d m an us cr ip t https://doi.org/10.1016/j.powtec.2013.05.009 https://doi.org/10.1016/J.powtec.2016.09.025 https://doi.org/10.1016/j.powtec.2008.08.005 http://dx.doi.org/10.2298/JSC0505775G https://doi.org/10.1021/ie00062a028