Analysis of the Inﬂuence of Process and Formulation Properties on the Drying Behavior of Pharmaceutical Granules in a Semi-Continuous Fluid Bed Drying System

: In the last decade, twin-screw wet granulation became an essential technology for continuous pharmaceutical tablet production. Consequently, interest in (semi-)continuous ﬂuidized bed drying systems as a subsequent processing unit has grown. In parallel, it has become pivotal to fully understand and control manufacturing processes in line with in the quality-by-design paradigm. Formulation-generic prediction models would enormously facilitate digitally enhanced process development and require dedicated experimental data collection and process knowledge. To obtain this knowledge, three experimental campaigns were performed in this work. Firstly, an investigation into the effect of dryer process settings on drying behavior is presented. Secondly, the effect of active pharmaceutical ingredient properties on drying was assessed by producing granules of similar particle size and porosity and evaluating their drying and breakage behavior. Finally, additional experiments with varying active pharmaceutical ingredients and drug load were conducted to increase the genericity of the data set. This knowledge can be used in mathematical process modelling.


Introduction
A commonly used technique for particle enlargement in pharmaceutical manufacturing is granulation. Based on the European Public Assessments Reports, Leane et al. [1] calculated that 55% of the formulations include a wet granulation step in the tablet manufacturing pathway. At the same time, continuous manufacturing technology has gained interest from pharmaceutical companies, as this production method is supported by the Food and Drug Administration. Since then, the redesign of production processes and the establishment of continuous wet granulation lines in the pharmaceutical industry has manifested itself [2]. However, the number of industrial applications is limited, which highlights the complexity of developing and applying the scientific understanding of this innovative method in actual pharmaceutical products [2,3].
During twin-screw wet granulation (TSWG), powder mixtures are agglomerated into semi-permanent aggregates by means of the addition of a binder and liquid to the processed material [4]. The granulation liquid, together with the binder, binds the powder particles together in the wet granule, through a combination of capillary and viscous forces. Subsequently, more permanent bonds are formed through binder solidification by a pivotal drying step [5]. Despite the requirement of both unit operations in a continuous process line, the scientific literature has mainly focused on the optimization of the granulation settings (liquid-to-solid (LS) ratio, mass flow rate (MFR), screw speed) [6][7][8] or screw configurations [9][10][11]. To the authors' knowledge, only the studies of De Leersnyder et al. [12] and Ryckaert et al. [13] have experimentally investigated the influence of process parameters on the breakage and drying behavior of pharmaceutical granules in semi-continuous fluid bed drying (FBD) processes after wet granulation.
However, both studies focused on a single formulation, making it difficult to generalize the results. To establish knowledge over a range of formulations, it is important to compare the breakage and drying behavior of multiple formulations. Especially, for the development and validation of generic FBD models, more accurate data, most notably on residual moisture content (RMC) and granule size, are necessary [14]. For this reason, three different experimental campaigns were executed in this study on a continuous TSWG manufacturing device. The first two campaigns had a different focus, allowing us to quantify the impact of, respectively, the process settings and active pharmaceutical ingredient (API) properties on the drying behavior of wet granules. The third campaign aimed to further generalize the data set through adding experiments with different formulations in order to increase the genericity of the data set.

