CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[snappyHexMesh] snappy hangs while layer generation

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 25, 2009, 06:05
Default snappy hangs while layer generation
  #1
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 530
Rep Power: 20
bastil is on a distinguished road
Hi group,

I have snappy hanging while layer generation quite regularely for our geometries at the following stage:

Setting up information for layer truncation

snappy is still using all the CPUs and memory but does not seem to work at all. How can I find out more about the problem? Thanks.

BastiL
bastil is offline   Reply With Quote

Old   November 27, 2009, 06:38
Default
  #2
Member
 
bruce
Join Date: May 2009
Location: Germany
Posts: 42
Rep Power: 17
bruce is on a distinguished road
Hi
what i understood when i had the same issue that the Layer parameters are inconsistant to get fast extrusion.

e.g too high layer size -> takes long time to push mesh
look at how many faces were stopped during extrusion? that will give some info.

Post your Layer.log file, small test case as well.

--
bruce is offline   Reply With Quote

Old   October 29, 2010, 06:00
Default process killed while layer generation
  #3
New Member
 
CFD user
Join Date: Apr 2010
Location: Germany
Posts: 28
Rep Power: 16
subhkirti is on a distinguished road
Hi bruce,

I am having the same problem as bastil. The entire process of mesh generation by sHM gets killed when it starts writing the layer information in the 3rd time directory.

Specifically,

Layer addition iteration 0
--------------------------
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.

Determining displacement for added points according to pointNormal ...
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.

shrinkMeshMedialDistance : Smoothing using Medial Axis ...
shrinkMeshMedialDistance : Reduce layer thickness at 1242 nodes where thickness to medial axis distance is large
shrinkMeshDistance : Removing isolated regions ...
Number isolated points extrusion stopped : 1
shrinkMeshDistance : Smoothing field ...
Iteration 0 residual 0.301
shrinkMeshMedialDistance : Moving mesh ...
Iteration 0
Moving mesh using diplacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 65 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.02 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.02 : 0
faces on cells with determinant < 0.001 : 0
shrinkMeshMedialDistance : Successfully moved mesh
shrinkMeshMedialDistance : Finished moving mesh ...
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.
truncateDisplacement : Unextruded 0 faces due to non-consecutive vertices being extruded.

Setting up information for layer truncation ...
Killed

I am not able to figure out the problem. Can you help me out with this?

The layer part of my sHMDict file is given below:

/ Settings for the layer addition.
addLayersControls
{

relativeSizes true;

// Per final patch (so not geometry!) the layer information
layers
{
gondel
{
nSurfaceLayers 5;
}
}

// Expansion factor for layer mesh
expansionRatio 1.15;

finalLayerThickness 0.2;

//- Wanted thickness of final added cell layer. If multiple layers
// is the
// thickness of the layer furthest away from the wall.
// Relative to undistorted size of cell outside layer.
finalLayerRatio 0.1;

//- Minimum thickness of cell layer. If for any reason layer
// cannot be above minThickness do not add layer.
// Relative to undistorted size of cell outside layer.
minThickness 0.01;

//- If points get not extruded do nGrow layers of connected faces that are
// also not grown. This helps convergence of the layer addition process
// close to features.
nGrow 1;


// Advanced settings

//- When not to extrude surface. 0 is flat surface, 90 is when two faces
// make straight angle.
featureAngle 90;

//- Maximum number of snapping relaxation iterations. Should stop
// before upon reaching a correct mesh.
nRelaxIter 10;

// Number of smoothing iterations of surface normals
nSmoothSurfaceNormals 1;

// Number of smoothing iterations of interior mesh movement direction
nSmoothNormals 3;

// Smooth layer thickness over surface patches
nSmoothThickness 10;

// Stop layer growth on highly warped cells
maxFaceThicknessRatio 0.5;

// Reduce layer growth where ratio thickness to medial
// distance is large
maxThicknessToMedialRatio 0.3;

// Angle used to pick up medial axis points
minMedianAxisAngle 130;

// Create buffer region for new layer terminations
nBufferCellsNoExtrude 0;

// Overall max number of layer addition iterations
nLayerIter 10;
}


Quote:
Originally Posted by bruce View Post
Hi
what i understood when i had the same issue that the Layer parameters are inconsistant to get fast extrusion.

e.g too high layer size -> takes long time to push mesh
look at how many faces were stopped during extrusion? that will give some info.

Post your Layer.log file, small test case as well.

--
__________________
Thanks,
subhkirti is offline   Reply With Quote

Old   October 29, 2010, 08:19
Default process killed while layer generation
  #4
