  In this paper, we develop and analyze a novel numerical scheme for the steady
incompressible Navier-Stokes equations by the weak Galerkin methods. The
divergence-preserving velocity reconstruction operator is employed in the
discretization of momentum equation. By employing the velocity construction
operator, our algorithm can achieve pressure-robust, which means, the velocity
error is independent of the pressure and the irrotational body force. Error
analysis is established to show the optimal rate of convergence. Numerical
experiments are presented to validate the theoretical conclusions.
