[santiagomarquezd]

Hi all FOAMers, I'm wondering why final phi at the end of each timestep in icoFoam is calculated as:

a)

phi -= pEqn.flux();

instead of doing:

b)

phi=linearInterpolate(U) & mesh.Sf();

particularly having in mind that at this point U is already corrected by PISO.

I've checked the values of both methods and they are slightly different, maybe with a difference of 10-30%. It's interesting that, if you are running continued, phi for the next timestep is calculated like a), but if you stop running and then restart phi for the next timestep then phi is calculated following b).


It should be equal? If not, What is correct and why? Why one is used for continued runnings and the other one for restarting?

Thx in advance.