# Maximum number of iterations exceeded - rhoSimpleFoam

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

 June 5, 2016, 17:01 Maximum number of iterations exceeded - rhoSimpleFoam #1 New Member   A Z Join Date: Dec 2015 Posts: 16 Rep Power: 7 Hi everyone, I'm try to simulate a Pipe using rhoSimpleFoam. I get after a few teps the "Maximum number of iterations exceeded" error. I tried some solutions I found here but without results. (changing the relaxationFactors & schemes, setting TMin, Tmax ... , checking and changing Boundary Conditions) The Mesh looks fine. Below the logfile: Code: ```Create time Create mesh for time = 0 SIMPLE: convergence criteria field p tolerance 0.001 field U tolerance 0.001 field h tolerance 0.001 field "(k|epsilon|omega)" tolerance 0.001 Reading thermophysical properties Selecting thermodynamics package { type hePsiThermo; mixture pureMixture; transport const; thermo hConst; equationOfState perfectGas; specie specie; energy sensibleEnthalpy; } Reading field U Reading/calculating face flux field phi Creating turbulence model Selecting RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; C3 -0.33; sigmak 1; sigmaEps 1.3; } No MRF models present No finite volume options present Starting time loop Time = 1e-05 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.022547925, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.055572477, No Iterations 3 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.05735183, No Iterations 3 smoothSolver: Solving for h, Initial residual = 0.99560741, Final residual = 0.00054553968, No Iterations 1 GAMG: Solving for p, Initial residual = 1, Final residual = 6.8005097e-09, No Iterations 35 GAMG: Solving for p, Initial residual = 0.029825433, Final residual = 7.8671775e-09, No Iterations 17 GAMG: Solving for p, Initial residual = 0.0016980265, Final residual = 8.9338721e-09, No Iterations 13 time step continuity errors : sum local = 7.1954987e-10, global = -9.7011997e-12, cumulative = -9.7011997e-12 rho max/min : 1.2639489 1.2014065 smoothSolver: Solving for epsilon, Initial residual = 0.18539279, Final residual = 0.018263051, No Iterations 2 bounding epsilon, min: -1.0623939 max: 113.81697 average: 1.8930911 smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.097817915, No Iterations 2 ExecutionTime = 251.38 s ClockTime = 251 s Time = 2e-05 smoothSolver: Solving for Ux, Initial residual = 0.028023654, Final residual = 0.00072812086, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.36155382, Final residual = 0.033268639, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.33901192, Final residual = 0.0069851927, No Iterations 1 smoothSolver: Solving for h, Initial residual = 2.885836e-06, Final residual = 6.91766e-09, No Iterations 1 GAMG: Solving for p, Initial residual = 0.021754383, Final residual = 5.2425385e-09, No Iterations 31 GAMG: Solving for p, Initial residual = 0.88325208, Final residual = 8.0467665e-09, No Iterations 26 GAMG: Solving for p, Initial residual = 0.20759689, Final residual = 8.2151254e-09, No Iterations 22 time step continuity errors : sum local = 2.9362986e-08, global = 2.2426915e-10, cumulative = 2.1456795e-10 rho max/min : 1.4111591 1.1097508 smoothSolver: Solving for epsilon, Initial residual = 0.99460634, Final residual = 0.087019508, No Iterations 3 bounding epsilon, min: -4188.6099 max: 58007085 average: 1010.3673 smoothSolver: Solving for k, Initial residual = 0.99980048, Final residual = 0.072515354, No Iterations 3 bounding k, min: -56.453502 max: 8.5420722e+08 average: 820.94691 ExecutionTime = 503.97 s ClockTime = 504 s Time = 3e-05 smoothSolver: Solving for Ux, Initial residual = 0.72842759, Final residual = 0.01241609, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.63153151, Final residual = 0.028186555, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.72692847, Final residual = 0.056586838, No Iterations 2 smoothSolver: Solving for h, Initial residual = 1, Final residual = 0.0014389051, No Iterations 1 GAMG: Solving for p, Initial residual = 0.98094699, Final residual = 8.5390336e-09, No Iterations 54 GAMG: Solving for p, Initial residual = 6.8142427e-06, Final residual = 7.5921063e-09, No Iterations 14 GAMG: Solving for p, Initial residual = 5.3619733e-06, Final residual = 7.6303402e-09, No Iterations 5 time step continuity errors : sum local = 0.0049250584, global = -9.5693406e-05, cumulative = -9.5693192e-05 rho max/min : 1.5289273 1.0278007 smoothSolver: Solving for epsilon, Initial residual = 0.83663086, Final residual = 0.080689562, No Iterations 2 bounding epsilon, min: -2182.0455 max: 1.0398029e+16 average: 9.9707159e+09 smoothSolver: Solving for k, Initial residual = 0.99736547, Final residual = 0.069018573, No Iterations 3 bounding k, min: -199.49558 max: 9.9985416e+13 average: 1.6617361e+08 ExecutionTime = 745.56 s ClockTime = 747 s Time = 4e-05 smoothSolver: Solving for Ux, Initial residual = 0.43323326, Final residual = 0.011644508, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.25795325, Final residual = 0.010723061, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.29775014, Final residual = 0.0085255883, No Iterations 1 smoothSolver: Solving for h, Initial residual = 0.089014636, Final residual = 0.00040056378, No Iterations 1 --> FOAM FATAL ERROR: Maximum number of iterations exceeded From function Foam::scalar Foam::species::thermo::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo::*)(Foam::scalar)const) const [with Thermo = Foam::hConstThermo >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo = Foam::species::thermo >, Foam::sensibleEnthalpy>] in file /home/openfoam/OpenFOAM/OpenFOAM-3.0.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::species::thermo >, Foam::sensibleEnthalpy>::THs(double, double, double) const at ??:? #3 Foam::hePsiThermo >, Foam::sensibleEnthalpy> > > >::calculate() at ??:? #4 Foam::hePsiThermo >, Foam::sensibleEnthalpy> > > >::correct() at ??:? #5 ? at ??:? #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #7 ? at ??:? Abort (core dumped)``` CheckMesh Code: ```Create time Create polyMesh for time = 0 Time = 0 Mesh stats points: 3151471 faces: 8646263 internal faces: 8378086 cells: 2760772 faces per cell: 6.1665176 boundary patches: 4 point zones: 0 face zones: 0 cell zones: 0 Overall number of cells of each type: hexahedra: 2297184 prisms: 73432 wedges: 0 pyramids: 0 tet wedges: 800 tetrahedra: 2 polyhedra: 389354 Breakdown of polyhedra by number of faces: faces number of cells 4 24349 5 16295 6 21202 7 206911 8 47896 9 50514 10 1042 11 352 12 15902 13 32 14 70 15 4696 16 1 17 1 18 91 Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology inlet 4097 4165 ok (non-closed singly connected) outlet_beifahrer 6283 7226 ok (non-closed singly connected) outlet_fahrer 6262 7211 ok (non-closed singly connected) wall 251535 276060 ok (non-closed singly connected) Checking geometry... Overall domain bounding box (1.315 -0.38629527 -0.10150999) (3.3232987 0.38629978 0.042797109) Mesh has 3 geometric (non-empty/wedge) directions (1 1 1) Mesh has 3 solution (non-empty) directions (1 1 1) Boundary openness (1.0739519e-16 1.9630194e-16 1.0030706e-16) OK. Max cell openness = 1.6036898e-15 OK. Max aspect ratio = 43.785327 OK. Minimum face area = 3.457976e-10. Maximum face area = 1.8677262e-05. Face area magnitudes OK. Min volume = 8.0322233e-13. Max volume = 8.0108398e-08. Total volume = 0.012347979. Cell volumes OK. Mesh non-orthogonality Max: 64.69395 average: 8.0658553 Non-orthogonality check OK. Face pyramids OK. ***Max skewness = 10.13589, 75 highly skew faces detected which may impair the quality of the results <

 July 6, 2016, 06:24 #2 Member   carno Join Date: Mar 2009 Posts: 70 Rep Power: 14 Are you using simpleFoam solution to initialize this? Or using simpleFoam files?

 March 29, 2017, 07:52 #3 Senior Member     Attesz Join Date: Mar 2009 Location: Munich Posts: 367 Rep Power: 14 Hi don6alfonso, did you find a solution to your problem at the end? __________________ I am doing CFD Consulting Services. Ich biete CFD Strömungssimulationen an.

March 29, 2017, 09:21
#4
New Member

A Z
Join Date: Dec 2015
Posts: 16
Rep Power: 7
Quote:
 Originally Posted by Attesz Hi don6alfonso, did you find a solution to your problem at the end?
Yes, I initialized the internalField in the 0/U and 0/p with the solution from potentialFoam. And I change the relaxationFactors for rho to 0.05.

 Tags rhosimplefoam