[Pervispasco]

Hi everyone,

I'm having troubles running a simple 3D RANS flow over a cube.

It works perfectly fine with uniform inlet velocity, but as soon as I use the atmBoundarylayerinlet for U and Epsilon, everything stop working.

I simply copied the two files of the turbine tutorial in simpleFoam directory.

These are mine U and Epsilon files:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/*  Windows port by CFD support (www.cfdsupport.com) [based on Symscape]     *\
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0.01";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    buildings
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    ground
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    frontAndBack
    {
        type            symmetry;
    }

    inlet
    {
            type            atmBoundaryLayerInletVelocity;
                Uref                 10.0;
                Zref                 20;
                zDir                 (0 0 1);
                flowDir              (1 0 0);
                z0                   uniform 0.1;
                zGround              uniform 935.0;
                value                $internalField;

    }
    outlet
    {
        type            zeroGradient;
    }
}


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

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/*  Windows port by CFD support (www.cfdsupport.com) [based on Symscape]     *\
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0.01";
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 9.98e-06;

boundaryField
{

    buildings
    {
        type            epsilonWallFunction;
        value           uniform 9.98e-06;
    }
    ground
    {
        type            epsilonWallFunction;
        value           uniform 9.98e-06;
    }
    frontAndBack
    {
        type            symmetry;
    }
    inlet
    {
     type            atmBoundaryLayerInletEpsilon;
        Uref                 10.0;
        Zref                 20;
        zDir                 (0 0 1);
        flowDir              (1 0 0);   
        z0                   uniform 0.1;
        zGround              uniform 935.0;
        value                $internalField;
    }
    outlet
    {
        type            zeroGradient;
    }
}


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

And this is what I get from OpenFoam when I try to run the case:

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
smoothSolver: Solving for Uz, Initial residual = nan, Final residual = nan, No Iterations 1000
GAMG: Solving for p, Initial residual = nan, Final residual = nan, No Iterations 1000
time step continuity errors : sum local = nan, global = nan, cumulative = nan
smoothSolver: Solving for epsilon, Initial residual = nan, Final residual = nan, No Iterations 1000
smoothSolver: Solving for k, Initial residual = nan, Final residual = nan, No Iterations 1000
ExecutionTime = 386.812 s ClockTime = 386 s

forces forces output:
sum of forces:
pressure : (nan nan nan)
viscous : (nan nan nan)
porous : (0 0 0)
sum of moments:
pressure : (nan nan nan)
viscous : (nan nan nan)
porous : (0 0 0)

forceCoeffs forceCoeffs output:
Cm = nan
Cd = nan
Cl = nan
Cl(f) = nan
Cl(r) = nan

fieldAverage fieldAverage1 output:
Calculating averages

Time = 5


What could be the issue?

Thanks