[ss32]

I have flow over a cylinder, simple case to play with BCs and turbulence models. Everything is good for laminar and solving turbulence with SpalartAllmaras but trying to implement K-Epsilon diverges when trying to solve for epsilon. At the first time step it diverges and I get a floating point exception. Is there an issue with my boundary conditions for epsilon?

Create time

Create mesh for time = 0


SIMPLE: convergence criteria
    field p     tolerance 1e-05
    field U     tolerance 1e-05
    field nuTilda     tolerance 1e-05

Reading field p

Reading field U

Reading/calculating face flux field phi

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

No MRF models present

No finite volume options present


Starting time loop

Time = 0.001

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 5.8534e-09, No Iterations 50
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 9.37619e-09, No Iterations 44
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0899236, No Iterations 7
GAMG:  Solving for p, Initial residual = 0.098076, Final residual = 0.00642081, No Iterations 1
time step continuity errors : sum local = 0.00354511, global = -1.46424e-05, cumulative = -1.46424e-05
smoothSolver:  Solving for epsilon, Initial residual = 0.995688, Final residual = nan, No Iterations 1000
smoothSolver:  Solving for k, Initial residual = nan, Final residual = nan, No Iterations 1000
ExecutionTime = 1.83 s  ClockTime = 2 s

Time = 0.002

smoothSolver:  Solving for Ux, Initial residual = nan, Final residual = nan, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = nan, Final residual = nan, No Iterations 1000
GAMG:  Solving for p, Initial residual = nan, Final residual = nan, No Iterations 1000
GAMG:  Solving for p, Initial residual = nan, Final residual = nan, No Iterations 1000


--> FOAM FATAL IO ERROR: 
wrong token type - expected Scalar, found on line 0 the word 'nan'

file: /home/scott/OpenFOAM/scott-v3.0+/run/tutorials/incompressible/simpleFoam/cylinder/system/data.solverPerformance.p at line 0.

    From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::doubleScalar&)
    in file lnInclude/Scalar.C at line 93.

FOAM exiting

k:

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0.008;

boundaryField
{
    frontback
    {
        type            empty;
    }
    inlet
    {
        type            fixedValue;
        value           uniform 0.008;
    }
    outlet
    {
        type            zeroGradient;
    }
    topbottom
    {
        type            zeroGradient;
    }
    walls
    {
        type            kqRWallFunction;
        value           uniform 0.008;
    }
}


// ************************************************************************* //

epsilon - my calculated epsilon is much lower but reading through old posts some people suggested raising it to fix the error but no such luck for me.

FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -3 0 0 0 0];

internalField   uniform 0.5;

boundaryField
{
    frontback
    {
        type            empty;
    }
    inlet
    {
        type            fixedValue;
        value           uniform 0.5;
    }
    outlet
    {
        type            zeroGradient;
    }
    topbottom
    {
        type            zeroGradient;
    }
    walls
    {
        type            epsilonWallFunction;
        value           uniform 0.5;
    }
}


// ************************************************************************* //