Next Article in Journal
Photochemical Internalization with Fimaporfin: Enhanced Bleomycin Treatment for Head and Neck Cancer
Next Article in Special Issue
Integrated Continuous Wet Granulation and Drying: Process Evaluation and Comparison with Batch Processing
Previous Article in Journal
Nanoparticle-Mediated Therapy with miR-198 Sensitizes Pancreatic Cancer to Gemcitabine Treatment through Downregulation of VCP-Mediated Autophagy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Material Transport Characteristics in Planetary Roller Melt Granulation

1
Laboratory of Solids Process Engineering, Department of Biochemical and Chemical Engineering, TU Dortmund University, 44227 Dortmund, Germany
2
MeltPrep GmbH, 8020 Graz, Austria
*
Author to whom correspondence should be addressed.
Pharmaceutics 2023, 15(8), 2039; https://doi.org/10.3390/pharmaceutics15082039
Submission received: 15 June 2023 / Revised: 20 July 2023 / Accepted: 27 July 2023 / Published: 28 July 2023
(This article belongs to the Special Issue Pharmaceutical Continuous Manufacturing: Then and Now)

Abstract

:
Melt granulation for improving material handling by modifying particle size distribution offers significant advantages compared to the standard methods of dry and wet granulation in dust reduction, obviating a subsequent drying step. Furthermore, current research in pharmaceutical technology aims for continuous methods, as these have an enhanced potential to reduce product quality fluctuations. Concerning both aspects, the use of a planetary roller granulator is consequential. The process control with these machines benefits from the enhanced ratio of heated surface to processed volume, compared to the usually-applied twin-screw systems. This is related to the unique concept of planetary spindles flowing around a central spindle in a roller cylinder. Herein, the movement pattern defines the transport characteristics, which determine the energy input and overall processing conditions. The aim of this study is to investigate the residence time distribution in planetary roller melt granulation (PRMG) as an indicator for the material transport. By altering feed rate and rotation speed, the fill level in the granulator is adjusted, which directly affects the average transport velocity and mixing volume. The two-compartment model was utilized to reflect these coherences, as the model parameters symbolize the sub-processes of axial material transport and mixing.

1. Introduction

Within the last decade, research in the field of pharmaceutical technology has focused on maximizing the potential of a systematic process design, following the Quality-by-Design approach as proclaimed in the ICH Q8 (R2) guideline [1] by establishing methods for continuous operation [2,3]. As well as a reduced machinery footprint, among others [4], a central advantage is the potential to reduce product quality fluctuations in comparison to batch processes. This is based on processing in a dynamic steady state, which is beneficial for implementing process control strategies [5], as here the central task is to maintain the process in balance.
This paradigm shift in the preferred manufacturing type also includes the unit operation granulation [6,7]. In comparison to dry [8,9] and wet granulation [10,11], as the most commonly applied methods in pharmaceutical technology, melt granulation offers the advantage of reduced fines and dust production, as well as the saving of a subsequent drying step [12]. This is especially valuable in terms of sustainability. In addition, this method enables the handling of moisture-sensitive materials [13] or the processing of high-dose formulations with poor compactibility [14]. To date, continuous melt granulation is typically executed on twin-screw machines [15,16,17]. However, the use of a planetary roller granulator, which belongs to the class of multiple screw machines which have already been applied for wet granulation [18,19,20], is a promising alternative [21] due to the basic process concept (Figure 1).
A firmly locked roller cylinder surrounds a motor-driven central spindle. Both parts have a 45° toothing and individual heating systems to regulate the energy input over the heated surface. The implemented planetary spindles determine the free processing volume. Within the processing section these are free-flowing and have to be distributed symmetrically across the radial cross-sectional area due to stability reasons. The minimum number is three in any case, while the maximum number depends on the apparatus size. During processing, the movement of the central spindle is superimposed onto the planetary spindles due to the toothing matching the other main parts. The resulting overall direction of circulation of the planetary spindles is concurrent to the central spindle, while its self-rotation is counter current. This radial movement pattern, in addition to the axial conveying related to the 45° toothing, defines the material transport characteristics during processing. The enhanced ratio of heated surface to free processing volume is beneficial in terms of temperature control. This is relevant for melt granulation, which is based on the softening or melting of the utilized binder component and the processing of temperature-sensitive materials.
During planetary roller melt granulation (PRMG), the product temperature and overall thermal stress correspond to the energy input, which is determined by the applied specific mechanical energy, the specific thermal energy, and the exposure time to the corresponding mechanisms. These can be characterized by the residence time distribution (RTD), an indirect process parameter frequently applied to evaluate and control various methods of pharmaceutical technology covering, e.g., mixing/blending [22,23] or tableting [24,25]. This, of course, also includes wet [26,27] or melt granulation [28] executed on twin-screw machines.
The characteristics of an RTD (Figure 2) provide information on the sub-processes occurring during melt granulation, including the melting of the binder and the granulation mechanisms [29]. The on-set is the point in time of the first discharged material after entering the process at t = 0. Therefore, the time interval in between is the minimum duration of the material in the process. Thus, a temperature above the binder’s melting point has to be obtained in order to start the granulation. The off-set symbolizes the latest discharge of material after entering the process at t = 0. Therefore, the time interval in between is the maximum duration of applied shear and thermal stress to the material. This has to be considered in the context of breakage mechanisms during granulation [29] or thermal material degradation. The width as the time span between the on-set and off-set of the distribution reflects the overall back-mixing of the material during processing. This is crucial for the fundamental particle growth mechanisms [29]. Furthermore, the kinetic between peak and off-set is usually referred to as material “wash-out” and prolonged in comparison to the “loading” kinetic between on-set and peak. The median is representative for the overall RTD and average transport velocity.
The execution of a pulse experiment is the main procedure for determining the RTD in the form of the density function (E) [30,31,32]. Therefore, a tracer is added to the material stream at the inlet in order to measure the response signal in the outlet stream. This density function represents the probability of an individual entity exiting the process at a specific time (t). Via a step signal experiment, which utilizes a shift from a standard formulation to a tracer formulation or vice versa, the cumulative function (F) can be determined. This has a limit value of one (Equation (1)) as it is the sum function of the density distribution. In any case, a representative parameter for the RTD is the median (t50) representing the average transportation time (Equation (2)).
F   t = 0 t E ( t ) d t = 1
F   t 50 = 0 t 50 E ( t ) d t = 0.5 .
The aim of this study is to identify the material transport and mixing characteristics during PRMG. Therefore, the RTD is applied as a key indicator for the processing conditions. In this respect, an on-line residence time determination method is evaluated in terms of repeatability as well as the suitability of a model approach from the literature to capture the experimental data. Furthermore, the impact of the process parameter variations on the machine fill level is discussed. These insights serve as basis for the characterization of the blending characteristics during processing, which refers to axial back-mixing or radial material distribution. Thus, the identified correlations between model parameters and processing conditions are highlighted.

