Influence of Thermal Conditions on Particle Properties in Fluidized Bed Layering Granulation

Fluidized bed layering granulation is frequently used to formulate particles of high quality. From previous studies, it is well known that the dynamic behavior of the process, as well as the product properties depend on operating parameters. The process is characterized by heat and mass transfer between fluidized particles and the surrounding fluidization medium. To investigate the mutual influence between particle phase and fluidization medium, a dynamic model is introduced. The model comprises two parts: a population balance model to describe the evolution of the particle sizes and a system of ordinary differential equations to account for thermal conditions. For the first time, the dynamic model considers the bidirectional coupling of particles and fluidization medium in fluidized bed layering granulation. By means of simulations, it is shown that the derived model is capable of reproducing the experimental findings.


Introduction
Fluidized bed layering granulation (FBLG) is a widely-applied particle formulation process.For instance, pharmaceutical, food, and fertilizer industries utilize FBLG to produce granules of high quality (Mörl et al. [1]).Therein, a solid-containing liquid, for instance a suspension or solution, is sprayed onto fluidized particles.Since the fluidization medium is heated, the liquid content of the spray evaporates, while the solid remains on particles surface, leading to a layer-wise growth, as illustrated in Figure 1 (Tsotsas [2]).As a result, the solid product is more durable and easier to handle than its liquid equivalent.Since the produced granules are often an intermediate product in a process chain, subsequent process steps require specific particle properties.For example, more than a minimum particle strength prevents breakage and dust formation during transport (Antonyuk et al. [3]), whereas solubility is a key feature if the granules have to be dissolved at a later stage.
To achieve high throughput, FBLG is operated in continuous mode.As particles are removed from the process continuously, nuclei have to be supplied permanently.Depending on the configuration of the FBLG, seed nuclei have to be provided by internal sub-processes such as overspray (Grünewald et al. [4]) or, as presented in Figure 2, by comminution of oversized particles (Cotabarren et al. [5]).process time t As Schmidt et al. [6][7][8] proved by experiments, operating parameters influence the dynamic behavior of the continuous FBLG significantly: while one parameter set Π 1 leads to a stable steady state, characterized by a constant particle size distribution over time, a second parameter set Π 2 may induce self-sustained oscillations.Thus, the process is classified as stable for Π 1 and unstable for Π 2 .
Besides the dynamic behavior, process conditions affect the solidification process and thereby the product properties.The influence of selected process conditions on product porosity as a key property was investigated by Rieck et al. [9] and Hoffmann et al. [10] for the formation of sodium benzoate granules.The underlying solidification process was studied in great detail by Sondej et al. [11].Therein, the influence of different operation parameters on the structure of the deposits from drying a sessile sodium benzoate droplet was studied.It was shown that the initial concentration of the solid in the droplet, as well as the thermal conditions of the surrounding fluidization medium influence the porosity and footprint of the remaining salt deposit.To account for thermal conditions within the FBLG, Rieck et al. [9] introduced the drying potential η.For granulation of sodium benzoate, a linear relation between the drying potential η and the porosity of the particle shell shell was found experimentally by Rieck et al. [9] and Hoffmann et al. [10] for FBLG with top spray configuration.Diez et al. [12] showed that a linear correlation between shell and η also holds for a bottom spray configuration.Since the particle strength, and thereby breakage behavior of particles, depends on porosity, it is an important property.Further relevant properties are the apparent particle density ρ p and average moisture content X of product particles.
As a supplement to experiments, simulations are a promising way to investigate FBLG and improve the process understanding.As Ramkrishna [13] showed, population balance models are a suitable framework to derive dynamic models of particulate processes.In previous studies, numerous authors established a variety of population balance models (PBMs) for different process configurations and purposes.For instance, Vreman et al. [14] and Heinrich et al. [15] studied the dynamic behavior of FBLG.While Vreman et al. [14] investigated the dynamics of FBLG with internal nucleation by thermal overspray, the start-up behavior of an FBLG equipped with an external screen-mill-cycle was in the focus of Heinrich et al. [15].Based on the PBM presented by Heinrich et al. [15], Radichkov et al. [16] investigated the influence of operation parameters on the qualitative behavior of the process by means of a rigorous stability analysis.So far, the dynamic models assumed an equal wetting of particle surface, leading to a uniform particle growth.From the nature of the process, as shown in Figure 2, it becomes clear that not all particles are sprinkled at a time.Therefore, Hampel et al. [17], Dreyschultze et al. [18], and Neugebauer et al. [19] introduced a dynamic model with two interconnected functional zones: wet and drying zone.In the first, particles are sprinkled with the injected liquid, while in the second, drying of the particles takes place.Based on the derived models, Dreyschultze et al. [18], Bück et al. [20], and Neugebauer et al. [19] investigated the influence of zone formation and process configurations on process behavior by means of a rigorous stability analysis.However, thermal conditions were not taken explicitly into account.Further, it turned out that the effect of zone formation on overall dynamics is often small and is therefore neglected throughout this contribution.
Models that consider the thermal condition were presented by Drechsler et al. [21] and Dosta et al. [22].There, the impact of the particulate phase on the thermal states was studied.On the contrary, the models did not consider the influence of thermal conditions on particle properties.Thus, the models presented so far are not capable of predicting any product property besides particle size L.
In the present contribution, the interplay of thermal conditions and product properties, as well as the dynamic behavior of the FBLG depicted in Figure 2 is studied.Therefore, the population balance model presented by Radichkov et al. [16] is extended by a system of ordinary differential equations to account for thermal conditions.The correlation of η and shell , as presented by Rieck et al. [9] and Hoffmann et al. [10], is utilized to link thermal conditions and particle growth.The developed dynamic model is presented in Section 2. In Section 3, the influence of operating parameters on product properties, as well as the interplay of particulate phase and thermal conditions are investigated by means of process simulations for granulation of sodium benzoate.Finally, conclusions and an outlook on future work are given in Section 4.

