CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Bugs

SnappyHexMesh and planar STL surfaces

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

LinkBack Thread Tools Display Modes
Old   August 1, 2008, 08:23
Default Description: When planar STL
Hannes Kröger
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 97
Rep Power: 11
hannes is on a distinguished road
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
#1 Foam::error::abort() in
Foam::indexedOctree<foam::treedatatrisurface>::div ide(Foam::treeBoundBox
const&, Foam::DynamicList<foam::list<int>, 0u, 2u, 1u>&, int) const in
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/"
#4 Foam::triSurfaceMesh::tree() const in
#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
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/"
#7 Foam::meshRefinement::updateIntersections(Foam::Li st<int> const&) in
#8 Foam::meshRefinement::meshRefinement(Foam::fvMesh& , double,
Foam::refinementSurfaces const&, Foam::shellSurfaces const&) in
#9 main in
"/home/cluster/josef/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/snap pyHexMesh"
#10 __libc_start_main in "/lib64/"
#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:

bb.min()[0] >= bb.max()[0]
|| bb.min()[1] >= bb.max()[1]
|| bb.min()[2] >= bb.max()[2]
<< "Badly formed bounding box:" << bb
<< abort(FatalError);

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


Source file:

Linux x86_64, Fedora 9

silentdynamics GmbH -
open source CAE software solutions & support
hannes is offline   Reply With Quote

Old   August 1, 2008, 08:47
Default Thanks for report. The indexed
Super Moderator
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,419
Rep Power: 18
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


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 16: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 17: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 16:09.