Modelling of pharmaceutical tablet swelling and dissolution using discrete element method

Abstract

This work presents a novel use of the Discrete Element Method (DEM) combined with inter-particle mass transfer in order to simulate polymer swelling and dissolution. Each particle can absorb water and swell, pushing on its neighbours and causing an overall expansion. Once the disentanglement threshold is reached, the polymer dissolves and the particle reduces in size. This paper applies DEM to simulate the radial swelling and dissolution of cylindrical tablets. The method was validated against exact numerical solution of the same system to assess the accuracy of the DEM simulations for different DEM particle sizes. Parametric studies were done to assess the impact of physical parameters – namely the concentration-dependent diffusion coefficient of water through the polymer, the dissolution rate constant of the polymer and the disentanglement threshold of the polymer – on the radial expansion of the tablet. It was found that different settings of the concentration-dependent water diffusion coefficient function could produce similar radial expansion curves but with different internal concentration profiles. Increasing the dissolution rate constant or decreasing the disentanglement threshold of the polymer caused a reduction in the maximum radius of tablet. Lastly, ATR-FTIR spectroscopic imaging was used to obtain chemical images of a pure hydroxy-propyl methylcellulose (HPMC) tablet swelling and dissolving. The model was optimised to match both the HPMC tablet radius and the concentration profiles over time.