Mathematical Model
The coupling of particles and the surrounding fluidization medium is characteristic for FBLG.Both particles and fluidization medium are characterized by various properties.A particle for instance can be described by the corresponding individual size, temperature, and moisture content, while the state of the fluidization medium is determined by its temperature and moisture.From the operator's point of view, the size of particles is often a key property in particulate processes such as FBLG.To gain detailed information about the size of all particles within the bed, the particle size distribution n (t, L) is introduced.It represents the number of particles with size L at specific time t.To describe the dynamic behavior of n (t, L), a population balance model (PBM) is introduced.Further important process variables are the temperature, moisture content, and dry mass of particles and fluidization medium, respectively.It is assumed that both particles and fluidization medium are well mixed.Therefore, the temperature of particles θ p (t) and fluid θ f (t), the moisture content of particles X (t) and fluid Y (t), and the dry mass of particles m p,dry (t) and fluid m f ,dry (t) can be modeled as lumped states.The dynamic behavior of those states is described by ordinary differential equations.In the remainder of this section, the dynamic model of the FBLG is introduced.The proposed model describes continuously-operated FBLG with the external screen-mill-cycle only.However, to account for further process configurations, e.g., a batch-wise operation, only slight adaptations of the dynamic model are required.As depicted in Figure 3, the model consists of two subsystems.The first subsystem is the PBM describing the dynamics of n (t, L) over time t and particle size L. The second subsystem, a system of ordinary differential equations (ODEs), describes the dynamics of the thermal conditions within the FBLG.It is derived by utilizing mass and energy balances.As will be shown, the two sub-systems are interconnected.

