Fluidized Bed Co-Melt Granulation: New Insights in the Influence of Process Variables and Validation of Regime Map Theory

: The understanding of granule growth mechanisms and the effects of formulation and operating conditions over product quality and process performance in ﬂuidized bed co-melt granulation is nowadays of great interest. In this sense, this work systematically studies the combined effects of binder content (W PEG ) and ﬂuidization air ﬂowrate (F A ) and temperature (T A ) on granules’ quality and process-related variables (product mass (M P ), elutriated ﬁnes (M f ), mass stuck on walls (M W )) by using a Box–Behnken-type design of experiments (DoE), as it is a statistical tool suggested by the Quality by Design (QbD) initiative. It was found that the granules’ size and powder ﬂowability are signiﬁcantly affected by W PEG (higher W PEG , higher granule size and better ﬂowability). Interestingly, T A is the process variable that signiﬁcantly affects M P , enhancing process performance at high temperature values. Regarding F A , it signiﬁcantly affects d 10 , promoting the formation of small particles due to breakage at high ﬂowrates and the presence of non-elutriated powder at low ﬂowrates. As a consequence, intermediate F A is the optimum for obtaining higher M P . Regarding regime map studies, most runs experienced a rapid growth regime, which is in accordance with the granules’ high pore saturation. This result agrees with the observed high increment in particle size and the morphology of the ﬁnal granules, allowing researchers to validate and extend existing previous maps.


Introduction
Granulation is an important process step applied in the manufacturing of numerous pharmaceutical products because it improves material handling and flowability, allows obtaining a desired particle size, enhances product appearance, controls particle moisture content, and it improves compressibility, content uniformity and, in some cases, the dissolution rate of the materials once they are formulated into a final solid dosage form (e.g., tablets) [1][2][3].Typically, three components are needed to produce granules: initial seeds or nuclei, mixing, and a binder.The seeds are always agitated to achieve a good binder distribution; therefore, granulators are often classified into mechanical (e.g., pans, drums, high shear granulators, twin-screw) or pneumatic (fluidized-bed granulators) agitated, according to the mixing principle.Fluidized-bed granulators offer some advantages with respect to other granulation systems, such as the simultaneous spraying, granulation, drying and/or cooling stages, and control, within certain limits, of the granule physical properties by the manipulation of some operating variables [4].
According to the binder nature, granulation processes are usually classified as wet, dry or melt.In melt granulation, particles are enlarged using binders that are in their solid states at room temperature, but they soften and flow at temperatures above their glass transition temperature or melting point [4].The absence of water or other solvents makes this approach quite valuable for moisture-sensitive ingredients.It completely avoids the use of solvents and the disadvantages associated with their recovery and final disposal, while it minimizes the energy cost related to solvent evaporation [5].Furthermore, the binders are added to the system either as: (1) powders that melt during the granulation process (co-melt or in situ granulation) or (2) atomized molten liquids (spray-on melt granulation) [3].
In this sense, fluidized bed melt granulation (FBMG) has received wide interest in the pharmaceutical research field, being the most common application the preparation of pharmaceutical granules for tableting [6].Although granulation processes involving the incorporation of the binder as droplets have received significantly more attention from a scientific point of view, both for wet and melt granulation, it is notable that co-melt processes are much more practical from an industrial point of view, since they avoid the spraying stage with hot-melt flows, leading to a simpler production process [7].To our knowledge, for co-melt granulation, the Walker et al. have been the first to extensively study this process.Walker et al. [8][9][10] and Zhai et al. [11,12] explored the effect of granulation time, binder/seed mass ratio, binder viscosity, and size of the seeds on the granule's final size (Particle Size Distribution, PSD) and growth mechanisms on systems composed of lactose or glass ballotini as seeds and different size fractions of PEG (of different grades according to desired variation in viscosity) and Poloxamer.Additionally, Mangwandi et al. [13] investigated the effects of fluidization air velocity, fluidization air temperature, and granulation time on granule mean size, PSD, granule homogeneity, and strength with lactose as seeds and PEG 1500 as solid binder.Ansari and Stepanek [14] used co-melt FBMG to study the growth mechanisms in hollow granules produced from mannitol granulated with PEG flakes of different initial size.Kidokoro et al. [15] prepared granules of lactose and PEG of different grades to study the binder/seed mass ratio and fluidization air temperature on granule size and flow properties.Passerini et al. [16] compared the effects of high-shear and FBMG techniques on granule size, flow properties, granulation yield (fraction between 100 and 1400 µm) and other variables related to tablet dissolution and stability for lactose granulated with PEG 6000 and ibuprofen.Prado et al. [17] studied the ability of different formulations of Gelucire to generate granules of lactose and ibuprofen, which were able to be tableted, and their effect on granules' PSD, flowability, tablet drug release, and granulation yield.Pauli Bruns et al. [18] used co-melt FBMG to prepare sustained release matrix granules of theophylline, using talc as a filler and wax as a binder, exploring the effect of formulation variables and different size fractions of wax on granules' PSD, sphericity, process yield, granules porosity, crushing resistance, and drug dissolution time.
As it can be seen, there is a great interest in understanding the granule growth mechanisms and effects of process conditions over granules' quality.However, still there is a lack of knowledge about the predictability and modeling of FBMG.In this sense, the US Food and Drug Administration (FDA) Quality by Design (QbD) initiative encourages the use of statistical tools for improving the development of high-quality pharmaceuticals.One of these tools is the design of experiments (DoE), which allows a better understanding of how formulation and process variables can influence product quality by defining the design space.From the experimental results of a response surface design, a polynomial model, describing the relation between a response and considered factors, is built.Afterwards, graphical and/or statistical analysis is carried out to determine the optimal operating conditions in order to obtain the desired formulation properties [7].In this sense, a reduced group of authors used DoE for elucidating FBMG under a more systematic approach.Masic et al. [19] studied the effect of binder/seed mass ratio, binder particle size and type of binder, fluidization air flowrate, and granulation time on granules' PSD and morphology, flow properties, and the tablet compression and drug release of a system of lactose and paracetamol granulated with PEG 2000 and Precirol ATO 5 as binders, using a 2 3 full factorial design.The same group [20,21] studied the granulation of paracetamol and lactose with Gelucire of different size fractions through a 2 3 full factorial design, Box-Behnken, and neural networks to evaluate the effect of binder/seed mass ratio, binder particle size, and granulation time on the obtained granules, tablets properties and granulation yield (fraction between 125 and 1250 µm).Mangwandi et al. [13] also performed studies with a DoE approach (a two-level factorial design).Kukec et al. [7] granulated Carvediol and lactose with a poloxamer to study the binder/seed mass ratio, fluidization air temperature, and endpoint temperature on granules' PSD, morphology, flow properties, and drug release through a Circumscribed Central Composite Design.More recently, Korteby et al. [6] related granules' PSD, aspect ratio, porosity, and granulation yield with the binder/seed mass ratio, viscosity and size of the binder, and outlet temperature through another empirical technique as Artificial Neural Networks (ANNs).
Also, to find an optimal operation window from a more mechanistic approach, the socalled regime maps have been developed in order to predict granulation and enable process design at all scales based on dimensionless numbers.In comparison to the DoE approach, once the map has been developed, a reduced number of experiments must be carried out [5].Each area of the regime map represents a granulation mechanism where the controlling parameters influence to different extents.However, a wide range of formulations, process conditions and granulators are necessary for the confirmation and validation of a regime map.In this sense, a turning point contribution is the one developed by Iveson and Litster to predict wet granulation regimes in drum granulators [22], which are extensively used to predict granulation also in high-shear granulators [23].However, regarding FBMG, scarcely few contributions can be found related to regime map development.Indeed, the same regime map was used by the group of Walker in some of the previously mentioned contributions [9,24] to predict granulation mechanisms in FBMG, even if the type or granulation unit was different (i.e., fluidized bed).
Overall, all these studies have announced a significant potential of FBMG for the preparation of granules with excellent properties for compression into tablets, relating operating conditions to granules properties through different approaches ranging from empirical but systematic, like DoE, to more predictive over different scales, such as regime map development.However, even if the studied process and formulation variables are of great importance in the FBMG process, no incidence of the combined effect of binder/seed mass ratio, fluidization air flowrate, and fluidization air temperature is assayed given that they are critical for binder melting and the distribution process and, therefore, granules' final properties.Furthermore, other critical process performance variables such as the quantity of fines collected on the FB granulator filters and/or the granulated particles that are lost due to stickiness on the equipment wall were not studied.Therefore, it is the aim of this work to systematically study the combined effect of (1) binder/seed mass ratio, (2) fluidization air flow rate, and (3) fluidization air temperature on granules' PSD, morphology, flow properties, and process-related variables (product collected in the granulation chamber, fines generation, loss due to wall stickiness) through a DoE approach by a Box-Behnken design.Furthermore, the Walker's regime map, which has been developed for a small number of experiments at similar conditions [5], is tested to identify the growth regimes of new experimental data obtained at conditions that incorporate variations in the fluidization air flowrate and temperature.Thus, it is expected to validate the regime map theory and expand its predictability to co-melt granulation in fluidized beds.

