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

buoyantPimpleFoam and incompressible fluid (temperature dependent properties)

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

Like Tree1Likes
  • 1 Post By agustinvo

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 18, 2015, 11:36
Default buoyantPimpleFoam and incompressible fluid (temperature dependent properties)
  #1
Senior Member
 
Agustín Villa
Join Date: Apr 2013
Location: Fuenlabrada
Posts: 257
Rep Power: 11
agustinvo is on a distinguished road
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:

Code:
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:
Code:
--> 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 is offline   Reply With Quote

Old   September 22, 2015, 11:56
Default
  #2
Senior Member
 
Agustín Villa
Join Date: Apr 2013
Location: Fuenlabrada
Posts: 257
Rep Power: 11
agustinvo is on a distinguished road
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

Code:
dpdt

I will update if there is something new.
Ramzy1990 likes this.
agustinvo is offline   Reply With Quote

Old   May 26, 2016, 18:33
Default Numerical difficulties with incompressible liquid and buoyantPimpleFoam
  #3
New Member
 
Vilma62
Join Date: Nov 2014
Posts: 9
Rep Power: 8
Vilma62 is on a distinguished road
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.
Vilma62 is offline   Reply With Quote

Reply

Tags
buoyantpimplefoam, dpdt, incompressible, polynomial, thermophysicalproperties

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
1.7.x -> buoyantPimpleFoam -> hRhoThermo -> incompressible and icoPoly3ThermoPhysics? will.logie OpenFOAM Programming & Development 1 February 16, 2011 21:52
1.7.x -> buoyantPimpleFoam -> hRhoThermo -> incompressible and icoPoly3ThermoPhysics? will.logie OpenFOAM 0 December 16, 2010 08:08


All times are GMT -4. The time now is 18:48.