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

[snappyHexMesh] snappyHexMesh does not create any mesh except one for the reference cell

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

Like Tree1Likes
  • 1 Post By Arman_N

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 16, 2019, 14:55
Default snappyHexMesh does not create any mesh except one for the reference cell
  #1
New Member
 
Arman N
Join Date: Mar 2019
Posts: 13
Rep Power: 7
Arman_N is on a distinguished road
Hello to all,

I am somewhat new to CFD modeling and mesh generation.

Basically, I am trying to generate a mesh for a simple porous geometry and have created an .stl file for it (file attached below). I have checked it with surfaceCheck and it says that the surface is closed.

However, when I run the snappyHexMesh command, it removes all the background mesh except for the cell which contains the point specified by locationInMesh in snappyHexMeshDict (figure attached below). It seems that somehow the rest of the sections of the mesh are unreachable and hence it removes them all.

I would really appreciate if anyone could help me find out what I am doing wrong and how to solve this since it seems to be a trivial problem.

Many thanks,
Arman


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

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

castellatedMesh true;
snap            false;
addLayers       false;

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

castellatedMeshControls
{
    maxLocalCells 500000;
    maxGlobalCells 2000000;
    minRefinementCells 0;
    maxLoadUnbalance 0.10;
    nCellsBetweenLevels 1;

    features
    (
        //{
        //    file "all.eMesh";
        //    level 0;
        //}
    );

    refinementSurfaces
    {
        all
        {
            level (0 0);
        }
    }

    resolveFeatureAngle 30;

    refinementRegions
    {	
    }

    locationInMesh (-0.00123 0.00003 0);
    allowFreeStandingZoneFaces true;
}

snapControls
{
    nSmoothPatch 3;
    tolerance 1.0;
    nSolveIter 300;
    nRelaxIter 5;
        nFeatureSnapIter 10;
        implicitFeatureSnap false;
        explicitFeatureSnap true;
        multiRegionFeatureSnap false;
}

addLayersControls
{
    relativeSizes true;

    layers
    {
        //"flange_.*"
        //{
        //    nSurfaceLayers 1;
        //}
    }

    expansionRatio 1.0;
    finalLayerThickness 0.3;
    minThickness 0.1;
    nGrow 0;
    featureAngle 0;
    nRelaxIter 5;
    nSmoothSurfaceNormals 1;
    nSmoothNormals 3;
    nSmoothThickness 10;
    maxFaceThicknessRatio 0.5;
    maxThicknessToMedialRatio 0.3;
    minMedianAxisAngle 90;
    nBufferCellsNoExtrude 0;
    nLayerIter 50;
    nRelaxedIter 20;
}

meshQualityControls
{
    #include "meshQualityDict"

    relaxed
    {
        maxNonOrtho 80;
    }

    maxBoundarySkewness 2;
    maxInternalSkewness 2;
    maxSkewness 2;
    maxConcave 80;
    minFlatness 0.5;
    minVol 0;
    minTetQuality 1e-10;
    minArea 0;
    minTwist 0.5;
    minDeterminant 0.25;
    minFaceWeight 0.2;
    minVolRatio 0.5;
    minTriangleTwist -1;
    nSmoothScale 4;
    errorReduction 0.75;
}

writeFlags
(
    scalarLevels
    layerSets
    layerFields
);

debug 0;

mergeTolerance 1E-6;

// ************************************************************************* //

snappyHexMesh output:
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1806                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v1806
Arch   : "LSB;label=32;scalar=64"
Exec   : snappyHexMesh
Date   : Apr 16 2019
Time   : 17:50:34
Host   : "Arman"
PID    : 16369
I/O    : uncollated
Case   : /media/arman/EDU/ProSoft/Installed/OpenFOAM/GeoChemFoam/run/validationCase/mesh_all
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

--> FOAM Warning : 
    From function static Foam::IOstreamOption::compressionType Foam::IOstreamOption::compressionEnum(const Foam::word&)
    in file db/IOstreams/IOstreams/IOstreamOption.C at line 86
    Unknown compression specifier 'uncompressed', assuming no compression
Create mesh for time = 0

Read mesh in = 0.03 s

Overall mesh bounding box  : (-0.0018 -0.00055 -2.5e-05) (0.0018 0.00055 2.5e-05)
Relative tolerance         : 1e-06
Absolute matching distance : 3.76464e-09

Reading refinement surfaces.
Read refinement surfaces in = 0.01 s

Reading refinement shells.
Read refinement shells in = 0 s

Setting refinement level of surface to be consistent with shells.
For geometry all.stl detected 0 uncached triangles out of 5000
Checked shell refinement in = 0.01 s

Reading features.
Read features in = 0 s