Materials
The meltable binder used in this study was polyethylene glycol 6000 (PEG 6000) purchased from Sistemas Analiticos S.A. (Buenos Aires, Argentina).It has a melting point around 60 • C [25].The binder was pre-conditioned by milling the commercial flakes using a laboratory hammer mill (Laboratory Mill 120, Perten, Stockholm, Sweden) equipped with a 2 mm hole-size sieve.After milling, PEG was sieved, and the fraction within 0.250-0.105mm was saved to be used in the granulation process.Lactose monohydrate (HMS (NZ)) was purchased from Chutrau S.A.C.I.F.(DFE PHARMA, Buenos Aires, Argentina) and used as provided.Several key specifications of lactose are shown in Table 1.The lactose particle size distribution (PSD) was measured by a light-scattering analyzer (LA 950 V2, Horiba, Kyoto, Japan) giving the diameters reported in Table 1 (i.e., d 10 , d 50 and d 90 ).The bulk density (ρ B ), tap density (ρ T ), Carr Index (CI), and angle of repose (AR) were determined as specified in Section 2.4.2.The solid density (ρ s ) of PEG 6000 and lactose was determined by nitrogen adsorption (NOVA 1200e, Quantachrome Instruments, Boyton Beach, FL, USA).This technique provides a density value that does not consider the volume of pores connected to the surface of the particle, thus giving a value very close to the solid density.The fraction of PEG 6000 used as a meltable binder has a ρ s = 1.163 g/mL.The experiments were performed in a fluidized bed granulator of ad hoc design (CON-ICET, Bahía Blanca, Argentina).A schematic diagram together with the chambers' geometrical parameters is shown in Figure 1.The equipment comprises a stainless-steel bottom conical vessel (1) and a cylindrical column (8) on top of it.The air distributor is a stainless-steel perforated plate (2).The fluidization air was supplied by a centrifugal blower (3) with a bypass system that can be actioned by means of a manual butterfly valve (6).This bypass allows better control of the air flow within the chamber.Before the bypass, the total fluidization air flowrate was measured by an orifice flowmeter (4) and preheated by an electrical heater (5) to maintain the air temperature at the desired level.The bypass system also presents an orifice flowmeter (7); therefore, the air flow entering the granulation chamber is inferred as the difference between the flow measured by ( 4) and ( 7) flowmeters.The elutriated fine particles were collected by a set of three filter bags located at the top of the fluidized bed freeboard (9).These filters were periodically blown back by air pulses to disengage the particulate matter.To monitor granulation temperature, the chamber was equipped with 3 thermocouples located at different distances from the air distributor (10).A Programmable Logic Control system (PLC) was used to register and control the process temperatures.

Experimental Procedure
For each run, the granulator chamber was loaded with 150 g of solid lactose monohydrate and PEG 6000 in different mass proportions (WPEG) according to the experience.The solid mixtures were previously heated to 50 °C in an oven to avoid a drop in temperature once charged to the unit.Furthermore, the granulator chamber was also previously

