Mass flux correction in SIMPLE
 

Hi, as my unanswered previous posts indicate, I have a problem with my SIMPLE code while I use variable density. As an everyday routine (for a month already) I was reading books, googling and studying my code then I came across this detail in OpenFOAM wiki and also in CFD-online wiki:

1. Set the boundary conditions.
2. Solve the discretized momentum equation to compute the intermediate velocity field.
3. Compute the mass fluxes at the cells faces.
4. Solve the pressure equation and apply under-relaxation.
5. Correct the mass fluxes at the cell faces.
6. Correct the velocities on the basis of the new pressure field.
7. Update the boundary conditions.
8. Repeat till convergence.

I don't correct mass fluxes after solving pressure correction equation however, my code works fine when I use constant density, it even works fine when I use variable density for all subroutines except pressure correction. Mass flux correction is not mentioned in Vertseeg book nor in patankar's. Where can I learn more about this step and can someone confirm if this is the cause of larger errors when I use variable density?