CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Bugs

problems with midPoint set sampling

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By grandgo

Reply
 
LinkBack Thread Tools Display Modes
Old   July 4, 2011, 05:45
Default problems with midPoint set sampling
  #1
Member
 
Join Date: Oct 2010
Location: Stuttgart
Posts: 35
Rep Power: 6
grandgo is on a distinguished road
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:
Sample bugg with midPoint cell interpolationScheme

best regards,
grandgo
bmikuz likes this.
grandgo is offline   Reply With Quote

Old   October 3, 2011, 18:49
Default
  #2
Member
 
ubald's Avatar
 
Nicolas Lussier Clément
Join Date: Apr 2009
Location: Montréal, Qc, Canada
Posts: 47
Rep Power: 8
ubald is on a distinguished road
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.
ubald is offline   Reply With Quote

Old   December 12, 2011, 14:16
Default
  #3
Senior Member
 
musaddeque hossein
Join Date: Mar 2009
Posts: 307
Rep Power: 9
musahossein is on a distinguished road
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.
musahossein is offline   Reply With Quote

Old   December 12, 2011, 14:24
Default
  #4
Member
 
ubald's Avatar
 
Nicolas Lussier Clément
Join Date: Apr 2009
Location: Montréal, Qc, Canada
Posts: 47
Rep Power: 8
ubald is on a distinguished road
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.
ubald is offline   Reply With Quote

Old   December 13, 2011, 04:30
Default reporting bugs
  #5
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 359
Rep Power: 11
romant is on a distinguished road
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.
__________________
~roman
romant is offline   Reply With Quote

Old   December 13, 2011, 09:58
Default
  #6
Senior Member
 
musaddeque hossein
Join Date: Mar 2009
Posts: 307
Rep Power: 9
musahossein is on a distinguished road
Quote:
Originally Posted by romant View Post
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?
musahossein is offline   Reply With Quote

Old   December 13, 2011, 10:34
Default
  #7
Senior Member
 
romant's Avatar
 
Roman Thiele
Join Date: Aug 2009
Location: Stockholm, Sweden
Posts: 359
Rep Power: 11
romant is on a distinguished road
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.
__________________
~roman
romant is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Questions for a species transport problems (snapshots attached) aleisia FLUENT 2 October 9, 2011 04:40
OpenFOAM 1.7.1 installation problem on Fedora 14 armonica OpenFOAM Installation 16 March 31, 2011 13:16
Compiling OpenFOAM13 on AMD64 with Redhat Enterprise mbeaudoin OpenFOAM Installation 20 June 17, 2008 06:43
Problems about multi swirl tubes simulation using Jason FLUENT 4 May 9, 2008 00:33
Problem in installation of OpenFOAM sachin OpenFOAM Installation 7 January 22, 2008 02:40


All times are GMT -4. The time now is 01:30.