The Stokes Number Approach to Support Scale-Up and Technology Transfer of a Mixing Process

doi:10.1208/s12249-012-9818-z

AAPS PharmSciTech. 2012 September; 13(3): 928–933.

Published online 2012 June 26. doi: 10.1208/s12249-012-9818-z

PMCID: PMC3429663

The Stokes Number Approach to Support Scale-Up and Technology Transfer of a Mixing Process

Tofan A. Willemsz, Ricardo Hooijmaijers, Carina M. Rubingh, Henderik W. Frijlink, Herman Vromans, and Kees van der Voort Maarschalk

Department of Pharmaceutical Technology and Biopharmacy, University of Groningen, A. Deusinglaan 1, 9713 AV Groningen, The Netherlands

Pharmaceutical Sciences and Clinical Supply, Merck MSD, PO Box 20, 5340 BH Oss, The Netherlands

Quality of Life, TNO, PO Box 360, 3700 AJ Zeist, The Netherlands

Department of Pharmaceutics, Utrecht Institute for Pharmaceutical Sciences (UIPS) Utrecht University, PO Box 80082, 3508 TB Utrecht, The Netherlands

Competence Center Powders and Formulations, Purac Biochem, PO Box 21, 4200 AA Gorinchem, The Netherlands

Kees van der Voort Maarschalk, Phone: +31-412-662342, Fax: +31-412-662617, Email: tofan.willemsz/at/merck.com.

Corresponding author.

Received February 8, 2012; Accepted June 8, 2012.

Abstract

Transferring processes between different scales and types of mixers is a common operation in industry. Challenges within this operation include the existence of considerable differences in blending conditions between mixer scales and types. Obtaining the correct blending conditions is crucial for the ability to break up agglomerates in order to achieve the desired blend uniformity. Agglomerate break up is often an abrasion process. In this study, the abrasion rate potential of agglomerates is described by the Stokes abrasion (St_Abr) number of the system. The St_Abr number equals the ratio between the kinetic energy density of the moving powder bed and the work of fracture of the agglomerate. In this study, the St_Abr approach demonstrates to be a useful tool to predict the abrasion of agglomerates during blending when technology is transferred between mixer scales/types. Applying the St_Abr approach revealed a transition point between parameters that determined agglomerate abrasion. This study gave evidence that (1) below this transition point, agglomerate abrasion is determined by a combination of impeller effects and by the kinetic energy density of the powder blend, whereas (2) above this transition point, agglomerate abrasion is mainly determined by the kinetic energy density of the powder blend.

Key words: dry mixing, scale-up, stokes number

INTRODUCTION

One of the challenges during scale-up and technology transfer of a mixing process is that there usually are considerable differences in blending conditions between different mixer scales and types (1). Finding the blending conditions that lead to the desired blend uniformity is not a trivial operation (2,3). For this reason fill volume, impeller rotational speeds and blender geometry are preferably kept constant as much as possible during such a technology transfer (4,5). In practice this is not always possible. Moreover, such an approach does not necessarily mean that the blend conditions are always similar (6). Still, obtaining the correct blend conditions is crucial because only this safeguards formation of a sufficiently uniform blend. A specific, but frequently occurring case is a blend that contains cohesive components that tend to form agglomerates. Removal of these agglomerates and prevention of the formation of new agglomerates is often the critical step in the assessment of the uniform blend (7,8).

Removal of agglomerates in a dry mixing system predominantly occurs via abrasion (8,9). The abrasion process is typically characterized by a high frequency of impacts by filler particles on the surfaces of the agglomerates. Mechanical properties of the agglomerates and product and process-related parameters like filler particle size and the rotational rate of the impeller were found to affect the size reduction rate of the agglomerates (8). It appeared to be possible to describe the abrasion process via definition of a Stokes number, the Stokes abrasion number (St_Abr). This number is the ratio of the kinetic energy density of the powder bed to the work of fracture of the agglomerate (9). Results so far showed that it is possible to predict effects of filler particle size and impeller rotation rate on abrasion rate with reasonable accuracy. The purpose of this study is to check the validity of the St_Abr approach and apply it to technology transfer situations where mixers of the same working mechanisms, but different production scales and different geometries are applied. Additionally, this paper identifies limits of applicability of the St_abr approach.

