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

I have pEqnflux question

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 18, 2006, 15:56
Default For a Steady State SIMPLE solv
  #1
New Member
 
Daniel Mills
Join Date: Mar 2009
Posts: 13
Rep Power: 8
daniel_mills is on a distinguished road
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?
daniel_mills is offline   Reply With Quote

Old   October 18, 2006, 16:00
Default Very easy: just swap the state
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Very easy: just swap the statements.

Unfortunately, it would also be not too clever: if you update the fluxes after pressure under-relaxation, you will violate continuity.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   October 18, 2006, 16:13
Default Thank you for your quick respo
  #3
New Member
 
Daniel Mills
Join Date: Mar 2009
Posts: 13
Rep Power: 8
daniel_mills is on a distinguished road
Thank you for your quick response. I am a CFD newbie. I was just going off the cfd-online 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 under-relaxation 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 under-realxing it. I was just wondering if this had somthing to do with the fluxes of phi?
daniel_mills is offline   Reply With Quote

Old   October 18, 2006, 16:14
Default I might as well include the li
  #4
New Member
 
Daniel Mills
Join Date: Mar 2009
Posts: 13
Rep Power: 8
daniel_mills is on a distinguished road
I might as well include the link.

http://www.cfd-online.com/Wiki/SIMPLE_algorithm
daniel_mills is offline   Reply With Quote

Old   October 18, 2006, 17:23
Default Are steps 6,7, and 8 for SIMPL
  #5
New Member
 
Daniel Mills
Join Date: Mar 2009
Posts: 13
Rep Power: 8
daniel_mills is on a distinguished road
Are steps 6,7, and 8 for SIMPLE_algorithm on the cfd-online Wiki in the wrong order or am I just mis-interpreting the way it is laid out?
daniel_mills is offline   Reply With Quote

Old   October 20, 2006, 10:20
Default Can someone please answer a qu
  #6
New Member
 
Daniel Mills
Join Date: Mar 2009
Posts: 13
Rep Power: 8
daniel_mills is on a distinguished road
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.
daniel_mills is offline   Reply With Quote

Old   October 21, 2006, 01:37
Default 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
msrinath80 is on a distinguished road
Normally step 6 (Update the pressure field: where urf is the under-relaxation 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).
msrinath80 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
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


All times are GMT -4. The time now is 19:51.