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

From function specieThermo<thermo>::T(scalar f, scalar T0, scalar ... Error

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

Like Tree1Likes
  • 1 Post By kalle

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 21, 2011, 08:24
Default From function specieThermo<thermo>::T(scalar f, scalar T0, scalar ... Error
  #1
Member
 
Join Date: Nov 2010
Posts: 86
Rep Power: 15
abminternet is on a distinguished road
Hi foamers,

Has anyone encountered this error before and succesfully solved it??

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 /opt/software/openfoam/extend/OpenFOAM-1.6-ext/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 73
Maximum number of iterations exceeded. Rescue by HJ

I am getting it in two different cases, one steady state and one transient. I have read some threads and they say lowering some relaxation factors could solve it, so I tried it for the steady state case, and it just goes on for a couple more iterations but then it fails again. For the transient case I get it once and the solver stops. Any suggestions??? Thanks in advance.
abminternet is offline   Reply With Quote

Old   March 8, 2011, 14:41
Default
  #2
New Member
 
Miguel Suarez
Join Date: Aug 2010
Location: El Paso Texas
Posts: 28
Rep Power: 15
masuarez is on a distinguished road
I have the same problem... i'm using chtMultiRegionFoam. If you find out anything, please let me know. I'll keep looking and see what I can find out.

Miguel
masuarez is offline   Reply With Quote

Old   March 29, 2012, 05:57
Default
  #3
New Member
 
Peter Bishop
Join Date: Jan 2012
Posts: 20
Rep Power: 14
PeterBishop is on a distinguished road
I'm experiencing the same problem, any news?
PeterBishop is offline   Reply With Quote

Old   December 12, 2012, 02:30
Default Maximum number of iterations exceeded
  #4
Member
 
Farshad
Join Date: Oct 2010
Posts: 76
Rep Power: 15
Farshad_Noravesh is on a distinguished road
Hi,

I have the same problem in rhoCentralFoam.

Farshad
Farshad_Noravesh is offline   Reply With Quote

Old   December 12, 2012, 05:34
Default
  #5
Senior Member
 
Karl-Johan Nogenmyr
Join Date: Mar 2009
Location: Linköping
Posts: 279
Rep Power: 21
kalle is on a distinguished road
This thread may provide solution(s):

http://www.cfd-online.com/Forums/ope...ombustion.html

If specieThermo cannot invert h(T) to give T iteratively, usually two things can be wron: Either is the h(T) function erroneous, like in the thread above, or your enthalpy equation is diverging in your solver, giving extreme values of h. In the former case, review the thermo input data. In the latter case, you should investigate why you get such values for h.

K
calf.Z likes this.
kalle is offline   Reply With Quote

Old   December 15, 2012, 15:50
Default
  #6
Senior Member
 
Join Date: Nov 2012
Posts: 171
Rep Power: 13
hz283 is on a distinguished road
Hi,

I also had this problem last week. The density is negative, which may be caused by wrong calculation of pressure or enthalpy (energy). But I think the true reason is diverse and case by case. For my simulation, I modified the mesh (remove some cell of bad quality), and the simulation seems good now. Maybe this issue can be caused by not suitable BC specification, and so on.
hz283 is offline   Reply With Quote

Old   January 3, 2019, 23:38
Default Maximum number of iterations exceeded: 1000
  #7
Senior Member
 
Jianrui Zeng
Join Date: May 2018
Location: China
Posts: 157
Rep Power: 7
calf.Z is on a distinguished road
Quote:
Originally Posted by kalle View Post
This thread may provide solution(s):

http://www.cfd-online.com/Forums/ope...ombustion.html

If specieThermo cannot invert h(T) to give T iteratively, usually two things can be wron: Either is the h(T) function erroneous, like in the thread above, or your enthalpy equation is diverging in your solver, giving extreme values of h. In the former case, review the thermo input data. In the latter case, you should investigate why you get such values for h.

K
Thank you for giving me some hints.

According to my simulation, the error appeared during the first iteration:

PIMPLE: iteration 1
DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 0.0720087, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.063316, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 0.0684356, No Iterations 1
DILUPBiCGStab: Solving for h, Initial residual = 1, Final residual = 0.0749308, No Iterations 2

--> FOAM FATAL ERROR:
Maximum number of iterations exceeded: 1000
......
FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::heRhoTabularThermo<Foam::rhoThermo, Foam:ureMixture<Foam::tabularTransport<Foam::spe cies::thermo<Foam::hTabularThermo<Foam::tabularEOS <Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#3 Foam::heRhoTabularThermo<Foam::rhoThermo, Foam:ureMixture<Foam::tabularTransport<Foam::spe cies::thermo<Foam::hTabularThermo<Foam::tabularEOS <Foam::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?


I think this error is related to h and T, but I have no idea where I should check. Maybe there are some mistakes with my code or setting.

Any suggestion is highly appreciated.
calf.Z 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
CGNS lib and Fortran compiler manaliac Main CFD Forum 2 November 29, 2010 07:25
[blockMesh] BlockMesh FOAM warning gaottino OpenFOAM Meshing & Mesh Conversion 7 July 19, 2010 15:11
POSDAT problem piotka STAR-CD 4 June 12, 2009 09:43
[blockMesh] Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Meshing & Mesh Conversion 10 April 2, 2007 15:00
Installation problem with GCC Norma McKee (Mckee) OpenFOAM Installation 10 March 4, 2007 08:09


All times are GMT -4. The time now is 19:54.