2. Materials and Methods

2.1. Melt Granulation

Melt granulation was executed on a lab-scale planetary roller granulator (PWE 30, Entex Rust & Mitschke GmbH, Bochum, Germany). The processing section of the machine consisted of a single module unit with a total length of 222 mm and an inner diameter of 35.4 mm of the roller cylinder. This was equipped with five standard planetary spindles with a diameter of 8.8 mm and a central spindle with a diameter of 17.7 mm. The length of these was equal to the processing section. Therefore, the overall free processing section volume ( V f r e e )  was 96 mL, while the ratio of heated surface to free process volume was 0.58 m−1. An open discharge nozzle with a large diameter for pressure-free release of the granules (Figure 1) was used at the outlet. A gravimetric feeder (DDW-M-DS(R) 28, Brabender GmbH & Co. KG, Duisburg, Germany) was used for dosing of the model formulation, which was blended before processing with a wheel mixer (RRM ELTE 650, J. Engelsmann AG, Lugwigshafen, Germany) for 10 min. Lactose monohydrate (Lactose 310, Foremost Farms USA, Baraboo, WI, USA) served as model compound (90 wt%) and hydroxypropylcellulose (Klucel EXF Pharm, Ashland Inc., Covington, GA, USA) was utilized as melt binder (10 wt%). The set temperature for both heat systems was 150 °C for all experiments. Internal heat exchangers were used with water as tempering media with a volume flow for the central spindle of 7 L min−1 and 17 L min−1 for the roller cylinder heating. Product temperature was determined on-line with an IR-camera (Test 875, Testo SE & Co. KGaA, Titisee-Neustadt, Germany).

2.2. Experimental Design Space

A design space for the execution of the residence time experiments was determined a priori. The set points within the experimental study cover four equidistant levels, each for the investigated process parameter, including the feed rate ( m ˙ ) from 0.3 to 1.2  k g   h 1  and rotation speed ( n ) from 60 to 240  m i n 1 . In combination, these direct process parameters define the specific feed load ( S F L ), which is a typical characteristic parameter for twin-screw processes. This approach is adapted for planetary roller processing (Equation (3)) and considers the material volume flow, defined by the ratio of feed rate to untapped bulk density ( ρ b u l k ), in relation to the potential transport volume flow per one revolution of the central spindle. This correlates to the rotation speed, the pitch of the central spindle ( h c s ), and the free cross-sectional area ( A f r e e ), which in combination with the length of the central spindle ( L c s ) determines the free processing volume ( V f r e e ).
S F L = m ˙ n   ρ b u l k   h c s   A f r e e = m ˙ n   ρ b u l k   h c s L c s   V f r e e .

2.3. Residence Time Characterization

