[Houthuys]

Dear all,

To the best of my knowldege, in PISO algorithm we do several pressure correction cycles to ensure that the flux is divergence free. At the end of each cycle we calculate the corrected velocity based on the updated pressure field,

U = HbyA - rAU*fvc::grad(p);

I can see that this new velocities will be used to update H(U) operator.

But why we also calculate phi here:

phi = phiHbyA - pEqn.flux();

I have checked that the value of rAU will be constant for each time step (correct me if that's wrong). I also checked that when I removed this line the results will be the same ( I've tested cavity case for with icoFoam).

So, for which purpose we calculate the new flux at the end of each PISO loop in icoFOAM?

Best regards