A mathematical model was established to characterize the freeze concentration process based on the volume increase of materials, which was correlated with the crystallinity of the solution. The relationship between the height of liquid column and the solute concentration was measured on a multi-pass freeze concentrator. The results were compared with the model prediction and were used for model modification and the goodness of fit was high. The crystallization rate error was within 5%.