Multigrid Methods for Nonlinear PDE Systems with Applications in Phase-Field Models