[srinath]

Hello

I find the following quite puzzling.
In icoFoam.C, i find the foll. lines

solve(UEqn == -fvc::grad(p));

// --- PISO loop

for (int corr=0; corr<nCorr; corr++)
{
volScalarField rUA = 1.0/UEqn.A();

U = rUA*UEqn.H();
......
}

Solving UEqn presumably changes the value of U. But what is the point if couple of lines later, we store rUA*Ueqn.H() in U

If solving UEqn is a predictor for a U value, what is the point if we never use the value and overwrite it?

I would appreciate it if someone could clarify this as i find it quite puzzling.

Regards
Srinath
p.s. commenting out the solve UEqn line, seems to give visually acceptable results for the cavity case, but diff says the results are different. So clearly it is not a bug, but is doing something quite subtle