A high accuracy hybrid method for two-dimensional Navier–Stokes equations