We study extensions of an earlier developed energy and helicity preserving scheme for the 3D Navier-Stokes equations and apply them to a more general class of problems. The scheme is studied together with stabilizations of grad-div type in order to mitigate the effect of the Bernoulli pressure error on the velocity error. We prove stability, convergence, discuss conservation properties, and present numerical experiments that demonstrate the advantages of the scheme
