|
[Sponsors] |
June 21, 2010, 02:26 |
cyclics and snappyHexMesh
|
#1 |
Member
Fabian Peng Karrholm
Join Date: Mar 2009
Posts: 61
Rep Power: 17 |
Hello,
I'm not sure if this is classified as a bug or a case of user error, but I've found out that when using snappyHexMesh on a block with cyclics, creation of layers from the stl crashes if you have more than one patch of cyclics. It does not crash if you collect all the cyclics into one patch with createPatch however. I've used the latest (as of last week) version of OpenFOAM 1.6.x. I only tried it with one stl surface, but the output from my debug compilation of OpenFOAM is: patch faces layers avg thickness[m] near-wall overall ----- ----- ------ --------- ------- stlfile_mypatch 78324 1 0.0208 0.0208 medialAxisSmoothingInfo : Calculate distance to Medial Axis ... shrinkMeshDistance : Smoothing normals ... Iteration 0 residual 0.00948 --> FOAM FATAL ERROR: index 2 out of range 0 ... 1 From function UList<T>::checkIndex(const label) in file ~/OpenFOAM/OpenFOAM-1.6.x/src/OpenFOAM/lnInclude/UListI.H at line 110. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) in "~/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPDebug/libOpenFOAM.so" #1 Foam::error::abort() at ~/OpenFOAM/OpenFOAM-1.6.x/src/OpenFOAM/lnInclude/error.C:250 #2 Foam::Ostream& Foam::operator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) at ~/OpenFOAM/OpenFOAM-1.6.x/src/OpenFOAM/lnInclude/errorManip.H:87 #3 Foam::UList<Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> >*>::checkIndex(int) const at ~/OpenFOAM/OpenFOAM-1.6.x/src/OpenFOAM/lnInclude/UListI.H:114 #4 Foam::UList<Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> >*>::operator[](int) at ~/OpenFOAM/OpenFOAM-1.6.x/src/OpenFOAM/lnInclude/UListI.H:141 #5 Foam::PtrList<Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> > >::set(int, Foam::PrimitivePatch<Foam::face, Foam::SubList, Foam::Field<Foam::Vector<double> > const&, Foam::Vector<double> >*) at ~/OpenFOAM/OpenFOAM-1.6.x/src/OpenFOAM/lnInclude/PtrListI.H:65 #6 Foam::PointEdgeWave<Foam::pointData>::calcCyclicAd dressing() at ~/OpenFOAM/OpenFOAM-1.6.x/src/meshTools/lnInclude/PointEdgeWave.C:80 #7 Foam::PointEdgeWave<Foam::pointData>::PointEdgeWav e(Foam::polyMesh const&, Foam::List<int> const&, Foam::List<Foam::pointData> const&, Foam::List<Foam::pointData>&, Foam::List<Foam::pointData>&, int) at ~/OpenFOAM/OpenFOAM-1.6.x/src/meshTools/lnInclude/PointEdgeWave.C:748 #8 Foam::autoLayerDriver::medialAxisSmoothingInfo(Foa m::motionSmoother const&, int, int, double, Foam::GeometricField<Foam::Vector<double>, Foam::pointPatchField, Foam::pointMesh>&, Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>&, Foam::GeometricField<double, Foam::pointPatchField, Foam::pointMesh>&) const at ~/OpenFOAM/OpenFOAM-1.6.x/src/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriverShrink.C:793 #9 Foam::autoLayerDriver::addLayers(Foam::layerParame ters const&, Foam::dictionary const&, Foam::List<int> const&, int, Foam::decompositionMethod&, Foam::fvMeshDistribute&) at ~/OpenFOAM/OpenFOAM-1.6.x/src/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C:2846 #10 Foam::autoLayerDriver::doLayers(Foam::dictionary const&, Foam::dictionary const&, Foam::layerParameters const&, bool, Foam::decompositionMethod&, Foam::fvMeshDistribute&) at ~/OpenFOAM/OpenFOAM-1.6.x/src/autoMesh/autoHexMesh/autoHexMeshDriver/autoLayerDriver.C:3349 #11 main at ~/OpenFOAM/OpenFOAM-1.6.x/applications/utilities/mesh/generation/snappyHexMesh/snappyHexMesh.C:470 #12 __libc_start_main in "/lib64/libc.so.6" #13 _start in "~/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPDebug/snappyHexMesh" /Fabian |
|
June 22, 2010, 12:51 |
|
#2 |
Senior Member
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,419
Rep Power: 26 |
Thanks, was a real bug - used bool instead of label. Fixed in 16x.
Thanks, Mattijs |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[snappyHexMesh] SnappyHexMeshDict | fabrizio | OpenFOAM Meshing & Mesh Conversion | 25 | October 14, 2011 07:00 |