CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

missing p.relax() in pisoFOAM

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree2Likes
  • 1 Post By akidess
  • 1 Post By akidess

Reply
 
LinkBack Thread Tools Display Modes
Old   December 10, 2010, 12:00
Default missing p.relax() in pisoFOAM
  #1
New Member
 
Join Date: Nov 2010
Posts: 2
Rep Power: 0
vgalindo is on a distinguished road
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.
vgalindo is offline   Reply With Quote

Old   December 11, 2010, 01:59
Default
  #2
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to 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
nimasam is online now   Reply With Quote

Old   December 13, 2010, 06:04
Default
  #3
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 919
Rep Power: 17
akidess will become famous soon enough
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.

Quote:
Originally Posted by nimasam View Post
piso Foam is a transient solution so relaxation can not be true for a transient solution
sharonyue likes this.
akidess is offline   Reply With Quote

Old   September 2, 2013, 23:47
Default
  #4
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 669
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by akidess View Post
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:
Quote:
As before, under - relaxation is required with the above procedure to stabilise the caculation process.
sharonyue is offline   Reply With Quote

Old   September 3, 2013, 04:02
Default
  #5
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 919
Rep Power: 17
akidess will become famous soon enough
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.

Quote:
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 likes this.
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
*Help define the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oam-technology
akidess is offline   Reply With Quote

Old   June 1, 2015, 05:50
Default
  #6
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 669
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by sharonyue View Post
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:


Code:
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.
__________________
Im the translator of OpenFOAM User Chinese Edition!
But always newbie on CFD. Never too old to learn CFD. Worship these CFD experts.
sharonyue is offline   Reply With Quote

Reply

Tags
p.relax(), pisofoam

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Differences in solution method for pisoFoam and buoyantBoussinesqPisoFoam mchurchf OpenFOAM 6 September 10, 2013 23:24
Buoundary Condition for PisoFoam vishal OpenFOAM 4 October 21, 2010 09:17
the result of pitzDaily using simpleFoam and pisoFoam Kr_kim OpenFOAM Running, Solving & CFD 1 March 16, 2010 14:38
Problem with paraFoam ata OpenFOAM 14 November 30, 2009 05:23
Missing tets along line when meshing with ICEM CFD Georges P. Côté CFX 6 March 23, 2006 01:34


All times are GMT -4. The time now is 12:35.