# buoyantBoussinesqSimpleFoam convergence problem.

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

 September 25, 2012, 10:56 buoyantBoussinesqSimpleFoam convergence problem. #1 New Member   Eskil Aursand Join Date: Sep 2012 Posts: 7 Rep Power: 5 I am trying to get a buoyantBoussinesqSimpleFoam simulation with a k-eps turbulence model to converge, but it never does. The output states that the residuals are smaller than the SIMPLE tolerances, but it still doesn't stop. Here is an example output: Code: ```DILUPBiCG: Solving for Ux, Initial residual = 0.00776199, Final residual = 3.05559e-05, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.0033988, Final residual = 1.59953e-05, No Iterations 1 DILUPBiCG: Solving for T, Initial residual = 0.000345836, Final residual = 1.40283e-05, No Iterations 1 DICPCG: Solving for p_rgh, Initial residual = 0.0300916, Final residual = 0.0002878, No Iterations 15 time step continuity errors : sum local = 7.9083e-05, global = -1.09848e-05, cumulative = 0.00106109 DILUPBiCG: Solving for epsilon, Initial residual = 0.000168969, Final residual = 2.31597e-06, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.001979, Final residual = 3.89602e-05, No Iterations 1``` Here is my fvSolution: Code: ```solvers { p_rgh { solver PCG; preconditioner DIC; tolerance 1e-08; relTol 0.01; } "(U|T|k|epsilon|R)" { solver PBiCG; preconditioner DILU; tolerance 1e-05; relTol 0.1; } } SIMPLE { nNonOrthogonalCorrectors 0; pRefCell 0; pRefValue 0; residualControl { p_rgh 1e-2; U 1e-4; T 1e-2; // possibly check turbulence fields "(k|epsilon|omega)" 1e-3; } } relaxationFactors { fields { p_rgh 0.5; } equations { U 0.2; T 0.5; "(k|epsilon|R)" 0.5; } }``` As seen the final residuals are smaller than the tolerances set in residualControl. If i set the SIMPLE tolerances to 1e-1 it converges however, but I don't understand why they have to be so low. If it matters: The case is 2D with a cold circle in the center, hot walls on the sides and below, and an insulating boundary above. Gravity points down. Solver is buoyantBoussinesqSimpleFoam. Any help and tips will be appreciated!

August 5, 2014, 01:25
Exact same problem
#2
Senior Member

Srivathsan N
Join Date: Jan 2013
Location: India
Posts: 101
Rep Power: 4
Hi all,

I am facing a similar problem. My last time step residuals are
Code:
```DILUPBiCG:  Solving for Ux, Initial residual = 1.13167e-09, Final residual = 1.13167e-09, No Iterations 0
DILUPBiCG:  Solving for Uy, Initial residual = 3.7544e-10, Final residual = 3.7544e-10, No Iterations 0
DILUPBiCG:  Solving for Uz, Initial residual = 0.000492656, Final residual = 4.17232e-06, No Iterations 1
DILUPBiCG:  Solving for T, Initial residual = 9.9924e-06, Final residual = 9.9924e-06, No Iterations 0
DICPCG:  Solving for p_rgh, Initial residual = 1.05994e-08, Final residual = 7.33015e-09, No Iterations 1
time step continuity errors : sum local = 3.02676e-11, global = 9.96442e-20, cumulative = -1.45381e-16```
However my fvSolution file reads as follows,

Code:
```solvers
{
p_rgh
{
solver          PCG;
preconditioner  DIC;
tolerance       1e-08;
relTol          0.01;
}

"(U|T|k|epsilon|R)"
{
solver          PBiCG;
preconditioner  DILU;
tolerance       1e-05;
relTol          0.1;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
pRefCell        0;
pRefValue       0;

residualControl
{
p_rgh           1e-4;
U               1e-5;
T               1e-2;
}```
The residuals also seem to fluctuate very much except temperature that seems to flatten out.
Q: Why hasn't the solution converged even though the SIMPLE tolerances have been met?

Any help would be great!
Attached Images
 Residuals.jpg (38.7 KB, 16 views)
__________________
Regards,

Srivaths

Last edited by Sherlock_1812; August 20, 2014 at 08:30.

 August 20, 2014, 08:59 #3 Member   Cong Join Date: Nov 2012 Posts: 77 Rep Power: 4 Hi, I remember I was faced with the same problems with the buoyantSimpleFoam. The residual never goes down below 1e-6. I had no idea why.

 August 20, 2014, 09:16 #4 Senior Member   Srivathsan N Join Date: Jan 2013 Location: India Posts: 101 Rep Power: 4 I have tried increasing the pseudo time step in the controlDict. There are no oscillation if I do that but the solution diverges after a few iterations. Again, have no idea why. I'm using OF 2.2.2, by the way. __________________ Regards, Srivaths

 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 Forrest_Lei OpenFOAM 3 July 19, 2011 06:00 commonyue Main CFD Forum 1 December 1, 2009 04:54 nasdak CFX 2 June 29, 2009 01:17 Emily FLUENT 2 March 21, 2007 23:18 Balraj Main CFD Forum 3 December 9, 2004 01:24

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