Experimental Procedure
For each run, the granulator chamber was loaded with 150 g of solid lactose monohydrate and PEG 6000 in different mass proportions (W PEG ) according to the experience.The solid mixtures were previously heated to 50 • C in an oven to avoid a drop in temperature once charged to the unit.Furthermore, the granulator chamber was also previously heated by means of the fluidization air at the desired process temperature (T A ). Once the solids were charged in the unit, fluidization air was again turned on at a set flowrate (F A ), which was determined by the aperture of the bypass valve (V A ).A controller manipulates the fluidization air inlet temperature (T A ) to maintain the temperature at the specified value.Once the thermal regime of the bed was established, the granulation was carried out for 20 min.After that, the fluidization air temperature was turned down, allowing the granules to cool.Once the granules' temperature reached 30 • C, the fluidization air was turned off, and the run stopped.Afterwards, the material deposited in the different parts of the granulation unit-granulation chamber (M P ), walls (M W ), and filters (M f )-were collected and separately saved in hermetically sealed plastic bags for characterization.The granules collected in the granulation chamber (M P ) were considered as each run product.

Experimental Design
To evaluate the influence of the main operating parameters on the properties of the obtained granules and process conditions, we used a Box-Behnken incomplete factorial design consisting of three factors with three central points: % aperture of bypass valve (X 1 ), binder mass content (X 2 ), and fluidization air temperature (X 3 ).These were selected according to similar works in which these variables had significant effects, and the settings were determined based on the equipment practical considerations and the authors' initial experimentation and prior experience.Since Box-Behnken designs are specially made to require only three levels, coded as −1, 0, and +1, the factor levels chosen in terms of actual values were X 1 = 0, 50, 100%; X 2 = 10, 15, 20 wt%; and X 3 = 80, 90, 100 • C. In the case of X 1 , the variable that is manipulated at the end is the fluidization air flowrate that, as explained in Section 2.2.1, is determined by the difference between the total fluidization air flowrate and the one derived to the bypass according to the valve % aperture setting; the higher the aperture, the lower the fluidization air flowrate for the granulation chamber.
As responses, several variables were studied.Regarding the produced granules, the following variables were analyzed: size distribution (by estimation of representative granule diameters d 10 , d 50 , d 90 ), angle of repose (AR), mass of granular product collected in the chamber (M P ), mass of fines (M f ), and mass on the walls (M W ).
A matrix design was performed using the Design Expert 10.0 software (Trial version, Stat-Ease, Inc., Minneapolis, MN, USA), which consisted of 15 runs with the combination of factors and levels shown in Table 2.The software was used to find suitable fitting models in which the coefficients were calculated based on the experimental responses.Multiple quadratic regression models were obtained from the full or reduced following equation: where ŷ is the response studied, β 0 is the intercept, β i is the main effect coefficient for the ith factor, β ii is the quadratic effect coefficient for the ith factor, β ij is the model coefficient for the interaction between factors ith and jth, and x i and x j are each of the considered factors.ANOVA was used to determine the statistical significance of the model factors.
Model terms with p-values less than 0.05 were considered significant, and non-significant terms were deleted from the equation except when they increased the quality of the fit.The adequacy of the developed regression models was evaluated by calculation of the modeling statistical coefficient R-squared (R 2 ) and the adequate precision (AP) which is a measure of the signal-to-noise ratio.An AP value greater than 4 is desirable to use the model to navigate the design space.Last but not least, the lack of fit value (LOF) should not be significant (i.e., LOF p-values should be higher than 0.05 to indicate that the model fits the data well) [26].The PSD of the product granules was evaluated by sieving analysis using a vibratory sieve shaker (ZONYTEST EJR 2000, Zonytest, Argentina) and a series of standard sieves in a range from 3.36 to 0.053 mm following a 2 1/2 progression.Sieve analysis was performed by duplicate, with each sample weighing 30 g and 15 min of sieving time.d 10 , d 50 , and d 90 were determined from the mass cumulative distribution.
To further characterize the granules, the morphology of the external surface of some selected granules (G1, G3, and G5) was assessed in an EVO 40-XVP, LEO Scanning Electron Microscope (SEM).The samples were previously metalized with gold in a PELCO 91,000 sputter coater.

Bulk and Tapped Densities. Carr Index. Angle of Repose
Bulk densities (ρ B ) of granules and lactose powder were measured by pouring 1.5 g in a 10 mL graduated cylinder and calculated as the ratio between the weight of the sample and the occupied volume.The tapped density (ρ B ) of the same samples was measured by tapping the cylinder 100 times and registering the final volume.
For the same purpose, the angle of repose (AR) was determined by pouring 1.5 g of granules or lactose powder through a funnel located at a fixed height on a flat horizontal surface of fixed radius (r) and measuring the height (h) of the conical pile formed.The tangent of the angle of repose is given by the h/r ratio [27].
Powder compressibility was evaluated using Carr's compressibility index (CI), which was estimated by Equation (2) [28]: All samples were analyzed at least in triplicate.

Granule Porosity
To determine the granules' porosity (ε), the envelope density (ρ E ) was tested by using a pycnometer and ethanol as the displacement fluid, considering both lactose and PEG 6000 insolubility or very low solubility in this substance.Afterwards, the samples' solid density was estimated from the lactose, and the PEG 6000 solid density was determined by nitrogen adsorption according to the corresponding mass fraction for each run.Then, the porosity of granules was calculated by Equation (3).

Granulation Performance
The granulation performance was evaluated in terms of the mass of the granular product collected in the chamber (M P ).The process losses were estimated by weighting the material collected in the filters (M f ) and on the chamber walls (M w ).

