# convergence in pimpleFoam

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

 February 10, 2012, 06:29 convergence in pimpleFoam #1 Senior Member   Andrea Ferrari Join Date: Dec 2010 Posts: 275 Rep Power: 7 Hi all, I normally work with multiphase flow and interFoam, but i need some single phase simulations of flow in a channel with obstacles. My geometry is very simple: inlet, outlet and walls (2D). The simulation type is laminar. I'm using fixedVelocity at the inlet, fixedPressure at the outlet and zero velocity on the walls, the others BC are set to zeroGradient. I found less convergence to a stable solution when the injection velocity is lower. For u=1e-2 m/s and u=1e-3 the solution becomes stable in a short time and the residuals are good, but for u=1e-4 m/s the pressure solution is oscillating and residual doesn't converge. Any ideas of why? at most I would expect the opposite. regards andrea

 February 10, 2012, 07:57 #2 Senior Member     A_R Join Date: Jun 2009 Posts: 118 Rep Power: 8 Dear Andrea you decrease your velocity at inlet so the field of velocity and pressure, which you have solved, have lower value. you should decrease you residual limit one or two orders. check your fvsolution.

 February 10, 2012, 10:59 #3 Senior Member   Andrea Ferrari Join Date: Dec 2010 Posts: 275 Rep Power: 7 Hi and thanks for reply. this is my fvSolution for u=1e-4 m/s U { solver PBiCG; preconditioner DILU; tolerance 1e-9; relTol 0.001; } UFinal { solver PBiCG; preconditioner DILU; tolerance 1e-9; relTol 0; } p { solver GAMG; tolerance 1e-9; relTol 0.001; smoother GaussSeidel; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } pFinal { solver GAMG; tolerance 1e-9; relTol 0; smoother GaussSeidel; cacheAgglomeration true; nCellsInCoarsestLevel 10; agglomerator faceAreaPair; mergeLevels 1; } } PIMPLE { nOuterCorrectors 2; nCorrectors 2; nNonOrthogonalCorrectors 0; pRefCell 0; pRefValue 0; } I've obtained stable solution for u=1e-2 m/s after 0.06s with tolerance set to 1e-5. I would expect almost immediate convergence with u=1e-4, or not? lower tolerance has the effect of making the pressure solution more stable (smaller or no oscillation) but after 5s of simulations the resisual are still high. What is your opinion? best andrea

 February 10, 2012, 14:20 #4 Senior Member     A_R Join Date: Jun 2009 Posts: 118 Rep Power: 8 Dear Andrea if you found your stable solution when you set residuals in 10e-9 for u=10e-2 and 10e-3. it is better to reduce your residual to 10e-11 or -12 for u=10e-4. try it

 February 10, 2012, 15:23 #5 Senior Member   Andrea Ferrari Join Date: Dec 2010 Posts: 275 Rep Power: 7 Dear niaz, sorry i have been not clear. I found stable solution for u=1e-2/1e-3 setting residuals to 1e-5, so i guess 1e-9 is ok for u=1e-4. btw the solution doesn't converge even with 1e-12. it's really strange! regards andrea

 February 11, 2012, 03:32 #6 Senior Member     A_R Join Date: Jun 2009 Posts: 118 Rep Power: 8 Dear Andrea run your test case and save your residual report in a file them check that the residuals decrease to your settings or not? if not, change the value of loop corrector in the fvsolution.

 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 Centurion2011 FLUENT 24 May 9, 2015 08:02 colopolo CFX 13 October 4, 2011 22:03 nasdak CFX 2 June 29, 2009 01:17 ganesh Main CFD Forum 4 June 30, 2006 14:20 Chetan FLUENT 3 April 15, 2004 19:13

All times are GMT -4. The time now is 15:09.