2018 Volume 35 Pages 160-171
This study focuses on the periodic motion observed in the FT4 Powder Rheometer and aims to develop a new testing method to evaluate powder flowability. This new method is based on the autocorrelation analysis of the torque measurements. The test results of powders belonging to different Geldart’s groups show that the torque measurements are periodic and that the cycle time (period) depends on the flowability. When the powder cohesiveness increases, the oscillations’ amplitude and the cycle time increase. Conversely, a free-flowing powder exhibits a liquid-like behavior and shows almost no periodic motion.
Powders or bulk solids are widely used in the industries. A good powder flowability is of crucial importance since poor powder flow can lead to many problems during manufacturing such as inaccurate dosing and off-specification products (Valverde J.M., 2012). Especially during a number of processes such as filling small bags or tablet press dies, the powder materials are under very low stress conditions and their flow behavior is more difficult to predict (Zafar et al., 2015).
To replicate and characterize the powder flow at low stress conditions, a number of characterization techniques have been developed over the last 50 years (Jenike A.W., 1980; Prescott J.J. and Barnum R.A., 2000; Schwedes J., 2003; Zafar et al., 2015).
The shear cell test is most commonly used for the characterization of powder flowability. It gives the relationship between the shear stress (which leads to flow of the material) and the normal stress (Bek et al., 2016). The yield locus in the low stress region and the tensile strength are obtained indirectly by extrapolating to zero of the measured yield locus. However, the yield locus tends to curve downwards at low stresses for cohesive powders so that the results can be overestimated (Vasilenko et al., 2013).
Other techniques that allow powder flow characterization under low stress conditions are reported in the literature. The Sevilla Powder Tester (SPT) and the Raining Bed Method (RBM) are two similar test methods that can directly measure the tensile yield stress at low stress levels. Both methods use a controlled gas flow to introduce small stresses in the powder bed. The tensile yield stress is then calculated by measuring the pressure drop across the powder bed at failure (Formisani et al., 2002; Valverde et al., 2000).
Another recently developed test method, the ball indentation method (BIM), measures the resistance to plastic deformation of powder materials at very low stress to assess the flowability (Hassanpour and Ghadiri, 2007). However, this method was originally designed for continuum solids and the extension to the testing of cohesive powder beds has only recently been analyzed (Pasha et al., 2013).
Zafar et al. (2015) proposed a comparison of these four techniques (shear cell, SPT, RBM and BIM), and showed that:
Quintanilla et al. (2006) studied the behavior of a series of granular materials of different cohesiveness in a half-filled and slowly rotating horizontal drum. The authors used a CCD camera and digital image processing to detect and characterize the avalanches. By statistical analysis, they studied the periodicity of the signals as well as the probability distributions of the avalanche size and the maximum angle of stability. The results showed that the system cohesiveness was correlated to these probability distributions: the larger the cohesion, the larger the average size and maximum angle of stability of the avalanches.
More in-depth studies and new testing methods still need to be developed to better characterize powder flow properties under low stress conditions. This work aims to develop a new dynamic testing method to evaluate powder flowability under unconfined flow conditions. This new method is based on the autocorrelation analysis of the torque measurements obtained with the FT4 Powder Rheometer.
Six powders belonging to Geldart’s groups A, B and C were used in this study. The properties of these six powders are summarized in Table 1. The particle size distributions were measured using a laser diffraction size analyzer (Mastersizer 2000, Malvern Instruments) and the particle density by a helium pycnometer Accupyc 1330 (Micromeritics). In addition, a liquid honey was used in this study as a reference material to simulate an “ideal” flow behavior. As it will be shown later in this paper, the liquid honey is viscous enough to obtain significant torque measurements. At the same time, these torque measurements are close to those corresponding to a “perfectly flowing” granular material. As a consequence, it will be considered as a “reference material” and will be used to simulate an “ideal” flow behavior.
Physical properties of the powders.
powder | d[3, 2] (μm) | d[4, 3] (μm) | d[10] (μm) | d[50] (μm) | d[90] (μm) | span | ρ (kg/m3) | shape | Geldart’s classification |
---|---|---|---|---|---|---|---|---|---|
CaCO3 (a) | 3.4 ± 0.5 | 11.7 ± 1.7 | 1.2 ± 0.2 | 8.6 ± 1.4 | 26.1 ± 3.9 | 2.9 | 2750 | angular | C |
CaCO3 (b) | 4.0 ± 0.3 | 15.9 ± 1.4 | 1.5 ± 0.2 | 12.3 ± 1.1 | 35.2 ± 2.9 | 2.8 | 2750 | angular | C |
CaCO3 (c) | 8.4 ± 0.2 | 51.5 ± 1.2 | 4.9 ± 0.4 | 35.2 ± 0.3 | 104.6 ± 2.1 | 2.8 | 2750 | angular | C |
ZrO2 | 23.7 ± 0.5 | 60.1 ± 0.5 | 12.8 ± 0.2 | 55.1 ± 0.4 | 113.0 ± 1.2 | 1.7 | 5790 | angular | A |
SiC | 79.1 ± 0.7 | 66.0 ± 3.6 | 43.4 ± 0.9 | 73.7 ± 0.5 | 122.7 ± 0.4 | 1.1 | 3056 | angular | A |
sand | 342.2 ± 15.7 | 287.1 ± 8.1 | 186.3 ± 6.6 | 306.3 ± 11.4 | 517.6 ± 20.6 | 1.1 | 2546 | angular | B |
Values after “±” indicate the 95 % confidence intervals;
span = (d90 – d10)/d50
Tests on the FT4 powder rheometer are based on a sophisticated twisted blade which passes through the powder bed along a predetermined helical path (Fig. 1). When the blade moves clockwise (conditioning cycle), the powder is lifted up and conditioned to obtain a uniform and reproducible state. When the blade moves anticlockwise (test cycle), it creates a compaction zone which is a localized high-stress region in front of the blade and forces the powder to flow (Leturia M. et al., 2014).
Test cycle (left) and conditioning cycle (right) with the FT4 Powder Rheometer.
The FT4 standard “dynamic test” measures the evolution of the torque and force for a downward flow pattern in a conditioned controlled volume of powder. The energy required to induce this movement is then calculated (by integration of the force and torque measurements) to evaluate the flow properties of the powder (Conesa C. et al., 2004). This test is divided into two main parts.
The first part, called “stability test”, corresponds to seven identical complete cycles (conditioning + test cycles) at a constant blade tip speed of 100 mm/s. The objective is to detect any possible change in the powder flow properties after repeated tests (attrition, agglomeration, etc.). The stability test also allows a stable energy level to be achieved. The basic flow energy (BFE), as defined by Freeman Technology, corresponds to the stabilized flow energy. It represents the energy required to force a conditioned powder sample to flow during a downwards test:
(1) |
Usually, fine powders have relatively high interparticulate forces compared to gravitational forces. As a consequence, they behave in a cohesive way: Agglomerates are formed and a large amount of air is trapped inside the bulk powder (Freeman R., 2007). In this case, particles forced to flow at the blade face can be accommodated by the voids that exist between the agglomerates (high compressibility). The resulting flow or stress transmission zone is relatively localized. Thus the displaced volume of powder is small and the measured BFE value is particularly low. Powders made of larger particles often behave in a non-cohesive manner due to the minimal effect of cohesive forces compared to gravitational forces. Less air is trapped in the bulk powder resulting in a relatively lower compressibility. The flow zone is transmitted far ahead of the blade (Freeman R., 2007). In this case, the displaced volume of powder is large and the measured flow energy is high.
The second part, called “variable flow rate test”, corresponds to 4 complete cycles at decreasing blade tip speeds (100, 70, 40 and 10 mm/s). The objective is to evaluate the powder sensitivity to the flow rate. The flow rate index (FRI) corresponds to the ratio of the flow energy at 10 mm/s blade tip speed to the flow energy at 100 mm/s blade tip speed:
(2) |
Many factors can affect the value of the FRI index, e.g.:
Although the combination of the values of the BFE and FRI helps to identify the powder flow properties, the standard dynamic test has many limitations for the classification of powder flowability.
Fig. 2 and Table 2 present the standard dynamic test results and the corresponding BFE and FRI values of four powder samples:
Results of the FT4 standard dynamic test.
Results of the FT4 standard dynamic test.
powder | BFE (mJ) | FRI |
---|---|---|
CaCO3 (a) | 715.4 | 2.0 |
CaCO3 (b) | 816.6 | 2.2 |
CaCO3 (c) | 1696.2 | 1.9 |
ZrO2 | 2426.1 | 2.0 |
In both cases, it is difficult to make a conclusion about which powder has better flowability. Therefore, the standard dynamic test on the FT4 cannot be used to compare two different powders and it can be difficult to distinguish between very similar materials.
3.2 Observation of a periodic motionIn order to optimize the FT4 standard procedures, the powder flow behavior during testing was closely observed (Fig. 3) and the resulting raw data were analyzed (normal force and torque measurements before integration for the energy calculation). For cohesive powders, the blade leaves an empty space behind it when passing through the powder bed. This empty space grows in size as the blade moves forward and collapses when it reaches a maximum volume. It then grows in size and collapses again, following a cyclic pattern. It was also observed that this empty zone had a smaller maximum volume at a lower blade tip speed. For non-cohesive powders, the blade does not leave any significant empty space behind it and particles flow around the blade, almost like a liquid. The cyclic pattern is also reflected in the torque measurements where oscillations can be observed (Fig. 4). These oscillations are much more significant in a cohesive sample compared to a free-flowing powder.
Flow behavior of non-cohesive powders (sand, left) and cohesive powders (CaCO3 (a), right).
Torque variation during a standard dynamic test (blade tip speed = 70 mm/s).
A “verification test” was performed to reveal the causes of this cyclic pattern. In this test, only the surface of the powder bed was forced to flow so that the flow behavior could easily be observed. A 160-ml conditioned sample of CaCO3 (a) was placed in the 50-mm ID cylinder vessel. The blade was moved just below the powder surface (70 mm from the bottom) and rotated for 60 seconds at 4 different blade tip speeds (100 mm/s, 70 mm/s, 40 mm/s and 10 mm/s). At the same time, a video camera recorded the powder motion around the blade. Since the blade stayed at the same height, there was no significant normal force. Under these conditions, the total energy needed to move the blade through the powder is proportional to the torque acting on the blade. The torque oscillations were evaluated through an autocorrelation analysis, which is a mathematics tool commonly used in signal processing and time-series analysis to investigate repeating patterns (Dunn and Patrick F., 2005)
In statistics, the autocorrelation function for a continuous-time process, Xt, is the correlation of the process with itself at a later time. It represents the similarity between observations as a function of the time lag (Box G.E.P. et al., 1994):
(3) |
where “E” is the expected value operator. For a stationary process, the mean μ and the variance σ2 are time-independent and the autocorrelation depends only on the time lag τ. In other words, the correlation depends only on the time distance between the two chosen values but not on their position on the time axis. Particularly, the autocorrelation function of a periodic function has the same period as the original function, which is the key to investigating the repeating patterns. The values of the autocorrelation function (ACF) should be in the range [−1, 1], with 1 indicating a perfect correlation (the signals exactly overlap when time is shifted by τ) and −1 indicating a perfect anti-correlation (Priestley M.B., 1982). For a discrete process with n observations such as the measured torque series in this test, an estimate of the autocorrelation may be obtained as (Chatfield C., 2004):
(4) |
Fig. 5 presents the time evolution of the torque measurements and the corresponding autocorrelation functions (ACF) for the “verification test”. On the one hand, a series of oscillations in the torque measurements are observed and the autocorrelation function brings to light the existence of a cyclic pattern. At high blade tip speeds, the CaCO3 (a) material has a period corresponding to one rotation of the blade. It should be noted that small peaks (small “events”) can be found at each half-rotation of the blade. As the blade tip speed decreases, it can be observed that:
Torque measurements and autocorrelation function (ACF) of CaCO3 (a) for the verification test.
The test video was analyzed frame by frame. The images clearly show (Fig. 6) that a “block” of powder moves in front of the blade while a large empty space is left behind it.
Image of the verification test (the blade rotates in anticlockwise direction).
During the “verification test”, the “powder block” in front of the blade was observed in the form of several agglomerates (Fig. 7). Whenever some agglomerates fall off the blade (as the big one in the red frame in Fig. 7-A), a sudden drop is observed in the torque measurements (Fig. 7-B). The torque goes up again when this agglomerate is picked up by the other side of the blade in Fig. 7-C, D. Since the agglomerates are not exactly of the same size, the oscillations do not have the same amplitude. In Fig. 7-E, F, G, H, torque oscillations can also be observed as several small agglomerates fall off the blade one after another.
Torque variation and agglomerates behavior (blade tip speed = 70 mm/s).
Indeed, while the blade pushes the powder to flow, an unstable “compacted powder zone” builds up in front of each side of the blade, forming two “powder blocks” and two empty spaces behind them. Since these powder blocks are unstable, they may eventually collapse to fill the empty spaces. The fallen agglomerates stay in the empty space for a short time and are picked up by the other side of the blade. The second powder block then increases in volume and becomes less stable. It will finally collapse like the first powder block and this flow pattern is repeated continuously.
The flow speed has an important impact on this periodic motion:
This behavior could be explained by the fact that the collapsing process is not immediate as a given powder has a certain response time to plastic flow deformation. As a consequence, there may be a critical flow speed above which the response time is longer than a 1/2 rotation, resulting in a cycle time of one complete rotation of the blade.
A specific testing procedure was developed on the FT4 Powder Rheometer. First, a 160-ml conditioned sample was placed in the 50-mm ID cylinder vessel. Then the twisted blade was set to rotate in the powder bed at a height of 30 mm from the bottom with a blade tip speed of 100 mm/s for 200 seconds. The blade tip speed was then gradually decreased from 100 mm/s to 10 mm/s (by 10 mm/s steps for each speed reduction), and the rotation was repeated for 200 seconds at each step. A conditioning cycle was performed between each step.
The torque needed to induce flow was measured continuously. To ensure a stable powder flow (steady state), only the data measured during the last 50 seconds were analyzed for each step. The autocorrelation function was then used to determine the cycle time and correlation level. By this means, the periodic motion was characterized and the powders could then be classified depending on their behavior.
4.1 Results obtained with the proposed testing methodThe liquid honey was first used to simulate a regular flow behavior (reference material). Indeed, it was viscous enough to obtain significant torque measurements while behaving like an ideal material under the testing conditions. Fig. 8 (top) presents the variations of the required torque to force the liquid honey to flow at different flow rates. The measured torque is very stable (small oscillations) and the corresponding ACF values (Fig. 8 bottom) are always close to 0. This means that the torque needed to force a viscous liquid to flow at a given speed is a constant and no periodic motion can be detected with the FT4. From this point of view, the liquid honey can be regarded as an “ideal” material (in this flow rate range). On the other hand, it can also be noted that the average torque (Fig. 8, center) is proportional to the flow speed. This corresponds to the flow behavior of a Newtonian fluid.
Torque variation as a function of the blade tip speed for the liquid honey and autocorrelation analysis.
The CaCO3 (a), SiC and sand powders were first characterized by this method, since each of them belongs to a different Geldart’s group. Visual observations during the tests showed that:
Fig. 9 presents the relative torque measurements of the 3 powder specimens (+ honey) as a function of the number of rotations of the blade (at a blade tip speed of 100 mm/s). The relative torque is defined as:
Relative torque evolution ( (T–
(5) |
where
The results show that the torque measurements can be correlated to the system cohesiveness:
These results can be related to those obtained by Quintanilla et al. (2006), who used a half-filled and slowly rotating horizontal drum. Indeed, the collapse of the powder blocks can be assimilated to the avalanches in the rotating drum. In both cases, a periodic behavior is observed and the periodicity of the signals can be correlated to the system cohesiveness. Also, in both studies, it is shown that a higher cohesion leads to larger events (collapses or avalanches).
4.2 Periodic motion analysis and powder flowability classificationDuring the tests, the torque was recorded every 0.04 seconds and the 50-second measurements could then be treated as a discrete process with 1250 observations for the autocorrelation analysis. The autocorrelation functions of the three powders are presented in Fig. 10. These graphs indicate that the flow within the powder bed exhibits the same cyclic pattern as the flow on the surface (Fig. 5). However, the periodic motion on the surface shows a lower correlation level, which might be caused by the surface effects.
Autocorrelation function at different blade tip speeds. a) blade tip speed = 100 mm/s; b) blade tip speed = 70 mm/s; c) blade tip speed = 40 mm/s; d) blade tip speed = 10 mm/s.
These ACF curves show that:
The cycle times (expressed in the number of rotations of the blade) of the three samples at different flow speeds are summarized in Fig. 11. Based on previous results, a threshold value of 0.2 was used to differentiate periodic from non-periodic flow patterns. All the processes with a correlation level lower than 0.2 were ignored and considered as non-periodic flow patterns. For convenience, the cycle time of such processes was then considered as 0.
Cycle time at different blade tip speeds.
The results of the autocorrelation analysis indicate that:
These results suggest that for any powder material, there are three flow behaviors:
These results suggest the existence of two flow transition speeds (Fig. 12): TS[1, 2] is the transition speed from flow type 1 to flow type 2 and TS[2, 3] is the transition speed from flow type 2 to flow type 3. For a given powder, the values of TS[1, 2] and TS[2, 3] depend on its response time to plastic flow deformation. Therefore, the values of the two transition speeds can be used to identify the powder flow properties:
Other tests were also performed to verify these conclusions. As explained in section 3.1, the standard dynamic test cannot discriminate the flowabilities between CaCO3 (a), CaCO3 (b), CaCO3 (c) and ZrO2. The results of the new testing method clearly show these differences (Fig. 13):
Transition speeds and powder flow behaviors in the FT4 Powder Rheometer.
Test results of the new testing method: autocorrelation function at 100 mm/s (left); cycle time at different blade tip speeds (right).
In this study, a new flowability testing method is presented. This method is based on the FT4 standard dynamic test and quantifies the powder periodic motion under forced flow conditions. The test results prove that the dependence of the torque cycle time on the flow rate can be linked to the powder flow properties: when the powder cohesiveness increases, the autocorrelation level and cycle time increase; conversely, a free-flowing powder tends to act like a liquid and shows no periodic motion.
This work shows that this new method has a good capability to identify and classify powder flow properties. The flowabilities of different powders which cannot be discriminated with the standard dynamic method of the FT4 Powder Rheometer are easily identified with this new method
Further studies still need to be done to optimize this method such as the influence of the blade height (30 mm in this study), the variation of the periodic motion in a wider flow speed range, as well as the influence of the blade shape.
Basic Flow Energy, [mJ]
d[10]10 % of the population lies below the d[10], [m]
d[3, 2]Mean surface-volume diameter (Sauter mean diameter), [m]
d[4, 3]Volume or mass moment mean (De Brouckere mean diameter), [m]
d[50]Mass median diameter, [m]
d[90]90 % of the population lies below the d[90], [m]
EExpected value operator
FRIFlow Rate Index, [−]
kObservation number in the discrete time process, [−]
nTotal observations of the discrete process
ACFAutocorrelation function, [−]
Autocorrelation function of a discrete process, [−]
span(d[90]-d[10])/d[50], [−]
tTime, [s]
TTorque, [mN m]
Average torque, [mN m]
TS[1, 2]Transition speed between flow type 1 and flow type 2, [mm/s]
TS[2, 3]Transition speed between flow type 2 and flow type 3, [mm/s]
XSignal of a continuous-time process
μMean
ρParticle density, [kg/m3]
σStandard deviation
τTime lag, [s]
Ming Li
Ming Li is a Ph.D. student in the Chemical Engineering Department at Compiègne University of Technology (Sorbonne Universities). His research is focused on the comparison of powder characterization methods and the improvement of testing protocols for a better evaluation of powder flow properties.
Mikel Leturia
Mikel Leturia is a chemical engineer from Compiègne University of Technology (2009). He prepared a Ph.D. thesis in the Chemical Engineering Department at Compiègne University of Technology on the experimental study and modeling of fluidizedbed reactors using cohesive powders (2013). He is currently associate professor in the same university and his research focuses on powder technology including fluidization, caking and flowability testing methods.
Khashayar Saleh
Khashayar Saleh prepared a Ph.D. thesis on the coating of fine powders in the Chemical Engineering Laboratory of Toulouse and obtained his Ph.D. degree in 1998 from the Institute National Poly-technique de Toulouse (France). Professor Khashayar Saleh is currently the director of the Chemical Engineering Department at Compiègne University of Technology (Sorbonne Universities). His work is focused on powder technology including size enlargement technology, powder electrification and powder characterisation methods.