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

Problem using PISO

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 6, 2009, 13:20
Default Problem using PISO
  #1
Senior Member
 
Antonio Martins
Join Date: Mar 2009
Location: Porto, Porto, Portugal
Posts: 112
Rep Power: 8
titio is on a distinguished road
Send a message via MSN to titio Send a message via Skype™ to titio
Hi Foamers,

I am having some problems using PISO to simulate the flow of complex fluids, in particular polymers, where some extra stresses have to be considered. I can already simulate the flow in a simple system, infinite parallel, but when I compare the analytical solution with the predictions of OpenFoam, I manage to obtain the correct velocity profile but the predicted pressure gradient is smaller than expected. On possible problem is the analytical solution is wrong, but I doube check it and it is correct.

I am thinking on the the possibility of doing more PISO steps, to improve the calculation of the pressure field, or changing the PISO algorithm implemented in icoFoam, that I am using.

Has anybody any ideas of what is going on? Any help will be much appreciated.

Regards,

Titio
titio is offline   Reply With Quote

Old   July 7, 2009, 02:26
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hi Antonio,

first, an elementary question: is the pressure gradient wrong or is the value of pressure simply different? I ask because if you use an incompressible solver in OF, the actual value of the pressure you find depends on the reference value. However in such a case the gradient should be OK.

About your specific case, without knowing more on the model you use to describe the non-Newtonian stress tensor, and how you implemented it, it is impossible to reply.

Best,
A.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   July 7, 2009, 05:48
Default
  #3
Senior Member
 
Antonio Martins
Join Date: Mar 2009
Location: Porto, Porto, Portugal
Posts: 112
Rep Power: 8
titio is on a distinguished road
Send a message via MSN to titio Send a message via Skype™ to titio
Quote:
Originally Posted by alberto View Post
Hi Antonio,

first, an elementary question: is the pressure gradient wrong or is the value of pressure simply different? I ask because if you use an incompressible solver in OF, the actual value of the pressure you find depends on the reference value. However in such a case the gradient should be OK.

About your specific case, without knowing more on the model you use to describe the non-Newtonian stress tensor, and how you implemented it, it is impossible to reply.

Best,
A.
Thanks Alberto,

The pressure value is simply different, by at must 10%, I tried to change the mesh and the problem remains. I am testing the possibility of iterating several times in a PISO step, this time including the momentum predictor. Do you have more ideas?

Regards,

Titio
titio is offline   Reply With Quote

Reply

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
Incoherent problem table in hollow-fiber spinning Gianni FLUENT 0 April 5, 2008 10:33
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13
Problem in Tutorial problem of fluent Phanindra FLUENT 5 April 17, 2007 09:57
GAMBIT meshing problem Gauthier Lambert Main CFD Forum 1 August 3, 2000 09:22
PISO vs. SIMPLEC Andreas Abdon Main CFD Forum 1 March 4, 2000 03:18


All times are GMT -4. The time now is 13:50.