[Hemmati]

Hello. i run blockMesh, and checkMesh is ok. but after run the snappyHexMesh -overwrite the che checkMesh is not ok, and 5 error is taken. please help me,
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 135138
faces: 330617
internal faces: 283767
cells: 98704
faces per cell: 6.22451
boundary patches: 7
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 87765
prisms: 1008
wedges: 0
pyramids: 0
tet wedges: 54
tetrahedra: 0
polyhedra: 9877
Breakdown of polyhedra by number of faces:
faces number of cells
5 172
6 2256
7 1530
8 120
9 4541
10 1
11 1
12 1114
15 114
18 28

Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
w-inlet 258 348 ok (non-closed singly connected)
outlet 258 348 ok (non-closed singly connected)
atmosphere 399 536 ok (non-closed singly connected)
walls 325 441 ok (non-closed singly connected)
Front 15118 16452 ok (non-closed singly connected)
Back 13924 15041 ok (non-closed singly connected)
Spillway_profile3 16568 18353 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (510.5 853.32 0) (514.5 856.5 0.3)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-3.8561e-17 -1.06038e-16 1.37107e-14) OK.
***Open cells found, max cell openness: 0.5, number of open cells 36
<<Writing 36 non closed cells to set nonClosedCells
<<Writing 9110 cells with high aspect ratio to set highAspectRatioCells
***Zero or negative face area detected. Minimum area: 0
<<Writing 18868 zero area faces to set zeroAreaFaces
Min volume = 2e-300. Max volume = 0.000120907. Total volume = 3.22843. Cell volumes OK.
Mesh non-orthogonality Max: 180 average: 28.9768
*Number of severely non-orthogonal faces: 3599.
***Number of non-orthogonality errors: 21337.
<<Writing 24936 non-orthogonal faces to set nonOrthoFaces
***Error in face pyramids: 93 faces are incorrectly oriented.
<<Writing 93 faces with incorrect orientation to set wrongOrientedFaces
***Max skewness = 1.67617e+298, 1477 highly skew faces detected which may impair the quality of the results
<<Writing 1477 skew faces to set skewFaces
Coupled point location match (average 0) OK.

Failed 5 mesh checks.

End

[arsalan.dryi]

Hi Mousa,
It seems there is something wrong with your snappyHexMeshDict, please provide the output of snappy command ( snappyHexMesh &> log.snappy) and your snappyHexMeshDict.
I think the mesh quality controls option is set to relax in your snappyHexMeshDict.

Regards,
Arsalan.

[Hemmati]

Hi Arasalan. Tank's for answer. snappyHexMesh &> log.snappy as fallows:

Create time

Create mesh for time = 0

Read mesh in = 1.34 s

Overall mesh bounding box : (510.5 853.32 0) (514.5 856.5 0.3)
Relative tolerance : 1e-06
Absolute matching distance : 5.11883e-06

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.
Checked shell refinement in = 0 s

Reading features.
Read features in = 0 s


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

Edge intersection testing:
Number of edges : 330617
Number of edges to retest : 330617
Number of intersected edges : 10308
Calculated surface intersections in = 2.02 s

Initial mesh : cells:98704 faces:330617 points:135138
Cells per refinement level:
0 27792
1 4445
2 30290
3 36177

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

Patch Type Region
----- ---- ------
Spillway:

6 wall Spillway_profile3

Added patches in = 0.01 s

Selecting decompositionMethod none

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

Found point (511 855 0.1) in cell 13570 on processor 0

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

Marked for refinement due to surface intersection : 0 cells.
Marked for refinement due to curvature/regions : 256 cells.
Determined cells to refine in = 0.64 s
Selected for refinement : 256 cells (out of 98704)
hexRef8 : Dumping cell as obj to "/home/mosa/Desktop/phdt1/cell_79157.obj"


--> FOAM FATAL ERROR:
cell 79157 of level 2 uses more than 8 points of equal or lower level
Points so far:8(43971 51028 53963 59172 115160 115161 115163 115165)

From function hexRef8::setRefinement(const labelList&, polyTopoChange&)
in file polyTopoChange/polyTopoChange/hexRef8.C at line 3743.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam220/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Foam::hexRef8::setRefinement(Foam::List<int> const&, Foam:olyTopoChange&) in "/opt/openfoam220/platforms/linuxGccDPOpt/lib/libdynamicMesh.so"
#3 Foam::meshRefinement::refine(Foam::List<int> const&) in "/opt/openfoam220/platforms/linuxGccDPOpt/lib/libautoMesh.so"
#4 Foam::meshRefinement::refineAndBalance(Foam::strin g const&, Foam::decompositionMethod&, Foam::fvMeshDistribute&, Foam::List<int> const&, double) in "/opt/openfoam220/platforms/linuxGccDPOpt/lib/libautoMesh.so"
#5 Foam::autoRefineDriver::surfaceOnlyRefine(Foam::re finementParameters const&, int) in "/opt/openfoam220/platforms/linuxGccDPOpt/lib/libautoMesh.so"
#6 Foam::autoRefineDriver::doRefine(Foam::dictionary const&, Foam::refinementParameters const&, bool, Foam::dictionary const&) in "/opt/openfoam220/platforms/linuxGccDPOpt/lib/libautoMesh.so"
#7
in "/opt/openfoam220/platforms/linuxGccDPOpt/bin/snappyHexMesh"
#8 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#9
in "/opt/openfoam220/platforms/linuxGccDPOpt/bin/snappyHexMesh"

[arsalan.dryi]

well, snappyHexMesh has failed and nothing could be said with these low details.
Please provide your sHMDict and a whole description of geometry, also with blockMeshDict.

Regards,
Arsalan.

[Hemmati]

Hi.
Geometry is steeped spillway. I guess the snappyHexMeshDict isnot set god. then i offered snappyHexMeshDict, BlockMeh, stl file, and iamage paraview.
best wishes

[Hemmati]

Hi.
Fortunately, the error was resolved. Stl files had not plotted in the origin.