Regime Map
To elaborate Walker's regime map, it is necessary to estimate for the system the Stokes deformation number (St de f ) and the maximum pore saturation (s max ) given by the dimensionless groups presented in the following equations: where ρ g is the granules' density, U the inter-granule collision velocity, Y g is the granule dynamic yield stress, w is the mass ratio of liquid to solid, ρ SS is the solid seeds density, ρ L is the liquid binder density, and ε is the granules' porosity.For s max , w was determined from the binder/solid ratio for each run, ρ SS and ρ L were the ones measured by nitrogen adsorption, and ε represents the granules' porosity estimated as explained in Section 2.4.3.Regarding the Stokes deformation number, ρ g was estimated as ρ E (Section 2.4.3), and the granule collision velocity for simplicity was approximated as the superficial gas velocity, while Y g was more challenging to determine.According to Walker et al. [9], the granule dynamic yield stress should be replaced by the characteristic flow stress of the particle-binder slurry τ g , whose results, of course, are quite complex to determine.The same authors proposed, in another contribution, a novel technique to measure τ g , based on a compression test at temperatures similar to the ones experimented by the granules during the granulation process [10].Based on this, an ad hoc procedure was implemented in this contribution that consisted of measuring granules' compression breakage strength using a Texture Analyzer (TA Plus Lloyd Instruments, England) equipped with a 50 N cylindrical load cell.The deformation behavior and fracture strength of individual granules corresponding to the selected representative runs (G1, G3, G5 and G8) were measured by compression between a stainless-steel base and a cylindrical probe of 5 mm diameter.The probe descended at a constant velocity of 1 cm/min.The test was carried out at least for 10 granules from each run, which were manually selected from the fraction retained between 1.19 and 1.68 mm sieves.The granules were heated on a hotplate, for at least 5 min prior to compression testing, to a temperature equal to the average chamber temperature (T TEST ) of the corresponding test.In this way, the first appreciable fall of the force-deformation curve was considered as the granule breakage point.The particle-binder slurry flow stress, τ g , was taken to be proportional to the breakage force divided by the cross-sectional area of the granules.Image analysis of the granules selected for the compression test was carried out to determine each particle's cross-sectional area.  2 presents the combination of factors and levels for the 15 runs together with the corresponding experimental results obtained for responses d 10 , d 50 , d 90 , and AR.Also, the responses M P , M f , and M w are shown, which will be discussed in Section 3.1.2.Furthermore, Table 3 presents the analysis of variance (ANOVA) for the same responses, which were analyzed by the application of response surface methodology.The p-values for each model term of the fitted responses, as well as the corresponding statistical coefficients, are presented.The corresponding regression equations are presented in the Appendix A. As shown, the binder mass content results showed significant factors for all the analyzed responses.The fluidization air temperature significantly affects d 10 and d 50 , while the air flowrate (% aperture of bypass valve) only presents a significant term for d 10 .It is also important to note that the statistical coefficients values were reasonable for all responses, with R 2 greater than 0.75 for all cases, except AR (R 2 = 0.62).Furthermore, the AP is higher than 4 in all cases, indicating a good signal-to-noise ratio, and the LOF is not significant.As it can be seen in Figure 2a, d 10 increases as W PEG increases.The same trend regarding binder mass content was observed by several authors that studied FBMG with or without DoE [7][8][9]19].Regarding air temperature, it significantly affects the fines in the collected product; as it can be seen in Figure 2a, the higher the temperature, the higher the d 10 value, indicating that a decrease in binder viscosity favors the particle growth of small granules.A similar effect was observed by Walker et al. [8], who found bigger granules for lower viscosity binders.Interestingly, the quadratic term of fluidization air flowrate inversely significantly affects d 10 ; this would mean that as the fluidization air flowrate increases (0% aperture of bypass valve), particles with a smaller size are retained in the granulator product, which can be observed in Figure 2b (for W PEG and T A at intermediate levels).However, d 10 also tends to decrease for the lowest fluidization air flowrate (100% aperture of bypass valve).This indicates that at low air velocities, more fines are retained in the product, which is reasonable since they cannot be elutriated.As the air velocity increases, d 10 starts to increase, having its maximum at intermediate air flowrates.High air velocities tend to decrease d 10 again, which is probably due to the breakage or damage suffered by the formed agglomerates during the process, which experiences more intense impacts.Fluidization air velocity has been previously studied only by Masic et al. [19] and Mangwandi et al. [13].Masic et al. [19] found that fluidization air velocity alone (i.e., not in combination with air temperature) is not a significant factor for their system.On the other hand, Mangwandi et al. [13] found that air velocity strongly affects drug content and granule strength when studied in combination with air temperature.Regarding d50, Figure 2c shows an increasing trend with WPEG.As expected, the more binder in the system, the more the particles grow.Air temperature is the other significant factor for this variable; however, the tendency is the opposite to the one observed for d10.In fact, according to the corresponding model coefficients (Table 3), d10 presents a coefficient for X3 = 0.0637 (positive), while for d50, it is X3 = −0.2412(negative), indicating a direct relation for the first and an inverse one for the second.As shown in Figure 2c, bigger granules are obtained for lower air temperature.Indeed, temperature affects binder viscosity, which, in turn, influences the balance between the inertial forces (generated by the collisions between particles) and the cohesive forces (viscous forces of the formed liquid bonds between particles).The observed behavior suggests that at low temperatures, the binder's higher viscosity favors the dominance of viscous forces that tend to keep the particles stuck together, independently of its size.At higher temperatures, the size of the interme- Regarding d 50 , Figure 2c shows an increasing trend with W PEG .As expected, the more binder in the system, the more the particles grow.Air temperature is the other significant factor for this variable; however, the tendency is the opposite to the one observed for d 10 .In fact, according to the corresponding model coefficients (Table 3), d 10 presents a coefficient for X 3 = 0.0637 (positive), while for d 50 , it is X 3 = −0.2412(negative), indicating a direct relation for the first and an inverse one for the second.As shown in Figure 2c, bigger granules are obtained for lower air temperature.Indeed, temperature affects binder viscosity, which, in turn, influences the balance between the inertial forces (generated by the collisions between particles) and the cohesive forces (viscous forces of the formed liquid bonds between particles).The observed behavior suggests that at low temperatures, the binder's higher viscosity favors the dominance of viscous forces that tend to keep the particles stuck together, independently of its size.At higher temperatures, the size of the intermediate particles generates inertial forces higher than the viscous ones; therefore, the bonds are not resistant enough, and the granules tend to separate, decreasing d 50 .For small particles, viscous forces are always higher than inertial ones, independently of binder viscosity variation with temperature, which favors that d 10 always increases with temperature.This effect was also proposed by Mangwandi et al. [13], who stated that it is logical to expect a high rate of granule breakage at higher temperatures, since decreasing the binder viscosity reduces the strength.
For d 90 , the only variable that significantly affects the response is binder mass content.As it can be seen in Figure 2d, higher binder contents favor bigger granules (i.e., higher d 90 ), although d 90 tends to increase, in a lesser extent, for lower binder mass contents, as its quadratic term is significant as well.
The flow properties of powders can be described according to the obtained repose angle values.For AR, flow between 25 and 30 • is excellent, among 31-35 • the flow is good, within the range 36-40 • the flow is fair, and for values higher than 41 • , the powder has bad flow properties [29].As shown in Table 2, AR values are smaller than 30 • for all runs, which means excellent flow, except for one experience where this value falls within good flow.Regarding the statistical analysis for this response, Table 3 shows that AR is significantly affected only by binder mass content.As W PEG increases, the angle of repose decreases, becoming a better flowing material.In fact, this model presents a coefficient for X 2 = −3.660,indicating an inverse response with respect to the factor; that is, the higher the binder mass content, the lower the angle of repose.Indeed, it is well known that bigger particles favor powders flow.As previously seen, a higher binder mass content increases the granules size (higher d 10 , d 50 , and d 90 ), giving as a result smaller AR values.To complement this result, Figure 3 presents the measured CI for all runs.Carr s compressibility scale is more specific than AR, indicating excellent flow for values of 10%, good flowability between 11% and 15%, fair flow between 16% and 20%, acceptable flow between 21% and 25%, and poor flow between 26% and 31% [29].As it can be seen, the granules produced in 11 out of the 15 runs have excellent flow, while 1 run gives a product with good flow, 2 runs give a product with fair flow, and only 1 run gives a product with poor flow.Again, it is possible to observe that a higher binder mass content improves flowability (lower CI), while no specific trend with fluidization air flowrate (% aperture of bypass valve) or air temperature is seen.Both Masic et al. [19] and Kukec et al. [7] reported similar findings, demonstrating that a higher binder content led to better flowability of the granules due to bigger particles and narrower particle size distributions obtained under those conditions.
Powders 2023, 3, FOR PEER REVIEW 11 particles, viscous forces are always higher than inertial ones, independently of binder viscosity variation with temperature, which favors that d10 always increases with temperature.This effect was also proposed by Mangwandi et al. [13], who stated that it is logical to expect a high rate of granule breakage at higher temperatures, since decreasing the binder viscosity reduces the strength.For d90, the only variable that significantly affects the response is binder mass content.As it can be seen in Figure 2d, higher binder contents favor bigger granules (i.e., higher d90), although d90 tends to increase, in a lesser extent, for lower binder mass contents, as its quadratic term is significant as well.
The flow properties of powders can be described according to the obtained repose angle values.For AR, flow between 25 and 30° is excellent, among 31-35° the flow is good, within the range 36-40° the flow is fair, and for values higher than 41°, the powder has bad flow properties [29].As shown in Table 2, AR values are smaller than 30° for all runs, which means excellent flow, except for one experience where this value falls within good flow.Regarding the statistical analysis for this response, Table 3 shows that AR is significantly affected only by binder mass content.As WPEG increases, the angle of repose decreases, becoming a better flowing material.In fact, this model presents a coefficient for X2 = −3.660,indicating an inverse response with respect to the factor; that is, the higher the binder mass content, the lower the angle of repose.Indeed, it is well known that bigger particles favor powders flow.As previously seen, a higher binder mass content increases the granules size (higher d10, d50, and d90), giving as a result smaller AR values.To complement this result, Figure 3 presents the measured CI for all runs.Carr s compressibility scale is more specific than AR, indicating excellent flow for values of 10%, good flowability between 11% and 15%, fair flow between 16% and 20%, acceptable flow between 21% and 25%, and poor flow between 26% and 31% [29].As it can be seen, the granules produced in 11 out of the 15 runs have excellent flow, while 1 run gives a product with good flow, 2 runs give a product with fair flow, and only 1 run gives a product with poor flow.Again, it is possible to observe that a higher binder mass content improves flowability (lower CI), while no specific trend with fluidization air flowrate (% aperture of bypass valve) or air temperature is seen.Both Masic et al. [19] and Kukec et al. [7] reported similar findings, demonstrating that a higher binder content led to better flowability of the granules due to bigger particles and narrower particle size distributions obtained under those conditions.Powders flow is also enhanced by particles that present high sphericity and smooth surfaces.To illustrate the difference in granules surface as the amount of binder is modified,

