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

cyclics and snappyHexMesh

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

LinkBack Thread Tools Display Modes
Old   June 21, 2010, 02:26
Default cyclics and snappyHexMesh
Fabian Peng Karrholm
Join Date: Mar 2009
Posts: 61
Rep Power: 9
fabianpk is on a distinguished road
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

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/"
#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/"
#13 _start in "~/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPDebug/snappyHexMesh"

fabianpk is offline   Reply With Quote

Old   June 22, 2010, 12:51
Super Moderator
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 17
mattijs is on a distinguished road
Thanks, was a real bug - used bool instead of label. Fixed in 16x.


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
SnappyHexMeshDict fabrizio OpenFOAM Mesh Utilities 25 October 14, 2011 07:00

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