Determining initial surface intersections
-----------------------------------------

Edge intersection testing:
    Number of edges             : 25532
    Number of edges to retest   : 25532
    Number of intersected edges : 8130
Calculated surface intersections in = 0.07 s

Initial mesh : cells:6336  faces:25532  points:13050
Cells per refinement level:
    0	6336

Adding patches for surface regions
----------------------------------

Patch Type                Region                        
----- ----                ------                        
all:

1     wall                all                           

Added patches in = 0 s

Edge intersection testing:
    Number of edges             : 25532
    Number of edges to retest   : 0
    Number of intersected edges : 8130
Selecting decompositionMethod none [1]

Refinement phase
----------------

Found point (-0.00123 3e-05 0) in cell 3334 on processor 0

Surface refinement iteration 0
------------------------------

Marked for refinement due to surface intersection          : 0 cells.
Determined cells to refine in = 0.03 s
Selected for refinement : 0 cells (out of 6336)
Stopping refining since too few cells selected.


Removing mesh beyond surface intersections
------------------------------------------

Setting cellZones according to locationsInMesh:
Location : (-0.00123 3e-05 0)
    cellZone : none

For cellZone none found point (-0.00123 3e-05 0) in global region 1 out of 17 regions.
Keeping all cells containing inside points
Selected for keeping : 6336 cells.
Edge intersection testing:
    Number of edges             : 25532
    Number of edges to retest   : 0
    Number of intersected edges : 8130
Setting cellZones according to locationsInMesh:
Location : (-0.00123 3e-05 0)
    cellZone : none

For cellZone none found point (-0.00123 3e-05 0) in global region 1 out of 17 regions.
Keeping all cells containing inside points
Selected for keeping : 6336 cells.
Edge intersection testing:
    Number of edges             : 25532
    Number of edges to retest   : 0
    Number of intersected edges : 8130

Shell refinement iteration 0
----------------------------

Marked for refinement due to distance to explicit features : 0 cells.
Marked for refinement due to refinement shells             : 0 cells.
Determined cells to refine in = 0.13 s
Selected for internal refinement : 0 cells (out of 6336)
Stopping refining since too few cells selected.


Dangling coarse cells refinement iteration 0
--------------------------------------------

Determined cells to refine in = 0 s
Selected for refinement : 0 cells (out of 6336)
Stopping refining since too few cells selected.


Dangling coarse cells refinement iteration 0
--------------------------------------------

Determined cells to refine in = 0 s
Selected for refinement : 0 cells (out of 6336)
Stopping refining since too few cells selected.


Refinement transition refinement iteration 0
--------------------------------------------

Determined cells to refine in = 0 s
Selected for refinement : 0 cells (out of 6336)
Stopping refining since too few cells selected.


Splitting mesh at surface intersections
---------------------------------------

Introducing baffles for 8130 faces that are intersected by the surface.

Setting cellZones according to locationsInMesh:
Location : (-0.00123 3e-05 0)
    cellZone : none

For cellZone none found point (-0.00123 3e-05 0) in global region 1 out of 17 regions.
Edge intersection testing:
    Number of edges             : 26454
    Number of edges to retest   : 18630
    Number of intersected edges : 9052
Created baffles in = 0.11 s


After introducing baffles : cells:6336  faces:26454  points:13050
Cells per refinement level:
    0	6336

Introducing baffles to block off problem cells
----------------------------------------------

markFacesOnProblemCells : marked 6752 additional internal faces to be converted into baffles.
Analyzed problem cells in = 0.01 s


Introducing baffles to delete problem cells.

Edge intersection testing:
    Number of edges             : 33206
    Number of edges to retest   : 21695
    Number of intersected edges : 9052
Created baffles in = 0.05 s


After introducing baffles : cells:6336  faces:33206  points:13050
Cells per refinement level:
    0	6336
Setting cellZones according to locationsInMesh:
Location : (-0.00123 3e-05 0)
    cellZone : none

For cellZone none found point (-0.00123 3e-05 0) in global region 1956 out of 3626 regions.
Edge intersection testing:
    Number of edges             : 33206
    Number of edges to retest   : 25376
    Number of intersected edges : 9052

Remove unreachable sections of mesh
-----------------------------------

Keeping all cells containing points 1((-0.00123 3e-05 0))
Selected for keeping : 1 cells.
Edge intersection testing:
    Number of edges             : 6
    Number of edges to retest   : 0
    Number of intersected edges : 2
Split mesh in = 0.14 s


After subsetting : cells:1  faces:6  points:8
Cells per refinement level:
    0	1

Introducing zones for interfaces
--------------------------------

