nutkRoughWallFunction ERROR
 

Hi Foamers!

I'm working on a simulation using pisoFoam, I'm using the nutkRoughWallFunction with the following parameters:

The roughness height ks is required as the input of most nut wall equations, to assure a rough wall log layer between the center of the first near-bed grid cell and the wall. In our case, the first near-bed grid cell is equal to 5-cm, ks is 10-cm, and zo is approximately 2*ks (20-cm).


FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{
bed
{
type nutURoughWallFunction;
Ks uniform 0.1;
Cs uniform 2;
Cmu 0.09;
kappa 0.41;
E 9.8;
roughnessFactor 1;
value $internalField;
}
topWall
{
type slip;
}
sideRight
{
type slip;
}
sideLeft
{
type slip;
}
inlet
{
type fixedValue;
value uniform 1e-06;
}
out
{
type zeroGradient;
}
}


AND I GOT THIS ERROR:


Create time

Create mesh for time = 0

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type LES
Selecting LES turbulence model SpalartAllmarasDDES
Selecting LES delta type smooth
Selecting LES delta type cubeRootVol
Selecting patchDistMethod meshWave
LES
{
model SpalartAllmarasDDES;
turbulence on;
printCoeffs on;
delta smooth;
dynamicKEqnCoeffs
{
filter simple;
}
cubeRootVolCoeffs
{
deltaCoeff 1;
}
PrandtlCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
maxDeltaRatio 1.1;
}
Cdelta 0.158;
}
vanDriestCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
maxDeltaRatio 1.1;
}
Aplus 26;
Cdelta 0.158;
}
smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}
maxDeltaRatio 1.1;
}
Ce 1.048;
sigmaNut 0.66666;
kappa 0.41;
Cb1 0.1355;
Cb2 0.622;
Cw2 0.3;
Cw3 2;
Cv1 7.1;
Cs 0.3;
CDES 0.65;
ck 0.07;
}

No MRF models present

No finite volume options present
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 logf32x in "/lib/x86_64-linux-gnu/libm.so.6"
#4 Foam::nutURoughWallFunctionFvPatchScalarField::yPl us(Foam::Field<double> const&) const at ??:?
#5 Foam::nutURoughWallFunctionFvPatchScalarField::nut () const at ??:?
#6 Foam::nutWallFunctionFvPatchScalarField::updateCoe ffs() at ??:?
#7 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/pisoFoam"
#8 Foam::LESModels::SpalartAllmarasDES<Foam::Incompre ssibleMomentumTransportModel<Foam::transportModel> >::correctNut() at ??:?
#9 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/pisoFoam"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/pisoFoam"
Floating point exception (core dumped)

Can anyone help me to understand and solve the error, please?

Thanks