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

Number of cells in mesh does not equal size of cell Level

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

Like Tree2Likes
  • 2 Post By ssss

Reply
 
LinkBack Thread Tools Display Modes
Old   January 14, 2015, 09:06
Default Number of cells in mesh does not equal size of cell Level
  #1
Member
 
Join Date: Dec 2012
Posts: 75
Rep Power: 6
pizzaspinate is on a distinguished road
hi,
I have the same error message when using sHM. I already tried the tips of an other thread posted here but it is still not running :-(
I have to say that Im new to OpenFOAM so it might be a really stupid mistake.
Thank you very much in advance for your help



-> FOAM FATAL ERROR:
Number of cells in mesh:512000 does not equal size of cellLevel:5829368
This might be because of a restart with inconsistent cellLevel.

From function hexRef8::getLevel0EdgeLength() const
in file polyTopoChange/polyTopoChange/hexRef8.C at line 358.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::hexRef8::getLevel0EdgeLength() const at ??:?
#3 Foam::hexRef8::hexRef8(FoamlyMesh const&, bool) at ??:?
#4 Foam::meshRefinement::meshRefinement(Foam::fvMesh& , double, bool, Foam::refinementSurfaces const&, Foam::refinementFeatures const&, Foam::shellSurfaces const&) at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
Aborted (core dumped)



I used the motorbike sHMDict file and changed it were needed (well at least the stl names).
It now looks like this:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object snappyHexMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

// Which of the steps to run
castellatedMesh true;
snap true;
addLayers false;



geometry
{
UCAV.stl
{
type triSurfaceMesh;
name UCAV;
}

// refinementBox
// {
// type searchableBox;
// min (-1.0 -0.7 0.0);
// max ( 8.0 0.7 2.5);
// }
};



// Settings for the castellatedMesh generation.
castellatedMeshControls
{

// Refinement parameters
// ~~~~~~~~~~~~~~~~~~~~~


maxLocalCells 100000;



maxGlobalCells 2000000;


minRefinementCells 10;


maxLoadUnbalance 0.10;



nCellsBetweenLevels 3;



// Explicit feature edge refinement
// ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

features
(
{
file "UCAV.eMesh";
level 6;
}
);



// Surface based refinement
// ~~~~~~~~~~~~~~~~~~~~~~~~



refinementSurfaces
{
motorBike
{

level (5 6);


patchInfo
{
type wall;
inGroups (motorBikeGroup);
}
}
}

// Resolve sharp angles
resolveFeatureAngle 30;


// Region-wise refinement
// ~~~~~~~~~~~~~~~~~~~~~~


refinementRegions
{
refinementBox
{
mode inside;
levels ((1E15 4));
}
}


// Mesh selection
// ~~~~~~~~~~~~~~


locationInMesh (0 0 0);



allowFreeStandingZoneFaces true;
}



// Settings for the snapping.
snapControls
{

nSmoothPatch 3;

/
tolerance 2.0;


nSolveIter 30;


nRelaxIter 5;

// Feature snapping


nFeatureSnapIter 10;


implicitFeatureSnap false;


explicitFeatureSnap true;


multiRegionFeatureSnap false;
}



// Settings for the layer addition.
addLayersControls
{

relativeSizes true;


layers
{
"(lowerWall|motorBike).*"
{
nSurfaceLayers 1;
}
}


expansionRatio 1.0;


finalLayerThickness 0.3;


minThickness 0.1;


nGrow 0;

// Advanced settings


featureAngle 60;


slipFeatureAngle 30;


nRelaxIter 3;


nSmoothSurfaceNormals 1;


nSmoothNormals 3;


nSmoothThickness 10;


maxFaceThicknessRatio 0.5;


maxThicknessToMedialRatio 0.3;


minMedianAxisAngle 90;



nBufferCellsNoExtrude 0;



nLayerIter 50;
}




meshQualityControls
{
#include "meshQualityDict"


// Advanced


nSmoothScale 4;

errorReduction 0.75;
}


// Advanced

/
writeFlags
(
layerFields
);


mergeTolerance 1e-6;
pizzaspinate is offline   Reply With Quote

Old   January 14, 2015, 14:27
Default
  #2
Senior Member
 
anonymous
Join Date: Aug 2014
Posts: 198
Rep Power: 4
ssss is on a distinguished road
Are you running 2 consecutive times snappyHexMesh?

Try removing all the files in constant/polyMesh/ with *level* or *Level* in its filename
sina_mech and BlnPhoenix like this.
ssss is offline   Reply With Quote

Old   September 28, 2016, 23:22
Default
  #3
Member
 
Anonymouse
Join Date: Dec 2015
Posts: 71
Rep Power: 3
KarenRei is on a distinguished road
Just wanted to add (for posterity) that this worked for me when I encountered this problem
KarenRei 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
No layers in a small gap bobburnquist OpenFOAM Native Meshers: snappyHexMesh and Others 6 August 26, 2015 09:38
foam-extend_3.1 decompose and pyfoam warning shipman OpenFOAM 3 July 24, 2014 08:14
killed "snappyHexMesh" parkh32 OpenFOAM Pre-Processing 2 April 8, 2012 17:12
snappyHexMesh aborting Tobi OpenFOAM Native Meshers: snappyHexMesh and Others 0 November 10, 2010 04:23
Problems in compiling paraview in Suse 10.3 platform chiven OpenFOAM Installation 3 December 1, 2009 08:21


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