CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

FATAL ERROR: Max number of iterations exceeded

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 5, 2016, 16:55
Default FATAL ERROR: Max number of iterations exceeded
  #1
New Member
 
Henry
Join Date: May 2016
Posts: 15
Rep Power: 9
hbrist7 is on a distinguished road
I've been having a lot of trouble tracing the source of this error. I know it happens in the thermoI.H file, but I can't trace what is calling that, or figure out any reason why this is happening.

The case i am running is a simple wedge (meshed with snappy hex mesh) inside a blockmesh region using rhoCentral foam.

any help would be really appreciated

this is the error log file

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 3.0.1-119cac7e8750
Exec   : rhoCentralFoam -parallel
Date   : Jul 01 2016
Time   : 12:54:38
Host   : "cfd01"
PID    : 30359
Case   : /----/wedge15maxIterError
nProcs : 36
Slaves : 
35
(
"cfd01.30360"
"cfd01.30361"
"cfd01.30362"
"cfd01.30363"
"cfd01.30364"
"cfd01.30365"
"cfd01.30366"
"cfd01.30367"
"cfd01.30368"
"cfd01.30369"
"cfd01.30370"
"cfd02.26571"
"cfd02.26572"
"cfd02.26573"
"cfd02.26574"
"cfd02.26575"
"cfd02.26576"
"cfd02.26577"
"cfd02.26578"
"cfd02.26579"
"cfd02.26580"
"cfd02.26581"
"cfd02.26582"
"cfd03.29133"
"cfd03.29134"
"cfd03.29135"
"cfd03.29136"
"cfd03.29137"
"cfd03.29138"
"cfd03.29139"
"cfd03.29140"
"cfd03.29141"
"cfd03.29142"
"cfd03.29143"
"cfd03.29144"
)
Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Creating turbulence model

Selecting turbulence model type laminar
fluxScheme: Kurganov

Starting time loop

Mean and max Courant Numbers = 0.00521015 0.119072
deltaT = 2.4e-06
Time = 2.4e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.999285, Final residual = 0.00389594, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.999999, Final residual = 8.69665e-17, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 8.81879e-17, No Iterations 5
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for e, Initial residual = 1.46977e-11, Final residual = 1.46977e-11, No Iterations 0
ExecutionTime = 0.68 s  ClockTime = 1 s

Mean and max Courant Numbers = 0.00625242 0.142579
deltaT = 2.88e-06
Time = 5.28e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 6.23273e-10, Final residual = 6.23273e-10, No Iterations 0
smoothSolver:  Solving for Uy, Initial residual = 1.24042e-05, Final residual = 6.5351e-18, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.00028774, Final residual = 4.85118e-17, No Iterations 5
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for e, Initial residual = 5.7298e-12, Final residual = 5.7298e-12, No Iterations 0
[0] [1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] Maximum number of iterations exceeded
[1] 
[1]     From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar)const) const [with Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
[1]     in file /home/openfoam/OpenFOAM/OpenFOAM-3.0.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.
[1] 
FOAM parallel run aborting
[1] 
[2] 
[2] 
[2] --> FOAM FATAL ERROR: [3] 
[3] 
[3] --> FOAM FATAL ERROR: 
[3] Maximum number of iterations exceeded
[3] 
[3]     From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar)const) const [with Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
[3]     in file /home/openfoam/OpenFOAM/OpenFOAM-3.0.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66[7] 
[7] 
[7] --> FOAM FATAL ERROR: 
[7] Maximum number of iterations exceeded
[7] 
.
for reference, my P/T/U files are here:

p:

Code:
dimensions      [1 -1 -2 0 0 0 0];
internalField   uniform 18822.7;
boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 18822.7;
    }
    outlet
    {
        type            zeroGradient;
    }
    airplane
    {
	type		zeroGradient;
    }
    bottom
    {
        type            zeroGradient;
    }
    top
    {
        type            zeroGradient;
    }
    symmetry
    {
	type		symmetryPlane;
    }
    left
    {
	type 		zeroGradient;
    }
    defaultFaces
    {
        type            zeroGradient;
    }
}
T:

