Numerical study of a multigrid method with four smoothing methods for the incompressible Navier-Stokes equations in general coordinates