New Member
 
CFD user
Join Date: Apr 2010
Location: Germany
Posts: 28
Rep Power: 16
subhkirti is on a distinguished road
To add to my problem that i have above, the sMH works fine when I have a grid size of 1m in each dimension. However, this particular problem arose, when I refined the grid size to 0.75m or any further above.


Quote:
Originally Posted by subhkirti View Post
Hi bruce,

I am having the same problem as bastil. The entire process of mesh generation by sHM gets killed when it starts writing the layer information in the 3rd time directory.

Specifically,

Layer addition iteration 0
--------------------------
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.

Determining displacement for added points according to pointNormal ...
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.

shrinkMeshMedialDistance : Smoothing using Medial Axis ...
shrinkMeshMedialDistance : Reduce layer thickness at 1242 nodes where thickness to medial axis distance is large
shrinkMeshDistance : Removing isolated regions ...
Number isolated points extrusion stopped : 1
shrinkMeshDistance : Smoothing field ...
Iteration 0 residual 0.301
shrinkMeshMedialDistance : Moving mesh ...
Iteration 0
Moving mesh using diplacement scaling : min:1 max:1
Checking faces in error :
non-orthogonality > 65 degrees : 0
faces with face pyramid volume < 1e-13 : 0
faces with concavity > 80 degrees : 0
faces with skewness > 4 (internal) or 20 (boundary) : 0
faces with interpolation weights (0..1) < 0.02 : 0
faces with volume ratio of neighbour cells < 0.01 : 0
faces with face twist < 0.02 : 0
faces on cells with determinant < 0.001 : 0
shrinkMeshMedialDistance : Successfully moved mesh
shrinkMeshMedialDistance : Finished moving mesh ...
Prevented extrusion on 0 coupled patch points during syncPatchDisplacement.
truncateDisplacement : Unextruded 0 faces due to non-consecutive vertices being extruded.

Setting up information for layer truncation ...
Killed

I am not able to figure out the problem. Can you help me out with this?

The layer part of my sHMDict file is given below:

/ Settings for the layer addition.
addLayersControls
{

relativeSizes true;

// Per final patch (so not geometry!) the layer information
layers
{
gondel
{
nSurfaceLayers 5;
}
}

// Expansion factor for layer mesh
expansionRatio 1.15;

finalLayerThickness 0.2;

//- Wanted thickness of final added cell layer. If multiple layers
// is the
// thickness of the layer furthest away from the wall.
// Relative to undistorted size of cell outside layer.
finalLayerRatio 0.1;

//- Minimum thickness of cell layer. If for any reason layer
// cannot be above minThickness do not add layer.
// Relative to undistorted size of cell outside layer.
minThickness 0.01;

//- If points get not extruded do nGrow layers of connected faces that are
// also not grown. This helps convergence of the layer addition process
// close to features.
nGrow 1;


// Advanced settings

//- When not to extrude surface. 0 is flat surface, 90 is when two faces
// make straight angle.
featureAngle 90;

//- Maximum number of snapping relaxation iterations. Should stop
// before upon reaching a correct mesh.
nRelaxIter 10;

// Number of smoothing iterations of surface normals
nSmoothSurfaceNormals 1;

// Number of smoothing iterations of interior mesh movement direction
nSmoothNormals 3;

// Smooth layer thickness over surface patches
nSmoothThickness 10;

// Stop layer growth on highly warped cells
maxFaceThicknessRatio 0.5;

// Reduce layer growth where ratio thickness to medial
// distance is large
maxThicknessToMedialRatio 0.3;

// Angle used to pick up medial axis points
minMedianAxisAngle 130;

// Create buffer region for new layer terminations
nBufferCellsNoExtrude 0;

// Overall max number of layer addition iterations
nLayerIter 10;
}
__________________
Thanks,
subhkirti is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] Help with Snappy: no layers growing GianF OpenFOAM Meshing & Mesh Conversion 2 September 23, 2020 09:26
[snappyHexMesh] OpenFOAM SnappyHexMesh |snapping and Boundary Layer generation calculations. Srinivasan OpenFOAM Meshing & Mesh Conversion 0 April 8, 2019 06:24
[snappyHexMesh] snappy Error dhruv OpenFOAM Meshing & Mesh Conversion 2 June 29, 2012 05:55
[snappyHexMesh] Boundary layer generation in Snappy hex Aniruddha OpenFOAM Meshing & Mesh Conversion 0 September 28, 2011 06:25
errors Fahad Main CFD Forum 0 March 23, 2004 14:20


All times are GMT -4. The time now is 17:51.