[Raji]

Good morning People

I am trying to mesh a geometry with snappyhexmesh. That worked great. Now want to run pisoFoam and decompose that after. But OpenFoam then just show me that error:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.1-bcfaaa7b8660
Exec : decomposePar
Date : Jun 12 2016
Time : 09:27:47
Host : "pk"
PID : 21298
Case : /home/diana/OpenFOAM/diana-2.3.1/neuerVersuch
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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



Decomposing mesh region0

Create mesh

Calculating distribution of cells
Selecting decompositionMethod scotch

Finished decomposition in 3.89 s

Calculating original mesh data

Distributing cells to processors

Distributing faces to processors

Distributing points to processors

Constructing processor meshes

Processor 0
Number of cells = 635542
Number of faces shared with processor 1 = 1054
Number of processor patches = 1
Number of processor faces = 1054
Number of boundary faces = 121437


--> FOAM FATAL ERROR:
Cell 919contains face labels out of range: 6(2997 2998 2999 3000 -1 2992) Max face index = 2055946

From function polyMesh:olyMesh
(
const IOobject&,
const Xfer<pointField>&,
const Xfer<faceList>&,
const Xfer<cellList>&
)

in file meshes/polyMesh/polyMesh.C at line 654.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam:olyMesh:olyMesh(Foam::IOobject const&, Foam::Xfer<Foam::Field<Foam::Vector<double> > > const&, Foam::Xfer<Foam::List<Foam::face> > const&, Foam::Xfer<Foam::List<Foam::cell> > const&, bool) at ??:?
#3
at ??:?
#4
at ??:?
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6
at ??:?
^CAborted (core dumped)




Does anyone has an Idea what I made wrong?

[pupo]

Could you please run the comand "checkMesh. It will print out a log with possible errors in the Mesh.

[Raji]

diana@pk:~/OpenFOAM/diana-2.3.1/neuerVersuch2$ checkMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.1-bcfaaa7b8660
Exec : checkMesh
Date : Jun 12 2016
Time : 12:34:44
Host : "pk"
PID : 24811
Case : /home/diana/OpenFOAM/diana-2.3.1/neuerVersuch2
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 762410
faces: 2108375
internal faces: 1939810
cells: 673627
faces per cell: 6.009535
boundary patches: 16
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 629558
prisms: 30328
wedges: 0
pyramids: 0
tet wedges: 17
tetrahedra: 0
polyhedra: 13724
Breakdown of polyhedra by number of faces:
faces number of cells
4 53
5 799
6 608
7 460
9 11198
12 606

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
maxY 0 0 ok (empty)
minX 2432 2541 ok (non-closed singly connected)
maxX 2432 2541 ok (non-closed singly connected)
minY 0 0 ok (empty)
minZ 0 0 ok (empty)
maxZ 0 0 ok (empty)
inlet_1 912 1071 ok (non-closed singly connected)
inlet_2 912 1143 ok (non-closed singly connected)
inlet_3 912 1143 ok (non-closed singly connected)
inlet_4 905 1069 ok (non-closed singly connected)
inlet_5 906 1136 ok (non-closed singly connected)
inlet_6 904 1069 ok (non-closed singly connected)
inlet_7 858 1076 ok (non-closed singly connected)
inlet_8 858 1008 ok (non-closed singly connected)
outlet 1520 1617 ok (non-closed singly connected)
casing 155014 156939 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (0.277 0.20299501 0.74999845) (0.527 0.83880103 1.7500003)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (-3.6109514e-16 4.3811878e-15 1.4820177e-16) OK.
Max cell openness = 3.1918271e-16 OK.
Max aspect ratio = 5.5901407 OK.
Minimum face area = 1.2105523e-07. Maximum face area = 0.00030942757. Face area magnitudes OK.
Min volume = 3.543095e-10. Max volume = 3.2021256e-06. Total volume = 0.13153312. Cell volumes OK.
Mesh non-orthogonality Max: 60.461273 average: 5.1332881
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 3.1192572 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

diana@pk:~/OpenFOAM/diana-2.3.1/neuerVersuch2$

[Raji]

That´s the result of checkMesh.

[Atzori]

Hi!

Maybe I had a similar problem in the past, but I cannot remember well: are you allowed to share the case, to reproduce the issue?

BTW, did you see the post #17 (October 24, 2015, 10:08) in this thread? http://www.cfd-online.com/Forums/ope...out-range.html

Best Regards!

Marco Atzori