[vgalindo]

One can find at the end of simpleFOAM the code line

/ / Explicitly relax pressure for momentum corrector
p.relax ();

but this line is missing in pisoFoam.C.

Why?

best regards,
V. G.

[nimasam]

piso Foam is a transient solution so relaxation can not be true for a transient solution look CFD book you can find more about difference between piso and simple algorithm

[akidess]

You can also use relaxation with transient solutions, see eg. PIMPLE or transient SIMPLE. However, the whole point about the PISO algorithm is to eliminate the need for underrelaxation. I do agree with you VG. should have a look into a CFD book on this though.

Originally Posted by nimasam

piso Foam is a transient solution so relaxation can not be true for a transient solution

[sharonyue]

Originally Posted by akidess

You can also use relaxation with transient solutions, see eg. PIMPLE or transient SIMPLE. However, the whole point about the PISO algorithm is to eliminate the need for underrelaxation. I do agree with you VG. should have a look into a CFD book on this though.

But in pisoFoam.C, we have UEqn.relax(), why UEqn is relaxed?
and whats the difference between p.relax() with PEqn.relax()?

BTW, in Verssteeg'book page196, it said:

As before, under - relaxation is required with the above procedure to stabilise the caculation process.

[akidess]

It doesn't really matter that much if you relax UEqn. In fact, you can throw out the solution UEqn completely by deactivating the momentum predictor.

As for your question about the two types of relaxation, this has been answered multiple times in the past, e.g. here.

BTW, in Verssteeg'book page196, it said

Page numbers vary between editions of a book, and the quote you give is so isolated it is of limited use to someone without access to the book.

[sharonyue]

Originally Posted by sharonyue

But in pisoFoam.C, we have UEqn.relax(), why UEqn is relaxed?
and whats the difference between p.relax() with PEqn.relax()?

BTW, in Verssteeg'book page196, it said:

Okay, after a long time, I got some documents(Chinese) to talk about this, this is what they said, I tried to translate it:

In SIMPLE, cuz there is an omission of velocity correction of neighbor cells, using this u^' to correct velocity is moderate. But for pressure, using this u^' to correct pressure is exaggerated. Cause we need to do explicit relaxation in pressure field to make pressure correction to be moderate.

For velocity, we need to use relaxation to make matrix solver robuster. Then we use Ueqn.relax.

Looks interesting, shed a light on this but not very strightforward.