The new numerical method has to obey only a mild CFL condition, which is based on the fluid velocity and not on the sound speed. This makes the scheme particularly interesting for low Mach number flows, because large time steps are permitted. Moreover, being locally and globally conservative, the new method behaves also very well in the presence of shock waves. The proposed algorithm is first validated against the exact solution of a large set of one-dimensional Riemann problems for inviscid flows with three different EOS: the ideal gas law, the van der Waals EOS and the Redlich−Kwong EOS. In the final part of the paper, the method is extended to the two-dimensional viscous case.