During solid–solid precipitation, interface free energy anisotropy is known to drive faceting of precipitates. In this paper, using a recently developed phase field formulation based on higher order tensor terms, we develop and implement a family of phase field models and indicate the parameter choices which lead to faceted precipitate morphologies. We also indicate how to choose the parameters given either the known precipitate morphology or the interfacial free energy anisotropy. Specifically, we study the faceting of precipitates in systems with cubic and hexagonal anisotropies; in 2 and 3D implementation of our phase field model, the precipitates do show facets in accordance with the Wulff plot – including cases where the Wulff plot predicts facets made up of more than one family of planes. We also indicate the possible extensions of our model to study other problems of interest.