Modelling diffusive mixing in antisolvent crystallization

Diffusion controls local concentration profiles at interfaces between segregated fluid elements during mixing processes. This is important for antisolvent crystallization, where it is intuitively argued that local concentration profiles at interfaces between solution and antisolvent fluid elements can result in significant supersaturation overshoots over and above that at the final mixture composition, leading to poorly controlled nucleation. Previous work on modeling diffusive mixing in antisolvent crystallization has relied on Fickian diffusion, where concentration gradients are the driving force for diffusion. This predicts large overshoots in the supersaturation at interfaces between solution and antisolvent, as is often intuitively expected. However, chemical potential gradients provide a more physically realistic driving force for diffusion, and in highly nonideal solutions, such as those in antisolvent crystallization, this leads to nonintuitive behavior. In particular, as solute diffusion toward antisolvent is severely hindered, it can diffuse against its concentration gradient away from antisolvent. We apply thermodynamically consistent diffusion model based on the multicomponent Maxwell–Stefan formulation to examine diffusive mixing in a nonideal antisolvent crystallization system. Large supersaturation overshoots above that at the final mixture composition are not found when a thermodynamically consistent approach is used, demonstrating that these overshoots are modeling artifacts and are not expected to be present in physical systems. In addition, for certain conditions, localized liquid–liquid spinodal demixing is predicted to occur during the diffusive mixing process, even when the final mixture composition is outside the liquid–liquid phase separation region. Intermittent spinodal demixing driven by diffusive mixing may provide a novel explanation for differences of nucleation behaviors among various antisolvents.