CI (%)
W PEG (wt%) Powders flow is also enhanced by particles that present high sphericity and smooth surfaces.To illustrate the difference in granules surface as the amount of binder is modified, SEM images of two selected runs, G1 and G3, are shown in Figure 4.Both runs present the same settings for aperture of bypass valve and fluidization air temperature, but they also have minimum and maximum values of binder mass content.Figure 4a,b correspond to two magnifications of a representative granule of run G1, which was obtained with 10 wt% of PEG.The surface in this case is irregular, which is probably due to granules of lactose and PEG that come into contact and stay as agglomerates.Indeed, G1 presents a CI of 17.4, indicating acceptable flow.In Figure 4c,d, a surface formed by smaller particles, and overall, smoother is found.This granule corresponds to a run with a high amount of binder (20 wt%) and a product that presents a CI of 7.69, which is well within the range of excellent flow.As more binder is present, more liquid remains available in the granules surface, allowing small ungranulated lactose particles to stick to it and giving a layering growth with particles engulfment.Aleksic et al. [21] and Masic et al.
[20] also found that when large binder particles (with respect to seeds, as in the present study) were used in combination with a higher binder content, larger granules with smooth surfaces (i.e., surfaces with more saturated binder) were formed.This led to excellent flow properties in the resulting granules.Similarly, Korteby et al. [6] observed that when a low binder content was used, fragile and low-density granules with a low degree of sphericity were formed.
SEM images of two selected runs, G1 and G3, are shown in Figure 4.Both runs present the same settings for aperture of bypass valve and fluidization air temperature, but they also have minimum and maximum values of binder mass content.Figure 4a,b correspond to two magnifications of a representative granule of run G1, which was obtained with 10 wt% of PEG.The surface in this case is irregular, which is probably due to granules of lactose and PEG that come into contact and stay as agglomerates.Indeed, G1 presents a CI of 17.4, indicating acceptable flow.In Figure 4c,d, a surface formed by smaller particles, and overall, smoother is found.This granule corresponds to a run with a high amount of binder (20 wt%) and a product that presents a CI of 7.69, which is well within the range of excellent flow.As more binder is present, more liquid remains available in the granules surface, allowing small ungranulated lactose particles to stick to it and giving a layering growth with particles engulfment.Aleksic et al. [21] and Masic et al.
[20] also found that when large binder particles (with respect to seeds, as in the present study) were used in combination with a higher binder content, larger granules with smooth surfaces (i.e., surfaces with more saturated binder) were formed.This led to excellent flow properties in the resulting granules.Similarly, Korteby et al. [6] observed that when a low binder content was used, fragile and low-density granules with a low degree of sphericity were formed.
These results indicate that granulation, with any combination of factors, greatly improves the initial lactose powders' flowability, either by increasing particle size or giving granules with smooth surfaces.Indeed, lactose has an AR of 35.4° and CI of 35.71%, corresponding to fair and poor flow, respectively.

