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

SnappyHexMesh and planar STL surfaces

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 1, 2008, 09:23
Default Description: When planar STL
  #1
Member
 
Hannes Kröger
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 88
Rep Power: 8
hannes is on a distinguished road
Description:
When planar STL surfaces which are aligned with coordinate planes are supplied to snappyHexMesh, an error message like to following occurs:

==============================================
Edge intersection testing:
Number of edges : 74800
Number of edges to retest : 74800


Badly formed bounding box-0.0720968 0.027 -0.0410269) (0.080055 0.027
0.0410628)#0 Foam::error::printStack(Foam:stream&) in
"/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in
"/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#2
Foam::indexedOctree<foam::treedatatrisurface>::div ide(Foam::treeBoundBox
const&, Foam::DynamicList<foam::list<int>, 0u, 2u, 1u>&, int) const in
"/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libmeshTools.so"
#3
Foam::indexedOctree<foam::treedatatrisurface>::ind exedOctree(Foam::treeDataTriSu rface const&, Foam::treeBoundBox const&, int, double, double) in "/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libmeshTools.so"
#4 Foam::triSurfaceMesh::tree() const in
"/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libmeshTools.so"
#5 Foam::triSurfaceMesh::findLineAny(Foam::Field<foam ::vector<double> >
const&, Foam::Field<foam::vector<double> > const&,
Foam::List<foam::pointindexhit<foam::vector<double > > >&) const in
"/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libmeshTools.so"
#6
Foam::refinementSurfaces::findHigherIntersection(F oam::Field<foam::vector<double > > const&, Foam::Field<foam::vector<double> > const&, Foam::List<int> const&, Foam::List<int>&, Foam::List<int>&) const in "/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libautoMesh.so"
#7 Foam::meshRefinement::updateIntersections(Foam::Li st<int> const&) in
"/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libautoMesh.so"
#8 Foam::meshRefinement::meshRefinement(Foam::fvMesh& , double,
Foam::refinementSurfaces const&, Foam::shellSurfaces const&) in
"/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libautoMesh.so"
#9 main in
"/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/snap pyHexMesh"
#10 __libc_start_main in "/lib64/libc.so.6"
#11 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat,
Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in
"/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/snap pyHexMesh"


From function indexedOctree<type>::divide
in file lnInclude/indexedOctree.C at line 213.

FOAM aborting

==============================================

The statement after line 205 in $FOAM_SRC/meshTools/lnInclude/indexedOctree.C causes the error:

if
(
bb.min()[0] >= bb.max()[0]
|| bb.min()[1] >= bb.max()[1]
|| bb.min()[2] >= bb.max()[2]
)
{
FatalErrorIn("indexedOctree<type>::divide")
<< "Badly formed bounding box:" << bb
<< abort(FatalError);
}

If this check is removed and libmeshTools recompiled, the code runs fine without errors.


Solver/Application:
snappyHexMesh

Source file:
$FOAM_SRC/meshTools/lnInclude/indexedOctree.C

Platform:
Linux x86_64, Fedora 9

Version:
1.5
hannes is offline   Reply With Quote

Old   August 1, 2008, 09:47
Default Thanks for report. The indexed
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 15
mattijs is on a distinguished road
Thanks for report. The indexedOctree is for 3D geometry - am surprised it actually produces something useful for 2D geometry. You might get problems with truncation error. Is it possible to get the surface, blockMeshDict and snappyHexMeshDict? (m dot janssens)
mattijs 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
planar diffuser simulation in Fluent droper FLUENT 0 November 1, 2008 17:29
non-planar faces due to the deformation of grid gonski Main CFD Forum 5 November 12, 2007 16:25
gambit non-planar trim ahmet FLUENT 2 October 16, 2005 18:47
planar mixing layers Jitendra Main CFD Forum 5 March 19, 2002 08:59
planar conduction Ashutosh Joshi FLUENT 0 December 21, 2000 10:09


All times are GMT -4. The time now is 20:12.