Population Balance Model
The presented PBM is based on the dynamic model of FBLG presented by Radichkov et al. [16].To focus on characteristic process behavior, only particle growth, particle withdrawal, screening, and milling are considered.Breakage and aggregation can be neglected due to a proper choice of operation conditions in practice (Ennis et al. [23]).Therefore, the PBM reads: To determine the particle growth, the model presented by Mörl et al. [1] is utilized.Therein, a uniform particle growth is assumed depending on the total surface of the particle bed A p (t) = π ∞ 0 L 2 n (t, L) dL.To account for the porosity of the particle shell shell , the growth model is extended: with the total spray rate ṁinj , the solid fraction of the spray x inj , and the density of the solid fraction of the spray ρ inj,s .The shell porosity shell correlates, as Rieck et al. [9] and Hoffmann et al. [10] showed for the layering granulation of sodium benzoate, linearly to the drying potential η: The drying potential η, a measure of the efficiency of the drying process, is defined as the ratio between the actual drying rate and the drying rate at saturation: It depends on the moisture content of the fluidization medium inside the granulation chamber Y, the saturation moisture Y sat , and the moisture of the fluidization medium at inlet Y f ,in .The corresponding parameters ∆ shell and shell,0 are presented in the experimental study Rieck et al. [9].The apparent average porosity of particles is defined as: Since the investigated FBLG is operated continuously, product particles have to be withdrawn from the process, while nuclei have to be supplied.To remove particles from the FBLG, the non-classifying particle withdrawal ṅout (t, L) = K n (t, L) is applied.The outlet gain is determined such that the bed mass m p is constant over time: The withdrawn particles are classified by screening, according to: into oversized ṅoversized , product ṅproduct , and fines fraction ṅ f ines .The separation function T screen,i is described by: While the product fraction is removed from the process, the oversized fraction is milled and, together with fines fraction, re-fed to the granulation chamber, as illustrated in Figure 3: Thereby, the milled particles ṅmill serve as nuclei for the granulation process.In general, milling of granules is a complex process.For instance, the configuration of the mill itself, as well as particle properties, e.g., porosity, shape, and size, have a significant influence on the breakage behavior of particles (Antonyuk et al. [3]).Since this contribution focuses on the influence of thermal conditions on product properties in FBLG, an idealized model of the grinding process is utilized.It is assumed that the mill imprints the normalized particle size distribution q mill (L) on the oversized particles ṅoversized : where q mill (L) is considered as normally distributed with mean diameter µ mill and variance σ 2 mill .