Granulation Performance
Table 4 shows the ANOVA for MP, Mf, and Mw.The corresponding regression equations are presented in the Appendix section.As shown, the fluidization air temperature was shown to be a significant factor for all analyzed responses, while the binder mass content significantly affected the mass collected on the walls as a linear term, and in the filters, in the term interacting with air temperature.The statistical coefficients values were reasonable for all responses, with R 2 greater than 0.71 for all cases, an AP higher than 4, and a not significant LOF.These results indicate that granulation, with any combination of factors, greatly improves the initial lactose powders' flowability, either by increasing particle size or giving granules with smooth surfaces.Indeed, lactose has an AR of 35.4 • and CI of 35.71%, corresponding to fair and poor flow, respectively.

Granulation Performance
Table 4 shows the ANOVA for M P , M f , and M w .The corresponding regression equations are presented in the Appendix A. As shown, the fluidization air temperature was shown to be a significant factor for all analyzed responses, while the binder mass content significantly affected the mass collected on the walls as a linear term, and in the filters, in the term interacting with air temperature.The statistical coefficients values were reasonable for all responses, with R 2 greater than 0.71 for all cases, an AP higher than 4, and a not significant LOF.As shown in Figure 5a, the mass collected in the chamber increases with air temperature, which could mean that settings with low binder viscosity (i.e., high operation temperature) favor the formation of granules that remain in the chamber.This is expected since having a lower viscosity binder favors its spreading on the particles surface, which in turn allows granules' growth.
For the mass collected in the filters, again, it significantly depends on air temperature as well as its interaction with W PEG .As shown in the response surface of Figure 5d, M f increases with W PEG for low temperatures, while the opposite trend, for much lower values of M f , is observed at high temperatures, indicating an interaction between both factors.This is consistent with the findings previously observed: that at low T A , less product is collected in the granulator chamber.Low temperatures favor a high viscosity of the binder, which hinders its spreading on the initial seeds' surface, leaving more un-granulated particles able to be elutriated.In fact, it is possible that a portion of the binder droplets remain as such and solidify without even reaching a seed particle.The more binder there is in the system, the higher the probability of having this droplet solidification, which can explain the increment of fines at low temperatures and high W PEG .
Regarding the mass collected on the walls (Figure 5c), it is significantly affected by both binder mass content and fluidization air temperature.Interestingly, low W PEG and low T A favor materials' stickiness to the walls.Again, low temperatures correspond to the most unfavorable setting for M P , since both M W and M f are maximum, confirming mass balance closure.This can be attributed again to the effect of binder viscosity that generates stronger bonds at high viscosities (low temperatures); therefore, if granules coated with liquid binder reach the walls, the particle-wall liquid bridges are strong enough to remain stuck.
Regarding fluidization air velocity, even if it is not a statistically significant factor for any of the three analyzed responses, it is worth studying its effect, as it helps with understanding the particle's growth mechanisms.Figure 5d-f show the variation of M P , M f , and M W with the aperture of bypass valve for all the other factors on an intermediate level.As it can be seen, when the bypass valve is 0% open, which means high air velocity, M P is minimum, indicating that either ungranulated particles are able to be elutriated or that the particles that reach the walls remain adhered to them.Indeed, high velocities favor the entrainment of particles, which corresponds to the high mass on filters observed for 0% bypass valve aperture.At the same time, high air flowrates tend to generate higher temperatures in the granulation chamber, since there is more hot mass circulating, which favors binder melting.Furthermore, particle-wall collisions are also more frequent at high air velocities.If more particles covered with melted binder reach the wall, there is more probability for them to stay adhered to the chamber, which reflects the maximum of M P at this same setting of bypass valve aperture.On the contrary, lower fluidization air flowrates (100% bypass valve aperture) lead to lower temperatures within the chamber that can hinder granulation, leaving more ungranulated powder able to be elutriated.This is probably the reason why M f increases again for 100% bypass valve aperture, while M P is minimum at this setting.For the mass collected in the filters, again, it significantly depends on air temperature as well as its interaction with WPEG.As shown in the response surface of Figure 5d, Mf increases with WPEG for low temperatures, while the opposite trend, for much lower values of Mf, is observed at high temperatures, indicating an interaction between both factors.This is consistent with the findings previously observed: that at low TA, less product is collected in the granulator chamber.Low temperatures favor a high viscosity of the Overall, the combination of variables that most favor granulation with respect to process yield and less material losses is high temperature with high binder mass content at intermediate or low fluidization air velocities.This has not been previously reported and is an important result when defining the overall process performance.

