CFD Online Discussion Forums - Boussinesq equation of state

Hi foamers. I'm trying to simulate a natural convection of air over a hot cylinder in different conditions. For this specific condition, I want to use the Boussinesq equation of state in thermophysicalProperties. However, when I run the buoyantPimpleFoam solver an error appears.

(OBS: do not mind the weird value for mu, it's a specific value for this case.)

thermophysicalProperties:

Code:

thermoType

{

    type            heRhoThermo;

    mixture         pureMixture;

    transport       const;

    thermo          hConst;

    equationOfState Boussinesq;

    specie          specie;

    energy          sensibleEnthalpy;

}



pRef            101325;



mixture

{

    specie

    {

        nMoles          1;

        molWeight       28.9;

    }

     equationOfState

    {

        rho0            1.13644;

        T0              307.5;

        beta            3.252e-03;

    }

    thermodynamics

    {

        Cp              1007.3;

        Hf              0;

    }

    transport

    {

        mu              0.001571;

        Pr              0.70595;

    }

}

p_rgh:

Code:

internalField   uniform 101325;



boundaryField

{

    face_superior

    {

        type            totalPressure;

        p0              $internalField;

        U               U;

        phi             phi;

        rho             rho_0;

        psi             none;

        gamma           1.4;

        value           $internalField;

    }



    face_inferior

    {

        type            totalPressure;

        p0              $internalField;

        U               U;

        phi             phi;

        rho             rho_0;

        psi             none;

        gamma           1.4;

        value           $internalField;

    }



    entrada

    {

        type            totalPressure;

        p0              $internalField

        U               U;

        phi             phi;

        rho             rho_0;

        psi             none;

        gamma           1.4;

        value           $internalField;

    }



    saida

    {

        type            fixedFluxPressure;

        rho             rho_0;

        value           $internalField;

    }



    parede_cilindro

    {

        type            fixedFluxPressure;

        rho             rho_0;

        value           $internalField;

    }



    frente_traseira

    {

        type            empty;

    }

 }

The error that happens is as follows:
[/CODE]
--> FOAM FATAL ERROR: Maximum number of iterations exceeded 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::Boussinesq<Foam::specie> >; Type = Foam::sensibleEnthalpy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::Bou ssinesq<Foam::specie> >, Foam::sensibleEnthalpy>] in file /build/openfoam-8NJ_iy/openfoam-4.1+dfsg1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.FOAM aborting#0 Foam::error::printStack(Foam::Ostream&) at ??:?#1 Foam::error::abort() at ??:?#2 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam::Boussinesq<Foa m::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?#3 Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam::Boussinesq<Foa m::specie> >, Foam::sensibleEnthalpy> > > >::correct() at ??:?#4 ? in "/usr/bin/buoyantPimpleFoam"#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"#6 ? in "/usr/bin/buoyantPimpleFoam"Aborted (core dumped)
[/CODE]