Materials and Methods
A ConsiGma™-25 system (GEA engineering, Wommelgem, Belgium) was used to produce and dry the granular material, as described by Fonteyne et al. [15]. This continuous manufacturing equipment is composed of a twin-screw granulator, a six-segmented fluidbed dryer, and a product control hopper.
All experiments were performed with a fixed excipient base: 30% w/w Pharmatose 200M (DFE Pharma, Goch, Germany), 15% w/w Avicel PH-101 (DuPont, Wilmington, DE, USA), 5% w/w Methocel E15 (DuPont, Wilmington, DE, USA) and 50% w/w API. The API component is the only formulation fraction that cannot be substituted during formulation optimization. To quantify the effect of the API properties on the overall drying behavior, formulations with different components were compared in this study. Furthermore, none of the API components under study exhibit degradation behavior within the temperature range tested.
Before granulation, a batch blending procedure at 25 rpm for 15 min with a 20 L tumbling blender (Inversina-Bioengineering, Basel, Switzerland) was performed with the aforementioned powder mixtures. The load of each batch never exceeded 5 kg, to ensure sufficient free volume during blending.
The preblend was fed into the granulator using a twin-screw loss-in-weight feeder (KT20, K-Tron Soder, Niederlenz, Switzerland). The MFR was evaluated during the J1F50 experiments. However, no major effect was observed with respect to the MFR. Consequently, a fixed MFR value of 20 kg/h was used for the remaining experiments. The temperature of the cooling jacket was kept constant at 25°C with an active cooling system. The granulator is composed of two co-rotating self wiping screws of 25 mm diameter, with a length-to-diameter ratio of 20:1. The configuration of each screw was composed of four different zones or compartments: an initial zone of conveying elements, followed by two kneading compartments (six kneading elements at 60°forward angle) and a size control compartment (three size control elements at 60°forward angle) right before the granulator outlet. Each compartment was separated by conveying elements of the same length, equal to 1.5 times the diameter of each kneading compartment [16]. The screw speed was kept constant at 675 rpm. The LS ratio was varied according to the research focus of each campaign (see Section 2.1). Distilled water was used as granulation liquid for all experiments since the preblend already contained a dry-added binder component.
After granulation, the wet granules were pneumatically transferred to the six-segmented fluid bed dryer. The inlet air temperature was varied between 40°C and 60°C, while the inlet airflow was kept constant at 400 m 3 /h. A fixed nominal value for the inlet airflow was used during this study to limit the experimental workload. However, during each experiment, the fluidization quality was checked and found to be acceptable by visual assessment. Additionally, before each experimental run, a preheating period in an empty dryer was performed to ensure a stable environment in terms of inlet air temperature and velocity from the beginning of the run. Furthermore, each first-filled cell after start-up or a change in granulator settings was considered as a waste cell on which no samples were collected.
After the desired drying time had elapsed, the dried granules were removed from the dryer and pneumatically transferred to the product control hopper, where they were collected for further granule characterization. The contents of one drying cell were considered to be one population, representative of the drying process.

Experimental Design
In this study, three experimental campaigns, each with a different focus, were conducted. Different active pharmaceutical ingredients (APIs) with different material properties were used to perform the experiments. The APIs used are denoted by J1, J2, J3, U1, U2 or G1. An overview of the API properties, scaled between zero and one, is given in Table A1. To further distinguish clearly the powder mixtures used, the filling grade (semi-fine SF, fine F, regular R, or dense D) and API loading (mass percentage) were also added to the formulation code per experiment.

Impact of Process Parameters
The first experimental campaign was performed to unravel the effect of the different process settings on the overall drying behavior of pharmaceutical granules. The selected APIs, J1F50, and U1R50 are both highly soluble in water, as this property facilitates the production of granules by the TSWG process.
A two-level full factorial design for two factors (2 2 ) was used for this set of drying experiments. The factors were the inlet air temperature and the fill level of the drying cells. Furthermore, three centerpoint runs were included to measure process reproducibility and inherent variability. This design of experiments (DoE)DoE was performed for both formulations at two LS ratios close to the limits of the formulation-specific workable LS ratio. Table A2 summarizes the performed experiments for the first campaign. Furthermore, the order of execution was randomized to avoid bias.
To evaluate the effect of drying time on the RMC, granules were sampled at four different drying times. To obtain an estimation of the granule moisture content at the beginning of the drying process, the shortest drying time was set as 20 s after the filling time. The filling time is the time necessary to fill one drying cell and can be calculated by dividing the desired fill level with the MFR. The longest drying time was equal to six times the filling time due to the design and working principle of the equipment. However, to allow for the continuous flow of material, the total drying time had to be slightly reduced. Evaporation is an endothermic process, resulting in a decrease of both the bed and outlet air temperature when wet granules are being dried. After evaporation, the granules will heat up, because the drying rate on the surface of the granules is faster than the diffusion of water from inside the granules to the surface. Both, the second and third data points are taken before the bed temperature in the cell started to increase. This way, it could be assumed that the granules still contained moisture after sampling. Thus, taking samples of the granules at these drying times provided relevant information about the drying process. The granule temperature profile of the longest drying experiment was used to define these two sampling times (Figure 1). A summary of the first campaign experiments with their corresponding drying times is provided in Table A2

