nlc |
November 12, 2010 16:11 |
Sample bugg with midPoint cell interpolationScheme
I don't know if its a known bugg I just thought I should report it !!
using
interpolationScheme cell with midPoint for set will bugg with OpenFOAM 1.7.1 in a certain case !
Watt append hear is that the line I chose to sample on append to fit exactly with the faces of my structured mesh.
My 2D test case mesh is 1000m long and 100m high with 300 x 50 cells so at 100m I'm exactly on the face sheared by 30 and 31 column of cell !
the SampleDict is:
Code:
*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object sampleDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
setFormat raw;
surfaceFormat raw;
interpolationScheme cell;
fields
(
U
p
epsilon
k
);
sets
(
x0
{
type midPoint;
axis z;
start ( 0 0 100);
end ( 0 0 0 );
}
x100
{
type midPoint;
axis z;
start ( 100 0 100);
end ( 100 0 0 );
}
);
The error is :
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.1 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.7.1-03e7e056c215
Exec : sample -time 20000
Date : Nov 12 2010
Time : 15:40:56
Host : ubi
PID : 26371
Case : /home/niko/OpenFOAM/niko-1.7.1/run/AEMN/case/ABL/ABLzMesh
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 20000
--> FOAM FATAL ERROR:
Found cell 14730 using face 60080. But cell does not contain point (100 1.27012e-14 98.1575)
From function sampledSet::getCell(const label, const point&)
in file sampledSet/sampledSet/sampledSet.C at line 75.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccSPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccSPOpt/libOpenFOAM.so"
#2 Foam::sampledSet::getCell(int, Foam::Vector<float> const&) const in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccSPOpt/libsampling.so"
#3 Foam::midPointSet::genSamples() in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccSPOpt/libsampling.so"
#4 Foam::midPointSet::midPointSet(Foam::word const&, Foam::polyMesh const&, Foam::meshSearch&, Foam::dictionary const&) in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccSPOpt/libsampling.so"
#5 Foam::sampledSet::addwordConstructorToTable<Foam::midPointSet>::New(Foam::word const&, Foam::polyMesh const&, Foam::meshSearch&, Foam::dictionary const&) in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccSPOpt/libsampling.so"
#6 Foam::sampledSet::New(Foam::word const&, Foam::polyMesh const&, Foam::meshSearch&, Foam::dictionary const&) in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccSPOpt/libsampling.so"
#7 void Foam::PtrList<Foam::sampledSet>::read<Foam::sampledSet::iNew>(Foam::Istream&, Foam::sampledSet::iNew const&) in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccSPOpt/libsampling.so"
#8 Foam::sampledSets::read(Foam::dictionary const&) in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccSPOpt/libsampling.so"
#9 Foam::sampledSets::sampledSets(Foam::word const&, Foam::objectRegistry const&, Foam::dictionary const&, bool) in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccSPOpt/libsampling.so"
#10
at sample.C:0
#11
in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccSPOpt/sample"
#12 __libc_start_main in "/lib/libc.so.6"
#13
in "/home/niko/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccSPOpt/sample"
Abandon
If I run the same sampleDict with:
Code:
x100
{
type midPoint;
axis z;
start ( 101 0 100);
end ( 101 0 0 );
}
no error and every thing is great !
Hope it help to get this wonderfully well programmed C++ library, OpenFOAM, even better.
As now I have no clue on how this bugg could be fixed. :)
Happy openfoaming to every one !
Nicolas Lussier
|