Error for multiphaseInterFOAM (RASModel)
 

Dear OF_Users!

Has anyone setup the dambrake4phase tutorial as RAS (kepsilon) instead of Laminar? I tried but failed.

I did follow actions:

- Change in turbulanceProperties to RASModel
- Created a RASProperties-file with kepsilon
- Introduced BC for k, epsilon and nut

I get following messages:

Create time

Create mesh for time = 0


Reading g
Reading field p
Reading field U
Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
#0 Foam::error::printStack(Foam::Ostream&) in "/home/leorrc/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/leorrc/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 Foam::incompressible::RASModels::nutWallFunctionFv PatchScalarField::calcNut() const in
"/home/leorrc/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#4 Foam::incompressible::RASModels::nutWallFunctionFv PatchScalarField::updateCoeffs() in
"/home/leorrc/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#5 Foam::fvPatchField<double>::evaluate(Foam::Pstream ::commsTypes) in
"/home/leorrc/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/multiphaseInterFoam"
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in
"/home/leorrc/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/multiphaseInterFoam"
#7 Foam::incompressible::RASModels::kEpsilon::kEpsilo n(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField,
Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&)
in "/home/leorrc/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#8 Foam::incompressible::RASModel::adddictionaryConst ructorToTable<Foam::incompressible::RASModels::kEp silon>
::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,
Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in
"/home/leorrc/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#9 Foam::incompressible::RASModel::New(Foam::Geometri cField<Foam::Vector<double>, Foam::fvPatchField,
Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&,
Foam::transportModel&) in "/home/leorrc/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#10 Foam::incompressible::turbulenceModel::addturbulen ceModelConstructorToTable<Foam::incompressible::RA SModel>
::NewturbulenceModel(Foam::GeometricField<Foam::Ve ctor<double>, Foam::fvPatchField, Foam::volMesh> const&,
Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&) in
"/home/leorrc/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleRASModels.so"
#11 Foam::incompressible::turbulenceModel::New(Foam::G eometricField<Foam::Vector<double>, Foam::fvPatchField,
Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&,
Foam::transportModel&) in "/home/leorrc/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libincompressibleTurbulenceModel.so"
#12 main in "/home/leorrc/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/multiphaseInterFoam"
#13 __libc_start_main in "/lib64/libc.so.6"
#14 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception


Could you please advise what to do?

Thanks,
Mike

I am willing to bet quite alot that the initial value of epsilon (and/or bc) is zero.

Hello Niklas!

I read this possible failure in an other threat. Therefore I checked k and epsilon with paraview and they are definitely not zero.

nice greetings

Mike

Hello OF_User,
Paraview is usually printing node values. What I think Niklas meant by check your k/epsilon values is to look inside the boundary condition dictionaries and see if the values are different from zero.

Dragos

Hi Dragos!

Sorry, if i was not so clear in my statement. I checked the bc-files and did a cross check in paraview, if all the parameters are ok. But I found no mistake.

Mike

Hello Mike,
Could you post the k/epsilon dictionaries?

Dragos

Hello Dragos!

Here are my k/epsilon dictionaries.

Thanks for your help

Mike

k:

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

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

internalField uniform 0.01;

boundaryField
{
leftWall
{
type kqRWallFunction;
value uniform 0.01;
}
rightWall
{
type kqRWallFunction;
value uniform 0.01;
}
lowerWall
{
type kqRWallFunction;
value uniform 0.01;
}
atmosphere
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}
defaultFaces
{
type empty;
}
}

epsilon:

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

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

internalField uniform 0.01;

boundaryField
{
leftWall
{
type epsilonWallFunction;
value uniform 0.01;
}
rightWall
{
type epsilonWallFunction;
value uniform 0.01;
}
lowerWall
{
type epsilonWallFunction;
value uniform 0.01;
}
atmosphere
{
type inletOutlet;
inletValue uniform 0.1;
value uniform 0.1;
}
defaultFaces
{
type empty;
}
}

Hello Mike,
Indeed you were right, and there are no zero values set, but at least the epsilon values look peculiar. An estimation of both k and epsilon is presented in the documentation http://www.opencfd.co.uk/openfoam/do...tml#x5-40002.1 (eq. 2.8 and 2.9) as well as in any CFD book.
Another thing that cought my attention was the "empty" condition, which means that you have a 3D domain with only one cell thickness. Is this true?

Dragos

Yes, Dragos, you are right, I use a 2D-case.

The starting conditions for k and epsilon are difficult. Theoretically in this case they should be 0, because at the beginning, there is no movement. But since I read about the problems, if they are 0, I used the mentioned values.

Mike

Hmm, and it does that from the first iteration.
Well, the only suggestion I have is to follow equation 2.9 in the documentation and specify a value of for epsilon. Where ...sometimes it helps

Dragos

I also have ever tried to run it in RAS model, and met the same problem.

Best regards,
Jiejin Cai