Multigrid algorithm for three-dimensional incompressible high-Reynolds number turbulent flows