Regime Map
As described in the Introduction section, a semi-empirical approach to understanding the granulation process has been the well-known regime map study.This approach requires deeply understanding the effect of influencing parameters, such as formulation properties, process parameters and equipment parameters to determine convenient dimensionless groups that allow localizing a particular operation regime.Once the limits between the different growth regimes have been established, it is possible to define a design space in order to keep the dimensionless parameters constant or within the same operating window across different scales [5].
Even if there exist different approaches for regimes' maps development, the one proposed by Iveson and Litster to predict wet granulation regimes in drum and high-shear granulators [22,23] is, to our knowledge, the only one that has been extended to predict granules' growth in FBMG, which was completed by Walker s group [9,24].In Iveson and Litster's map, the dominant granulation mechanism is characterized by calculating the St de f (Equation ( 4)), which represents the ratio of externally applied kinetic energy and the energy required for deformation, and the s max (Equation ( 5)), that reflects the system liquid content.
In general, these growth regime studies have proposed for FBMG that, initially, granules are mainly obtained by two different nucleation mechanisms: distribution (molten binder is distributed on the surface of the solid particles) and immersion mechanisms (solid particles are immersed into the surface of the molten binder droplets).In some previous studies on FBMG [30], it was found that the prevalence of one or the other mechanism is mainly governed by the ratio between the size of molten binder particles (droplets once they melt) and the size of solid particles.Nucleation by immersion and subsequent layering with small solid particles is promoted by larger molten binder droplets, while nucleation by distribution and further growth by coalescence or agglomeration occurs when the molten binder droplets are not greater than the solid particles.
Once the nucleation phase is completed, two main types of growth can be distinguished: steady growth and induction-type growth.Steady growth occurs when granules deform and coalesce easily (high St de f ); the collision facilitates liquid binder to squeeze out to the granules' surface, generating new areas of contact where new bonds can be formed.If these bonds are resistant enough, larger granules are formed.In this regime, the growth rate increases when more binder is present in the system (high s max ).In contrast, induction behavior occurs when the granules possess low deformability (low St de f ) and is characterized by the deficient formation of contact area.Clumps of granules can break easily as the bonds are weak.The duration of the induction period decreases with increasing liquid level [5].At high liquid contents (s max ~100%), rapid growth occurs since sufficient liquid arises in the surface that facilitates strong bond formations.Further liquid addition (s max > 100%) causes a transition to an over-wet mass [22].If the system is exposed to high impact velocities (very high St de f ), no permanent granules are generated, as they are too weak to resist any separating forces, and a crumb behavior is exhibited [5].
With the aim of extending the FBMG regime map developed by Walker and coworkers, the dimensionless numbers s max and St de f presented in Section 2.6 need to be calculated.To this end, and as it was mentioned in the Materials and Methods section, it was necessary to measure the characteristic flow stress of the particle-binder slurry (τ g ) by the compression tests at granulation temperatures described in Section 2.6.Figure 6 shows a typical stress-deformation profile obtained in these tests.The arrow indicates the value taken as the stress of the particle-binder slurry, considering that it is estimated as the force value registered during the compression tests divided by the projected area of the granule obtained by image analysis prior to the assay.Table 5 summarizes the average τ g values and deviations for all the tests.As it can be seen, the flow stress presents variations between the different samples, although no direct relation with operating conditions can be established.Consequently, it was decided to take for all runs an average value for τ g of 0.045 MPa.Using this value, St de f values for all runs were calculated.Table 6 presents the measured values for each run envelope density ( ), particle density ( ), granules' porosity () and the corresponding dimensionless numbers  and . As it can be seen, porosities range between 0.018 and 0.216; 0.107 is the average value for all runs.In general, granules' porosity is rarely reported in the literature, which is mainly because it is a variable quite difficult to measure considering that it must consider the pores that are formed between the particles and binder during the granulation process.In this case, lower porosities than those reported by Walker et al. [9] are found, although these authors used an average porosity of 0.17 for  calculation and stated that it was independent of the mass of binder.In our case, porosity tends to decrease with temperature, while as observed by Walker, no clear trend with WPEG and fluidization air velocity is found.Indeed, high temperatures favor binder spreading due to lower viscosities, filling all the internal pores within the granules.Considering that most of the obtained  are higher than 1, it is possible to infer that the process involves the formation of granules where binder can fill most of the internal voids.Regarding  , the values are within 0.0187 and 0.0726, widening the spectra of deformation values tested by Walker et al. [9], taking into account that in the present contribution, variations in fluidization air flowrate are considered.Table 6 presents the measured values for each run envelope density (ρ E ), particle density (ρ s ), granules' porosity (ε) and the corresponding dimensionless numbers St de f and s max .As it can be seen, porosities range between 0.018 and 0.216; 0.107 is the average value for all runs.In general, granules' porosity is rarely reported in the literature, which is mainly because it is a variable quite difficult to measure considering that it must consider the pores that are formed between the particles and binder during the granulation process.In this case, lower porosities than those reported by Walker et al. [9] are found, although these authors used an average porosity of 0.17 for s max calculation and stated that it was independent of the mass of binder.In our case, porosity tends to decrease with temperature, while as observed by Walker, no clear trend with W PEG and fluidization air velocity is found.Indeed, high temperatures favor binder spreading due to lower viscosities, filling all the internal pores within the granules.Considering that most of the obtained s max are higher than 1, it is possible to infer that the process involves the formation of granules where binder can fill most of the internal voids.Regarding St de f , the values are within 0.0187 and 0.0726, widening the spectra of deformation values tested by Walker et al. [9], taking into account that in the present contribution, variations in fluidization air flowrate are considered.Overall, Figure 7 depicts the obtained dimensionless numbers on the regime map proposed by Iveson and Litster [22,23] and later extended to predict granules' growth in FBMG by Walker s group [9,24].G7 was excluded from the map to better visualize the runs' spread.As it can be seen, most of the runs fall within the rapid growth area.That is, final granules are formed by the agglomeration of multiple smaller granules, giving a broad PSD that is substantially bigger than the initial seeds.Indeed, this is indicated by the values of mean diameters reported in Table 2, which are one order of magnitude bigger than that corresponding to lactose (d 50 = 0.038 mm).As an example, the SEM images of selected representative runs (G1 and G5) are shown in Figure 8.As it can be seen, as the pore saturation increases (i.e., greater than 1), granules are formed by the agglomeration of smaller granules (G1).Furthermore, when the binder content is high and it is able to wet the external surface, small lactose particles adhere to the granules, promoting a growth similar to coating.Similar results were observed by Korteby et al. [6], who observed that a large binder particle size, such as the one used in this study, together with high binder viscosity (G5 presents low T A ) promotes the successful attachment of filler particles to the surface of the binder.The runs that present higher fluidization air flowrates (red symbols) are within high deformation numbers.As the collision velocity between granules increases, it is expected to favor granules' deformation.The opposite occurs for the lower fluidization air flowrates (gray symbols).Still, fluidization air velocities are low enough to avoid a crumb type of growth.In fact, all granules seem well consolidated, suggesting that breakage is not the dominating process.Regarding T A , it is possible to observe that at higher fluidization air temperature, there is higher pore saturation (i.e., solid symbols −80 • C are toward the left of the map while stripped −90 • C and empty −100 • C are more to the right), meaning that the melted binder can fill and saturate the granules' pores.As previously explained, this can be related to the decrease in viscosity at higher temperatures.Finally, for W PEG , even if squares (20 wt%), circles (15 wt%), and triangles (10 wt%) are found across the regime map, there is a greater concentration of intermediate and high binder mass content toward the right, indicating that more PEG in the system favors the formation of granules with saturated pores.
Regarding growth regimes, most of the runs are categorized as rapid growth, there are two runs with steady growth (G1 and G6) and one in the nucleation area (G2).Particularly, these three runs present the worst flowability (Table 2 and Figure 3), indicating that the combination of operating variables hinders the production of big granules with smooth surfaces.Regarding growth regimes, most of the runs are categorized as rapid growth, there are two runs with steady growth (G1 and G6) and one in the nucleation area (G2).Particularly, these three runs present the worst flowability (Table 2 and Figure 3), indicating that the combination of operating variables hinders the production of big granules with smooth surfaces.

