  We consider a pressure correction temporal discretization for the
incompressible Navier-Stokes equations in EMAC form. We prove stability and
error estimates for the case of mixed finite element spatial discretization,
and in particular that the Gronwall constant's exponential dependence on the
Reynolds number is removed (for sufficiently smooth true solutions) or at least
significantly reduced compared to the commonly used skew-symmetric formulation.
We also show the method preserves momentum and angular momentum, and while it
does not preserve energy it does admit an energy inequality. Several numerical
tests show the advantages EMAC can have over other commonly used formulations
of the nonlinearity. Additionally, we discuss extensions of the results to the
usual Crank-Nicolson temporal discretization.