EXPERIMENTAL

Materials

The materials used were microcrystalline cellulose (Avicel PH-101, FMC, Philadelphia, USA) and α-lactose monohydrate (Pharmatose® 100M from DMV Fonterra Excipients, Goch, Germany, with a bulk density of 750 kg/m³).

Methods

Model Agglomerates (Brittle Calibrated Test Particles)

The model agglomerates or spherical brittle calibrated test particles (bCTPs) were prepared as described before by Willemsz et al. (8). The porosities of all bCTPs produced were measured from the diameters and the weights of the bCTPs. The true density of the MCC was determined using a pycnometer (AccuPyc 1330, Micromeritics, Norcross, U.S.A.) using nitrogen as test gas and was found to be 1,600 kg/m³. The mechanical properties of the bCTPs have been described in Willemsz et al. (9).

Blending Tests

The blending experiments reported in this study were performed using convective mixers with bowl volumes of 25 l (Fukae Powtec model FS-GS-25J, Japan, bottom-driven impeller) and 50 l (Glatt model VG50, Germany equipped with a top-driven impeller). The chopper was never installed. Table I lists the parameters that were varied. A test was started by adding selected test particles to a powder sample of lactose 100 M. This mixture was placed in the blender. After a given blending time, the blend was sieved over a 500-μm sieve to collect the test particles. The weights and dimensions of the bCTPs were determined as a function of blending time as mentioned in the previous communication (8).

Table I

Geometry of the High Shear Mixers and Experimental Set-Up for the blending Tests

The Froude number (Fr) was calculated as follows:

With N the impeller rotational rate, D bowl diameter, and g is the acceleration of gravity.

Powder Surface Velocimetry

To collect data for powder surface velocimetry analysis, a plexiglass lid was placed on the 25-L granulator. The powder surface velocimetry data from the 50-L high-shear mixer were collected by recording through the watch glass of the apparatus. The powder flow was recorded using a high-speed video camera (Casio-EX-F1, Casio computer co., LTD, Tokyo, Japan) operating at a speed of 600 frames per second. The data were analyzed according to Willemsz et al. (10). The camera was placed perpendicular to the bowl and such that about 50% of the total powder surface was visible for the 25-L high-shear mixer and 25% of the surface in the 50-L high-shear mixer.

Statistical Analyses

The outlier diagnostics, standard deviation (SD) and 95% confidence interval calculations described in this paper were performed using SAS V9.1 software (SAS institute Inc., North Carolina, USA).

RESULTS AND DISCUSSION

Abrasion Rate Constant (ξm) Measurements in the 25- and 50-L High Shear-Mixer

In this study, two vertical axis high-shear mixers of different geometries (as described in Table II) were used to assess how the abrasion rate constants (ξ_m) of brittle agglomerates scale with process variables in high-shear mixers. The bCTPs mass reduction (M_rel) over time was determined. It obeys apparent first order kinetics, with mass reduction rate constant as described by Willemsz et al. (8):

with M(t) as the mass after blending time t and M₀ as the initial mass.

Table II

Fit Parameters of Mechanical Properties Using equation M3

with X

Y or X

σ_c and X_0.55

σ_c or Y at ε

