High-order relaxation methods for incompressible Navier-Stokes equations