Determination of the residence time distribution density function at different process parameter sets was executed with an on-line camera system (ExtruVis, MeltPrep GmbH, Graz, Austria) via pulse signal experiments (Figure 3), which have already been applied successfully to characterize the RTD in twin-screw processes [33,34]. During the steady stage of the process, E104 or Ponceau 4R was added as a tracer to the feed formulation through the hopper, with a ratio of feed rate to tracer mass at 20.83 per second. The corresponding signal was measured at the discharge blend with a frequency of 36 to 51 Hz. To characterize the residence time distribution, a pre-data treatment of the experimental data was executed to derive the E-function. This procedure included a baseline correction as well as a normalization by the overall integral value, which was calculated via the trapezoidal rule. The two-compartment model [35,36] was fitted to the experimental data with the least square method. This model (Equation (4)) considers a material plug flow for the axial transport superimposed by the mixing function of a continuously stirred tank reactor.
E t = k 0 c 0   0.5   e 0.5   t m i x 3 σ 2 t m i x 1 t t d s   e r f c t m i x 1 σ 2 t t d s 2 σ .
The dead time ( t d s ) of the process is the transport time based on the dead process volume in relation to the material flow, whereby the standard deviation ( σ ) characterizes the axial material distribution of the transport function. For the mixing process, the ratio of the mixing volume to the material volume flow describes the mixing duration ( t m i x ) within the reactor and therefore reflects the time extension of the corresponding mechanisms. Furthermore, the initial tracer concentration signal ( c 0 ) is the ratio of the scaling parameter ( k 0 ) to the mixing parameter.

3. Results and Discussion

3.1. Determination and Characterization of the Residence Time

The repeatability concerning the residence time determination was tested at four different set points. These cover the combination of the minimum and maximum levels for each of the adjusted direct process parameters. In steady-state processing, signal impulse experiments were executed in triplicate. The various symbols represent the different runs for each parameter set (Figure 4). The residence time distribution (RTD) is expressed by the density function. At the same time, the suitability of the two-compartment model for the application in planetary roller melt granulation (PRMG) is evaluated by this experimental study. The model fit to the data is displayed in solid lines.
In each case, the characteristic of the RTDs is similar, including a single on-set, peak, and off-set, as well as a prolonged kinetic between peak and off-set in comparison to the kinetic between on-set and peak. This profile matches the common RTD profiles in twin-screw granulation (Figure 2). Furthermore, a sufficient agreement between model fit and experimental data is found for all runs at each fixed parameter set. Therefore, the two-compartment model is classified as suitable for representing the RTD in PRMG.
At the same time, the displayed data also imply a sufficient agreement among the individual repetitions of the residence time determination at a fixed parameter set in a steady state. The visual impression is thereby supported by the fitted model parameters. For the three repetitions, the coefficient of variation (CV) of the individual factors as a ratio of standard deviation to mean value is generally in the range of 10% or less (Table 1). This is sufficient and highlights the repeatability of the RTD determination method.
For the parameter set combining the highest rotation speed at 240  m i n 1  and the lowest mass flow at 0.30  k g   h 1 , the density of the material flow at the outlet is minimized. This can cause some irregularity in the form of pulsation, as in run 1 (Figure 4, bottom right). The corresponding impact on the RTD measurement and determined profile is limited to the “outwash” kinetic of the RTD. This is the part of the function related to the model mixing parameter ( t m i x ). Due to the coupling for the applied RTD model (Equation (4)),  k 0  is also affected. Overall, this is not considered critical, since the described effect has only been observed as an exception for this parameter set in all experimental investigations. However, the potential fluctuation of the model parameters  t m i x  and k0 for this single parameter set has to be considered during the evaluation of the transport characteristics in PRMG.

3.2. Process Parameter Impact on the RTD in PRMG

The determined RTD profiles for the 16 parameter sets within the design space of the experimental study are given in Figure 5. The data are staggered for the clarity of the depiction in four categories based on a constant rotation speed. The connected model fit is represented by the lines, which are overall in sufficient agreement with the experiments.
For a constant rotation speed, an increase of the feed rate (encoded by the symbol grey scale) generally results in a shift of the peak as well as of the  t 50  (Table 2) to shorter exit times, while the absolute width of the RTD decreases. Therefore, the peak height increases as well. In comparison, increasing the rotation speed (encoded by the symbols) for a constant feed rate has the same effect on the RTD profile. These findings are similar to those in twin-screw processes and the impact of direct process parameter variations on the RTD [37]. Furthermore, the other determined characteristic model parameters (Equation (4)) match the visual findings concerning the process parameter impact.
First of all, the  c 0  values in the range of one indicate a successful normalization of the experimental data to obtain the fundamental density function. The minor deviations are related to the fit of a model, assuming an ideal RTD profile to an experimental data set underlying a method-specific signal-to-noise ratio. In this respect,  k 0  is nearly the reciprocal value of the mixing time.
The dead time represents an exit time for a plug flow in model theory [35], which is not affected by any superimposed material blending or mixing.
Therefore,  t d s  is connected to the part of the RTD profile before the peak, including the on-set, and decreases for larger feed rates and rotation speeds. The parameter σ is subjected to the same impact, since this represents the absolute width of the exit time distribution of the plug flow function. Instead, the correlations for the mixing time are ambiguous. For a constant rotation speed, the values of  t m i x  decrease with increasing throughput, while an alternating impact is found for an increase of the rotation speed at a constant feed rate. This implies a mixing volume variability, which could originate from an affected material transition pattern between the planetary and central spindle.

