CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (http://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   problems with midPoint set sampling (http://www.cfd-online.com/Forums/openfoam-bugs/90196-problems-midpoint-set-sampling.html)

grandgo July 4, 2011 05:45

problems with midPoint set sampling
 
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:

Code:

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

Code:

--> 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:

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

ubald October 3, 2011 18:49

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 December 12, 2011 14:16

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.

ubald December 12, 2011 14:24

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 December 13, 2011 04:30

reporting bugs
 
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 December 13, 2011 09:58

Quote:

Originally Posted by romant (Post 335634)
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 December 13, 2011 10:34

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.


All times are GMT -4. The time now is 22:43.