CFD Online Discussion Forums - --> FOAM FATAL ERROR: Maximum number of iterations exceeded

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - --> FOAM FATAL ERROR: Maximum number of iterations exceeded (https://www.cfd-online.com/Forums/openfoam-solving/105561-foam-fatal-error-maximum-number-iterations-exceeded.html)

Quote:

Originally Posted by slash89 (Post 522649)

Hi all,
i got the same problem. I am using the buoyantSimpleRadiationFoam. The problems is always at the second time step, when solving the G file. Any suggestions to fix this problem?

Thank you,

Best regards

Greetings slash89,

please attach your log file, in this way we can help you much more easily

Here is the log file. I am sorry but i could not upload the file and i don't know why!

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.1.1-221db2718bbb
Exec : buoyantSimpleRadiationFoam -parallel
Date : Dec 05 2014
Time : 18:58:35
Host : "millegradi-nb"
PID : 4112
Case : /home/bolzo/TERMIGNONI/run/prove_solver/prova_5
nProcs : 4
Slaves :
3
(
"millegradi-nb.4113"
"millegradi-nb.4114"
"millegradi-nb.4115"
)

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading g
Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
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;
Prt 1;
}

Calculating field g.h

Reading field p_rgh

Selecting radiationModel P1
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting scatterModel constantScatter

SIMPLE: convergence criteria
field p_rgh tolerance 0.01
field U tolerance 0.001
field h tolerance 0.001
field G tolerance 0.001
field "(k|epsilon|omega)" tolerance 0.001

Starting time loop

Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.00766817, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.00749132, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.00791072, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.00766938, No Iterations 1
DICPCG: Solving for G, Initial residual = 1, Final residual = 0.096472, No Iterations 335
DICPCG: Solving for p_rgh, Initial residual = 0.999948, Final residual = 0.00829887, No Iterations 450
time step continuity errors : sum local = 0.145697, global = -0.0018876, cumulative = -0.0018876
rho max/min : 79.6338 1.1739
DILUPBiCG: Solving for epsilon, Initial residual = 0.12008, Final residual = 0.00390459, No Iterations 1
bounding epsilon, min: -5.64769 max: 1881.23 average: 13.0219
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.0865008, No Iterations 1
ExecutionTime = 23.58 s ClockTime = 23 s

Time = 2

DILUPBiCG: Solving for Ux, Initial residual = 0.998573, Final residual = 0.0137559, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.999799, Final residual = 0.0141418, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.999881, Final residual = 0.0139621, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.013323, No Iterations 1

Thak you!

Are you sure about the BCs you exerted on your geometry? and I wanted the log file not the fist two iteration log ;)

it seems that your cumulative time step error will grow up and the range of rho variation I think is very much!!!

Code:

time step continuity errors : sum local = 0.145697, global = -0.0018876, cumulative = -0.0018876

rho max/min : 79.6338 1.1739

if you're sure about the BCs, then you should check your schemes.

Sorry but i don't understand which file you need. Do you need the G file?

Best regards

just run your case and get the log file:

Code:

solverName >log

and the errors appeared in your terminal.
such as: http://www.cfd-online.com/Forums/ope...tml#post375086

Code:

--> FOAM FATAL ERROR: 

[2] Maximum number of iterations exceeded

[2] 

[2]     From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const

[2]     in file /home/opencfd/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.

[2] 

FOAM parallel run aborting

[2] 

[2] #0  Foam::error::printStack(Foam::Ostream&)[0] 

[0] 

[0] --> FOAM FATAL ERROR: 

[0] Maximum number of iterations exceeded

[0] 

[0]     From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const

[0]     in file /home/opencfd/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.

[0] 

FOAM parallel run aborting

[0] 

[0] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

[2] #1  Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

[0] #1  Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

[2] #2  Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::T(double, double, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

[0] #2  Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::T(double, double, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const, double (Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> >::*)(double) const) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"

[2] #3  Foam::hPsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"

[0] #3  Foam::hPsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::calculate() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"

[0] #4   in "/opt/openfoam2Foam::hPsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct()11/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"

[2] #4  Foam::hPsiThermo<Foam::pureMixture<Foam::constTransport<Foam::specieThermo<Foam::hConstThermo<Foam::perfectGas> > > > >::correct() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"

[0] #5   in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"

[2] #5  



[0]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantSimpleRadiationFoam"

[0] #6  __libc_start_main[2]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantSimpleRadiationFoam"

[2] #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"

[0] #7   in "/lib/x86_64-linux-gnu/libc.so.6"

[2] #7  



[0]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantSimpleRadiationFoam"

--------------------------------------------------------------------------

MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD 

with errorcode 1.



NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.

You may or may not see output from other processes, depending on

exactly when Open MPI kills them.

--------------------------------------------------------------------------

[2]  in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantSimpleRadiationFoam"

--------------------------------------------------------------------------

mpirun has exited due to process rank 0 with PID 2812 on

node millegradi-nb exiting without calling "finalize". This may

have caused other processes in the application to be

terminated by signals sent by mpirun (as reported here).

--------------------------------------------------------------------------

[millegradi-nb:02811] 1 more process has sent help message help-mpi-api.txt / mpi-abort

[millegradi-nb:02811] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

Hope this is what we need! :)

The problem is probably related to the pressure field. I guess you have negative pressure, so the thermophysical library crashes.

check the following:

check your grids, checkmesh
use 1st order scheme for temporal discretization
modify the tolerance tol_ from 1.0e-4 to something higher e.g. 1.0e-3, but this is not a good advice. thanks to dmoroian
at last, modify the maxIter_ from 100 to something larger. you can find the implementation process in http://www.cfd-online.com/Forums/ope...tml#post179437

hope these help you

It still does not work. At this point the problem is in the BCs. Can I upload them here?

Quote:

Originally Posted by slash89 (Post 522710)

It still does not work. At this point the problem is in the BCs. Can I upload them here?

it's the best thing you can do! ;)

Attachment 35793

Thak you!

Maximum number of iterations exceeded

hi

i have this problem
can anyone help me?
if i change the mesh ,solve?
or if i remove the energy equation

best regards

masoud

--> FOAM FATAL ERROR:
Maximum number of iterations exceeded

From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const
in file /home/opencfd/OpenFOAM/OpenFOAM-2.0.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 67.

Need a lot more info

hi Derek

what do you want to know?
did you see this problem later?

best regards
masoud

Hi

I got the sam problem,Can anyone solve it?
I have this problem in my project,I do anything such as mesh,Bc , ... but doesn't work.
if I find anything tell here ,please help me if the problem solve

Best Regards
Masoud

Quote:

Originally Posted by masoudsh (Post 596939)

Hi Masoud,
Take a look at this:
http://openfoam.ir/questions/questio...AE%D8%B7%D8%A7

Best.