3.3. Machine Fill Level during PRMG

A key parameter influencing the material transport and mixing during PRMG is the fill level of the granulator during processing. This is defined as the ratio of the processing volume, where the material is actually flowing through to the free processing volume in the roller cylinder, which is not occupied by the central or planetary spindles. For a maximum fill level, the hydrodynamic residence time ( τ ) represents the average transport due to plug flow. This parameter (Equation (5)) correlates to the granulator load.
τ = V f r e e   V ˙   = V f r e e   m ˙ ρ b u l k     = 1   n   h c s L c s   S F L     .
For the overall PRMG process, the median of the RTD represents the average axial transport time based on the covered free processing volume within the machine by the material stream. Therefore, the ratio of  t 50  to  τ  can be applied as a surrogate for the granulator fill level (Figure 6).
In general, the surrogate for the fill level is constantly below 100% (Table 2) and increases for a higher SFL, and the corresponding correlation fits a power model approach. Therefore, the impact of the investigated process parameter variations and the fill level adjustment in steady state depends on the balance of the acting forces related to the material transport in PRMG.
On the one hand, this includes the applied specific mechanical energy, which is transmitted at the contact surface of planetary spindles and central spindle, respectively, and the roller cylinder. Due to the 45° toothing and the orbital motion of the planetary spindles, the transferred force ( F t r a n s f e r ) is converted into axial drive and radial distribution of the bulk material. Thus, the material is sheared by the moving mechanical parts. Consequently, the overall transferred forces rely on the shear rate depending on the rotation speed, viscosity ( η ), and the covered contact surface. This must directly correlate to the covered cross-sectional area by the material stream ( A f r e e , c o v ) , which represents the granulator fill level. On the other hand, the axial drive of the material resulting in a transport velocity is limited due to friction and dissipation. In fluid dynamics, this is usually referred to as pressure loss. The corresponding force ( F f r i c t i o n ) depends on the feed rate ratio in the power of two to the bulk density and the covered cross-sectional area by the material stream, which represents the granulator fill level.
The overall force balance for the material transport can be summarized (Equation (6)). However, the material-related parameters can be neglected for the conducted experimental investigations, as the material formulation was constant.
F t r a n s f e r r e d n ,   η ,   A f r e e , c o v = F f r i c t i o n m ˙ 2 ,   ρ b u l k 1 ,   A f r e e , c o v 1 .
For a constant rotation speed, a higher SFL refers to a higher feed rate. The corresponding amplifying impact on the friction force is balanced by increasing the flow through the cross-sectional area. This reduces the flow velocity and, thereby, the friction force and enhances the transferred energy. For a constant feed rate, a higher SFL refers to a lower rotation speed. The corresponding attenuating impact on the transferred force is balanced by increasing the flow through the cross-sectional area. This reduces the friction force, even while the throughput is maintained.
However, the efficiency of the mechanical energy transfer decreases with higher fill levels since the material is initially distributed in the contact area of the mechanical parts. Later on, the material fills up the void spaces between the planetary spindles, which barely contribute to the energy input. Furthermore, based on the determined particle size distributions of the granules in one study [21], a higher fill level contributes to the mechanisms of nucleation and agglomeration [29] and leads to an enlarged particle size.

3.4. Material Transport and Mixing Characteristics in PRMG

With respect to the applied RTD model, the axial material transport is characterized by  t d s , while axial back-mixing is related to  σ  and radial material distribution by  t m i x . Due to the machine concept, including the 45° toothing of the mechanical parts as well as the orbital motion of the planetary spindles around the central spindle, all mentioned parameters and the corresponding mechanisms are coupled. Moreover, these are related concerning the model theory to the granulator fill level, which correlates to the SFL (Figure 6). In this context, a transformation of the time parameters (Equation (5)) is considered to display the data (Figure 7).
In general, all transformed parameters correlate linearly to the SFL, which justifies the applied parameter transformation. The values for σ are the highest at each set point ((σ n)−1 = 2.21 SFL + 0.0104), as the width of the axial material transport is absolutely lower than the duration of the material transport (Table 2). Furthermore, transport and mixing time are in the same order of magnitude ((tds n)−1 = 0.711 SFL + 0.00196, (tmix n)−1 = 0.530 SFL + 0.00175), since the corresponding mechanisms refer to the same fill volume inside the granulator. The higher values for  t m i x  result from the superimposition of the orbital planetary spindle movement and their self-rotation on the radial distribution by the part toothing. This leads to an apparent increase in the mixing volume in comparison, due to a circulation within and exchange in between distinct material portions continuously during processing along the granulator.
Furthermore, the impact of the process parameter adjustments on the transformed model parameters is directly related to the granulator fill level. This is also manipulated by a variation of the feed rate or the rotation speed (Figure 7), while the other one is constant. However, the relative factor a parameter adaption leading to a fill level increase is exceeding the relative factor of the corresponding fill volume increase. Therefore, the transformed surrogates decrease while, in turn, the model parameters representing a duration as ratio of a volume to feed rate increase. This effect is due to the balance regarding the transferred shear forces by the mechanical parts to the material and the friction force limiting the material discharge (Equation (6)). Here, the fill level acts as an equalizer having a multiple influence on the acting forces. Furthermore, the manipulation of the mechanical energy transferred into the material is less effective with an increasing load, as the material in the void space is barely taking part in the interaction in the contact area of the mechanical parts. At the same time, the mixing volume is related to the fill volume as well as the number of material circulations. These are driven by the self-rotation of the planetary spindles and the transition patterns between planetary and central spindles. Therefore, a reduction in the rotation speed has an increasing effect on the fill level, while the reduced circulation of the fill volume has a decreasing effect on the mixing volume. This adds to the limited effectiveness of the parameter adaption to manipulate the fill level.