Code:
dimensions      [0 0 0 1 0 0 0];
internalField   uniform 216.65;
boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 216.65;
    }
    outlet
    {
        type            zeroGradient;
    }
    airplane
    {
	type		zeroGradient;
    }
    bottom
    {
        type            zeroGradient;
    }
    top
    {
        type            zeroGradient;
    }
    left
    {
	type		zeroGradient;
    }
    symmetry
    {
	type		symmetryPlane;
    }
    defaultFaces
    {
        type            zeroGradient;
    }
}
U:

Code:
dimensions      [0 1 -1 0 0 0 0];
internalField   uniform (400 0 0); //(246.933 0 0);
boundaryField
{
    inlet
    {
        type		fixedValue;
        value           uniform (400 0 0); //(246.933 0 0);
    }
    outlet
    {
        type            zeroGradient;
    }
    airplane
    {
	type		slip;
    }
    bottom
    {
        type            zeroGradient; 
    }
    top
    {
        type            zeroGradient; 
    }
    left
    {
	type		zeroGradient;
    }
    symmetry
    {
	type		symmetryPlane;
    }
    defaultFaces
    {
        type            zeroGradient;
    }
}
hbrist7 is offline   Reply With Quote

Old   July 6, 2016, 17:02
Default
  #2
New Member
 
Henry
Join Date: May 2016
Posts: 15
Rep Power: 9
hbrist7 is on a distinguished road
Also, I changed the boundary conditions to supersonicFreestream and waveTransmissive in hopes that that would help, but i am getting the same problem.

I have my entire post-run case directory on dropbox, with the error located in the log.rho file

https://www.dropbox.com/sh/7fvzzazmx...liXvZj0ga?dl=0

Does anyone have any ideas what the problem is?
hbrist7 is offline   Reply With Quote

Old   March 15, 2019, 09:15
Default
  #3
New Member
 
Hao Li
Join Date: May 2017
Posts: 7
Rep Power: 8
haolee0928 is on a distinguished road
Quote:
Originally Posted by hbrist7 View Post
I've been having a lot of trouble tracing the source of this error. I know it happens in the thermoI.H file, but I can't trace what is calling that, or figure out any reason why this is happening.

The case i am running is a simple wedge (meshed with snappy hex mesh) inside a blockmesh region using rhoCentral foam.

any help would be really appreciated

this is the error log file

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 3.0.1-119cac7e8750
Exec   : rhoCentralFoam -parallel
Date   : Jul 01 2016
Time   : 12:54:38
Host   : "cfd01"
PID    : 30359
Case   : /----/wedge15maxIterError
nProcs : 36
Slaves : 
35
(
"cfd01.30360"
"cfd01.30361"
"cfd01.30362"
"cfd01.30363"
"cfd01.30364"
"cfd01.30365"
"cfd01.30366"
"cfd01.30367"
"cfd01.30368"
"cfd01.30369"
"cfd01.30370"
"cfd02.26571"
"cfd02.26572"
"cfd02.26573"
"cfd02.26574"
"cfd02.26575"
"cfd02.26576"
"cfd02.26577"
"cfd02.26578"
"cfd02.26579"
"cfd02.26580"
"cfd02.26581"
"cfd02.26582"
"cfd03.29133"
"cfd03.29134"
"cfd03.29135"
"cfd03.29136"
"cfd03.29137"
"cfd03.29138"
"cfd03.29139"
"cfd03.29140"
"cfd03.29141"
"cfd03.29142"
"cfd03.29143"
"cfd03.29144"
)
Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Creating turbulence model

Selecting turbulence model type laminar
fluxScheme: Kurganov

Starting time loop

Mean and max Courant Numbers = 0.00521015 0.119072
deltaT = 2.4e-06
Time = 2.4e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.999285, Final residual = 0.00389594, No Iterations 5
smoothSolver:  Solving for Uy, Initial residual = 0.999999, Final residual = 8.69665e-17, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 8.81879e-17, No Iterations 5
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for e, Initial residual = 1.46977e-11, Final residual = 1.46977e-11, No Iterations 0
ExecutionTime = 0.68 s  ClockTime = 1 s