Setting cellZones according to locationsInMesh:
Location : (-0.00123 3e-05 0)
    cellZone : none

For cellZone none found point (-0.00123 3e-05 0) in global region 0 out of 1 regions.
Setting faceZones according to neighbouring cellZones:

Edge intersection testing:
    Number of edges             : 6
    Number of edges to retest   : 0
    Number of intersected edges : 2

Handling cells with snap problems
---------------------------------

Introducing baffles for 2 faces that are intersected by the surface.

Setting cellZones according to locationsInMesh:
Location : (-0.00123 3e-05 0)
    cellZone : none

For cellZone none found point (-0.00123 3e-05 0) in global region 0 out of 1 regions.
Edge intersection testing:
    Number of edges             : 6
    Number of edges to retest   : 6
    Number of intersected edges : 2
Created baffles in = 0 s


After introducing baffles : cells:1  faces:6  points:8
Cells per refinement level:
    0	1

Introducing baffles to block off problem cells
----------------------------------------------

markFacesOnProblemCells : marked 0 additional internal faces to be converted into baffles.
Analyzed problem cells in = 0 s


Introducing baffles to delete problem cells.

Created baffles in = 0 s


After introducing baffles : cells:1  faces:6  points:8
Cells per refinement level:
    0	1
Setting cellZones according to locationsInMesh:
Location : (-0.00123 3e-05 0)
    cellZone : none

For cellZone none found point (-0.00123 3e-05 0) in global region 0 out of 1 regions.
Edge intersection testing:
    Number of edges             : 6
    Number of edges to retest   : 6
    Number of intersected edges : 2

Remove unreachable sections of mesh
-----------------------------------

Split mesh in = 0 s


After subsetting : cells:1  faces:6  points:8
Cells per refinement level:
    0	1

Merge free-standing baffles
---------------------------

freeStandingBaffles : detected 0 free-standing baffles out of 0

Detected free-standing baffles : 0
Merged free-standing baffles in = 0 s


dupNonManifoldPoints : Found : 0 non-manifold points (out of 8)
Detected unsplittable baffles : 0

Merge refined boundary faces
----------------------------

Merging 0 sets of faces.

No faces merged ...

Merging all points on surface that
- are used by only two boundary faces and
- make an angle with a cosine of more than 0.707107.

No straight edges simplified and no points removed ...
Refined mesh : cells:1  faces:6  points:8
Cells per refinement level:
    0	1
Writing mesh to time 1
Wrote mesh in = 0 s.
Mesh refined in = 0.46 s.
Checking final mesh ...
Checking faces in error :
    non-orthogonality > 65  degrees                        : 0
    faces with face pyramid volume < 0                     : 0
    faces with face-decomposition tet quality < 1e-10      : 0
    faces with concavity > 80  degrees                     : 0
    faces with area < 0     m^2                            : 0
    faces with skewness > 2   (internal) or 2   (boundary) : 0
    faces with interpolation weights (0..1)  < 0.2         : 0
    faces with volume ratio of neighbour cells < 0.5       : 0
    faces with face twist < 0.5                            : 0
    faces on cells with determinant < 0.25                 : 0
Finished meshing without any errors
Finished meshing in = 0.59 s.
End
stl file.zip

figure.jpg

Last edited by Arman_N; April 17, 2019 at 02:49.
Arman_N is offline   Reply With Quote

Old   May 20, 2019, 18:16
Default Solved!
  #2
New Member
 
Arman N
Join Date: Mar 2019
Posts: 13
Rep Power: 7
Arman_N is on a distinguished road
UPDATE:

I managed to solve the problem by simply inverting the .stl file (attached below).

It seems that snappyHexMesh has trouble meshing a geometry with holes in it (i.e. porous media). So instead, the holes should be created into a .stl file and then removed from the background mesh using snappyHexMesh.

Cheers.
Attached Images
File Type: png Figure.png (59.3 KB, 93 views)
Kavi likes this.
Arman_N is offline   Reply With Quote

Reply

Tags
mesh generation, porous, snappyhexmesh, stl

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
LEMOS InflowGenerator r_gordon OpenFOAM Running, Solving & CFD 103 December 18, 2018 01:58
[OpenFOAM] snappyHexMesh - can't open mesh in ParaView tfigueiro ParaView 0 June 11, 2015 20:05
[snappyHexMesh] snappyHexMesh & Mesh around a missile sasanghomi OpenFOAM Meshing & Mesh Conversion 2 October 15, 2014 20:54
Where's the singularity/mesh flaw? audrich FLUENT 3 August 4, 2009 02:07
Error with Wmake skabilan OpenFOAM Installation 3 July 28, 2009 01:35


All times are GMT -4. The time now is 01:02.