Numerical Discretization of Variational Phase Field Model for Phase Transitions in Ferroelectric Thin Films. Commun. Comput. Phys. 2021;29(3):729-746. doi:10.4208/cicp.OA-2020-0118