[yaqb]

Hi All,

I am trying to run reactingfoam ( a bit modified - with moving mesh). I try to change the reaction scheme from single stage to multi-stage adding chemikin database to the calculations. Everything goes fine until around 0.035 sek when the temperature reach around 1030 K, than the simulation crashes giving the following error massage :

--> 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.1.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 69.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::specieThermo<Foam::janafThermo<Foam:erfect Gas> >::T(double, double, double (Foam::specieThermo<Foam::janafThermo<Foam:erfec tGas> >::*)(double) const, double (Foam::specieThermo<Foam::janafThermo<Foam:erfec tGas> >::*)(double) const, double (Foam::specieThermo<Foam::janafThermo<Foam:erfec tGas> >::*)(double) const) const at ??:?
#3 Foam::ODEChemistryModel<Foam:siChemistryModel, Foam::sutherlandTransport<Foam::specieThermo<Foam: :janafThermo<Foam:erfectGas> > > >::solve(double, double) at ??:?
#4 Foam::combustionModels::PaSR<Foam::combustionModel s:siChemistryCombustionModel>::correct() at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?

I tired to increase the number of allowed iterations in specieThermo.C and even decrease the tolerance limit (to 1e-3) no results.
Does anyone have idea how to deal with this problem?

[Max1234]

Hi,

did you finally find a solution for your problem?

I get the same error. I am using a custom solver with implemented YEqn and EEqn of the reactingFoam Solver of Version 2.2, based on the PaSR Concept.
The Setup is a closed volume with a primixed air - fuel mixture using a quasi-global mechanism (12 species). Turbulence is calculated using the oneEqEddy LES Model. The mixture is ignited using a scalarSemiImplicitSource (fvOptions).

It works fine, and suddenly:

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

From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
in file /home/openfoam/OpenFOAM/OpenFOAM-2.4.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting


Thanks for your help!

[Momo1805]

Hey max,

I'm having exactly the same problem with PaSR model and trying to simulate combustion with ignition. Did you finally solve the problem???

Would be so helpful for me!

Thanks and cheers

Moritz

[laurentD]

Maybe you can post a logFile which will let us to see what is exactly happening.
This is a problem i have faced many times but there is a lot of possible causes...

[Momo1805]

Hey laurentD,

enclosed the logFile of my case. Please not that I use the standard aachenBomb geometry so far. I validated the thermophysical data of my gases. Would be awesome if you could provide some hints where to look for the bug in my simulation.

Thanks mate

Moritz

[laurentD]

It's not sure but whan i look your logFile, i am curious about your boundary conditions ;-)

[hbrist7]

Hi laurentD,

I've been getting the same max number of iteration exceeded problem too with rhoCentralFoam

I think my boundary conditions are set correctly, but here's a dropbox link to my whole case post run. my log file is log.rho and is where im getting the error.

Edit: It's been solved - my units for specific heat were off by a magnitude of 1000. A stupid error, actually.

thanks,

henry

[Momo1805]

Hey Laurent,

I am simulating the injection, evaporation and combustion of an isobutane/Propylene mixture in a one stroke gas engine under ambient pressure. BC are only walls like in the aachenBomb tutorial. Just tried to get all the thermodynamics working so far.

Cheers

Moe