[grandgo]

hi foamers,

i don't know if this is a bug or not, but weird anyway.

i want to sample some points on the radius of a pipe-flow via sampleDict. the diameter of the pipe is 0.05.
by using this code in the sampleDict file:

setFormat raw;
surfaceFormat null;
interpolationScheme cell;
fields
(
p
U
);
sets
(
lineRadiusY4DmidPoint
{
        type        midPoint;
        axis        x;
        start       (0 0.2 0); 
        end         (0.025 0.2 0);
}
);
surfaces();

i get this error-message

--> FOAM FATAL ERROR: 
None of the neighbours of face 289377 contains point (-1.68086e-18 0.199999 6.11069e-07)

    From function sampledSet::getCell(const label, const point&)
    in file sampledSet/sampledSet/sampledSet.C at line 104.

FOAM aborting
#0  Foam::error::printStack(Foam::Ostream&) in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::error::abort() in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2  Foam::sampledSet::getCell(int, Foam::Vector<double> const&) const in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libsampling.so"
#3  Foam::midPointSet::genSamples() in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libsampling.so"
#4  Foam::midPointSet::midPointSet(Foam::word const&, Foam::polyMesh const&, Foam::meshSearch&, Foam::dictionary const&) in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libsampling.so"
#5  Foam::sampledSet::addwordConstructorToTable<Foam::midPointSet>::New(Foam::word const&, Foam::polyMesh const&, Foam::meshSearch&, Foam::dictionary const&) in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libsampling.so"
#6  Foam::sampledSet::New(Foam::word const&, Foam::polyMesh const&, Foam::meshSearch&, Foam::dictionary const&) in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libsampling.so"
#7  void Foam::PtrList<Foam::sampledSet>::read<Foam::sampledSet::iNew>(Foam::Istream&, Foam::sampledSet::iNew const&) in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libsampling.so"
#8  Foam::sampledSets::read(Foam::dictionary const&) in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libsampling.so"
#9  Foam::sampledSets::sampledSets(Foam::word const&, Foam::objectRegistry const&, Foam::dictionary const&, bool) in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libsampling.so"
#10  Foam::OutputFilterFunctionObject<Foam::sampledSets>::allocateFilter() in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libsampling.so"
#11  Foam::OutputFilterFunctionObject<Foam::sampledSets>::start() in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libsampling.so"
#12  Foam::functionObjectList::read() in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#13  Foam::Time::operator++() in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#14  Foam::Time::loop() in "/sw/OpenFOAM/OpenFOAM-1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#15  
 in "/home/stss8/OpenFOAM/stss8-1.7.x/applications/bin/linux64GccDPOpt/pisoFoamTNuLESEff"
#16  __libc_start_main in "/lib64/libc.so.6"
#17  
 at /usr/src/packages/BUILD/glibc-2.11.2/csu/../sysdeps/x86_64/elf/start.S:116

But, if i change the sampleDict file a little bit to this code:

setFormat raw;
surfaceFormat null;
interpolationScheme cell;
fields
(
p
U
);
sets
(
lineRadiusY4DmidPoint
{
        type        midPoint;
        axis        x;
        start       (0.001 0.2 0); 
        end         (0.025 0.2 0);
}
);
surfaces();

everything works fine.

the problem is: with a finer mesh it doesn't work at all, no matter what what the starting point is and i get an error message similar to that one above. in this case i have to use 'uniform' instead of 'midPoint'.

i'm wondering if this is a bug or just incorrect usage.

same problem here:
http://www.cfd-online.com/Forums/ope...ionscheme.html

best regards,
grandgo

[nlc]

Did you find a solution ?? A have this problem with a fine mesh and it dos not work what ever I do !!

Regards
Nicolas L.

[musahossein]

I think the founders of OpenFOAM have not looked at the bug part of the forum for a while and may not be doing so in the future. There are a lot of unanswered questions. So dont hold your breath for a solution. Just a suggestion. I hope I am wrong.

[nlc]

I understand.

I live with the error for now as It dos not bug me enough for now to give me the will of finding the error in the C++ code. In any case I'm using now my own tool to sample and post process the data for now.

Thank you for you r interest musahossein.

[romant]

Hej,

for bug reports you should use http://www.openfoam.com/mantisbt/view_all_bug_page.php which is the official bug tracker for openfoam.

[musahossein]

Originally Posted by romant

Hej,

for bug reports you should use http://www.openfoam.com/mantisbt/view_all_bug_page.php which is the official bug tracker for openfoam.

Is this also located in the users forum?

[romant]

no. in order to file a bug report you have to create an account on the bugtracker website. this way the OF developers get it faster and can assign the right ppl to it.