Impact of API Properties
To unravel the effect of API properties on the drying behavior of wet granules, a second experimental campaign with varied drug characteristics was performed. The U2R50 and U2SF50 blends were more hydrophilic compared to the hydrophobic G1R50 formulation. U2R50 and U2SF50 were blends with the same API but varied milling grades; in this way, differences in flowability could be observed.
To isolate the effect of API, granules with similar properties needed to be created with the various blends. Preliminary granulation experiments, with a step-wise increase of the LS ratio at increments of 0.5% w/w , were conducted to determine those granulation settings at which granules with similar granule properties could be obtained (e.g., granule size distribution (GSD), granule porosity and initial moisture content). The LS ratio was varied between 18 and 25% w/w and within this range granules could be produced for all formulations in study. Next, the volume-based GSDs were subjected to the kernel principle component analysis (KPCA) technique, by which the distributions were transformed into a single point in a dimensional space described by principal components [17]. The results are visualized in Figure A1.
The maximum mean discrepancy (MMD) is a distance-measure between distributions, which is defined as the squared distance between their embeddings in a dimensional space [18]. The more similar the distributions are, the lower the MMD value is. The MMD was calculated between the distributions obtained with the same granulation settings, Figure A2. Those granulation settings for which the MMD is the lowest have the most similar GSD [17]. Based on this analysis, the LS ratios of 18.5% w/w , 19.5% w/w and 21.5% w/w were chosen to perform the drying experiments.
It is assumed that the same LS ratio results in a similar initial moisture content. Furthermore, based on expert knowledge, it was assumed that using the same LS ratio also resulted in similar granule porosity values after granulation. This was additionally measured and confirmed using He-pycnometry (Accupyc 1330 pycnometer, Micromeritics Instruments Inc., Norcross, GA, USA).
The experimental plan of the second campaign investigated trends in the drying behavior of the granules under investigation, by varying the LS ratio, fill level and inlet air temperature. The inlet air temperature was varied from 40°C to 60°C and the filling degree from 500 g to 1000 g. Each process setting was varied one by one while the other process settings were kept at the values of the baseline experiment. An LS ratio of 21.5% w/w and inlet air temperature of 40°C were chosen for this experiment. To minimize material consumption, a fill level of 500 g was used in the baseline experiment. Three repetitions were conducted of this experiment in order to have a measure of the experimental variability. Table A3 gives an overview of the conducted experiments during the second experimental campaign.
To analyze the drying curves, different samples were taken at different drying times during each experiment. These sample points were determined similarly to the method used in the first experimental campaign. However, in this case, only the baseline experiments were used to determine the sample times. Subsequently, the same time steps were used, making a one-to-one comparison between the different experiments possible.

Generalisation of the Data Set
A third experimental campaign was conducted to further increase the generic nature of the data set with the aim to enable the development of a mathematical model that predicts the drying process based on the API properties. Additionally, these experiments could also be used to confirm previously observed phenomena or to validate drying model predictions. Experiments for three additional blends were performed: J1R50, J2F5, and J3R50. Within the previously tested process settings ranges, three additional experiments were defined for each blend. One experiment where a slow drying rate was expected, one where this rate was expected to be fast, and an intermediate situation. The J2F5 experiments were performed with a 5% drug load to diversify the overall data set. An overview of the conducted experiments is given in Table A4.

