In this study, we carried out a para-xylene crystallization experiment at constant temperature and concentration levels. Throughout the process, the kinetics of nucleation, growth, breakage, and aggregation of para-xylene particles were measured and built using a morphological approach. An additional a three-stage continuous suspension crystallization separation experiment was carried out, the process for which was simulated using the population balance model based on correlated kinetic equations. The population balance equation was solved using an extended moment of classes algorithm, and the solving process was implemented in MATLAB. In this case, the predicted particle size distribution of the products matched well with the experiment. In order to provide references for the optimization of the industrial para-xylene crystallization process, a three-stage suspension crystallization separation experiment was designed and conducted, in which each crystallizer had a distinct operating temperature and mean residence time. The effects of operating parameters on the final product were investigated further. The proposed models and algorithms can also be applied in other cases and provide an alternative approach for optimizing continuous crystallization processes.