[agustinvo]

Hi

I would like to know if someone of you have already tried to run an incompressible simulation with buoyantPimpleFoam for an incompressible fluid (now I am working with air, but later on I will use other kind of fluids).

I have temperature dependence in all the variables, so my thermophysicalProperties looks like:

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       polynomial;
    thermo          hPolynomial;
    equationOfState  icoPolynomial ;
    specie          specie;
    energy          sensibleEnthalpy;
}
specie
    {
        nMoles          1;
        molWeight       28.85;
    }
    equationOfState
    {
        
        rhoCoeffs<8>    ( 4.0097 -0.016954 3.3057e-05 -3.0042e-08 1.0286e-11 0 0 0 );
    }
    thermodynamics
    {
        Tlow         300;
        Thigh        500;
        Tcommon      300; 
        Hf              0;
        Sf              0;
        CpCoeffs<8>     ( 948.76 0.39171 -0.00095999 1.393e-06 -6.2029e-10 0 0 0 );
    }
    transport
    {
        muCoeffs<8>     ( 1.5061e-06 6.16e-08 -1.819e-11 0 0 0 0 0 );
        kappaCoeffs<8>  ( 0.0025219 8.506e-05 -1.312e-08 0 0 0 0 0 );
    }

It runs ok, but suddenly crashes:

--> 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 /software/alternate/openfoam/fc20/2.3.1/OpenFOAM-2.3.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::THs(double, double, double) const at ??:?
#3  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > >::calculate() at ??:?
#4  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > >::correct() at ??:?
#5  
 at ??:?
#6  __libc_start_main at ??:?
#7  
 at ??:?
Aborted

Do you know what could be happening? The temperature doesn't go to unlogical values, and I have polynomials expressions.

[agustinvo]

Hi everyone,

I am able to run it! I paid attention to some threads here (I will edit with their links), and it seems the problem was with the

dpdt

I will update if there is something new.

[Vilma62]

Hello agustinvo,

I am trying to simulate buoyancy driven oil (incompressible liquid) flow and heating in the closed cavity by using buoyantPimpleFoam solver, and run into numerical difficulties (--> FOAM FATAL ERROR: Maximum number of iterations exceeded).

Could you please elaborate in more detail how you succeed to stabilize buoyantPimpleFoam solver in your simulation (please note; I'm a total beginner)?

Thank you in advance.