[agustinvo]

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:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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