ABSTRACT In the inviscid limit, the energy of a velocity field satisfying the incompressible Navier–Stokes equations is conserved. Non‐dissipative numerical methods that discretely mimic this energy conservation feature have been demonstrated in the literature to be extremely valuable for robust and accurate large‐eddy simulations of high Reynolds number incompressible turbulent flows. For complex geometries, such numerical methods have been traditionally developed using the finite volume framework and they have b…