Due to the small neutron absorption cross section and excellent thermal creep performance, zirconium alloy is one of the most important cladding materials for fuel rods in commercial fission reactors. However, quantitative analysis of the effects of temperature and grain boundaries on the corrosion microstructure evolution of zirconium alloys is still needed. The establishing of a phase field simulation for the corrosion process of polycrystalline zirconium alloy and the systematical investigating of the thermodynamic influence are both very important. In this study, the phase field model of the corrosion process in zirconium alloys is developed by combining corrosion electrochemistry through calculating the interfacial energy at the metal-oxide and oxide-fluid boundaries. Then the model is used to investigate the uniform corrosion behavior on the surface of Zr-2.5Sn alloy, which demonstrates that the corrosion kinetic curve follows a cubic rule. Subsequently, the influence of temperature on the corrosion thickening curve of zirconium alloy is examined, and good agreement between simulation and experimental results is achieved. It is observed that during early stage of oxide layer formation, there is a high growth rate with minimal temperature dependence; however, as the oxide layer thickness increases, temperature becomes a significant factor affecting its growth rate, with higher temperatures resulting in faster corrosion rates. Furthermore, the effect of polycrystalline zirconium alloy matrices on corrosion rate is investigated, revealing that the grain boundaries accelerate oxide layer thickening due to enhanced oxygen diffusion rates. At metal-oxide interface, O ^2–bands are formed in areas with higher O ^2–concentration along these grain boundaries towards the metal matrix, which mainly influences oxidation-corrosion rate during the initial oxidation stage.