# I have pEqnflux question

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

 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 under-relaxation, 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 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?

 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 I might as well include the link. http://www.cfd-online.com/Wiki/SIMPLE_algorithm

 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 cfd-online Wiki in the wrong order or am I just mis-interpreting 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 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).

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post dbxmcf OpenFOAM Running, Solving & CFD 1 July 11, 2015 15:21 massume Main CFD Forum 0 August 24, 2007 05:37 DAK565656 CFX 0 August 6, 2007 00:06 Kerem FLUENT 0 April 3, 2006 19:24 galeazzo FLUENT 3 March 22, 2005 05:02

All times are GMT -4. The time now is 23:07.