The phase-field method has been widely employed recently for simulating grain growth. Phase-field grain growth models are classified into two types according to their conservation constraints for phase-field variables: the multi-phase-field model and the continuum-field model. In addition, within the multi-phase-field model framework, three models with different formulations exist. These models are reported to accurately simulate grain growth under conditions of isotropic or weakly anisotropic grain boundary energy and mobility. However, for cases of strongly anisotropic grain boundary properties, the accuracy of these models has not yet been examined in detail. In this study, using the continuum-field model and three different multi-phase-field models, systematic grain growth simulations with anisotropic grain boundary energies and mobilities are performed. Through the detailed investigation of the accuracy of the simulated results, the suitability of each model for anisotropic grain growth simulations is revealed. Furthermore, based on the higher-order terms, accuracy improvement of the phase-field models is attempted.