4. Conclusions

The residence time distribution is a key indicator for the processing conditions during planetary roller melt granulation, as this indirect process parameter reflects the material transport. This especially includes the actual fill level, which represents the covered fraction of the free processing volume by the material during processing. The free volume in a single module is determined by the type and number of planetary spindles considered in a configuration. During processing, the fill level itself is adjustable by the direct process parameters, which simultaneously affect the axial material transport and the radial material distribution. These, in turn, determine the transport velocity of the material and the mixing time. Both aspects contribute to the granulation performance during planetary roller melt granulation. However, in the common setup, which was also used for these experimental investigations, the pre-mix of binder and model component resulted in a superimposition of granulation and melting during processing. Therefore, the residence time distribution reflects both mechanisms and cannot be directly correlated to the product particle size distribution.

Author Contributions

Conceptualization, M.T. and J.B.; investigation, T.L. and J.B.; methodology, T.L. and J.B.; resources, A.B.; supervision, J.B.; visualization, J.B.; writing—original draft, J.B.; writing—review & editing, A.B. and M.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors would like to thank Merck KgaA (Darmstadt, Germany) for providing lactose as a model compound and Ashland Industries Deutschland GmbH (Düsseldorf, Germany) for providing the binder used in the melt granulation, as well as Entex Rust & Mitschke GmbH (Bochum, Germany) for the granulator availability and affiliated equipment.

Conflicts of Interest

The authors declare no conflict of interest. The company had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Nomenclature

A f r e e Free area of the planetary roller processing cross-section[m2]
A f r e e , c o v Free area of the planetary roller processing cross-section covered with material during processing[m2]
c 0 Signal strength of initial tracer concentration[*]
CVCoefficient of variation[-]
E Residence time density function[s−1]
ηDynamic viscosity[Pa s]
F Cumulative residence time function[-]
F f r i c t i o n   Friction force[N]
F t r a n s f e r r e d Transferred force via rotating parts [N]
h c s Pitch of central spindle[mm]
k 0 Scaling parameter[s−1]
l c s Length of central spindle[mm]
m ˙ Feed rate of material pre-mix[kg h−1]
nRotation speed of central spindle [min−1]
ρ b u l k Untapped bulk density [kg m−3]
SFLSpecific feed load[-]
σTransport function width[s]
tTime[s]
t 50 Median of residence time distribution[s]
t d e a d Process dead time[s]
t m i x Process mixing time[s]
τHydrodynamic mean residence time[s]
V ˙ Process volume flow[m3 h−1]
V f r e e Free volume of the planetary roller processing section[m3]

