[dav.dap83]

Hi everyone,

I am trying to simulate convection into a water-filled horizontal cylinder. One end is kept at a higher temperature than the other boundaries. The solver gives the following error:

--> 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::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::specieThermo<Foam::hPolynomialThermo<Foam::icoPolynomial<8>, 8> >::T(double, double, double (Foam::specieThermo<Foam::hPolynomialThermo<Foam::icoPolynomial<8>, 8> >::*)(double) const, double (Foam::specieThermo<Foam::hPolynomialThermo<Foam::icoPolynomial<8>, 8> >::*)(double) const, double (Foam::specieThermo<Foam::hPolynomialThermo<Foam::icoPolynomial<8>, 8> >::*)(double) const) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#3  Foam::hRhoThermo<Foam::pureMixture<Foam::polynomialTransport<Foam::specieThermo<Foam::hPolynomialThermo<Foam::icoPolynomial<8>, 8> >, 8> > >::calculate() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#4  Foam::hRhoThermo<Foam::pureMixture<Foam::polynomialTransport<Foam::specieThermo<Foam::hPolynomialThermo<Foam::icoPolynomial<8>, 8> >, 8> > >::correct() in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libbasicThermophysicalModels.so"
#5  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantPimpleFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  
 in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/buoyantPimpleFoam"
Aborted (core dumped)

My thermophysicalProperties file looks like:

thermoType hRhoThermo<pureMixture<icoPoly8ThermoPhysics>>;

pRef            100000;

mixture
{
    specie
    {
        nMoles          1;
        molWeight       1;
    }

    equationOfState
    {
        rhoCoeffs<8>   ( 1000 0 0 0 0 0 0 0);
    }
    thermodynamics
    {
        Hf              0;
        Sf              0;
        CpCoeffs<8>     ( 4186 0 0 0 0 0 0 0);
    }
    transport
    {
        Pr              7.01;
        muCoeffs<8>
           ( 1.002e-3 0 0 0 0 0 0 0);
        kappaCoeffs<8>
           ( 0.596 0 0 0 0 0 0 0);
    }
}

Thank you in advance.