Digital Design of Batch Cooling Crystallization Processes: Computational Fluid Dynamics Methodology for Modeling Free-Surface Hydrodynamics in Agitated Crystallizers

A framework for the digital design of batch cooling crystallization processes is presented comprising three stages, which are based on different levels of process complexity, integrating crystallizer hydrodynamics with crystallization kinetics and consequently with expected crystal size distribution. In the first stage of the framework, a computational fluid dynamics methodology is developed to accurately assess hydrodynamics in a typical batch crystallizer configuration, comprising a 20 L scale dish-bottom vessel with a single beavertail baffle agitated by a retreat curve impeller, used in the pharmaceutical as well as in the fine chemical industries. The hydrodynamics of crystallizers with such configurations is characterized by vortex formation on the free liquid surface. It is therefore important to model the free surface using the Volume-of-Fluid (VoF) method. Comparison of the predicted mean velocity components with experimental measurements using laser Doppler anemometry reveals that improved predictions are obtained using a differential Reynolds-stress transport model for turbulence coupled with the VoF for modeling the gas-liquid interface compared with those using the Shear-stress transport model and with a flat liquid surface. This study demonstrates that an accurate treatment of the liquid free surface for capturing vortex formation is essential for reliable predictions of the crystallizer’s flow field. While the vortex depth is predicted to increase with increasing impeller Reynolds number, the dependence of hydrodynamic macroparameters, including power number, impeller flow number, and secondary circulation flow number, on Reynolds number reveals that they are essentially constant within the turbulent regime but fluctuate when the flow is in the transitional and laminar regimes as fluid viscosity increases.