Heat and Mass Transfer
The dynamic behavior of thermal conditions within FBLG is described by six states: temperature of particles θ p (t) and fluid θ f (t), moisture content of particles X (t) and fluid Y (t), the mass of the dry hold-up of particles m p,dry (t), and fluidization medium m f ,dry (t).To derive the corresponding state equations, mass and enthalpy balances are used, assuming an ideally-mixed system.The moisture contents of particle and fluid phase are defined as ratios of solvent to dry mass such that X (t) = m p,solvent/m p,dry and Y (t) = m f ,solvent/m f ,dry .To describe the temporal behavior, the time derivatives of X (t) and Y (t): are derived.The required time derivatives of the solvent and the dry mass of the particle and fluid phase result from balancing the mass flows: ṁp,dry (t) = x inj ṁinj + ṁp,recycle,dry − ṁp,out,dry ṁp,solvent (t) = 1 − x inj ṁinj + X recycle ṁp,recycle,dry − X ṁp,out,dry − ṁevap (15) While x inj , ṁinj , mass flow ṁ f ,dry,in , and moisture of fluidization medium at inlet Y f ,in , as well as the mass flow of fluidization medium at outlet ṁ f ,dry,out are operation parameters, the mass flows of the recycle ṁp,recycle,dry and of the outlet ṁp,out,dry depend on the PBM presented in Section 2.1.They are determined by ṁp,i,dry (t The mass flow rate ṁevap denotes the mass transfer of evaporated solvent from particle phase to fluidization medium.Since the evaporation decreases the moisture of particles X, drying takes place.
In general, drying describes the exchange of mass and heat between a solid and a fluid due to evaporation of a liquid.According to van Meel [24], the drying process consists of three periods.As illustrated in Figure 4, the first drying period is characterized by X ≥ X crit , where X crit denotes the critical moisture of particles.In this period, the evaporation of the solvent is limited by the state of the fluidization medium only.In the second period, X crit ≤ X ≤ X eq , the evaporation rate decreases due to particle-side resistances.In the third drying period, the particle moisture X is, however, below the thermodynamic adsorption equilibrium moisture X eq such that X < X eq is yielded.In this period, no further drying takes place.
X eq X crit 0 To address the decrease of the evaporation rate ṁevap over the drying periods, van Meel [24] introduced the normalized drying velocity ν, defined as the ratio of the actual evaporation rate to the maximal evaporation rate that occurs in first drying period.Therefore, ν is unity in the first period, while it becomes zero in third period.In the second drying period, ν depends on the particle moisture X such as specific material properties.To address the influence of material properties on drying characteristics, the parameter p is introduced.Thus, the normalized drying velocity can be written as: where δ describes the normalized moisture content, which is determined by δ (X) = (X−Xeq)/ (Xcrit−Xeq).
The influence of the drying period and parameter p on ν is depicted in Figure 4.The difference between saturation moisture Y sat and Y drives the evaporation.Therefore, with the mass transfer coefficient β f p , calculated by the correlation according to Gnielinski [25], and the density of the dry fraction of the fluidization medium ρ f ,dry yielded.The saturation moisture Y sat is defined as the maximum amount of vapor the fluidization medium can carry.It is determined by utilizing the temperature θ f ,in and moisture Y f ,in .The balance equations of particle and fluid temperature θ p (t) and θ f (t) are derived from the total enthalpy of particles H p = m p,dry c p,p + c p,l X θ p and fluidization medium The expressions of the enthalpies are derived with respect to time and, afterwards, reformulated such that: and are yielded.By balancing the enthalpy flows, the time derivatives of particle enthalpy Ḣp and fluid enthalpy Ḣf are determined: It becomes evident that particle phase and fluidization medium are coupled by heat transfer: depending on the heat transfer coefficient α f p presented in Gnielinski [25] and enthalpy transfer caused by the evaporated solvent Ḣevap .The latter correlates to the evaporation rate ṁevap introduced in Equation ( 19) such that: with specific evaporation enthalpy ∆h evap at θ = 0 • C yielded.The remaining time derivatives of the enthalpies Ḣi result from the mass transfer of the injection: the particle transfers of recycle and withdrawal: Ḣout (t) = ṁp,out,dry c p,p + c p,l X θ p , (28) and the enthalpy transfer induced by the fluidization: Ḣf ,out = ṁ f ,dry,out c p,g In the remainder, the focus is on granulation of sodium benzoate, whereby water is used as the solvent.The specific heat capacities of liquid water c p,l , water vapor c p,v , dry air c p,g , as well as the evaporation enthalpy ∆h evap depend on the respective temperature.For the calculation, the correlations presented in Glück [26] are utilized.On the contrary, the specific heat capacity of sodium benzoate c p,p depends hardly on temperature under the investigated conditions (Kong et al. [27]).Therefore, the parameter value given in Appendix A is used for c p,p .
Examining the structure of the derived dynamic model, it becomes clear that the two sub-systems, presented in Sections 2.1 and 2.2, are interconnected.For instance, as depicted in Figure 5, the particle size distribution n (t, L) determines the total surface of the particle bed A p .According to (19) and ( 24), the evaporation rate ṁevap and the heat transfer Q f p are correlated to A p .Heat transfer has an impact on the enthalpy flows ( 22) and ( 23) and therefore on temperatures θ p and θ f by Equations ( 20) and ( 21).In addition, the moisture contents X and Y depend on the evaporation rate by Equations ( 13), (16), and (17).Thus, the population balance model affects the thermal conditions.Otherwise, the drying potential η depends on Y.As Equation (3) states, η affects the shell porosity shell and thus the particle growth G presented in (2).In conclusion, the thermal conditions have an impact on n (t, L); thus, a bidirectional coupling between particulate and gas phase is present.
PBM n (t, L)  system of ODEs • three-period evaporation model [? ] heat and mass transfer drying potential . Bidirectional coupling of particulate phase and thermal conditions represented by PBM n (t, L) and the system of ODEs, respectively.

Results and Discussion
The influence of operating parameters on product properties was investigated by means of simulations.For this purpose, the model presented in Section 2 was implemented in MATLAB (2018a, MathWorks, Natick, MA, USA, 2018) applying a method of lines approach.The partial differential Equation (1) was discretized using a finite volume method, in particular a first order upwind scheme with 200 equidistant grid points was used.The resulting system of ordinary differential equations was solved with the MATLAB built-in solver ode15s.
To validate the dynamic model, simulation results were compared to the outcome of the theoretical studies of Radichkov et al. [16] and Dreyschultze et al. [18], as well as to the experimental findings of Schmidt et al. [7], Rieck et al. [9], and Hoffmann et al. [10].Table 1 presents the simulation scenarios.While the first simulation scenario, inspired by [7,16,18], put the focus on the dynamic behavior with respect to the size of milled particles, the second and third scenario investigated the influence of thermal conditions on particle porosity p .Rieck et al. [9] proved by means of batch experiments, illustrated in Figure 6 and Table 2, a correlation of p with operation parameters: an increase of the injection rate ṁinj and/or a reduction of the temperature of the fluidization medium θ f ,in leads to a more porous particle surface.In the experiments, Rieck et al. [9] coated glass particles with sodium benzoate.However, similar behavior was observed by Hoffmann et al. [10] for continuous granulation of pure sodium benzoate particles.During long-term FBLG, the inlet moisture content of the fluidization medium Y f ,in may vary over time, thus the fourth simulation scenario presents the response to this typical process disturbance.
increase injection rate reduce temperature Figure 6.Influence of drying conditions on layer porosity as presented in Rieck et al. [9]: glass particles coated with sodium benzoate according to the operating parameters presented in Table 2.For all simulations, the steady state of the dynamic system with the default parameter set Π 0 , presented in Appendix A, is used as the initial state.At simulation time t = 2 h, a single parameter is manipulated according to the scenarios given in Table 1.After 15 h of simulation time, the parameter is reset to the default value.

Variations of the Size of Milled Particles
The experimental study Schmidt et al. [7] proved the significant impact of the milling process on the qualitative behavior of FBLG with an external screen-mill-cycle.While coarse milling leads to the formation of a stable steady state, characterized by a constant particle size distribution n (t, L), fine milling induces self-sustained oscillations of n (t, L).Those oscillations lead to variations of the particle size distribution and mass flows.In the worst case, they may cause a breakdown of the process.Therefore, they are highly undesired.To predict and avoid instabilities in FBLG processes, the dynamic model, presented in Section 2, has to be capable of reproducing these different forms of dynamic behavior.To prove this, the first simulation scenario is utilized.
As depicted in Figure 7a, the reduction of the diameter of milled particles from µ mill = 0.8 mm to 0.7 mm at t 1 = 2 h induces self-sustained oscillations of the normalized particle size distribution q (t, L).As a consequence, the Sauter diameter d 32 , defined as: and presented in Figure 8, varies with time, as well.By means of Equation ( 6), a constant bed mass is ensured.If, in accordance with the simulation results shown in Figure 8, the variation of particle porosity p is neglected, the overall bed volume V p can be considered as constant.Therefore, the oscillations of d 32 are motivated by the variations of the total bed surface A p . Figure 7. Normalized particle size distribution q 0 (t, L) related to the first, second, and third simulation scenario.
. Simulation results of the first scenario: the setpoint switch of mean size of milled particles µ mill from 0.8 mm to 0.7 mm at t 1 = 2 h leads to the arising of self-sustained oscillations.After the reset of µ mill at t 2 , the oscillations decay and the dynamic system reaches steady state again.
Alterations of A p change the heat transfer between fluidization medium and particles, Q f p .Thus, the temperatures of particles θ p and fluidization medium θ f vary, leading to a change of the drying process.As a consequence, the moisture contents X and Y change, inducing a variation of the drying potential η.As shell depends on η, the porosity of the shell changes, as well.This leads throughout simulation time to a change of the apparent particle porosity p and, as a consequence, to a variation Figure 8. Simulation of the first scenario: the setpoint switch of mean size of milled particles µ mill from 0.8 mm to 0.7 mm at t 1 = 2 h leads to the arising of self-sustained oscillations.After the reset of µ mill at t 2 , the oscillations decay and the dynamic system reaches steady state again.
Alterations of A p change the heat transfer between fluidization medium and particles, Q f p .Thus, the temperatures of particles θ p and fluidization medium θ f vary, leading to a change of the drying process.As a consequence, the moisture contents X and Y change, inducing a variation of the drying potential η.As shell depends on η, the porosity of the shell changes, as well.This leads throughout simulation time to a change of the apparent particle porosity p and, as a consequence, to a variation of the particles density ρ p .
At t 2 = 15 h, the parameter µ mill is reset to 0.8 mm.As a consequence, the oscillations decay such that the system settles, after sufficient long simulation time, at steady state again.Thereby, the steady state and the initial state share the same values.

Variations of Injection Rate
As depicted in Figure 6, an increased injection rate ṁinj leads to a more porous particle, surface, resulting in a higher shell porosity shell .To reproduce this experimental observation by simulations, the second scenario was introduced: the injection rate ṁinj was increased from 40 kg /h to 50 kg /h at t 1 = 2 h and, at t 2 = 15 h, reset to the default value.
The increase of the injection rate ṁinj at t 1 raised the amount of supplied salt and solvent.While the first had an impact on the particle growth and thus on normalized particle size distribution q (t, L), as well as Sauter diameter d 32 , hardly observed in simulation results shown in Figures 7b and 9, the second lead to significant variations of the thermal conditions.Since more solvent was injected, the particles moisture X increased from 47.8 g wet/kg dry to 49.1 g wet/kg dry .Considering Equation ( 18), this led, as long as the second period drying proceeds, to an enhanced evaporation rate ṁevap .As a consequence, the moisture of the fluidization medium Y was raised from 23.3 g wet/kg dry to 27.6 g wet/kg dry .The resulting drop of the drying potential η led, in accordance with Equation (3), to a higher shell porosity shell .As can be seen in Figure 9, the thermal conditions discussed so far are characterized by a fast transition from the initial to the new steady state.The transition of the apparent particle porosity p to the new steady state is much slower.This can be explained by the nature of p : the apparent particle porosity is an integral value.Therefore, all particles with an apparent porosity p differing from current shell porosity shell have to be removed from the process.After that, the apparent particle porosity p and the porosity of the shell shell share the same value.
Simulation results of the second scenario: an increase of injection rate ṁinj from 40 kg /h to 50 kg /h at t 1 induces variations of thermal conditions, resulting in an increased particle porosity p .After resetting ṁinj to the default value, the dynamic states and particle properties settle at the (initial) steady state again.Regard the different time scales of the step response at t 1 and t 2 illustrated in Figure 10.
At t 2 , the injection rate ṁinj was decreased to the default value.Again, normalized particle size distribution q (t, L) and Sauter diameter d 32 varied only slightly.The thermal states X, Y, θ p , and θ f reached the new steady state rapidly.Due to the decreasing moisture of the fluidization medium Y, the drying potential η increased.This led to a reduced shell porosity shell and, as a consequence, to a decreased particle porosity p .The transition of p to the corresponding steady state value, shell , required a long period of time again.A comparison of the transition of p surrounding the steps Figure 9. Simulation results of the second scenario: an increase of injection rate ṁinj from 40 kg /h to 50 kg /h at t 1 induces variations of thermal conditions, resulting in an increased particle porosity p .After resetting ṁinj to the default value, the dynamic states and particle properties settle at the (initial) steady state again.Regard the different time scales of the step response at t 1 and t 2 illustrated in Figure 10.
At t 2 , the injection rate ṁinj was decreased to the default value.Again, normalized particle size distribution q (t, L) and Sauter diameter d 32 varied only slightly.The thermal states X, Y, θ p , and θ f reached the new steady state rapidly.Due to the decreasing moisture of the fluidization medium Y, the drying potential η increased.This led to a reduced shell porosity shell and, as a consequence, to a decreased particle porosity p .The transition of p to the corresponding steady state value, shell , required a long period of time again.A comparison of the transition of p surrounding the steps at t 1 and t 2 is presented in Figure 10.Therefore, ∆t = t − t i and the deviation of particle porosity ∆ p = p (t) − p (t i ) were introduced for i = {1, 2}.It became evident that the time constant τ 1 , belonging to the step at t 1 , was smaller than time constant τ 2 .The simulation results of the second scenario are in qualitative agreement with the experimental observations of Rieck et al. [9] presented in Figure 6.The increased injection rate ṁinj led to a higher shell porosity shell and an increased particle porosity p .

Variations of the Temperature of the Fluidization Medium
In further experiments, depicted in Figure 6 and Table 2, Rieck et al. [9] showed that the temperature of the fluidization medium θ f ,in influences the shell porosity shell : a lowering of θ f ,in increased shell .To address these experiments, the third simulation scenario was introduced.There, the temperature θ f ,in was decreased from 95 • C to 70 • C at t 1 and, at t 2 , reset to the default value of 95 • C.
As presented in Figures 7c and 11, the manipulation of θ f ,in at t 1 has a slight influence on the normalized particle size distribution q (t, L) and Sauter diameter d 32 .Otherwise, the reduction of θ f ,in led to a shift of the saturation point.As a consequence, the saturation moisture Y sat descended from 32.0 g wet/kg dry to 23.3 g wet/kg dry .A reduced saturation moisture Y sat resulted, according to Equation (19), in a decreased evaporation rate ṁevap .Following Equation (3), a reduction of the saturation moisture Y sat increased the drying potential η.Since shell is negatively correlated with η, the porosity of the shell increased.This lead, over simulation time, to an increasing particle porosity p .At t 2 , the temperature of the fluidization medium was reset.As in the previous simulation scenario, the system settled at the steady state corresponding to the default parameter set.
Again, the simulation results and experimental findings of Rieck et al. [9] are in qualitative agreement: A reduction of θ f ,in leads to an increased shell and particle porosity.

Disturbance of the Inlet Moisture Content of the Fluidization Medium
In general, continuously-operated FBLG is a slow process.Therefore, the environmental conditions may alter over time.For instance, the temperature and moisture content of the fed fluidization medium, e.g., air, may change.While the temperature θ f ,in can be changed by heating the fluid to the desired value, the moisture of the fluidization medium Y f ,in can only be manipulated with much effort.By accepting variations of Y f ,in as process disturbances, one often forgoes controlling Y f ,in .
To study the influence of the described disturbance on the dynamic behavior, as well as the product properties, the four scenario was introduced: Y f ,in was increased from default value 6.0 g wet/kg dry to 15.0 g wet/kg dry at t 1 and reset at t 2 .Again, the simulation results and experimental findings of Rieck et al. [9] are in qualitative agreement: A reduction of θ f ,in leads to an increased shell and particle porosity.

Disturbance of the Inlet Moisture Content of the Fluidization Medium
In general, continuously-operated FBLG is a slow process.Therefore, the environmental conditions may alter over time.For instance, the temperature and moisture content of the fed fluidization medium, e.g., air, may change.While the temperature θ f ,in can be changed by heating the fluid to the desired value, the moisture of the fluidization medium Y f ,in can only be manipulated with much effort.By accepting variations of Y f ,in as process disturbances, one often forgoes controlling Y f ,in .
To study the influence of the described disturbance on the dynamic behavior, as well as the product properties, the four scenario was introduced: Y f ,in was increased from default value 6.0 g wet/kg dry to 15.0 g wet/kg dry at t 1 and reset at t 2 .
Referring to Figure 12, the disturbance has neither a significant influence on the qualitative behavior, nor on particle properties.While the raised Y f ,in at t 1 increased the moisture content Y from 23.3 g wet/kg dry to 32.2 g wet/kg dry , it also shifted the saturation point.As a result, Y sat was raised from 32.7 g wet/kg dry to 36.5 g wet/kg dry .Even if the changes of Y f ,in , Y, and Y sat were quite large, the drying potential η varied only slightly as Y Sat − Y f ,in remained almost constant.The small shift of η was reduced from 0.334 to 0.312, inducing slight variations of shell porosity shell and particle porosity p .Therefore, the normalized particle size distribution q (t, L) stayed almost constant, resulting in a uniform Sauter diameter d 32 .
Figure 11.Simulation results of the third scenario: the reduction of the temperature of the fluidization medium θ f ,in from 95 • C to 70 • C at t 1 shifts the saturation point.This induces variations of thermal conditions and thus to a change of shell shell and particle porosity p .After resetting θ f ,in to the default value at t 2 , the dynamic states and particle properties settle at the (initial) steady state again.
Referring to Figure 12, the disturbance has neither a significant influence on the qualitative behavior, nor on particle properties.While the raised Y f ,in at t 1 increased the moisture content Y from 23.3 g wet/kg dry to 32.2 g wet/kg dry , it also shifted the saturation point.As a result, Y sat was raised from 32.7 g wet/kg dry to 36.5 g wet/kg dry .Even if the changes of Y f ,in , Y, and Y sat were quite large, the drying potential η varied only slightly as Y Sat − Y f ,in remained almost constant.The small shift of η was reduced from 0.334 to 0.312, inducing slight variations of shell porosity shell and particle porosity p .Therefore, the normalized particle size distribution q (t, L) stayed almost constant, resulting in a uniform Sauter diameter d 32 .

Figure 4 .
Figure 4. Normalized drying velocity ν over moisture content of particles X as introduced by van Meel [24] dependent on different p.

t 2 Figure 10 .
Figure 10.Normalized particle porosity ∆ p surrounding the set point switches of ṁinj at t 1 and t 2 according to the second simulation scenario.

Table 1 .
Simulation scenarios according to the simulation results presented in Sections 3.1-3.4.
Simulation results of the third scenario: the reduction of the temperature of the fluidization medium θ f ,in from 95 • C to 70 • C at t 1 shifts the saturation point.This induces variations of thermal conditions and thus to a change of shell shell and particle porosity p .After resetting θ f ,in to the default value at t 2 , the dynamic states and particle properties settle at the (initial) steady state again.
Simulation results of fourth scenario: the disturbance of the moisture content of the fluid at inlet Y in at t 1 results in small variations of the thermal conditions and particle properties.