Numerical and experimental validation of a detailed non-isothermal CFD-DEM model of a pilot-scale Wurster coater

Improvement of a pharmaceutical powder mixing process in a tote blender via DEM simulations

An industrial-scale pharmaceutical powder blending process was studied via discrete element method (DEM) simulations. A DEM model of two active pharmaceutical ingredient (API) components and a combined excipient component was calibrated by matching the simulated response in a dynamic angle of repose tester to the experimentally observed response. A simulation of the 25-minute bin blending process predicted inhomogeneous API distributions along the rotation axis of the blending container. These concentration differences were confirmed experimentally in a production-scale mixing trial using high-performance liquid chromatography analysis of samples from various locations in the bin. Several strategies to improve the blend homogeneity were then studied using DEM simulations. Reversing the direction of rotation of the blender every minute was found to negligibly improve the blending performance. Introducing a baffle into the lid at a 45° angle to the rotation axis sped up the axial mixing and resulted in a better final blend uniformity. Alternatively, rotating the blending container 90° around the vertical axis five minutes prior to the process end was predicted to reduce axial segregation tendencies.

Modeling the coating layer thickness in a pharmaceutical coating process

Although mechanistic numerical simulations can offer great insights into a process, they are limited with respect to resolved process time. While statistical models provide long-term predictability, determining the underlying probability distributions is often challenging. In this work, detailed CFD-DEM simulations of a pharmaceutical Wurster coating process for microspheres are used to evaluate the input parameters for a novel Monte-Carlo simulation approach. The combined strengths of both modeling approaches make it possible to predict the coating mass and thickness distributions over the entire process time. It was observed that smaller beads receive a thicker coating layer since they pass the spray zone closer to the nozzle. Moreover, it was established that, in contrast to the airflow rate, the spray rate has a great impact on the inter-particle coating variability. A stochastic model was developed to investigate the relative contribution of coating layer variability and fill weight variability to the product non-uniformity in a capsule filling process of Multiple Unit Pellet Systems (MUPS).