[gushrqcosta]

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:

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:

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:rintStack(Foam::Ostream&) at ??:?#1 Foam::error::abort() at ??:?#2 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam::Boussinesq<Foa m::specie> >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?#3 Foam::heRhoThermo<Foam::rhoThermo, Foam:ureMixture<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]

[gushrqcosta]

UPDATE:
I did not manage to solve this error, what I did instead is just install and use OpenFOAM-10 since it's more developed and it worked just fine.