CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Non-linearity Pressure Equation -- PISO algorithm

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

LinkBack Thread Tools Display Modes
Old   August 23, 2012, 04:58
Default Non-linearity Pressure Equation -- PISO algorithm
New Member
Gilles De Neyer
Join Date: Nov 2011
Posts: 18
Rep Power: 5
gdeneyer is on a distinguished road
Dear Foamers,

I analysed precisely the code of IcoFoam after having read a lot of documentation about the PISO algorithm and the discretization procedure in OpenFOAM. The algorithm is completely known for me but I still have two questions concerning the implementation:

It is said in Mr. Jasak's thesis that the non linearity is treated during the discretization procedure with decomposition of the normal surface vector in two contributions: one orthogonal and another chosen by diverse maners, like Least Square, ... My expectations were that the same procedure was used for the discretization of the laplacian of u and the laplacian of p (the fact that the decomposition is done for dealing the non linearity, not the model of decomposition used)

However, during the PISO loop, there's a nonOrthogonal loop and, as I analyse that, during this loop only two operations are done
PEqn = ... // creating the equation for pressure with the last value of P
solve PEqn // solving for P ==> P has now a new value

I guess that the effect of non linearity kind of disappear while solving again and again the same equation but with the pressure changing. I would like to know why it is so ? Note : only the pressure changes from one nonOrth iteration to another. The fluxes are corrected only after nNonOrth is reached !!

My second question concerns the compressible PISO algorithm. I read the really interesting doc by Mr. Jasak :
but I still don't understand one thing : in this doc it is said that the momentum equation is solved first, then the PISO loop containing pressure equation and (not compulsory if the coupling density-pressure is much less important that coupling velocity-pressure) the energy equation. Then, after this is done, the density is given by : either solving the continuity equation with new velocity field or by using the law of perfect gazes . What I don't understand is that in the PISO code, both are done ! We first explicitely solve the continuity equation with old velocity field to find density, then we do the algorithm as explained and the density is again computed with the law of perfect gazes. My opinion is that the first operation solve(rhoEqn) is not needed as the density from the old timestep already satisfy the continuity equation with the old velocity field !

Thank you in advance for reading all of this :-),

gdeneyer is offline   Reply With Quote

Old   August 23, 2012, 05:19
New Member
Gilles De Neyer
Join Date: Nov 2011
Posts: 18
Rep Power: 5
gdeneyer is on a distinguished road
I saw that for my first question the loop was based on a deffered correction approach with the orthogonal part of the laplacian calculated implicitely and the non orthogonal part explicitely. Is this correct ?

Then two others questions :
- does someone have any documentation about the convergence/efficiancy of this method
- why we don't use it when solving the momentum predictor ? are we considering that the non linear effects can be treated only in the pressure equation and will suffice have a velocity non dependant of non-linear effect after velocity correction ?
gdeneyer is offline   Reply With Quote


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
Pressure Outlet Guage pressure Mohsin FLUENT 34 January 12, 2012 06:45
Pressure BC for combustion chamber Giuki FLUENT 1 July 19, 2011 11:35
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
SIMPLER Algorithm Bharath Main CFD Forum 2 July 10, 2009 10:53
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 22:02

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