Display omitted
•The coarse-grain model was applied to investigate mixing state in a spouted bed.•Mixing state well agreed between the original and coarse grain particle systems.•The coarse-grain ...model makes it possible to reduce the calculation cost drastically.•Adequacy of the coarse grain model was proven in solid mixing in the spouted bed.
We describe the effectiveness of the coarse-grain model of the discrete element method (DEM) in solid mixing in a spouted bed. The coarse-grain model was developed originally to simulate large-scale DEM simulations efficiently, where the coarse-grain particle represents a group of original particles. In previous studies, the adequacy of the coarse-grain model was proven in a pneumatic conveying system and bubbling fluidized beds through verification and validation tests. Namely, In this study, the coarse-grain model is applied to solid mixing in a spouted bed. Agreement of the mixing state is shown between an original particle system and coarse grain model systems. Subsequently, correlation between solid mixing and macroscopic behavior of the solid particles is examined in the present study. The macroscopic properties such as solid-particle spatial distribution, pressure drop and velocity distribution of the solid phase are confirmed to correspond quantitatively in the coarse-grain model and an original particle system. Finally, the calculation efficiency of the coarse-grain model is evaluated. Consequently, the coarse-grain model is shown to be able to efficiently investigate the solid mixing in a spouted bed.
Drug delivery from a capsule-based dry powder inhaler depends on the inhaler's design, the drug's formulation, and the inhalation maneuver. The latter affects both the air flow and the capsule motion ...in the inhaler. It is well known that patient-to-patient variability is a major challenge in the design of new inhaler types. Modeling and simulation are important tools for understanding such systems, yet quite complex. Simulation studies of capsule-based dry powder inhalers have disregarded the transient nature of the inhalation process, adopting a constant flow rate through the inhaler instead. In addition, either no capsules, a capsule in a fixed position, or a capsule rotating at a constant rate have been considered.
In this work, literature data for three inhalation flow profiles were incorporated into coupled simulations of the air flow and carrier particle motion through an Aerolizer® dry powder inhaler with a rotating capsule and compared to simulations at constant air flow rates. The results for the profile simulations indicated that the carrier powder experienced larger velocity fluctuations. Acceleration events were tracked as a measure of collision- and flow-induced dispersion. The majority of fast particle accelerations occurred when the particles collided with the swirl chamber walls. Of the two common particle dispersion metrics, only the peak particle force distribution appeared to be sensitive to the inhalation profiles, while the effect of the profiles on the cumulative impulse was small.
Display omitted
Protein crystals have many important applications in many fields, including pharmaceutics. Being more stable than other formulations, and having a high degree of purity and bioavailability, they are ...especially promising in the area of drug delivery. In this contribution, the development of a continuously operated tubular crystallizer for the production of protein crystals has been described. Using the model enzyme lysozyme, we successfully generated product particles ranging between 15 and 40 μm in size. At the reactor inlet, a protein solution was mixed with a crystallization agent solution to create high supersaturations required for nucleation. Along the tube, supersaturation was controlled using water baths that divided the crystallizer into a nucleation zone and a growth zone. Low flow rates minimized the effect of shear forces that may impede crystal growth. Simultaneously, a slug flow was implemented to ensure crystal transport through the reactor and to reduce the residence time distribution.
Coating thickness is a critical quality attribute of many coated tablets. Functional coatings ensure correct drug release kinetics or protection from light, while non-functional coatings are ...generally applied for cosmetic reasons. Traditionally, coating thickness is assessed indirectly via offline methods, such as weight gain or diameter growth. In the past decade, several methods, including optical coherence tomography (OCT) and Raman spectroscopy, have emerged to perform in-line measurements of various subclasses of coating formulations. However, there are some obstacles. For example, when using OCT, a major challenge is scattering pigments, such as titanium dioxide and iron oxide, which make the interface between the coating and the tablet core difficult to detect. This work explores novel OCT image evaluation techniques using unsupervised machine learning to compute image metrics. Certain image metrics of highly scattering coatings are correlated with the tablet thickness, and hence indirectly with the coating thickness. The method was demonstrated using a titanium dioxide rich coating formulation. The results are expected to be applicable to other scattering coatings and will significantly broaden the applicability of OCT to at-line and in-line coating thickness measurements of a much larger class of coating formulations.
Display omitted
Enteric coatings are designed to protect active pharmaceutical ingredients (APIs) against untimely release in the stomach. Acid protection of such coatings depends on the coating ...layer thickness and integrity, which must be determined in an accurate and reliable way to ensure the final product’s desired performance. Our work addresses the use of optical coherence tomography (OCT) for characterizing the coating thickness and variability of an enteric-coated drug product and linking them to resistance against gastric fluid. In this study, three batches of enteric-coated tablets drawn during the manufacturing process were investigated. An industrial OCT system was used to establish the coating thickness variability of single tablets (intra-tablet), all tablets in a batch (inter-tablet) and between the batches (inter-batch). Based on the large amount of OCT data, we calculated a critical coating thickness for the investigated film coating, which was found to be 27.4 µm. The corresponding distribution has a mean coating thickness of 44.3 µm ± 7.8 µm. The final coated product has a final mean coating thickness of 63.4 µm ± 8.7 µm, guaranteeing that all tablets meet the quality criterion (i.e., acid protection). Based on the measured thickness distributions, already known distribution functions were considered and an additional, new function was proposed for characterizing the coating thickness distributions in the early stages of industrial coating processes. The proposed approach can be transferred to in-line monitoring of the tablet coating processes, which could drastically improve the production efficiency by ultimately allowing real-time release testing (RTRT).
Display omitted
•High degree of batch intermixing inside the feeder poses challenge to material tracking.•DEM simulation of a full discharge process of a twin-screw feeder.•Particle-level material ...tracking over multiple feeder refills.•Reduced order models to lower computational cost.
Residence time distribution (RTD) modeling has proven to be a valuable tool for material tracking in continuous pharmaceutical processes. Refilling is thoroughly studied in the literature, but the main focus lies on the feed rate disturbances. The impact of the feeders themselves on intermixing of different material batches is often overlooked. Since the experimental methods to measure the RTD feeder discharging processes feeder are complex and material intensive, there is only limited experimental RTD data available in the literature. A DEM (discrete element method) simulation of a discharge of a twin-screw feeder shows that a large fraction of material that is moved and intermixed by the agitator. In addition to the intermixing, there is a tendency to discharge material located above the agitator early. In order to predict the behavior during multiple refill events, three models in order of increasing complexity are presented: (1) A simple exponential RTD assuming perfect intermixing of material batches; (2) a RTD model based on DEM results; (3) particle-level material tracking by extrapolation of the DEM results. All three of these models are able to predict the survival function of old material for late refills at low fill levels, however, earlier refills at high fill levels require more complex models to accurately represent the dynamics inside the hopper of the feeder.
The combination of Computational Fluid Dynamics (CFD) and Discrete Element Model (DEM) is a powerful tool for studying fluidized particulate systems and granular flows. In DEM, the individual ...interaction forces between particles are treated on a particle–particle pair basis, and therefore, this method is computational expensive. In addition, the CFD-calculation of the fluid flow increases the computational effort. Thus, current CFD–DEM simulations are limited to systems with particle numbers not exceeding 105. In order to simulate realistic systems, the recently available Compute Unified Device Architecture (CUDA) technology can be applied, which can perform massively-parallel DEM-simulations with several million particles on a single desk-side Graphics Processing Unit (GPU). The objective of this work is to present a new hybrid approach to solve CFD–DEM problems in gas–solid fluidized beds systems applying an efficient coupling method suitable for large-scale simulations. We are using the CUDA technology for the particle simulation and introducing a coupling methodology with a commercial CFD-code. The coupling method between a CFD-code, running on the CPU and our CUDA-based DEM-code running on the GPU, is introduced and discussed. The numerical results are compared to the CFD–DEM and the experimental results of Van Buijtenen et al. (2011). A good agreement was achieved. Finally, fluidized system simulations with up to 25 million particles are presented, which is an unprecented number.
•Analysis of three commonly used drag force models via simple fixed-bed simulation.•Introduction of hybrid CFD–DEM coupling method suitable for large-scale simulation.•Validation of the simulation and comparison of the results with literature data.•Discussion of the CFD–DEM coupling method advantages.•Evaluation of the computational performances up to 25 millions of particles.
A new framework for calibration of DEM models of cohesive materials is presented. DEM simulations are a widespread numerical tool for modelling granular assemblies and related processes, but require ...careful calibration to give realistic predictive results. To this aim routine small-scale laboratory tests are used to investigate powders rheology while limiting resource investment, in conjunction with a rational approach to reduce computational costs as well. A specific order in numerical tests and calibrated variables is defined, enabling the implementation of an automatic iterative routine to calibrate material parameters from a given dataset. This method is successfully applied to calibrate four common increasingly cohesive powders of pharmaceutical relevance. Our framework saves an important amount of computational expense when compared to traditional parameters design space approaches, and proved to be able to calibrate very cohesive powder exhibiting plastic bulk volume loss of up to 25% (after uniaxial consolidation up to 10kPa).
Display omitted
•New iterative fully-automatized calibration method for DEM models of cohesive materials is developed.•Physical considerations and model properties used alongside angle of repose, compression and shear tests.•Powder bed deformations of up to 25% during compression tests are modelled.•Increasingly cohesive pharmaceutical powders are successfully calibrated.
The dispersion of carrier-based formulations in capsule-based dry powder inhalers depends on several factors, including the patient's inhalation profile and the motion of capsule within the device. ...In the present study, coupled computational fluid dynamics and discrete element method simulations of a polydisperse cohesive lactose carrier in an Aerolizer® inhaler were conducted at a constant flow rate of 100 L/min and considering an inhalation profile of asthmatic children between 5 and 17 years approximated from literature data. In relevant high-speed photography experiments, it was observed that the powder was distributed to both capsule ends before being ejected from the capsule. Several methods of ensuring similar behavior in the simulations were presented. Both the constant flow rate simulation and the profile simulations showed a high powder retention in the capsule (7.37–19.00%). Although the inhaler retention was negligible in the constant flow rate simulation due to consistently high air velocities in the device, it reached values of around 7% in most of the profile simulations.
In all simulations, some of the carrier powder was ejected from the capsule as particle clusters. These clusters were larger in the profile simulation than in the constant flow rate simulation. Of the powder discharged from the capsule, a high percentage was bound in clusters in the profile simulation in the beginning and at the end of the inhalation profile while no more than 10% of the powder ejected from the capsule in the 100 L/min constant flow rate simulation were in clusters at any time.
The powder emission from the capsule was studied, indicating a strong dependency of the powder mass flow from the capsule on the angular capsule position. When the capsule holes face the inhaler's air inlets, the air flow into the capsule restricts the powder discharge.
The presented results provide a detailed view of some aspects of the powder flow and dispersion of a cohesive carrier in a capsule-based inhaler device. Furthermore, the importance of considering inhalation profiles in addition to conventional constant flow rate simulations was confirmed.
Display omitted
•Novel algorithm presented for contact detection in DEM between bi-convex tablets.•Especially suited for parallel GPU implementations due to higher execution convergence.•Comparison made with ...experimental data and glued-sphere approach, favoring the dedicated shape.•Performance numbers showed huge benefit over glued-sphere approach.
The coating of tablets to prevent product degradation or control dissolution is a typical process in its production. Coating uniformity is critical for the quality of final product and batch acceptance. Therefore, the coating process needs to be optimized in order to achieve the desired uniformity and reduce manufacturing costs. Thus, understanding how process parameters such as spray properties, equipment geometry and tablet shape influence the coating process is critical for process optimization and approval by regulatory bodies. However this is a non-trivial task as obtaining information about the detailed processes in a tablet coater via experimental means is limited. Thus, computational modeling is the most feasible option to obtain information about the physical processes affecting the performance of tablet coaters. The most widely used computational method for such numerical modelling is the Discrete Element Method (DEM) where individual particles (tablets) are simulated. However, the computational cost of representing the typical shape of tablets is high for industrially relevant simulations. Thus tablet shape is typically approximated by simpler shapes such as spheres or multi spheres. Even with such simplifications, typical simulations take months to complete making it unfeasible for process optimization and design. In the last decade, the Graphical Processor Unit (GPU) has enabled large-scale simulations of tens of millions of spheres and millions of shaped particles using the XPS code. In this paper, we present an algorithm for modeling accurate bi-convex tablets that is tailored to the GPU. We firstly validate the algorithm and implementation against a number of experiments. Finally we perform a simulation of 20 million tablets in a drum coater to illustrate the usefulness of GPU computing for industrial coating applications. We found that the proposed method yields a good match against the lab scale experiments. For the industrial simulation the proposed method gave a more accurate result compared to the multi sphere approach while being significantly faster.