Material dissolution is a critical attribute of many products in a wide variety of industries. The idealized view of dissolution through established prediction tools should be reconsidered because the number of new substances with low aqueous solubility is increasing. Due to this, a fundamental understanding of the dissolution process is desired. The aim of this study was to develop a tool to predict crystal dissolution performance based on experimentally measurable physical parameters. A numerical simulation, called the phase-field method, was used to simultaneously solve the time evolution of the phase and concentration fields of dissolving particles. This approach applies to diffusion-limited as well as surface reaction-limited systems. The numerical results were compared to analytical solutions, and the influence of particle shape and interparticle proximity on the dissolution process was numerically investigated. Dissolution behaviors of two different substances were modeled. A diffusion-limited model compound, xylitol, with a high aqueous solubility and a surface reaction-limited model compound, griseofulvin, with a low aqueous solubility were chosen. The results of the simulations demonstrated that phase-field modeling is a powerful approach for predicting the dissolution behaviors of pure crystalline substances.
Keywords: crystal; diffusion; dissolution process; interface velocity; numerical model; phase-field simulation; surface reaction.