0.55 (The Data Indicate the Average

95% Confidence (more ...)

The purpose of the current experiments was to investigate the effects of abrasion rate constants of agglomerates with different porosities when process variables are varied. The results are depicted in Fig. 1.

Fig. 1

The effect of process settings on the mass based abrasion rate constant for the 25-L (black square) and 50-L (white square) high-shear mixers. Fr represents the Froude number, the degree of fill of the equipment, ξ_m the abrasion rate (more ...)

Figure 1 shows that the abrasion rate constants (ξ_m) of the agglomerates in the 50-L mixer are always lower than those obtained in the 25-L mixer when Froude numbers are identical. Abrasion rates increase with Froude number but decrease with increasing fill levels. These results are in line with findings discussed in previous papers (e.g. [8, 11–13]).

It is reasonable to assume that a reduction in fill volume implies that the contribution of the impeller to the total rate of agglomerate abrasion will increase. To visualize this effect additional tests have been performed where the powder just covers the impeller. This corresponds with a relative fill level of 8% in the 25-L mixer. Figure 2 shows the abrasion rates (ξ_m) of the bCTPs at different fill levels for two different Fr numbers. Figure 2 clearly shows a considerable additional effect of the impeller at low fill levels on the abrasion rates of the particles. Moreover, there seems to be a step change in behavior: abrasions rates at fill levels above 16% are more or less in line, a low fill level gives much higher abrasion rates.

Fig. 2

The mass based abrasion rate constants (ξ_m) of agglomerates for the 25-L high-shear mixer at a (v/v) of 8% (white circle), 16% (black circle), 27% (black square), and 37% (black diamond)

Powder Surface Velocity and Abrasion

The powder surface velocity has been measured as previously described (10). The powder velocities (v_p) were determined at the conditions described in Table I, and Fig. 3 depicts the results.

Fig. 3

Powder surface velocities (v_p) at different fill volumes and Froude numbers (black diamond, white diamond: Fr

0.12; black square, white square: Fr

0.46; black triangle, white triangle: Fr

1.04) for (more ...)

Figure 3 shows decreasing powder surface velocity at increasing relative fill volume. This was observed in both mixers. The figure also shows that the powder velocities measured in the 50-L mixer are significantly lower compared to those in the 25-L mixer scale at comparable Fr numbers. Filler particle velocity is an important parameter in relation to the rate of abrasion of the agglomerates (8,9). The effects were correlated using the Stokes abrasion number (St_Abr). Applying St_Abr numbers gives the possibility to assess how the abrasion rate constant (ξ_m) scales with process variables in different types/scale of high-shear mixers.

The Stokes abrasion number (St_Abr) concept has been discussed in more detail earlier in our previous paper (9). St_Abr compares the energy density during blending ( equation M4

) with the work of fracture of an agglomerate ( equation M5

With ρ_b bulk density of the filler, v_p powder surface velocity, Y elastic moduli, and σ_c fracture stress.

The mechanical properties Y (elastic modulus) and σ_c (fracture stress) of the bCTPs have been calculated as previously described (9) and are based on the porosity (ε) values of the agglomerates. Table II summarizes the fit parameters for elastic moduli and fracture stresses after performing a least square fit analysis assuming exponential relationships (14,15).

Figure 4 shows the relationship between the abrasion rates and the Stokes abrasion numbers in the blending experiments at different working conditions. Visually, three distinct relationships can be seen: with the largest group of tests at fill levels larger than 16% and two groups that both describe relationships when fill level is low.

Fig. 4

The relationship between the abrasion rate constants (ξ_m) and the Stokes abrasion number (St_Abr) of various bCTPs at defined conditions. The two solid lines indicate the abrasion data deviating from the proposed regression model (dashed line) (more ...)

Model diagnostic plots of the data-set indicate that both variables (ξ_m and St_Abr) should be log-transformed before analysis to fulfill the statistical requirements for normal distribution of the values to identify outliers. From these data sets several outlier diagnostics (studentized residuals, DFFITTS, leverage, and DFbetas) were used to identify outliers in the data-set of Fig. 4.

From the five curves depicted in Fig. 4, the analyses marked three observations as real outliers. These outliers correspond with bCTPs collected during tests using the 50-L high-shear mixer. This mixer is larger which implicates that larger amounts of filler had to be sieved to collect the model agglomerates. It is likely that this introduces additional errors. This was the rationale to remove the outliers from the data set. These data points are not shown in Fig. 4. After removing these three outliers, different regression models correlating ξ_m and St_Abr have been produced:

Table III lists the models produced.

Table III

Regression Models Between ξ_m and St_Abr at Various Relative Fill Volumes for the Curves Depicted in Fig. 4

The results demonstrate a relationship between abrasion rate of agglomerates and the value of St_Abr. The R² presented in Table III indicates the extent that the values of St_Abr explain abrasion at various process conditions. Here, R² approaching 100% indicates that abrasion is fully explained by the parameters that describe the Stokes abrasion (St_Abr) number of the system. The R² presented in Table III shows that the St_Abr number is a reasonable way to predict agglomerate abrasion while there is no clear difference between the fits of the results between the different blenders. This makes it possible to combine these results into one model. This model includes 78% of the variance when the fill level exceeds 16%. The regression analysis in our previous study (9) included 84% of the variance using a smaller data set. The data-set in this study also covers the abrasion data for the 50 L high-shear mixer scale.

It is clear that a low fill level leads to much faster abrasion of the test particles (Fig. 4 and Table III). There is apparently a transition where the impeller starts to dominate the abrasion. To study the impact of fill level, regression analysis has been performed separately for all fill levels. These results are depicted in Table IV.

Table IV

Regression Models Between ξ_m and St_Abr at Various Fill Degree for Two Different High-Shear Mixer Scales Depicted in Fig. 4

To study the effect of the impeller, the fill degree of the blender has been defined relative to the impeller height, the relative fill height (Δh_powder):

With h_0,powder height of the stationary powder and h_impeller the impeller height (Table I). Figure 5 shows the relationship between the fit constants in Table III and the relative fill height of the powder in the blenders.

Fig. 5

Relationships between the variables slope (α, upper figure) and intercept (β, lower figure) with the relative distance between impeller and powder height ( [increment]

h_powder) for the 25 L (black circle) and 50 L (letter x (more ...)

The slope (α_i) of the fits is almost constant and has a value of around 1. This implicates that the relationships between abrasion rate and St_Abr are almost linear relationships. As a consequence, the intercept β describes the slope of the (almost) linear relationships in Fig. 4. The value of β increases drastically when the relative fill height is low.

The intercept between the dotted and solid lines shown in Fig. 5 has been calculated and gives a transition point at a [increment]

h_powder value of 3. This result shows that agglomerate abrasion is predominantly determined by the powder bed movements when the [increment]

h_powder value is larger than 3. Obviously, the presence of enough powder is a prerequisite for the applicability of the Stokes number approach. When insufficient powder is present, the impeller starts to dominate the process. Logically, direct contact between impeller and the bCTP’s yields a deviating abrasive phenomenon than the shear forces occurring when there is plenty of powder present.

CONCLUSION

The abrasion of agglomerates during dry mixing at different fill volumes, impeller rotational speeds, and two different high-shear mixer scales and types has been investigated. This study reveals that the St_Abr number is able to predict the abrasion potential of agglomerates at variable process conditions. This includes high-shear mixers that are geometrically different.

The study reveals a transition point between agglomerate abrasion completely dominated by the powder blend and where a combination of impeller and powder blend effects play a significant role.

The St_Abr number concept described in this study demonstrates to be a useful tool to predict the abrasion of agglomerates at conditions during process variation exercises in dry mixing such as transferring the same mixture composition (filler) to a different type and scale of mixer.

ACKNOWLEDGMENT

This study was performed within the framework of Top Institute Pharma project number D6-203. The authors would like to thank Clinton Rambanapasi and Jan van Dinther for their experimental assistance in this study.

Open Access

This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.

REFERENCES

1. Brone D, Alexander A, Muzzio FJ. Qtitative characterization of mixing of dry powders in V-Blenders. AICHE J. 1998;44(2):271–278. doi: 10.1002/aic.690440206. [Cross Ref]

2. Muzzio FJ, Shinbrot T, Glasser BJ. Powder technology in the pharmaceutical industry: the need to catch up fast. Powder Technol. 2002;124:1–7. doi: 10.1016/S0032-5910(01)00482-X. [Cross Ref]

3. Muzzio FJ, Alexander AW. Scale up of powder-blending operations. Pharm. Technol. 2005;34–42.

4. Landin M, York P, Cliff MJ, Rowe RC, Wigmore AJ. Scale-up of pharmaceutical granulation in a fixed bowl mixer-granulator. Int J Pharm. 1996;133:127–131. doi: 10.1016/0378-5173(95)04427-2. [Cross Ref]

5. Landin M, York P, Cliff MJ, Rowe RC, Wigmore AJ. The effect of batch size on scale-up of a pharmaceutical granulation in a fixed bowl mixer-granulator. Int J Pharm. 1996;134:243–246. doi: 10.1016/0378-5173(96)04461-4. [Cross Ref]

6. Litster JD, Hapgood KP, Michaels JN, Sims A, Roberts M, Kameneni SK. Scale-up of mixer granulators for effective liquid distribution. Powder Technol. 2002;124(3):272–280. doi: 10.1016/S0032-5910(02)00023-2. [Cross Ref]

7. Kuwagi K, Horio M. A numerical study on agglomerate formation in a fluidized bed of fine cohesive particles. Chem Eng Sci. 2002;57(8):4737–4744.

8. Willemsz TA, Oostra W, Hooijmaijers R, De Vegt O, Morad N, Vromans H, et al. Blending of agglomerates into powders 1: quantification of abrasion rate. Int J Pharm. 2010;387(1–2):87–92. doi: 10.1016/j.ijpharm.2009.12.006. [PubMed] [Cross Ref]

9. Willemsz TA, Hooijmaijers RAA, Rubingh CM, Tran TN, Frijlink HW, Vromans H, et al. Kinetic energy density and agglomerate abrasion rate during blending of agglomerates into powders. Eur J Pharm Sci. 2012;45:211–215. doi: 10.1016/j.ejps.2011.11.012. [PubMed] [Cross Ref]

10. Willemsz TA, Tran TN, van der Hoeven M, Hooijmaijers RAA, Frijlink HW, Vromans H, et al. A statistical method for velocity detection in moving powder beds using image analysis. AICHE J. 2011. doi:10.1002/aic.12624.

11. Loveday BK, Naidoo D. Rock abrasion in autogneous milling. Miner Eng. 1997;10:603–612. doi: 10.1016/S0892-6875(97)00039-3. [Cross Ref]

12. Flisyuk EV, Palechkin AV, Burakova MA, Markova LM, Sakanyan EI, Kotovskii BK. Modeling the process of abrasion of granules based on herbal extracts in fluidized bed. Pharm Chem J. 2005;39(7):54–56. doi: 10.1007/s11094-005-0164-2. [Cross Ref]

13. Khanal M, Morrison R. Discrete element method study of abrasion. Miner Eng. 2008;21:751–760. doi: 10.1016/j.mineng.2008.06.008. [Cross Ref]

14. Ryshkewitch E. Compression strength of porous sintered alumina and zirconia. J Am Ceram Soc. 1953;36:65–68. doi: 10.1111/j.1151-2916.1953.tb12837.x. [Cross Ref]

15. Spriggs RM. Expression for effect of porosity on elastic modulus of polycrystalline refractory materials, particularly aluminum oxide. J Am Ceram Soc. 1961;44:628–629. doi: 10.1111/j.1151-2916.1961.tb11671.x. [Cross Ref]

Articles from AAPS PharmSciTech are provided here courtesy of
American Association of Pharmaceutical Scientists