AbstractIn this study, we present a robust conservative time‐staggered scheme for variable density flow. This pressure correction scheme uses the compressible Navier–Stokes equations and is implemented in the collocated finite‐volume open‐source computational fluid dynamics solver code_saturne. The Helmholtz equation is solved for the pressure increment, taking the thermodynamic pressure into account and avoiding the acoustic time step limitation. The internal energy equation is used and completed by a source term derived from the discrete kinetic energy equation, thus enforcing total energy conservation and consistency for irregular solutions. A numerical analysis providing conditions ensuring the positivity of the thermodynamic variables is proposed. The scheme is verified and validated against analytical and experimental test cases. Its ability to reproduce the pressure variation while conserving the mass is demonstrated. Its conservative property and time convergence order are also verified. An irregular shock solution is studied, emphasizing the importance of the source term in the internal energy equation. Finally, the scheme is validated against reference numerical results on a two‐dimensional natural convection cavity and experimental data on a three‐dimensional ventilation test case. The comparison against experimental data is made using first‐and second‐order turbulent simulations.