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

Error with buoyantSimpleFoam, incompressible properties

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 10, 2015, 11:03
Default Error with buoyantSimpleFoam, incompressible properties
  #1
Senior Member
 
Agustín Villa
Join Date: Apr 2013
Location: Fuenlabrada
Posts: 290
Rep Power: 13
agustinvo is on a distinguished road
Hello

I want to simulate the benchmark of Ercoftac of a natural convection boundary layer (the one of Tsuji and Nagano 1988).

http://www.ercoftac.org/fileadmin/us...5.2/index.html

I am running this simulation with buoyantSimpleFoam, with the same boundaries as specified, by using the data they give as inlet values, except for them that are not available. I am running the case with polynomials for rho, Cp, kappa and viscosity.


When I run the case I get this log:

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.1-262087cdf8db
Exec   : buoyantSimpleFoam
Date   : Nov 10 2015
Time   : 16:52:40
Host   : "tulin23.private.vki.eu"
PID    : 26044
Case   : /nobackup/phd/villaortiz/PhD/benchmark/testCase_11
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Reading thermophysical properties

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       polynomial;
    thermo          hPolynomial;
    equationOfState icoPolynomial;
    specie          specie;
    energy          sensibleEnthalpy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
    Prt             1;
}

Calculating field g.h

Reading field p_rgh

No finite volume options present

Radiation model not active: radiationProperties not found
Selecting radiationModel none

SIMPLE: no convergence criteria found. Calculations will run for 12000 steps.


Starting time loop

wallShearStress wallShearStress1:
    processing wall patches: 
        heatedWall

Time = 1

DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 2.769719e-14, No Iterations 3
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 1.259212e-12, No Iterations 5
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 7.037925e-14, No Iterations 2
GAMG:  Solving for p_rgh, Initial residual = 0.9997015, Final residual = 0.00113553, No Iterations 10
time step continuity errors : sum local = 987.2223, global = -182.4999, cumulative = -182.4999
rho max/min : 1.22799 1.060341
DILUPBiCG:  Solving for epsilon, Initial residual = 0.9999087, Final residual = 1.175033e-11, No Iterations 8
DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 2.276777e-13, No Iterations 2
ExecutionTime = 5.75 s  ClockTime = 6 s

wallShearStress wallShearStress1 output:
    min/max(heatedWall) = (-1.905435 -0.01516225 0), (-0.000272035 395.6815 0)
wallShearStress wallShearStress1 output:
    writing field wallShearStress

calcFvcGrad wallGradT1 output:
    writing field wallGradT

Time = 2

DILUPBiCG:  Solving for Ux, Initial residual = 0.9982672, Final residual = 6.266588e-12, No Iterations 6
DILUPBiCG:  Solving for Uy, Initial residual = 0.03577351, Final residual = 6.744452e-12, No Iterations 5
DILUPBiCG:  Solving for h, Initial residual = 0.999311, Final residual = 2.060703e-12, No Iterations 6


--> 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
I know the error is related with my thermophysical model, but by using the same fluid properties in other cases, it runs. I attach the case in case you want to take a look.

Regards,
Agustin
Attached Files
File Type: gz testCase_11.tar.gz (93.6 KB, 22 views)
agustinvo is offline   Reply With Quote

Reply

Tags
buoyantsimplefoam, natural convection, polynomial properties

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
buoyantPimpleFoam and incompressible fluid (temperature dependent properties) agustinvo OpenFOAM Running, Solving & CFD 2 May 26, 2016 17:33
Effect of downstream flow properties on upstream in incompressible solvers Nick R CFX 0 December 3, 2013 20:49
incompressible thermo type (in buoyantSimpleFoam) andrea.pasquali OpenFOAM 8 November 21, 2013 15:10
depending properties to temprature in incompressible flow anijdon OpenFOAM 0 March 23, 2011 11:03
Turbulent flow of two incompressible fluids with the same physical properties ArneBab OpenFOAM Running, Solving & CFD 3 January 5, 2010 05:48


All times are GMT -4. The time now is 14:54.