References

  1. European Medicines Agency. ICH Q8 (R2) Pharmaceutical Development—Scientific Guideline. Available online: https://www.ema.europa.eu/en/ich-q8-r2-pharmaceutical-development-scientific-guideline (accessed on 19 July 2023).
  2. Kleinebudde, P.; Khinast, J.; Rantanen, J. Continuous Manufacturing of Pharmaceuticals, 1st ed.; Wiley: Hoboken, NJ, USA, 2017. [Google Scholar]
  3. Burcham, C.L.; Florence, A.J.; Johnson, M.D. Continuous manufacturing in pharmaceutical process development and manufacturing. Annu. Rev. Chem. Biomol. Eng. 2018, 9, 253–281. [Google Scholar] [CrossRef]
  4. Schlindwein, W.S.; Gibson, M. Pharmaceutical Quality by Design: A Practical Approach, 1st ed.; Wiley: Hoboken, NJ, USA, 2018. [Google Scholar]
  5. Jelsch, M.; Roggo, Y.; Kleinebudde, P.; Krumme, M. Model predictive control in pharmaceutical continuous manufacturing: A review from a user’s perspective. Eur. J. Pharm. Biopharm. 2021, 159, 137–142. [Google Scholar]
  6. Serno, P.; Kleinebudde, P.; Knop, K. Granulieren: Grundlagen, Verfahren, Formulierungen; ECV—Editio-Cantor-Verlag: Aulendorf, Germany, 2007. [Google Scholar]
  7. Salman, A.D.; Hounslow, M.J.; Seville, J.P.K. Granulation, 1st ed.; Elsevier: Amsterdam, The Netherlands; Boston, MA, USA, 2007. [Google Scholar]
  8. Kleinebudde, P. Roll compaction/dry granulation: Pharmaceutical applications. Eur. J. Pharm. Biopharm. 2004, 58, 317–326. [Google Scholar] [CrossRef]
  9. Sun, C.C.; Kleinebudde, P. Mini review: Mechanisms to the loss of tabletability by dry granulation. Eur. J. Pharm. Biopharm. 2016, 106, 9–14. [Google Scholar]
  10. Faure, A.; York, P.; Rowe, R.C. Process control and scale-up of pharmaceutical wet granulation processes: A review. Eur. J. Pharm. Biopharm. 2001, 52, 269–277. [Google Scholar]
  11. Liu, B.B.; Wang, J.M.; Zeng, J.; Zhao, L.J.; Wang, Y.J.; Feng, Y.; Du, R.F. A review of high shear wet granulation for better process understanding, control and product development. Powder Technol. 2021, 381, 204–223. [Google Scholar] [CrossRef]
  12. Parikh, D.M. Handbook of Pharmaceutical Granulation Technology, 4th ed.; CRC Press: Boca Raton, FL, USA, 2021. [Google Scholar]
  13. Weatherley, S.; Mu, B.O.; Thompson, M.R.; Sheskey, P.J.; O′Donnell, K.P. Hot-Melt Granulation in a Twin Screw Extruder: Effects of Processing on Formulations with Caffeine and Ibuprofen. J. Pharm. Sci. 2013, 102, 4330–4336. [Google Scholar] [PubMed]
  14. Lakshman, J.P.; Kowalski, J.; Vasanthavada, M.; Tong, W.Q.; Joshi, Y.M.; Serajuddin, A.T.M. Application of Melt Granulation Technology to Enhance Tabletting Properties of Poorly Compactible High-Dose Drugs. J. Pharm. Sci. 2011, 100, 1553–1565. [Google Scholar] [CrossRef] [PubMed]
  15. Van Melkebeke, B.; Vermeulen, B.; Vervaet, C.; Remon, J.P. Melt granulation using a twin-screw extruder: A case study. Int. J. Pharm. 2006, 326, 89–93. [Google Scholar] [CrossRef] [PubMed]
  16. Kittikunakorn, N.; Liu, T.Z.; Zhang, F. Twin-screw melt granulation: Current progress and challenges. Int. J. Pharm. 2020, 588, 119670. [Google Scholar] [CrossRef]
  17. Liu, T.Z.; Kittikunakorn, N.; Zhang, Y.; Zhang, F. Mechanisms of twin screw melt granulation. J. Drug Deliv. Sci. Technol. 2021, 61, 102150. [Google Scholar]
  18. Schroeder, R.; Steffens, K.-J. A new system for continuous wet granulation. Die Pharm. Ind. 2002, 64, 283–288. [Google Scholar]
  19. Schroeder, R.; Steffens, K.-J. Granulation of hydrophobic actives with the planetary roller extruder. Pharm. Ind. 2003, 65, 367–372. [Google Scholar]
  20. Schroeder, R.; Steffens, K.-J. Continuous granulation technologies. In Handbook of Pharmaceutical Granulation Technology, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2005; pp. 431–457. [Google Scholar]
  21. Nesges, D.; Lang, T.; Birr, T.; Thommes, M.; Bartsch, J. Planetary roller melt granulation (PRMG)—A new continuous method for powder processing. Powder Technol. 2023, 427, 118728. [Google Scholar] [CrossRef]
  22. Peterwitz, M.; Gerling, S.; Schembecker, G. Challenges in tracing material flow passing a loss-in-weight feeder in continuous manufacturing processes. Int. J. Pharm. 2022, 612, 121304. [Google Scholar]
  23. Tian, H.Y.; Bhalode, P.; Razavi, S.M.; Koolivand, A.; Muzzio, F.J.; Ierapetritou, M.G. Characterization and propagation of RTD uncertainty for continuous powder blending processes. Int. J. Pharm. 2022, 628, 122326. [Google Scholar] [CrossRef] [PubMed]
  24. Zimmermann, M.; Thommes, M. Residence time and mixing capacity of a rotary tablet press feed frame. Drug Dev. Ind. Pharm. 2021, 47, 790–798. [Google Scholar] [PubMed]
  25. Forgber, T.; Rehrl, J.; Matic, M.; Sibanc, R.; Sivanesapillai, R.; Khinast, J.G. Experimental and numerical investigations of the RTD in a GEA ConsiGma CTL25 tablet press. Powder Technol. 2022, 405, 117507. [Google Scholar] [CrossRef]
  26. Kumar, A.; Alakarjula, M.; Vanhoorne, V.; Toiviainen, M.; De Leersnyder, F.; Vercruysse, J.; Juuti, M.; Ketolainen, J.; Vervaet, C.; Remon, J.P.; et al. Linking granulation performance with residence time and granulation liquid distributions in twin-screw granulation: An experimental investigation. Eur. J. Pharm. Sci. 2016, 90, 25–37. [Google Scholar]
  27. Pauli, V.; Kleinebudde, P.; Krumme, M. From powder to tablets: Investigation of residence time distributions in a continuous manufacturing process train as basis for continuous process verification. Eur. J. Pharm. Biopharm. 2020, 153, 200–210. [Google Scholar]
  28. Mu, B.; Thompson, M.R. Examining the mechanics of granulation with a hot melt binder in a twin-screw extruder. Chem. Eng. Sci. 2012, 81, 46–56. [Google Scholar]
  29. Iveson, S.M.; Litster, J.D.; Hapgood, K.; Ennis, B.J. Nucleation, growth and breakage phenomena in agitated wet granulation processes: A review. Powder Technol. 2001, 117, 3–39. [Google Scholar]
  30. Wesholowski, J.; Berghaus, A.; Thommes, M. Inline Determination of Residence Time Distribution in Hot-Melt-Extrusion. Pharmaceutics 2018, 10, 49. [Google Scholar]
  31. Rahimi, S.K.; Paul, S.; Sun, C.C.; Zhang, F. The role of the screw profile on granular structure and mixing efficiency of a high-dose hydrophobic drug formulation during twin screw wet granulation. Int. J. Pharm. 2020, 575, 118958. [Google Scholar]
  32. Wesholowski, J.; Hoppe, K.; Nickel, K.; Muehlenfeld, C.; Thommes, M. Scale-Up of pharmaceutical Hot-Melt-Extrusion: Process optimization and transfer. Eur. J. Pharm. Biopharm. 2019, 142, 396–404. [Google Scholar] [PubMed]
  33. Mundozah, A.L.; Cartwright, J.J.; Tridon, C.C.; Hounslow, M.J.; Salman, A.D. Hydrophobic/hydrophilic powders: Practical implications of screw element type on the reduction of fines in twin screw granulation. Powder Technol. 2019, 341, 94–103. [Google Scholar]
  34. Meier, R.; Thommes, M.; Rasenack, N.; Moll, K.P.; Krumme, M.; Kleinebudde, P. Granule size distributions after twin-screw granulation—Do not forget the feeding systems. Eur. J. Pharm. Biopharm. 2016, 106, 59–69. [Google Scholar]
  35. Reitz, E.; Podhaisky, H.; Ely, D.; Thommes, M. Residence time modeling of hot melt extrusion processes. Eur. J. Pharm. Biopharm. 2013, 85, 1200–1205. [Google Scholar]
  36. Wesholowski, J.; Podhaisky, H.; Thommes, M. Comparison of residence time models for pharmaceutical twin-screw-extrusion processes. Powder Technol. 2019, 341, 85–93. [Google Scholar]
  37. Wesholowski, J.; Berghaus, A.; Thommes, M. Investigations Concerning the Residence Time Distribution of Twin-Screw-Extrusion Processes as Indicator for Inherent Mixing. Pharmaceutics 2018, 10, 207. [Google Scholar]
