RunTim Error for simpleFoam

Djub · April 9, 2013, 10:18

Hi,
I am trying to make a Foam model ,using simpleFoam, but Foam complains with an error that I don't understand. I would like to correct myself but the error message is too complicated for me...
Here is the error (I am running in parallel with 14 threads, on a machine named "titan"):

Code:

[12] #0  Foam::error::printStack(Foam::Ostream&)[4] #0[1] #0  Foam::error::printStack(Foam::Ostream&)  Foam::error::printStack(Foam::Ostream&)[0] #0  Foam::error::printStack(Foam::Ostream&)[8] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[12] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[4] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[4] #2   in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #2   in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2  DPOpt/lib/l in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #2  ibOpenFOAM.so"
[12] #2   in "/lib64/libc.so.6"
[4] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib64/libc.so.6"
[1] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib64/libc.so.6"
[0] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib64/libc.so.6"
[12] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib64/libc.so.6"
[8] #3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #4  void Foam::divide<Foam::fvPatchField>(Foam::FieldField<Foam::fvPatchField, double>&, Foam::FieldField<Foam::fvPatchField, double> const&, Foam::FieldField<Foam::fvPatchField, double> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[4] #4  void Foam::divide<Foam::fvPatchField>(Foam::FieldField<Foam::fvPatchField, double>&, Foam::FieldField<Foam::fvPatchField, double> const&, Foam::FieldField<Foam::fvPatchField, double> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[12] #4  void Foam::divide<Foam::fvPatchField>(Foam::FieldField<Foam::fvPatchField, double>&, Foam::FieldField<Foam::fvPatchField, double> const&, Foam::FieldField<Foam::fvPatchField, double> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[1] #4  void Foam::divide<Foam::fvPatchField>(Foam::FieldField<Foam::fvPatchField, double>&, Foam::FieldField<Foam::fvPatchField, double> const&, Foam::FieldField<Foam::fvPatchField, double> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[8] #4  void Foam::divide<Foam::fvPatchField>(Foam::FieldField<Foam::fvPatchField, double>&, Foam::FieldField<Foam::fvPatchField, double> const&, Foam::FieldField<Foam::fvPatchField, double> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[0] #5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[4] #5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[8] #5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[1] #5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[12] #5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[4] #6  Foam::incompressible::RASModels::kEpsilon::kEpsilon(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[0] #6  Foam::incompressible::RASModels::kEpsilon::kEpsilon(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[4] #7   in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
Foam::incompressible::RASModel::adddictionaryConstructorToTable<Foam::incompressible::RASModels::kEpsilon>::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&)[8] #6  Foam::incompressible::RASModels::kEpsilon::kEpsilon(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[0] #7  Foam::incompressible::RASModel::adddictionaryConstructorToTable<Foam::incompressible::RASModels::kEpsilon>::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[1] #6   in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[12] #6  Foam::incompressible::RASModels::kEpsilon::kEpsilon(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&)Foam::incompressible::RASModels::kEpsilon::kEpsilon(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[0] #8  Foam::incompressible::RASModel::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[12] #7  Foam::incompressible::RASModel::adddictionaryConstructorToTable<Foam::incompressible::RASModels::kEpsilon>::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[8] #7  Foam::incompressible::RASModel::adddictionaryConstructorToTable<Foam::incompressible::RASModels::kEpsilon>::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[4] #8  Foam::incompressible::RASModel::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[1] #7  Foam::incompressible::RASModel::adddictionaryConstructorToTable<Foam::incompressible::RASModels::kEpsilon>::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[12] #8  Foam::incompressible::RASModel::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[0] #9   in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[8] #8  Foam::incompressible::RASModel::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&)
 in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[12] #9   in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[1] #8  Foam::incompressible::RASModel::New(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::transportModel&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[4] #9   in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[8] #9  [0]  in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/simpleFoam"
[0] #10  __libc_start_main

 in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
[1] #9  [12]  in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/simpleFoam"
[12] #10  __libc_start_main[8]  in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/simpleFoam"
[8] #10  __libc_start_main[4]  in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/simpleFoam"
[4] #10  __libc_start_main in "/lib64/libc.so.6"
[0] #11  
 in "/lib64/libc.so.6"