Granule Size Measurements
Granule size was analyzed via dynamic image analysis (frame rate 70 Hz) using the QICPIC™ system (Sympatec, Clausthal-Zellerfeld, Germany) equipped with a vibrating dosing unit (VIBRI/L™) for gravimetrical feeding of the granules (feed rate 20%. The WINDOX software (Sympatec, Clausthal-Zellerfeld, Germany) was used, to report the measured GSD, visualizations were made in Python using the Seaborn package [19].

Residual Moisture Content
The RMC of the dried samples was determined using the loss-on-drying (LoD) method. Immediately after sampling at the hopper, the dried granules were sieved using a sieve tower (mesh sizes of 150 µm, 300 µm, 710 µm, 1000 µm, 1400 µm, 1700 µm, 2000 µm and 3150 µm). Subsequently, for each size class, a representative sample was taken and weighed (m 1 ). Then, the dishes were placed in an oven (Memmert N60, Memmert GmbH+Co.KG, Schwabach, Germany) for at least 24 h at 105°C after which they were reweighed (m 2 ). Of each measured weight, the mass of the aluminum dish (m 0 ) was subtracted after which the LoD (% w/w ) could be calculated according to Equation (1).

Flowability Test
The granule flowability was characterized using the GranuHeap™ (Granutools, Belgium), an automated heap shape measurement method. At the start of the measurement, a cylindrical tube was loaded with a volume of powder (60 mL) and lifted at a constant speed of 5 mm/s, resulting in a granule heap on the support plate. Next, an image recognition algorithm was used to determine the angle of the isosceles triangle or Angle of Repose (AoR) from 17 pictures, each taken with an angular difference of 10.5°. In addition, during the same measurement, the static cohesive index, a measure for the cohesive strength of the granules, was quantified. The static cohesive index indicates the deviation between the real heap interface and the isosceles triangular heap.

Friability Test
To quantify the mechanical stability, the friability (F) was measured by using a friabilator (Pharma Test Friabilator, PTFE, Hamburg, Germany) equipped with an abrasion drum. During this procedure, 10 g of granules (m F,1 ) with a granule size above 250 µm and 200 glass beads (diameter: 4 mm, mass: 8.5 mg) were rotated for 10 min at 25 rpm in the abrasion drum. Afterward, the glass beads and the fraction with diameter smaller than 250 µm were separated by sieving. The remaining granules were reweighed (m F,2 ). With m F,1 and m F,2 , F can be computed as expressed in Equation (2).

Results and Discussion
The experimental results of the first (Section 2.1.1) and second (Section 2.1.2) experimental campaigns are discussed in Sections 3.1-3.4. The experimental results are presented and discussed by data type, i.e., GSD, RMC, flowability, and friability, respectively. Subsequently, Section 3.5 discusses the experimental results of the third campaign (Section 2.1.3) in which the previously described consistent drying behavior of pharmaceutical granules was verified for three new formulations.

Granule Size Distribution
The GSDs were visualized for all experiments in function of the drying time. Note that the filling time was varied during the experiments as it is proportional to the desired fill level of the drying cell. To allow comparison between the experiments, all results were reported as a function of absolute drying time after cell filling.
Furthermore, for each experiment, a comparison was made with the oven dried GSD collected at the outlet of the granulator, the initial GSD. This GSD provides a good indication of the GSDs that were entering the ConsiGma™- 25   Furthermore, a modified screw configuration was used for the U1R50 experiments, which is different from the one discussed in Section 2. Although the screw configuration was very similar, the size-control compartment was excluded during the U1R50 experiments. The length and sequence of the other screw compartments remained unchanged concerning all other formulations. This adaption was necessary to obtain usable U1R50 granules for the corresponding drying experiments. Further investigation into the effect of formulation properties on the granulation behavior was considered to be beyond the scope of this study due to the focus on the FBD process.  For all experiments, even those with a short drying time, an increase in fine material was observed compared to the samples taken at the outlet of the granulator. This indicates that breakage and attrition occurs after the granulation process. However, this breakage and attrition cannot be directly ascribed to the drying stage only, since two pneumatic transfer lines were also involved prior to sampling. These are the wet and dry transfer lines before and after the fluidization process, respectively. If the breakage and attrition are caused by the drying process, the fines fraction should increase in function of drying time.
However, an increase of the fines fraction in function of process time is not observed. From the experiments with varied granulation settings ( Figure 2) it can be seen that breakage and attrition were much more pronounced for the low LS granules in comparison to the high LS granules. Badawy and Pandery [20] reported that wet granules are more prone to plastic deformation as opposed to breakage and attrition. This results in less breakage in the wet transfer line and early in the drying process. During the drying process, liquid bridges will be solidified, resulting in strong interparticle bridges. This process is more pronounced at high LS ratios, as more liquid bridges are present. The solidified interparticle bridges contribute to a higher granule strength at the later stage of the drying cycle or the dry transfer line. It seems that an optimal moisture content, to minimize breakage by pneumatic transport, could be reached. At that point, an optimal granule strength was obtained: a balance between weakly bound wet granules and overly brittle dry granules. This observation is in line with the work of Ryckaert et al. [13], who stated that the extent of breakage and attrition caused by the dry transfer line is positively correlated to the RMC of the granules above the optimal moisture content for granule strength. Furthermore, Ryckaert et al. [13] demonstrate that the dry transfer line generates extensive breakage, as the drying cells are forcefully emptied pneumatically.
In addition, Ryckaert et al. [13] described that higher drying rates will change the internal structure of granules, as a result of vapor accumulation in the granule pores. This accumulation result in an increase of the internal granule temperature and pressure, which can alter the microstructure making the granules more susceptible to mechanical stress. As a result, an effect of the inlet air temperature was expected. In contrast, based on Figures 2, A3, A5, A7, and A9, no clear trend was observed by considering the inlet air temperature in the collected generic data set. The results obtained are more in line with the observations of De Leersnyder et al. [12], which also did not observe any effect of the inlet air temperature on the GSDs after drying.
Leersnyder et al. [12] demonstrate that there is a positive link between inlet airflow and the yield fraction, at the expense of the oversized granules. The increase in interparticle and particle-wall collisions is considered to be the main cause for this correlation. Thus, a similar effect was expected at higher fill levels, due to the fact that increasing the amount of granular material also makes interparticle collisions more likely. In contrast, no clear effect of the fill level on the GSDs was observed for all formulations in the present work.
It must, however, be noted that the drying airflow rate was kept constant at 400 m 3 /s since the fluidization was acceptable across the different fill levels for all experiments. It might be possible that the fluidization was less optimal at higher fill levels, since the amount of material present was higher in the dryer cells. In that case, the granules would have had a lower velocity, resulting in fewer collisions. Further studies, with an adjusted experimental design, must be performed to prove whether the impact of filling degree on granule breakage and attrition behavior can be completely neglected during fluid bed drying.
The preliminary granulation experiments of the second experimental campaign resulted in granules of similar GSD for the U2SF50, U2D50, and G1R50 experiments at the beginning of the drying process. However, despite the similar GSD in the prior experiments, the drying process affected the different formulations differently. U2D50 ( Figure A8) showed the least and G1R50 ( Figure A10) the most granule degradation after drying. This observation could be attributed to the lower water solubility of G1R50 compared to U2SF50 and U2D50. Lower water solubility reduces the ability to form interparticle bridges. Fewer interparticle water bridges lead to more fragile G1R50 granules compared to the others. This phenomenon was most pronounced at the higher LS ratios, as more granulation liquid is available in that case.
In addition, a difference in granule degradation was observed between both API grades. It was expected that the grade with the smaller particles or higher contact area of the semi-fine formulation would permit more contact points or liquid bridges between the colliding particles. These bindings may then lead to stronger granules [21]. However, the opposite effect was observed-the U2D50 granules ( Figure A8) were stronger than the U2SF50 granules ( Figure A6). This is in line with the findings of Zuurman et al. [22], who stated that an increase in bulk density of lactose granules results in a decrease in its ability to consolidate, as the bulk is more porous and more difficult to pack.

Residual Moisture Content Distribution
The sampling times were chosen using the procedure described in Sections 2.1.1 and 2.1.2. Consequently, the drying times were not identical between the experiments of the different formulations which made direct comparison difficult. Therefore, the residual moisture content distributions were reported taking into account the drying time after filling. In this way a direct comparison between the different fill level experiments is possible.
Since the experimentally observed LoD behavior for all formulations was aligned with the behavior observed during the J1F50 experiments, solely the RMC data in the function of the process settings of the J1F50 experiments are reported in Figure 4. A detailed overview of the other collected data points can be found in Appendix D.2 ( Figures A11-A14).  The fact that the GSD were optimized before the drying experiment resulted in small sample sizes for the fines and oversized fractions. Consequently, those LoD measurements (standard deviation over repeated experiments, averaged over different formulations: 21%) were less reliable and more sensitive to variability compared to the yield fraction (standard deviation: 5%). Only qualitative conclusions could be made for the more extreme granule sizes, while the middle-sized classes (710-1400 µm) also allow a quantitative interpretation.
Based on the experimental data, it was concluded that granulation at low LS ratios results in shorter drying times compared to the higher LS ratios. This was expected, as the granules at low LS ratios contain less moisture. However, differences in initial moisture distributions in function of particle size were observed between high and low LS granulation settings. The distribution of the moisture across the sieve fractions seemed less uniform at low LS ratios. For the latter case, the larger granules contained more moisture compared to the smaller granules (e.g., Figure 4). This indicates less efficient distribution of the granulation liquid at low LS ratios, leading to wet and dry regions within the granulator. Furthermore, Figure 4 reveals a systematic increase of the moisture content with the granule size. However, the size-based RMC variation became narrower with increasing drying time. This observation indicates a proportional drying rate with granule size. However, this observation does not take into account the LoD measurements of the extreme size classes, which may give a biased result. Implementation and calibration of these phenomena in model structures contributes to the further knowledge development and quantification of the observed nonlinear drying processes.
The LoD measurements revealed that increasing the inlet air temperature resulted in a decreased necessary drying time for granule batches with similar properties. For example, the LoD values at the second drying time of experiments N1 and N3 were almost identical ( Figure 4). However, the second drying time of N3 was shorter than the second drying time of N1 with a lower inlet air temperature. This correlation was expected, since an increase in inlet air temperature results in higher drying rates. In other words, the RMC declined faster before the equilibrium moisture content was reached.
Considering the other drying parameters, it was observed that the LoD changes similarly at the shortest drying times for variable fill levels. However, after longer drying times, higher LoD values were observed at the high cell loads compared to the lower fill levels. As discussed in Section 3.1, there was probably less optimal fluidization apparent at the higher filling rates, which may have resulted in lower drying rates due to less contact between the drying air and the drying granules. In addition, at a higher fill level, there is more material in the dryer cell, so more liquid has to evaporate. The drying air in the cell will be more saturated, which may lead to a decrease in drying rate. However, by extending the total drying time, the equilibrium moisture content could be reached for all experiments. In conclusion, the rate at which the RMC was changing was dependent on the fill level, but this process variable is not a limiting factor. Analyzing the influence of the incoming air velocity was beyond the scope of this study.
In general, when comparing the G1R50, U2SF50 and U2D50 experiments, which were performed with the same granule porosity, GSD, and process settings, slight differences in drying rate were observed. For all sampled drying times, the G1R50 granules consistently contained the most moisture, while the U2D50 granules were the driest. However, this cannot be stated with certainty, since the differences were minimal and not greater than the standard deviation.
In any case, a consistent trend in RMC was observed with respect to the heat capacity (Cp) values or energy requirements for heating. G1R50, U2SF50 and U2D50 have a slightly lower Cp value, respectively. As such, more heat energy was required to heat up the solid fraction of G1R50 compared to the other formulations. This energy could no longer be used to evaporate moisture, resulting in lower drying rates or higher RMC values after a constant drying time.

Flowability
The flowability of the granules after drying was evaluated in this study, because one of the objectives of wet granulation is to improve this critical quality attribute (CGA). Poor flow characteristics can cause problems downstream during tableting. Both the AoR and the static cohesive index (CI) were measured to characterize the granule flowability.
The resulting AoRs following from the J1F50 experiments are reported in Figure 5. This data set contains the AoR in function of the LS ratio and drying time for the various drying scenarios. Figure 6 compares the U2SF50, U2D50, or G1R50 experiments, as they were conducted with similar granule properties. The AoR measurement is only reported for the samples with the longest drying time, as these drying times are most likely to be used for industrial applications. The effect of LS ratio, fill level (FL) and inlet air temperature were investigated. In addition, the flowability classification system described by Granutools™ [23] is given as a qualitative measure in Figures 5 and 6. Furthermore, the obtained CI results were in line with the measured AoR values, and can be found in Appendix D.3 ( Figures A17 and A18). Figure 5 shows that the variability in the AoR had no relation to the drying time. Furthermore, no clear difference was observed between the J1F50 experiments, provided that the granulator settings were the same. This suggested that the fill level and the inlet air temperature do not affect the AoR. Figure 6 confirms these observations for three additional formulations. This graph reveals no clear influence of the inlet air temperature on the granule flowability. This is in line with the observation that the inlet air temperature did not have an observable influence on the GSD. The granule size and shape are reported to be the granule characteristics that have the largest correlation to the flowability [24]. Results are compared between the low LS and high LS for air-dried granules, marker color corresponds to the drying times after the filling period. The x-axis corresponds with the experimental number as given in Table A2. With respect to the FL, it is difficult to make an unequivocal conclusion. The U2SF50 results showed a positive correlation between the AoR and the FL of the dryer cell, while this is not observed for the other formulations. This trend could not be explained so far as these values could not be related to the GSD results.
For both experimental campaigns, a clear difference in granule flowability was observed between the different granulator settings. The AoR was negatively impacted by the LS ratio, which corresponded to worse flowability. This can be fully attributed to the increase in fine material after granulation at lower LS ratios. The amount of fines has a negative impact on the observed AoR [24]. However, it is important to note that this difference was less visible in the samples taken at the end of the granulator before the drying process. So, it was concluded that the FBD unit induced the negative correlation between the LS ratio and the flowability. Once again, the reduced flowability could be attributed to breakage and attrition during fluidization.
Lastly, Figure 6 provides the opportunity to quantify the isolated blend effect, due to the fact that the granule porosity and GSD were similar across the tested formulations. G1R50 has the highest AoR values, while U2D50 showed the best flowability. This difference in flowability between G1R50 becomes more apparent with an increasing LS ratio. Again, the GSD is considered to be the best explanation of this phenomenon. As discussed in Section 3.1, the differences in GSDs could be attributed to the difference in solubility of both APIs. Section 3.1 also indicates that breakage and attrition were more pronounced for the U2SF50 granules compared to the U2D50 granules. An increased fine fraction and lower flowability for the U2SF50 granules was described. The observed difference in flowability between the two blends with different API grades decreases with the LS ratio. At lower LS ratios, less liquid is added, which reduces the formation of strong internal liquid bridges. This caused a decrease in strong granules and a more pronounced effect of degradation during fluidization or pneumatic transport. At the higher LS ratios, stronger granules are formed, which result consequently in a decreased difference of the flowability.

Friability
To avoid problems downstream in the tablet manufacturing process, the friability of the dried granules should be below 30% [25]. As such, another important CQA of wet granulation is friability. This property was measured for the U2SF50, U2D50, and G1R50 experiments on the samples with the longest drying time. These results are visualized in Figure 7. For all the applied settings and blends, the observed friability was below 30%, which was within the favourable range. Nevertheless, some changes in friability could be observed when the varying process phenomena were taken into account, Figure 7. The LS ratio has a negative correlation with the friability, when considering the results of the blends separately. This was in line with the conclusions of Kotamarthy et al. [26], who showed that an increased LS ratio leads to denser and harder granules and therefore a decreased friability of the granules. Lower LS ratios lead to fewer hydrogen bonds and fewer strong internal liquid bridges after drying, making the granules weaker. However, this trend is not proven to be significant based on these experimental data set, as the 95% confidence intervals within one blend overlapped.
Considering the inlet air temperature, no clear effect was observed as for all formulations an overlap of the confidence intervals was observed. Nevertheless, the lowest friability was always obtained at 50°C. Higher friability was expected at low inlet air temperatures after a constant drying time, since a higher moisture content leads to less internal solid bridges or weaker granules. However, in Section 3.1, it was discussed in detail that a high inlet air temperatures also leads to granule degradation during FBD. The decrease in granule plasticity at higher drying temperature, might have led to a similar friability as the experiments from 40°C. The granules dried at 50°C, have the best balance between plasticity and formed internal bridges, causing them to have the lowest friability in this data set.
Considering the fill level, no influence on the friability of the output granules was observed. With exception of the G1R50 experiment at the highest mass load, all 95% confidence intervals overlap across filling level as well as across formulation type. Indirectly, an effect of filling degree was expected through the fact that the moisture content decreases more slowly at higher cell loads. Wetter granules will show more degradation, so higher friability values were expected as well. On the other hand, when breakage occurs during the drying process, the weakest granules have already broken up and only the strongest granules remained after FBD.

Generalisation of the Data Set
This study showed consistent drying behavior given the process settings and API properties studied. However, in order to develop a general mathematical drying model, only a limited number of formulations were considered during the first two experimental campaigns. Therefore, during the third experimental campaign, both the GSD and RMC after drying were measured for three new formulations. For each formulation, three different combinations of process settings were considered, Table A4. The process settings were selected to evaluate a situation where fast, slow, or intermediate granule drying rates were expected. Only the J2F5 experiments are shown in this section, as the J1R50 and J3R50 results are consistent with the J2F5 results. A complete overview of the experimental data is in the Appendixes A-D. Figure 8 shows the granule size distributions of the J2F5 experiments in function of the LS ratio and drying time. The resulting GSDs shows a similar non-linear behavior as expected from the discussion in Section 3.1. Once more, it was evident that the drying granules reached an optimum strength point during the drying process. More granule breakage after transport through the dry transfer line occurred when the drying time was shorter or longer than the optimum time point, due to more brittle granules.   Figure 9 gives the RMC data in function of the granule size and drying time for the same J2F5 experiments as reported in Figure 8. Once again, it is concluded that the drying behavior of the investigated granules was aligned with the expectations discussed previously. Increasing the LS ratio or the degree of filling or decreasing the drying air temperature extends the required drying time to reach the equilibrium moisture content for the granules in all size classes. In general, the third experimental campaign provides a valuable source of additional data that can be used for the calibration and validation of generic predictive models for semi-continuous fluid bed drying (FBDs), as already proposed by Ghijs et al. [14] or Nitert and Hounslow [27]. Such mathematical models are tools able to capture the experimentally gained knowledge; they can be used as a virtual test environment to make predictions of the critical quality attributes at process conditions, which have not been tested before. Furthermore, such models can be used for the development of control systems, soft sensors, or a production device-specific digital twin. Especially, within the quality-by-design (QbD) framework, to which the pharmaceutical sector is transitioning, such models will provide tremendous benefits. From this perspective, this study provides an excellent experimental base for the further development of mathematical models that consider both process conditions and API properties.

Conclusions
This study elucidated the relationships between the general drying behavior in semicontinuous fluidized bed dryers and both the process settings and formulation properties were studied. It was found that, though API properties influence the drying and breakage behavior through water solubility, heat capacity, and powder particle size, the granulator LS ratio and dryer inlet air temperature dominated the process outcomes.
Similar trends concerning the residual moisture content were found across all experiments. The only process setting under study which showed a strong correlation with the drying rate was the inlet air temperature. Due to the sampling method used, this study revealed that the RMC increased with the granule size at the same drying conditions. Consequently, the mean LoD values could be biased. The width of these RMC distributions decreased with drying time. Next to that, the LS ratio was negatively impacting the breakage and attrition behavior. The largest granule size degradation was, however, not attributed to process settings but rather to the connection of the consecutive unit processes.
Additionally, different pharmaceutical blends with varying API characteristics were also considered, to increase the generic value of the obtained data set. Generic data sets are essential in the development of mathematical prediction models to make reliable predictions for blends with unknown processability. The water solubility and heat capacity could be identified as the major influencing formulation properties on the breakage and drying behavior during the fluidization process, respectively. Additionally, to a lesser extent, a negative effect of the bulk density was observed on the granule degradation behavior. Still, it was found that the effect of the process settings on the critical quality attributes was much greater than the API or blend properties. Nevertheless, it must be noted that the surface properties of the granules were not considered in this study. However, for followup studies, it might be useful to further investigate the influence of the granule surface properties on the drying behavior, for example with the BET (Brunauer-Emmett-Teller) analysis.
In addition to the breakage and evaporation behavior, the flowability and friability were also monitored during this study. Variations in both critical quality attributes (CQAs) have been fully explained by the observed differences in degradation and drying behavior. Consequently, the influential process settings and formulation properties were also found to be capable of describing the latter CQA.
Finally, this work provides an important source of data for the development of mathematical drying models for process optimization and process control, firstly through the unique generic nature of the data set and secondly through the recording of the RMC (after drying) as a function of granule size. Funding: This work is supported by the Consortium Agreement between Janssen Pharmaceutica, UCB and Ghent University.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: Appendix B Table A2. Summary of the process settings and sampling times for the first experimental campaign. From left to right, the tables gives respectively the number of the experiment, the formulation type, the LS ratio, the inlet air temperature, the mass flow rate, the total fill level of the dryer cell, the screw speed and the drying times for sampling.   Table A4. Summary of the process settings and sampling times for the third experimental campaign. From left to right, the tables contain respectively the number of the experiment, the formulation type, the liquid-to-solid (LS) ratio, the inlet air temperature, the mass flow rate, the total fill level of the dryer cell, the screw speed and the drying times for sampling.