We use geometrical arguments based on grain boundary symmetries to introduce crystalline interfacial energies for interfaces in polycrystalline thin films with a cubic lattice. These crystalline energies are incorporated into a multi-phase field model. Our aim is to apply the multi-phase field method to describe the evolution of faceted grain boundary triple junctions in epitaxially growing microstructures. In particular, we are interested in symmetry properties of triple junctions in tricrystalline thin films. Symmetries of triple junctions in tricrystalline films have been studied in experiments by Dahmen and Thangaraj.6,25 In accordance with their experiments, we find in numerical simulations that any two neighboring triple junctions belong to different symmetry classes. We introduce a local equilibrium condition at triple junctions which can be interpreted as a crystalline version of Young's law. The local equilibrium condition at triple junctions is purely determined by the grain boundary energies. In particular no triple junction energies are necessary to explain which triple junctions are possible. All triple junctions observed in the experiments as well as in the simulations fulfil the crystalline version of Young's law. Our approach is also capable of describing grain boundary motion in general polycrystalline thin films.