Mean and max Courant Numbers = 0.00625242 0.142579
deltaT = 2.88e-06
Time = 5.28e-06

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal:  Solving for rhoUz, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 6.23273e-10, Final residual = 6.23273e-10, No Iterations 0
smoothSolver:  Solving for Uy, Initial residual = 1.24042e-05, Final residual = 6.5351e-18, No Iterations 5
smoothSolver:  Solving for Uz, Initial residual = 0.00028774, Final residual = 4.85118e-17, No Iterations 5
diagonal:  Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for e, Initial residual = 5.7298e-12, Final residual = 5.7298e-12, No Iterations 0
[0] [1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] Maximum number of iterations exceeded
[1] 
[1]     From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar)const) const [with Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
[1]     in file /home/openfoam/OpenFOAM/OpenFOAM-3.0.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.
[1] 
FOAM parallel run aborting
[1] 
[2] 
[2] 
[2] --> FOAM FATAL ERROR: [3] 
[3] 
[3] --> FOAM FATAL ERROR: 
[3] Maximum number of iterations exceeded
[3] 
[3]     From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar)const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar)const) const [with Thermo = Foam::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
[3]     in file /home/openfoam/OpenFOAM/OpenFOAM-3.0.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66[7] 
[7] 
[7] --> FOAM FATAL ERROR: 
[7] Maximum number of iterations exceeded
[7] 
.
for reference, my P/T/U files are here:

p:

Code:
dimensions      [1 -1 -2 0 0 0 0];
internalField   uniform 18822.7;
boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 18822.7;
    }
    outlet
    {
        type            zeroGradient;
    }
    airplane
    {
	type		zeroGradient;
    }
    bottom
    {
        type            zeroGradient;
    }
    top
    {
        type            zeroGradient;
    }
    symmetry
    {
	type		symmetryPlane;
    }
    left
    {
	type 		zeroGradient;
    }
    defaultFaces
    {
        type            zeroGradient;
    }
}
T:

Code:
dimensions      [0 0 0 1 0 0 0];
internalField   uniform 216.65;
boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 216.65;
    }
    outlet
    {
        type            zeroGradient;
    }
    airplane
    {
	type		zeroGradient;
    }
    bottom
    {
        type            zeroGradient;
    }
    top
    {
        type            zeroGradient;
    }
    left
    {
	type		zeroGradient;
    }
    symmetry
    {
	type		symmetryPlane;
    }
    defaultFaces
    {
        type            zeroGradient;
    }
}
U:

Code:
dimensions      [0 1 -1 0 0 0 0];
internalField   uniform (400 0 0); //(246.933 0 0);
boundaryField
{
    inlet
    {
        type		fixedValue;
        value           uniform (400 0 0); //(246.933 0 0);
    }
    outlet
    {
        type            zeroGradient;
    }
    airplane
    {
	type		slip;
    }
    bottom
    {
        type            zeroGradient; 
    }
    top
    {
        type            zeroGradient; 
    }
    left
    {
	type		zeroGradient;
    }
    symmetry
    {
	type		symmetryPlane;
    }
    defaultFaces
    {
        type            zeroGradient;
    }
}
Hi, I met the nearly same trouble. Could you tell me how to fix it. Just a roadmap is ok
haolee0928 is offline   Reply With Quote

Old   May 13, 2022, 15:55
Default
  #4
New Member
 
Fuad Hasibul Hasan
Join Date: Feb 2022
Location: Dhaka, Bangladesh
Posts: 7
Rep Power: 4
fuad is on a distinguished road
Hi,
Have you solved this issue? I'm facing a similar error. I've been working on this error for more than a week but could not solve this yet. Please share if you did.
In my observation after trying a lot of things is that they are may be a unit problem of the initial and boundary conditions and thermophysical properties. The tutorial cases use the normalized parameters and as soon as I change them to S.I. units like here it causes this error.
fuad is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
chtMultiRegionSimpleFoam: maximum number of iterations excedeed. Nkl OpenFOAM Running, Solving & CFD 19 October 10, 2019 03:42
Floating point exception error lpz_michele OpenFOAM Running, Solving & CFD 53 October 19, 2015 03:50
[OpenFOAM.org] OF2.3.1 + OS13.2 - Trying to use the dummy Pstream library aylalisa OpenFOAM Installation 23 June 15, 2015 15:49
SLTS+rhoPisoFoam: what is rDeltaT??? nileshjrane OpenFOAM Running, Solving & CFD 4 February 25, 2013 05:13
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 12:16


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