Conclusions
There is a great interest in understanding the granule growth mechanisms and the effects of process conditions over granules' quality and the operation efficiency for the comelt fluidized bed granulation.In this sense, this work studies, by a Box-Behnken approach, the combined effect of binder mass content and fluidization air flow rate and temperature on granules' quality (particle size distribution, morphology, flow properties) and process-related variables (mass of granules collected as product in the granulator chamber, fines elutriation, and loss due to wall stickiness).It was found that binder content has a direct effect on the granules' size, and its increment also improves powder flowability.Furthermore, more granules are collected as product within the granulator for higher binder content.Regarding fluidization air temperature, it affects binder viscosity, impacting on the particle-particle liquid bridges' strength.Interestingly, for medium and coarsesize granules, higher temperatures tend to decrease their size.Furthermore, high air   Regarding growth regimes, most of the runs are categorized as rapid growth, there are two runs with steady growth (G1 and G6) and one in the nucleation area (G2).Particularly, these three runs present the worst flowability (Table 2 and Figure 3), indicating that the combination of operating variables hinders the production of big granules with smooth surfaces.

Conclusions
There is a great interest in understanding the granule growth mechanisms and the effects of process conditions over granules' quality and the operation efficiency for the comelt fluidized bed granulation.In this sense, this work studies, by a Box-Behnken approach, the combined effect of binder mass content and fluidization air flow rate and temperature on granules' quality (particle size distribution, morphology, flow properties) and process-related variables (mass of granules collected as product in the granulator chamber, fines elutriation, and loss due to wall stickiness).It was found that binder content has a direct effect on the granules' size, and its increment also improves powder flowability.Furthermore, more granules are collected as product within the granulator for higher binder content.Regarding fluidization air temperature, it affects binder viscosity, impacting on the particle-particle liquid bridges' strength.Interestingly, for medium and coarsesize granules, higher temperatures tend to decrease their size.Furthermore, high air

Conclusions
There is a great interest in understanding the granule growth mechanisms and the effects of process conditions over granules' quality and the operation efficiency for the co-melt fluidized bed granulation.In this sense, this work studies, by a Box-Behnken approach, the combined effect of binder mass content and fluidization air flow rate and temperature on granules' quality (particle size distribution, morphology, flow properties) and process-related variables (mass of granules collected as product in the granulator chamber, fines elutriation, and loss due to wall stickiness).It was found that binder content has a direct effect on the granules' size, and its increment also improves powder flowability.Furthermore, more granules are collected as product within the granulator for higher binder content.Regarding fluidization air temperature, it affects binder viscosity, impacting on the particle-particle liquid bridges' strength.Interestingly, for medium and coarse-size granules, higher temperatures tend to decrease their size.Furthermore, high air temperatures significantly enhance the process efficiency, since more granules are collected as product.Regarding the fluidization air flowrate, it affects the amount of particles smaller than d 10 in the system, promoting the formation of small particles due to breakage at high flowrate values as well as the presence of non-elutriated powder at low flowrates.Consequently, intermediate air flowrates are optimal for obtaining more mass collected as product.Regarding regime map studies, most runs were located in the rapid growth regime, according to the high pore saturation of the obtained granules.This result agrees with the high increment in particle size observed with respect to the initial seeds and the morphology of the final granules.Therefore, it was possible to extend Walker s regime map, identifying the growth regimes for runs that incorporate variations in fluidization air flowrate and temperature.

Figure 1 .
Figure 1.Schematic representation of the experimental set-up.

3 .
Results and Discussion 3.1.Experimental Design 3.1.1.Characterization of the Granules: Particle Size Distribution and Flow Properties Table

Figure 2 .
Figure 2. Contour plots for variations in binder mass content (W PEG ) and fluidization air temperature (T A ) for d 10 (a), d 50 (c), and d 90 (d) at 50% bypass valve apertures.d 10 for variations in bypass valve apertures (b) with other factors at intermediate level.

Powders 2023, 3 , 14 Figure 5 .
Figure 5. Contour plots and surface response for variations in binder mass content (WPEG) and fluidization air temperature (TA) for MP (a), Mf (b), and MW (c) at 50% bypass valve apertures.MP (d), Mf (e), and MW (f) for variations in bypass valve apertures with other factors at intermediate level.

Figure 5 .
Figure 5. Contour plots and surface response for variations in binder mass content (W PEG ) and fluidization air temperature (T A ) for M P (a), M f (b), and M W (c) at 50% bypass valve apertures.M P (d), M f (e), and M W (f) for variations in bypass valve apertures with other factors at intermediate level.

Table 2 .
Experimental design matrix used to evaluate the granulation process with factors as actual values and corresponding responses.

Table 3 .
Statistical values and coefficients of the mathematical models in terms of coded factors for d 10 , d 50 , d 90 and AR response variables.
a: X 1 , aperture of bypass valve; X 2 , binder mass content; X 3 , fluidization air temperature.b: intercept term.R 2 , R-squared; AP, adequate precision; LOF, lack of fit.Statistical significance of the model terms: significant p ≤ 0.05; non-significant, p greater than 0.05.

Table 4 .
Statistical values and coefficients of the mathematical models in terms of coded factors for MP, Mf, and MW response variables.

Table 4 .
Statistical values and coefficients of the mathematical models in terms of coded factors for M P , M f , and M W response variables.

Table 5 .
Average particle-binder slurry flow stress and deviation for the tested granules.

Table 6 .
Measured variables and dimensionless numbers for al runs.

Table 6 .
Measured variables and dimensionless numbers for al runs.