Classical finite volume schemes for compressible Euler system are not accurate on quadrangular mesh at low Mach number in the sense that they do not converge to the incompressible limit when the Mach number tends to zero [1]. The spurious mode that jeopardizes the convergence can be identified and corresponds to the long time limit of a first order wave system whose properties and discretization depend on the scheme used for the compressible system [2]. Then, the low Mach number accuracy can be analysed by studying the long time solution of the associated wave system. In this presentation, we will propose to enrich the velocity space approximation by adding a divergence free element in each cell. Using this new approximation space and a Godunov’ numerical flux, the long-time limit of the wave system discretization corresponds to the divergence free component of the initial condition (as expected). This divergence free component will be defined via a discrete Hodge-Helmoltz decomposition. As a consequence, for Euler equations, using this new approximation space and a numerical flux that degenerates in the low Mach number limit to a Godunov’ flux for the wave system, the numerical scheme provides pressure fluctuations of order Mach square and a divergence of the velocity field of order Mach in the low Mach number limit. This work can be seen as an extension to quadrangles of the properties obtained on triangles [3].