|
[Sponsors] |
August 1, 2008, 08:23 |
Description:
When planar STL
|
#1 |
Senior Member
Hannes Kröger
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 124
Rep Power: 18 |
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
__________________
silentdynamics GmbH - http://silentdynamics.de open source CAE software solutions & support |
|
August 1, 2008, 08:47 |
Thanks for report. The indexed
|
#2 |
Senior Member
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,419
Rep Power: 26 |
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)
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
planar diffuser simulation in Fluent | droper | FLUENT | 0 | November 1, 2008 15:29 |
non-planar faces due to the deformation of grid | gonski | Main CFD Forum | 5 | November 12, 2007 15:25 |
gambit non-planar trim | ahmet | FLUENT | 2 | October 16, 2005 17:47 |
planar mixing layers | Jitendra | Main CFD Forum | 5 | March 19, 2002 07:59 |
planar conduction | Ashutosh Joshi | FLUENT | 0 | December 21, 2000 09:09 |