The general mathematical model for batch cooling crystallization was established based on the population balance equation considering the change of slurry volume, and simulated with crystallization thermodynamics, kinetics and mass balance employing bed voidage. In the system of vitamin C-water-ethanol, reliability of this model was verified by comparison between simulation results and experimental data. The effects of operation parameters on product quality can be systematically investigated by modeling simulation.