Enhanced solution of 2D incompressible Navier–Stokes equations based on an immersed-boundary generalized harmonic polynomial cell method