“…To circumvent the limitations of standard time integration schemes, several general methods for conservative or dissipative ODEs such as orthogonal projection [29, Section IV.4] and relaxation [33,52,55] as well as more problem-dependent methods for dissipative ODEs such as artificial dissipation or filtering [28,40,63] have been proposed, mainly in the context of one-step methods. For fluid mechanics applications, orthogonal projection methods are not really suitable since they do not preserve linear invariants such as the total mass, [29]. In contrast, relaxation methods conserve all linear invariants and can still preserve the correct global entropy conservation/dissipation in time [33,51,52,55].…”