[booz]

Hi OpenFoam experts,

is it a difference for openFoam or get another results, if for example I write in simpleFoam-case solver
...
tmp<fvvectormatrix> UEqn
(
fvm::div(phi, U)
+ turbulence->divR(U)
+ fvc::grad(p)
);

UEqn().relax();
solve(UEqn() == 0);// change
...

instead of
...
tmp<fvvectormatrix> UEqn
(
fvm::div(phi, U)
+ turbulence->divR(U)
);

UEqn().relax();

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

I suppose i get another result, if i view the flows in paraFoam.

Regards,
Khaled

[grtabor]

Actually what you have done here is to change the SIMPLE algorithm very fundamentally. SIMPLE (and PISO) split the discretisation of the momentum equation into pressure and everything else; what you have done is to lump them all in together.

Gavin

[booz]

Hi Gavin,

that is it is a disadvantage with regards to the convergence?

Regards,
Khaled