CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Error using LaunderGibsonRSTM on SGI ALTIX 4700 (http://www.cfd-online.com/Forums/openfoam/60821-error-using-laundergibsonrstm-sgi-altix-4700-a.html)

jaswi March 19, 2008 10:46

Dear Forum Wish you all a v
 
Dear Forum

Wish you all a very nice day.

I have come across a strange error while using
LaunderGibsonRSTM on SGI ALTIX 4700. I am using dev version and it has been compiled with Intel 10.0. Below provided is the MPI traceback.

I request the experts to take a look at it. I do not understand most of the cryptic stuff the traceback has thrown out.

The case (with LaunderGibsonRSTM ) does not runs even with a single processor.It gives out a floating point exception error. Changing the turbulence model to kEpsilon or RNGkEpsilon leads to a successful run.

The same case work normally on a GNU compiled development version without any error.

MPI: --------stack traceback-------
Internal Error: Can't read/write file "/dev/mmtimer", (errno = 22)
MPI: Intel(R) Debugger for applications running on IA-64, Version 9.1-29 , Build 20070405
MPI: Reading symbolic information from /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/applications/bin/linuxIA64 I64DPOpt/InterRotati
onFoam...done
MPI: Attached to process id 1914 ....
MPI: stopped at [0xa000000000010621]
MPI: >0 0xa000000000010621
MPI: #1 0x2000000002c1b6f0 in __waitpid(...) in /lib/libc-2.4.so
MPI: #2 0x2000000003499340 in MPI_SGI_stacktraceback(...) in /usr/lib/libmpi.so
MPI: #3 0x200000000349ab60 in slave_sig_handler(...) in /usr/lib/libmpi.so
MPI: #4 0xa0000000000107e0
MPI: #5 0xa000000000010620
MPI: #6 0x2000000002b4c1c0 in raise(...) in /lib/libc-2.4.so
MPI: #7 0x2000000001e8f1a0 in _ZN4Foam6sigFpe13sigFpeHandlerEi(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linuxIA64
I64DPOpt/libOpenFOAM.so
MPI: #8 0xa0000000000107e0
MPI: #9 0x200000000049e9c0 in _ZN4Foam12FaceCellWaveINS_13wallPointDataINS_6Vect orIdEEEEEC9ERKNS_8polyMeshERKN S_4ListIiEERKNS9_IS4_EERNS_5
UListIS4_EESI_i(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linuxIA64I64DPOpt/libi ncompressibleTurbulenceModels.so
MPI: #10 0x200000000049d4c0 in _ZN4Foam8MeshWaveINS_13wallPointDataINS_6VectorIdE EEEEC9ERKNS_8polyMeshERKNS_4Li stIiEERKNS9_IS4_EEi(...) in
/home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linuxIA64I64DPOpt/libi ncompressibleTurbulenceModels.so
MPI: #11 0x200000000049b380 in _ZN4Foam13patchDataWaveINS_13wallPointDataINS_6Vec torIdEEEEE7correctEv(...) in /home/hlrb2/h0983/h0983af/Ope
nFOAM/OpenFOAM-1.4.1-dev/lib/linuxIA64I64DPOpt/libincompressibleTurbulenceModels .so
MPI: #12 0x2000000000499a70 in _ZN4Foam13patchDataWaveINS_13wallPointDataINS_6Vec torIdEEEEEC9ERKNS_8polyMeshERK NS_7HashSetIiNS_4HashIiEEEER
KNS_4ListIPNS_5FieldIS3_EEEEb(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linuxIA64I64DPOpt/libi ncompressibleTurbulen
ceModels.so
MPI: #13 0x2000000000498320 in _ZN4Foam12wallDistDataINS_13wallPointDataINS_6Vect orIdEEEEE7correctEv(...) in /home/hlrb2/h0983/h0983af/Open
FOAM/OpenFOAM-1.4.1-dev/lib/linuxIA64I64DPOpt/libincompressibleTurbulenceModels. so
MPI: #14 0x2000000000a78550 in _ZN4Foam12wallDistDataINS_13wallPointDataINS_6Vect orIdEEEEEC9ERKNS_6fvMeshERNS_1 4GeometricFieldIS3_NS_12fvPa
tchFieldENS_7volMeshEEEb(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linuxIA64I64DPOpt/libf initeVolume.so
MPI: #15 0x2000000000a66590 in _ZN4Foam18wallDistReflectionC9ERKNS_6fvMeshEb(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/
lib/linuxIA64I64DPOpt/libfiniteVolume.so
MPI: #16 0x2000000000a664f0 in _ZN4Foam18wallDistReflectionC1ERKNS_6fvMeshEb(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/
lib/linuxIA64I64DPOpt/libfiniteVolume.so
MPI: #17 0x200000000047ee80 in _ZN4Foam16turbulenceModels17LaunderGibsonRSTMC9ERK NS_14GeometricFieldINS_6Vector IdEENS_12fvPatchFieldENS_7vo
lMeshEEERKNS2_IdNS_13fvsPatchFieldENS_11surfaceMes hEEERNS_14transportModelE(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/l
ib/linuxIA64I64DPOpt/libincompressibleTurbulenceModels.so
MPI: #18 0x200000000047e7b0 in _ZN4Foam16turbulenceModels17LaunderGibsonRSTMC1ERK NS_14GeometricFieldINS_6Vector IdEENS_12fvPatchFieldENS_7vo
lMeshEEERKNS2_IdNS_13fvsPatchFieldENS_11surfaceMes hEEERNS_14transportModelE(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/l
ib/linuxIA64I64DPOpt/libincompressibleTurbulenceModels.so
MPI: #19 0x20000000004ab9b0 in _ZN4Foam15turbulenceModel31adddictionaryConstructo rToTableINS_16turbulenceModels 17LaunderGibsonRSTMEE3NewERK
NS_14GeometricFieldINS_6VectorIdEENS_12fvPatchFiel dENS_7volMeshEEERKNS5_IdNS_13f vsPatchFieldENS_11surfaceMeshEEERNS_14transportMod elE(...)
in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linuxIA64I64DPOpt/libi ncompressibleTurbulenceModels.so
MPI: #20 0x2000000000274f30 in _ZN4Foam15turbulenceModel3NewERKNS_14GeometricFiel dINS_6VectorIdEENS_12fvPatchFi eldENS_7volMeshEEERKNS1_IdNS_13fvsPatchFieldENS_11 surfaceMeshEEERNS_14transportM odelE(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/lib/linuxIA64I64DPOpt/libi ncompressibleTurbulenceModels.so
MPI: #21 0x40000000000299f0 in main(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/applications/bin/linuxIA64 I64DPOpt/InterRotationFoam
MPI: #22 0x2000000002b1fc20 in __libc_start_main(...) in /lib/libc-2.4.so
MPI: #23 0x40000000000274c0 in _start(...) in /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/applications/bin/linuxIA64 I64DPOpt/InterRotationFoam

MPI: -----stack traceback ends-----
MPI: On host a07, Program /home/hlrb2/h0983/h0983af/OpenFOAM/OpenFOAM-1.4.1-dev/applications/bin/linuxIA64 I64DPOpt/InterRotationFoam, Rank 14, Process 1914: Dumping core on signal SIGFPE(8) into directory /ptmp2/h0983/h0983af/OpenFOAM/run
MPI: MPI_COMM_WORLD rank 14 has terminated without calling MPI_Finalize()
MPI: aborting job
MPI: Received signal 8


Thanks in advance

With Best Regards
Jaswi

jaswi March 19, 2008 14:35

Dear Forum Good Evening.
 
Dear Forum

Good Evening.

Digging in for the source of error this is what I found:


LaunderGibsonRSTM has alongwith other objects:

(in LaunderGibsonRSTM.H file) wallDistReflection yr_;
(in LaunderGibsonRSTM.C file) yr_(mesh_) is initialized.

Next looking into wallDistReflection.H and .C, its
constructor is defined as:
-----------
Foam::wallDistReflection::wallDistReflection
(
const Foam::fvMesh& mesh,
const bool correctWalls
)
:
reflectionVectors(mesh),
wallDistData<wallpointdata<vector> >
(
mesh,
reflectionVectors::n_,
correctWalls
)
{}
------
reflectionVectors(mesh) shows no errors.

The source of error lies in wallDistData. Its constructor is defined as :
----------------------
template<class>
Foam::wallDistData<transfertype>::wallDistData
(
const Foam::fvMesh& mesh,
GeometricField<type,>& field,
const bool correctWalls
)
:
volScalarField
(
IOobject
(
"y",
mesh.time().timeName(),
mesh
),
mesh,
dimensionedScalar("y", dimLength, GREAT)
),
cellDistFuncs(mesh),
field_(field),
correctWalls_(correctWalls),
nUnset_(0)
{
correct();
}
----------------------
Everything is initialized but when the function correct() is called the simulation crashses with floating point exception error.

This is what the case run shows (with debug switches on) just before it crashes:

Constructing IOobject called turbulenceProperties of type IOobject
objectRegistry::checkIn(regIOobject&) : checking in turbulenceProperties
IOobject::readHeader(Istream&) : reading header for file "/ptmp2/h0983/h0983af/OpenFOAM/run/1st_case_370K_RSTM/constant/turbulencePropert ies"
.... read
Destroying regIOobject called turbulenceProperties of type regIOobject in directory "/ptmp2/h0983/h0983af/OpenFOAM/run/1st_case_370K_RSTM/constant"
objectRegistry::checkOut(regIOobject&) : checking out turbulenceProperties
Selecting turbulence model LaunderGibsonRSTM
Constructing IOobject called turbulenceProperties of type IOobject
objectRegistry::checkIn(regIOobject&) : checking in turbulenceProperties
IOobject::readHeader(Istream&) : reading header for file "/ptmp2/h0983/h0983af/OpenFOAM/run/1st_case_370K_RSTM/constant/turbulencePropert ies"
.... read
Calculating ypl
GeometricField<type,>::GeometricBoundaryField::Geo metricBoundaryField(const BoundaryMesh&, const Field<type>&, const word&)
fvPatchField<type>::New(const word&, const fvPatch&, const DimensionedField<type,>&) : constructing fvPatchField<type>
fvPatchField<type>::New(const word&, const fvPatch&, const DimensionedField<type,>&) : constructing fvPatchField<type>
Inside cellDistFuncs.C
PrimitivePatch<face,>::calcAddressing() : calculating patch addressing
PrimitivePatch<face,>::calcMeshData() : calculating mesh data in PrimitivePatch
PrimitivePatch<face,>::calcMeshData() : finished calculating mesh data in PrimitivePatch
PrimitivePatch<face,>::calcPointFaces() : calculating pointFaces
PrimitivePatch<face,>::calcPointFaces() finished calculating pointFaces
PrimitivePatch<face,>::calcAddressing() : finished calculating patch addressing
PrimitivePatch<face,>::calcAddressing() : calculating patch addressing
PrimitivePatch<face,>::calcMeshData() : calculating mesh data in PrimitivePatch
PrimitivePatch<face,>::calcMeshData() : finished calculating mesh data in PrimitivePatch
PrimitivePatch<face,>::calcPointFaces() : calculating pointFaces
PrimitivePatch<face,>::calcPointFaces() finished calculating pointFaces
PrimitivePatch<face,>::calcAddressing() : finished calculating patch addressing
Inside TurbulenceModel
Constructing IOobject called reflectionVectors of type IOobject
objectRegistry::checkIn(regIOobject&) : checking in reflectionVectors
GeometricField<type,>::GeometricBoundaryField::Geo metricBoundaryField(const BoundaryMesh&, const Field<type>&, const word&)
fvPatchField<type>::New(const word&, const fvPatch&, const DimensionedField<type,>&) : constructing fvPatchField<type>
fvPatchField<type>::New(const word&, const fvPatch&, const DimensionedField<type,>&) : constructing fvPatchField<type>
GeometricField<type,>::GeometricField : creating temporary
IOobject: volVectorField reflectionVectors "/ptmp2/h0983/h0983af/OpenFOAM/run/1st_case_370K_RSTM/0"

Inside reflectionVectors.C
Inside reflectionVectors.C --> Correct()
Constructing IOobject called y of type IOobject
objectRegistry::checkIn(regIOobject&) : checking in y
GeometricField<type,>::GeometricBoundaryField::Geo metricBoundaryField(const BoundaryMesh&, const Field<type>&, const word&)
fvPatchField<type>::New(const word&, const fvPatch&, const DimensionedField<type,>&) : constructing fvPatchField<type>
fvPatchField<type>::New(const word&, const fvPatch&, const DimensionedField<type,>&) : constructing fvPatchField<type>
GeometricField<type,>::GeometricField : creating temporary
IOobject: volScalarField y "/ptmp2/h0983/h0983af/OpenFOAM/run/1st_case_370K_RSTM/0"

Inside cellDistFuncs.C
Inside wallDistData.C
Inside cellDistFuncs.C
primitiveMesh::calcCells() : calculating cells
Floating point exception

Anybody who knows what is wrong here, Please help !.

Thanks in advance
Jaswi

lakeat April 29, 2008 10:54

Let's keep in touch, for I use
 
Let's keep in touch, for I use IRIX too, and I am suffering the same error as you.

Regards
Daniel


All times are GMT -4. The time now is 13:21.