
[Sponsors] 
October 18, 2006, 15:56 
For a Steady State SIMPLE solv

#1 
New Member
Daniel Mills
Join Date: Mar 2009
Posts: 13
Rep Power: 8 
For a Steady State SIMPLE solver how hard would it be to implicitly solve for p, then explicitly relax p, then update the fluxes based on the explicitly relaxed value of p, instead of implicitly solving, then updating fluxes and then explicitly relaxing?


October 18, 2006, 16:00 
Very easy: just swap the state

#2 
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21 
Very easy: just swap the statements.
Unfortunately, it would also be not too clever: if you update the fluxes after pressure underrelaxation, you will violate continuity. Hrv
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk 

October 18, 2006, 16:13 
Thank you for your quick respo

#3 
New Member
Daniel Mills
Join Date: Mar 2009
Posts: 13
Rep Power: 8 
Thank you for your quick response. I am a CFD newbie. I was just going off the cfdonline description of SIMPLE in which it states the following.
1.Set the boundary conditions. 2.Compute the gradients of velocity and pressure. 3.Solve the discretized momentum equation to compute the intermediate velocity field. 4.Compute the uncorrected mass fluxes at faces . 5.Solve the pressure correction equation to produce cell values of the pressure correction . 6.Update the pressure field: where urf is the underrelaxation factor for pressure. 7.Update the boundary pressure corrections . 8.Correct the face mass fluxes: 9.Correct the cell velocities: ; where is the gradient of the pressure corrections, is the vector of central coefficients for the discretized linear system representing the velocity equation and Vol is the cell volume. 10.Update density due to pressure changes. My problem is I actually get greater stability with the relaxation value for p set to 1 early in my start up. Than I do underrealxing it. I was just wondering if this had somthing to do with the fluxes of phi? 

October 18, 2006, 16:14 
I might as well include the li

#4 
New Member
Daniel Mills
Join Date: Mar 2009
Posts: 13
Rep Power: 8 

October 18, 2006, 17:23 
Are steps 6,7, and 8 for SIMPL

#5 
New Member
Daniel Mills
Join Date: Mar 2009
Posts: 13
Rep Power: 8 
Are steps 6,7, and 8 for SIMPLE_algorithm on the cfdonline Wiki in the wrong order or am I just misinterpreting the way it is laid out?


October 20, 2006, 10:20 
Can someone please answer a qu

#6 
New Member
Daniel Mills
Join Date: Mar 2009
Posts: 13
Rep Power: 8 
Can someone please answer a question for me?
When I am solving for an unsteady swirling flow using rhoSimpleFoam why does setting my relaxation for pressure to around 0.9 actually stabilize the unstable pressure oscillations that develop when the flow field starts to swirl? I would think the more I lower the pressure relaxation the more stable it would be, but the exact opposite is true. Why?? My relaxations are as follows: rho 0.1 p 0.9 U 0.7 h 0.7 This works well and gives me fairly similar numbers to fluent. 

October 21, 2006, 01:37 
Normally step 6 (Update the pr

#7 
Senior Member
Srinath Madhavan (a.k.a pUl)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 698
Rep Power: 12 
Normally step 6 (Update the pressure field: where urf is the underrelaxation factor for pressure.) is done at the last, but I guess it would not make any difference as the extrapolation of pressure correction at the boundary(s) or the updating of mass fluxes have nothing to do with the corrected pressure [at that point in time]. It is only in the next outer iteration that the corrected pressure itself is used for computation of the gradient (source term).


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
PEqnflux in compressible and incompressible PISO loop  dbxmcf  OpenFOAM Running, Solving & CFD  0  October 7, 2008 10:58 
question  massume  Main CFD Forum  0  August 24, 2007 05:37 
CFX 10 question  DAK565656  CFX  0  August 6, 2007 00:06 
UDF question  Kerem  FLUENT  0  April 3, 2006 19:24 
Question about y+  galeazzo  FLUENT  3  March 22, 2005 05:02 