Ferroelectric hafnium zirconium oxide (HZO) thin films show significant promise for applications in ferroelectric random-access memory, ferroelectric field-effect transistors, and ferroelectric tunneling junctions. However, there are shortcomings in understanding ferroelectric switching, which is crucial in the operation of these devices. Here a computational model based on phase field method is developed to simulate the switching behavior of polycrystalline HZO thin films. Furthermore, we introduce a novel