[12] #11  [0]  at /home/abuild/rpmbuild/BUILD/glibc-2.15/csu/../sysdeps/x86_64/elf/start.S:116
[titan:22422] *** Process received signal ***
[titan:22422] Signal: Floating point exception (8)
[titan:22422] Signal code:  (-6)
[titan:22422] Failing at address: 0x3ea00005796
 in "/lib64/libc.so.6"
[8] #11   in "/lib64/libc.so.6"
[4] #11  [1]  in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/bin/simpleFoam"
[1] #10  __libc_start_main[titan:22422] [ 0] /lib64/libc.so.6(+0x34da0) [0x7fd56d193da0]
[titan:22422] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x7fd56d193d25]
[titan:22422] [ 2] /lib64/libc.so.6(+0x34da0) [0x7fd56d193da0]
[titan:22422] [ 3] /opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam6divideERNS_5FieldIdEERKNS_5UListIdEES6_+0xd5) [0x7fd56e304395]
[titan:22422] [ 4] /opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam6divideINS_12fvPatchFieldEEEvRNS_10FieldFieldIT_dEERKS4_S7_+0x6c) [0x7fd57008be0c]
[titan:22422] [ 5] /opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4FoamdvINS_12fvPatchFieldENS_7volMeshEEENS_3tmpINS_14GeometricFieldIdT_T0_EEEERKS8_RKS7_+0x178) [0x7fd570098728]
[titan:22422] [ 6] /opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompressible9RASModels8kEpsilonC1ERKNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS3_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERNS_14transportModelERKNS_4wordESK_+0x8e9) [0x7fd570085199]
[titan:22422] [ 7] /opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompressible8RASModel31adddictionaryConstructorToTableINS0_9RASModels8kEpsilonEE3NewERKNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS6_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERNS_14transportModelERKNS_4wordE+0x58) [0x7fd570099398]
[titan:22422] [ 8] /opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so(_ZN4Foam14incompressible8RASModel3NewERKNS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldENS_7volMeshEEERKNS2_IdNS_13fvsPatchFieldENS_11surfaceMeshEEERNS_14transportModelERKNS_4wordE+0x199) [0x7fd570058399]
[titan:22422] [ 9] simpleFoam() [0x4218e1]
[titan:22422] [10] /lib64/libc.so.6(__libc_start_main+0xf5) [0x7fd56d180455]
[titan:22422] [11] simpleFoam() [0x4241fd]
[titan:22422] *** End of error message ***
 
 in "/lib64/libc.so.6"
[1] #11  --------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 22422 on node titan exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

(sorry for this ugly presentation, but at least it is complete!)

Here is the log file:

PHP Code:


			
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 RAS turbulence model kEpsilon
bounding k, min: 0 max: 0.72 average: 0.72
bounding epsilon, min: 0 max: 0.002 average: 0.002

I tried to initialize k uniform .72 and epsilon uniform .002 (to avoid bas 0 values).

I really don't understand what's happening...

I am using swak4foam for boundary conditions in my inlet. U, k and epsilon are type groovyBC, with expression as fallow:
U >> "vector(2.09*log(pos().z/2+5),0,0)";
k >> "1.05-.0033*pos().z";
epsilon >> "(.66-.002*pos().z)/((pos().z/2+5)*log(pos().z/2+5))";

Corresponding are:
p >> type zeroGradient;
nut >> type calculated

I don't understand what's wrong... I don't know if I gave you enough details to help me; but I can give you other dictionaries, parameters...

Thanks in advance.

Djub · April 12, 2013, 12:13

Hi,
It seems that the problem comes from groovyBC:
1/ the case does not work neither in serial mode
2/ replacing the groovyBC by a fixedValue BC, it works!
More precisely, here are the two answer of simpleFoam:
With fixedValue (working):

PHP Code:


			
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    sigmaEps        1.3;
}
No field sources present
 
SIMPLE: convergence criteria
    field p      tolerance 0.005
    field U      tolerance 0.005
    field "(k|epsilon|omega)"    tolerance 0.005
 
Starting time loop
Time = 1
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.0006714641, No Iterations 2
 
(...)

With groovyBC (not working):