Figure 1. Radial cross-section of planetary roller processing section, including central spindle, planetary spindle, and roller cylinder. The heating of the roller cylinder and central spindle is independent. The central spindle, in addition to the number of planetary spindles, determines the free radial cross-section and free processing volume in axial direction.
Figure 1. Radial cross-section of planetary roller processing section, including central spindle, planetary spindle, and roller cylinder. The heating of the roller cylinder and central spindle is independent. The central spindle, in addition to the number of planetary spindles, determines the free radial cross-section and free processing volume in axial direction.
Pharmaceutics 15 02039 g001
Figure 2. Typical profile of the residence time density function and corresponding cumulative form, including characteristic parameters.
Figure 2. Typical profile of the residence time density function and corresponding cumulative form, including characteristic parameters.
Pharmaceutics 15 02039 g002
Figure 3. Setup for non-destructive on-line determination of residence time distribution in planetary roller melt granulation experiments via ExtruVis system. Camera sensor placed in front of granulator discharge blend.
Figure 3. Setup for non-destructive on-line determination of residence time distribution in planetary roller melt granulation experiments via ExtruVis system. Camera sensor placed in front of granulator discharge blend.
Pharmaceutics 15 02039 g003
Figure 4. Experimentally determined residence time distribution in steady state at four fixed parameter sets, representing the investigated design space for the investigation of the transport characteristics in planetary roller melt granulation. For the clarity of visual data representation, only every 125th data point is plotted for each data set.
Figure 4. Experimentally determined residence time distribution in steady state at four fixed parameter sets, representing the investigated design space for the investigation of the transport characteristics in planetary roller melt granulation. For the clarity of visual data representation, only every 125th data point is plotted for each data set.
Pharmaceutics 15 02039 g004
Figure 5. Residence time distribution and model fit to the experimental data for each parameter set of the investigated design space. Granulator configuration, set temperature (central spindle, roller cylinder), and equipment were fixed for all experiments. The symbols equal the applied rotation speed; colour represents applied feed rates. For the clarity of visual data representation, only every 125th data point is plotted for each data set.
Figure 5. Residence time distribution and model fit to the experimental data for each parameter set of the investigated design space. Granulator configuration, set temperature (central spindle, roller cylinder), and equipment were fixed for all experiments. The symbols equal the applied rotation speed; colour represents applied feed rates. For the clarity of visual data representation, only every 125th data point is plotted for each data set.
Pharmaceutics 15 02039 g005
Figure 6. Ratio of RTD median to hydrodynamic residence time as surrogate for the granulator fill level parameter for each parameter set of the investigated design space (Figure 3). Granulator configuration, set temperature (central spindle, roller cylinder), and equipment were fixed for all experiments. The symbols equal the applied rotation speed, colour represents applied feed rates. Line represents the fit of a power model approach to the overall displayed data.
Figure 6. Ratio of RTD median to hydrodynamic residence time as surrogate for the granulator fill level parameter for each parameter set of the investigated design space (Figure 3). Granulator configuration, set temperature (central spindle, roller cylinder), and equipment were fixed for all experiments. The symbols equal the applied rotation speed, colour represents applied feed rates. Line represents the fit of a power model approach to the overall displayed data.
Pharmaceutics 15 02039 g006
Figure 7. Fitted and transformed transport and mixing-related parameters of the applied two-compartment model to the experimental data for each parameter set of the investigated design space. Granulator configuration, set temperature (central spindle, roller cylinder), and equipment were fixed for all experiments. The symbols equal the model parameter, colour represents applied feed rates. Rotation speeds decrease for higher SFL at a constant feed rate. Lines represent the fit of a linear model approach to the overall displayed data.
Figure 7. Fitted and transformed transport and mixing-related parameters of the applied two-compartment model to the experimental data for each parameter set of the investigated design space. Granulator configuration, set temperature (central spindle, roller cylinder), and equipment were fixed for all experiments. The symbols equal the model parameter, colour represents applied feed rates. Rotation speeds decrease for higher SFL at a constant feed rate. Lines represent the fit of a linear model approach to the overall displayed data.
Pharmaceutics 15 02039 g007
Table 1. Coefficient of variation (CV) of the determined RTD model parameters for the three repetitions at the different process parameter set points defined by the set rotation speed and feed rate.
Table 1. Coefficient of variation (CV) of the determined RTD model parameters for the three repetitions at the different process parameter set points defined by the set rotation speed and feed rate.
nset
[min−1]
m ˙ set
[kg h−1]
C V k 0
[%]
C V t d s
[%]
C V σ
[%]
C V t m i x
[%]
600.311.13.287.8310.1
601.22.591.343.692.44
2400.347.98.8111.041.3
2401.210.43.4611.66.49
Table 2. Parameter matrix for the determined RTD model fits to the experimental data.
Table 2. Parameter matrix for the determined RTD model fits to the experimental data.
nset
[min−1]
m ˙ set
[kg h−1]
SFL
[−]
t 50
[s]
c 0
[−]
k o
[−]
t d s
[s]
σ
[s]
t m i x
[s]
t 50 τ 1
[%]
600.30.00595269.91.010.00805172.449.6124.840.4
600.60.00300165.21.060.0097791.425.6108.434.6
600.90.00200129.71.050.0128270.419.081.929.8
601.20.00149107.71.040.0139653.917.874.227.1
1200.30.00191231.01.030.00937149.535.0110.449.4
1200.60.00601151.20.980.0104679.123.594.245.2
1200.90.00408118.81.020.0109954.114.492.846.1
1201.20.0029992.11.000.0162046.611.461.938.1
1800.30.01791199.30.970.00774109.417.8125.558.2
1800.60.00894154.20.960.0084971.719.9113.553.3
1800.90.00591107.41.000.0133052.715.075.648.2
1801.20.0044786.10.970.0161241.410.460.344.3
2400.30.02377181.31.080.0090294.524.0120.364.5
2400.60.01201127.31.000.0118965.217.784.555.1
2400.90.0078398.61.000.0136246.212.673.551.5
2401.20.0059783.50.980.0150036.69.865.650.0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lang, T.; Bramböck, A.; Thommes, M.; Bartsch, J. Material Transport Characteristics in Planetary Roller Melt Granulation. Pharmaceutics 2023, 15, 2039. https://doi.org/10.3390/pharmaceutics15082039

AMA Style

Lang T, Bramböck A, Thommes M, Bartsch J. Material Transport Characteristics in Planetary Roller Melt Granulation. Pharmaceutics. 2023; 15(8):2039. https://doi.org/10.3390/pharmaceutics15082039

Chicago/Turabian Style

Lang, Tom, Andreas Bramböck, Markus Thommes, and Jens Bartsch. 2023. "Material Transport Characteristics in Planetary Roller Melt Granulation" Pharmaceutics 15, no. 8: 2039. https://doi.org/10.3390/pharmaceutics15082039

APA Style

Lang, T., Bramböck, A., Thommes, M., & Bartsch, J. (2023). Material Transport Characteristics in Planetary Roller Melt Granulation. Pharmaceutics, 15(8), 2039. https://doi.org/10.3390/pharmaceutics15082039

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop