In this study, a comprehensive investigation of PET glycolysis has been performed. The research included the development of analytical techniques, experiments with different reactor systems, and the development of mathematical models to understand the kinetics of the process. Quantitative HPLC analysis was adopted and optimized for the detection of BHET, while the size-exclusion chromatography method was developed to determine the molecular weight distribution of solid PET residues. Over 33 experiments were performed with magnetically coupled shaft-stirred Amar reactors, resulting in more than 300 experimentally determined BHET concentration points at various reaction times, temperatures, catalyst concentrations, and ethylene glycol-to-PET ratios. Afterwards, a kinetic model was developed to describe the observed phenomena with a validation step.