PHP Code:


			
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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 RAS turbulence model kEpsilon
bounding k, min: 0 max: 0.72 average: 0.72
bounding epsilon, min: 0 max: 0.002 average: 0.002
[12] #0  Foam::error::printStack(Foam::Ostream&)[4] #0[1] #0  Foam::error::printStack(Foam::Ostream&)  Foam::error::printStack(Foam::Ostream&)[0] #0  Foam::error::printStack(Foam::Ostream&)[8] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[12] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[4] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
 
(...)

I tried both with
libs ( "libgroovyBC.so" ) ;
and
libs ( "libOpenFOAM.so" "libgroovyBC.so" ) ;
in my controldict, but with same lack of success.

I should precise:
swakVersion: 0.2.1 (Release date: 2012-10-18)

It tried both with pts().z and pos().z, but again with same lack of success.

So I think something's wrong with groovyBC. Any help would be very useful, thanks!

Djub · April 12, 2013, 12:51

Hi! I found one mistake: my calculated k value was incorrect (negative for high z values).
I changed :
k >> "pos().z>250 ? .225 : 1.05-.0033*pos().z";
epsilon >> "pos().z<250 ? (.66-.002*pos().z)/((pos().z/2+5)*log(pos().z/2+5)) : .141/((pos().z/2+5)*log(pos().z/2+5))";

But it still does not work. By the way, with:
k >> "pos().z>250 ? .225 : 1.05-.0033*pos().z";
epsilon >> type fixedValue value uniform 0.002)
it works!

So finally, I think there is no problem with groovyBC (my excuses for the author), but with the expression. Does anybody knows why it does not work?

Thanks in advance

April 12, 2013, 12:13	Problem linked with groovyBC ?	#2
Djub Senior Member Julien Join Date: Jun 2012 Location: France Posts: 152 Rep Power: 14	Hi, It seems that the problem comes from groovyBC: 1/ the case does not work neither in serial mode 2/ replacing the groovyBC by a fixedValue BC, it works! More precisely, here are the two answer of simpleFoam: With fixedValue (working): PHP Code: // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 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 RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; sigmaEps 1.3; } No field sources present SIMPLE: convergence criteria field p tolerance 0.005 field U tolerance 0.005 field "(k\|epsilon\|omega)" tolerance 0.005 Starting time loop Time = 1 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.0006714641, No Iterations 2 (...) With groovyBC (not working): PHP Code: // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 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 RAS turbulence model kEpsilon bounding k, min: 0 max: 0.72 average: 0.72 bounding epsilon, min: 0 max: 0.002 average: 0.002 [12] #0 Foam::error::printStack(Foam::Ostream&)[4] #0[1] #0 Foam::error::printStack(Foam::Ostream&) Foam::error::printStack(Foam::Ostream&)[0] #0 Foam::error::printStack(Foam::Ostream&)[8] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [12] #1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [4] #1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" (...) I tried both with libs ( "libgroovyBC.so" ) ; and libs ( "libOpenFOAM.so" "libgroovyBC.so" ) ; in my controldict, but with same lack of success. I should precise: swakVersion: 0.2.1 (Release date: 2012-10-18) It tried both with pts().z and pos().z, but again with same lack of success. So I think something's wrong with groovyBC. Any help would be very useful, thanks! Last edited by Djub; April 12, 2013 at 12:21. Reason: some precisions added

April 12, 2013, 12:51	error not groovyBC, but in the expression	#3
Djub Senior Member Julien Join Date: Jun 2012 Location: France Posts: 152 Rep Power: 14	Hi! I found one mistake: my calculated k value was incorrect (negative for high z values). I changed : k >> "pos().z>250 ? .225 : 1.05-.0033pos().z"; epsilon* >> "pos().z<250 ? (.66-.002pos().z)/((pos().z/2+5)log(pos().z/2+5)) : .141/((pos().z/2+5)log(pos().z/2+5))"; But it still does not work. By the way, with: k >> "pos().z>250 ? .225 : 1.05-.0033pos().z"; epsilon >> type fixedValue value uniform 0.002) it works! So finally, I think there is no problem with groovyBC (my excuses for the author), but with the expression. Does anybody knows why it does not work? Thanks in advance