[denis08]

Hello,


I am running a case using simpleFoam with a k-epsilon turbulence model, on a submarine stl geometry(named "hull" in my case). However, values of epsilon and k diverge quickly (after 5 iterations). I already check the mesh and tried to change the value of the time step but it doesnt solve my issue.



Has anyone an idea on how to solve my problem?


My files for k, epsilon, P, U and nut are:


k:

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

internalField uniform 0.00015;

boundaryField
{
inlet
{
type fixedValue;
value $internalField;
}
"(hull|top|bottom|side|side2)"
{
type kqRWallFunction;
value $internalField;
}

outlet
{
type zeroGradient;
}
}


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

internalField uniform 0.00025;

boundaryField
{
inlet
{
type fixedValue;
value $internalField;
}

"(hull|top|bottom|side|side2)"
{
type epsilonWallFunction;
value $internalField;
Cmu 0.09;
kappa 0.41;
E 9.8;
}


outlet
{
type zeroGradient;
}
}


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

internalField uniform 5e-07;

boundaryField
{
"(hull|top|bottom|side|side2)"
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value $internalField;
}


inlet
{
type calculated;
value $internalField;
}

outlet
{
type calculated;
value $internalField;
}
}


P:
dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
outlet
{

type fixedValue;
value $internalField;
}

inlet
{
type zeroGradient;
}

top
{
type zeroGradient;
}

bottom
{
type zeroGradient;
}

side
{
type zeroGradient;
}

side2
{
type zeroGradient;
}


hull
{
type zeroGradient;
}

}


U:
dimensions [0 1 -1 0 0 0 0];

internalField uniform (10 0 0);





boundaryField
{
inlet
{
type fixedValue;
value $internalField;
}


outlet
{
type pressureInletOutletVelocity;
value $internalField;
}

"(hull|top|bottom|side|side2)"
{
type fixedValue;
value uniform (0 0 0);
}

}






Thanks for your help