[hjasak]

For some reason, your patch merging has failed. Try switching on debug for:

polyTopoChanger
polyMesh
enrichedPatch
slidingInterface

and see what it says.

Hrv

[r2d2]

Sorry about the possibly stupid question...but how do I do that?
Radu

[hjasak]

Look at ~/.OpenFOAM-1.3/controlDict

If you haven't got that, it will live in the .OpenFOAM-1.3/controlDict directory of your installation.

Hrv

[r2d2]

Hrv,
Did what you said and the message is "richer" now, but still does not give me much of a clue:
blockMesh . filter
/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.3 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : blockMesh . filter
Date : Nov 21 2006
Time : 14:09:16
Host : marte
PID : 25722
Root : /home/radu/OpenFOAM/radu-1.3/run/Filtering
Case : filter
Nprocs : 1
Create time


Reading block mesh description dictionary

Creating block mesh

Creating blockCorners

Creating curved edges

Creating blocks

Creating patches

Creating block mesh topology

Default patch type set to empty
Constructing polyMesh from cell and boundary shapes.
--> FOAM Warning :
From function polyMesh::polyMesh(... construct from shapes...)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 573
Found 121 undefined faces in mesh; adding to default patch.
void polyMesh::initMesh(cellList& c) : calculating faceCells
--> FOAM Serious Error :
From function bool primitiveMesh::checkClosedCells(const bool report, labelHashSet*) const
in file meshes/primitiveMesh/primitiveMeshCheck.C at line 277
12 open cells found. Max cell openness: 2.8744
Number of non-orthogonality errors: 0. Number of severely non-orthogonal faces: 2.
--> FOAM Serious Error :
From function bool primitiveMesh::checkFacePyramids(const bool, const scalar, labelHashSet*) const
in file meshes/primitiveMesh/primitiveMeshCheck.C at line 704
Error in face pyramids: faces pointing the wrong way!
--> FOAM Warning :
From function primitiveMesh::checkFaceSkewness(const bool report, labelHashSet* setPtr) const
in file meshes/primitiveMesh/primitiveMeshCheck.C at line 838
Large face skewness detected. Max skewness = 405.041 percent.
This may impair the quality of the result.
2 highly skew faces detected.
Failed 3 mesh geometry checks.
Failed some mesh checks.
Mesh OK

Check block mesh topology

Basic statistics
Number of internal faces : 64
Number of boundary faces : 154
Number of defined boundary faces : 154
Number of undefined boundary faces : 0

Checking patch -> block consistency

Creating block offsets

Creating merge list .

Creating points

Creating cells

Creating patches

Creating mesh from block mesh

Default patch type set to empty
Constructing polyMesh from cell and boundary shapes.
void polyMesh::initMesh(cellList& c) : calculating faceCells
Number of non-orthogonality errors: 0. Number of severely non-orthogonal faces: 65.

Creating merge patch pairs

Adding point and face zones
Creating attachPolyTopoChanger
Sliding interface object couple0 :
master face zone: 0
slave face zone: 1
void Foam::slidingInterface::calcAttachedAddressing() const for object couple0 : Calculating zone face-cell addressing.
void Foam::slidingInterface::calcAttachedAddressing() const for object couple0 : Finished calculating zone face-cell addressing.
Sliding interface object couple1 :
master face zone: 3
slave face zone: 4
void Foam::slidingInterface::calcAttachedAddressing() const for object couple1 : Calculating zone face-cell addressing.
void Foam::slidingInterface::calcAttachedAddressing() const for object couple1 : Finished calculating zone face-cell addressing.
Sliding interface object couple2 :
master face zone: 6
slave face zone: 7
void Foam::slidingInterface::calcAttachedAddressing() const for object couple2 : Calculating zone face-cell addressing.
void Foam::slidingInterface::calcAttachedAddressing() const for object couple2 : Finished calculating zone face-cell addressing.
Sliding interface object couple3 :
master face zone: 9
slave face zone: 10
void Foam::slidingInterface::calcAttachedAddressing() const for object couple3 : Calculating zone face-cell addressing.
void Foam::slidingInterface::calcAttachedAddressing() const for object couple3 : Finished calculating zone face-cell addressing.
Sliding interface object couple4 :
master face zone: 12
slave face zone: 13
void Foam::slidingInterface::calcAttachedAddressing() const for object couple4 : Calculating zone face-cell addressing.
void Foam::slidingInterface::calcAttachedAddressing() const for object couple4 : Finished calculating zone face-cell addressing.
bool slidingInterface::changeTopology() const for object couple0 : Couple-decouple mode.
bool slidingInterface::projectPoints() : for object couple0 : Projecting slave points onto master surface.
Number of hits in point projection: 25 out of 26 points.
bool slidingInterface::projectPoints() for object couple0 : Adjusting point projection for integral match: a done.
Number of adjusted points in projection: 1
... projection OK.
... point merge OK.
Number of merged master points: 10
Number of adjusted slave points: 16
Processing slave edges
+ + + + + + + + + + + + .+ m.+ m.+ m.+ m.+ m.+ m..+ m..+ m..+ mm..+ mm..+ mm..+ mm...+ mm...+ mm..+ mm..+ mm.+ m.+ m..+ m..+ m.+ m.+ m.+ m+ .+ m
bool slidingInterface::projectPoints() for object couple0 : Finished projecting points. Topology = (Detached interface) changing.
bool slidingInterface::changeTopology() const for object couple1 : Couple-decouple mode.
bool slidingInterface::projectPoints() : for object couple1 : Projecting slave points onto master surface.
Number of hits in point projection: 37 out of 42 points.
bool slidingInterface::projectPoints() for object couple1 : Adjusting point projection for integral match: aaaaa done.
Number of adjusted points in projection: 5
... projection OK.
... point merge OK.
Number of merged master points: 18
Number of adjusted slave points: 24
Processing slave edges
+ + + + + + + + + + + + + + + + + + + + + + + + .+ .+ + + .+ .+ + + + + .+ .+ + + + + .+ .+ + + + + .+ .+ + + + + .+ .+ + + .+ .+ + + +
bool slidingInterface::projectPoints() for object couple1 : Finished projecting points. Topology = (Detached interface) changing.
bool slidingInterface::changeTopology() const for object couple2 : Couple-decouple mode.
bool slidingInterface::projectPoints() : for object couple2 : Projecting slave points onto master surface.
Number of hits in point projection: 18 out of 18 points.
bool slidingInterface::projectPoints() for object couple2 : Adjusting point projection for integral match: done.
Number of adjusted points in projection: 0
... projection OK.
... point merge OK.
Number of merged master points: 10
Number of adjusted slave points: 8
Processing slave edges
+ + + + + + + + + + + + .+ .+ + + .+ .+ + + .+ .+ + + +
bool slidingInterface::projectPoints() for object couple2 : Finished projecting points. Topology = (Detached interface) changing.
bool slidingInterface::changeTopology() const for object couple3 : Couple-decouple mode.
bool slidingInterface::projectPoints() : for object couple3 : Projecting slave points onto master surface.
Number of hits in point projection: 10 out of 10 points.
bool slidingInterface::projectPoints() for object couple3 : Adjusting point projection for integral match: done.
Number of adjusted points in projection: 0
... projection OK.
... point merge OK.
Number of merged master points: 6
Number of adjusted slave points: 4
Processing slave edges
+ + + + + + + + .+ .+ + + +
bool slidingInterface::projectPoints() for object couple3 : Finished projecting points. Topology = (Detached interface) changing.
bool slidingInterface::changeTopology() const for object couple4 : Couple-decouple mode.
bool slidingInterface::projectPoints() : for object couple4 : Projecting slave points onto master surface.
Number of hits in point projection: 6 out of 6 points.
bool slidingInterface::projectPoints() for object couple4 : Adjusting point projection for integral match: done.
Number of adjusted points in projection: 0
... projection OK.
... point merge OK.
Number of merged master points: 4
Number of adjusted slave points: 2
Processing slave edges
+ + + + + + +
bool slidingInterface::projectPoints() for object couple4 : Finished projecting points. Topology = (Detached interface) changing.
void slidingInterface::coupleInterface(polyTopoChange& ref) : Coupling sliding interface couple0
eeeeeeeeeeeeeeee
Processing slave edges
+ n-n-n-n-n-uun-n-n-
+ n-n-n-n-n-n-uun-n-
+ n-n-un-n-n-uuuu
+ n-n-n-n-n-n-uun-n-
+ n-n-n-n-n-n-uun-n-
+ uun-n-un-n-n-uu
+ n-n-uun-n-n-n-n-n-
+ n-n-uun-n-n-n-n-n-
+ n-n-uuuun-n-un-
+ n-n-uun-n-n-n-n-n-
+ n-n-n-un-un-n-n-n-
+ un-n-un-n-u
.+ uun-un-n-un-n-n-
.+ n-un-un-un-n-n-n-n-n-u
.+ n-n-un-n-n-n-un-un-n-n-
.+ n-n-un-n-un-n-n-un-n-n-n-n-n-
.+ n-n-n-n-un-un-n-n-un-un-un-
.+ n-n-un-un-n-un-n-un-un-n-n-
..+ n-n-n-uun-n-n-n-n-n-n-n-n-n-un-un-un-n-
..+ n-n-uun-n-n-n-n-n-un-un-un-n-n-n-n-n-n-
..+ n-n-n-n-uun-n-n-n-n-n-n-n-n-un-un-un-n-
..+ n-n-n-n-uun-n-n-n-n-n-un-un-un-n-n-n-n-
..+ n-un-n-n-n-uuun-n-n-n-n-n-n-n-un-un-u
..+ un-n-n-n-n-n-n-uuun-n-n-n-n-un-un-un-
...+ n-un-un-n-un-n-n-n-n-n-uuun-n-n-n-n-n-n-n-n-n-n-u
...+ un-un-un-n-n-n-n-n-n-un-un-un-n-n-n-n-n-n-n-n-un-
..+ n-n-n-n-n-un-un-n-un-n-n-n-n-uun-n-n-n-
..+ n-n-un-un-un-n-n-n-n-un-n-un-n-n-n-n-n-
.+ n-n-n-n-n-un-un-n-n-un-uun-
.+ n-n-un-un-un-n-n-n-un-un-n-
..+ n-n-n-n-n-n-n-n-un-n-n-un-n-uun-un-n-n-
..+ n-n-un-un-n-n-uun-un-n-n-n-n-n-n-
.+ n-n-n-n-n-n-n-un-n-n-uun-n-n-
.+ n-n-n-un-n-un-n-un-n-n-
.+ un-n-n-n-n-n-n-un-n-uu
+ un-n-n-un-u
.+ n-n-n-un-n-n-uuu
Enriched patch support OK. Slave faces: 12 Master faces: 20
local: 4(8 13 14 9) one side: 0 other side: 0
local: 4(14 70 68 9) one side: 0 other side: 1
Finished face 0
local: 4(68 70 15 10) one side: 1 other side: 1
local: 4(15 71 69 10) one side: 1 other side: 2
Finished face 1
local: 4(69 71 16 11) one side: 2 other side: 2
local: 4(16 17 12 11) one side: 2 other side: 3
Finished face 2
local: 4(12 17 30 24) one side: 3 other side: 4
local: 4(30 74 72 24) one side: 3 other side: 5
Finished face 3
local: 4(72 74 31 25) one side: 4 other side: 5
local: 4(31 32 26 25) one side: 4 other side: 6
local: 4(32 75 73 26) one side: 4 other side: 7
Finished face 4
local: 4(73 75 33 27) one side: 5 other side: 7
local: 4(33 34 28 27) one side: 5 other side: 8
local: 4(34 35 29 28) one side: 5 other side: 9
Finished face 5
local: 4(29 35 48 42) one side: 6 other side: 10
local: 4(48 49 43 42) one side: 6 other side: 11
local: 4(49 78 76 43) one side: 6 other side: 12
Finished face 6
local: 4(76 78 50 44) one side: 7 other side: 12
local: 4(50 51 45 44) one side: 7 other side: 13
local: 4(51 79 77 45) one side: 7 other side: 14
Finished face 7
local: 4(77 79 52 46) one side: 8 other side: 14
local: 4(52 53 47 46) one side: 8 other side: 15
Finished face 8
local: 4(47 53 64 60) one side: 9 other side: 16
local: 4(64 82 80 60) one side: 9 other side: 17
Finished face 9
local: 4(80 82 65 61) one side: 10 other side: 17
local: 4(65 83 81 61) one side: 10 other side: 18
Finished face 10
local: 4(81 83 66 62) one side: 11 other side: 18
local: 4(66 67 63 62) one side: 11 other side: 19
Finished face 11
Finished face 12
Finished face 13
Finished face 14
Finished face 15
Finished face 16
Finished face 17
Finished face 18
Finished face 19
Finished face 20
Finished face 21
Finished face 22
Finished face 23
Finished face 24
Finished face 25
Finished face 26
Finished face 27
Finished face 28
Finished face 29
Finished face 30
Finished face 31
Number of orphaned faces: master = 8 out of 20 slave = 12 out of 12
Retired 26 out of 26 points.
void slidingInterface::coupleInterface(polyTopoChange& ref) : Finished coupling sliding interface couple0
void slidingInterface::coupleInterface(polyTopoChange& ref) : Coupling sliding interface couple1
eeeeeeeeeeeeeeeeeeeeeeee
Processing slave edges
+ n-n-uun-n-n-
+ un-uun-uu
+ n-n-n-n-n-n-uun-n-
+ n-un-n-n-n-uuuu
+ n-n-n-n-n-n-uun-n-
+ n-n-n-n-n-n-uun-n-
+ n-n-un-n-n-uuuu
+ n-n-n-n-n-n-uun-n-
+ n-n-n-n-n-n-uun-n-
+ n-n-un-n-n-uuuu
+ n-n-n-n-n-n-uun-n-
+ n-n-n-n-n-n-uun-n-
+ n-n-un-n-n-uuuu
+ n-n-n-n-n-n-n-n-uu
+ n-n-n-n-n-n-n-n-uu
+ uun-n-un-n-n-uu
+ uun-n-n-n-n-n-n-n-
+ uuun-un-n-un-n-
+ n-n-un-un-n-
+ n-n-uuun-n-
+ n-n-n-uun-n-
+ n-n-un-un-n-
+ un-n-un-n-u
+ un-un-n-un-
.+ un-n-n-un-n-un-n-
.+ un-n-n-un-n-un-n-n-n-n-
+ n-un-n-n-n-n-un-u
+ n-un-n-n-n-un-un-
.+ n-un-n-n-n-n-n-un-un-n-
.+ n-un-n-n-n-n-n-un-un-n-n-n-n-
+ n-n-n-n-n-n-n-un-n-
+ n-n-n-n-n-n-un-n-n-
+ n-n-un-n-n-n-un-u
+ n-n-un-n-n-un-un-
.+ n-n-un-n-n-n-n-n-n-n-un-un-n-
.+ n-n-un-n-n-n-n-un-un-n-n-n-n-
+ n-n-n-n-n-n-n-un-n-
+ n-n-n-n-n-n-un-n-n-
+ n-n-un-n-n-n-un-u
+ n-n-un-n-n-un-un-
.+ n-n-n-un-n-n-n-n-n-n-un-un-n-
.+ n-n-n-n-un-n-n-n-n-un-un-n-n-
+ n-n-n-n-n-n-n-un-n-
+ n-n-n-n-n-n-un-n-n-
+ n-n-un-n-n-n-un-u
+ n-n-un-n-n-un-un-
.+ n-n-n-n-n-un-n-n-n-n-n-n-un-u
.+ n-n-n-n-n-n-n-un-n-n-n-un-un-
+ n-n-n-n-n-n-n-n-n-u
+ n-n-n-n-n-n-n-n-un-
+ n-un-n-un-n-n-n-u
+ un-n-n-un-n-n-un-
.+ n-un-n-n-n-n-n-un-n-n-n-n-n-u
.+ un-n-n-n-n-n-un-n-n-un-
+ n-un-n-un-n-un-n-
+ un-un-n-n-n-un-n-
.+ n-un-n-n-un-n-un-n-n-n-
.+ un-un-n-n-n-un-n-
+ n-n-n-uun-n-
+ n-n-uuun-n-
+ n-n-uun-n-n-
Enriched patch support OK. Slave faces: 20 Master faces: 8
local: 4(42 45 62 60) one side: 0 other side: 0
Finished face 0
local: 4(60 62 46 43) one side: 1 other side: 0
Finished face 1
local: 4(43 46 63 61) one side: 2 other side: 1
Finished face 2
local: 4(61 63 47 44) one side: 3 other side: 1
Finished face 3
local: 4(44 47 68 64) one side: 4 other side: 2
Finished face 4
local: 4(64 68 69 65) one side: 5 other side: 2
Finished face 5
local: 4(65 69 50 48) one side: 6 other side: 2
Finished face 6
local: 4(48 50 70 66) one side: 7 other side: 3
Finished face 7
local: 4(66 70 71 67) one side: 8 other side: 3
Finished face 8
local: 4(67 71 51 49) one side: 9 other side: 3
Finished face 9
local: 4(49 51 76 72) one side: 10 other side: 4
Finished face 10
local: 4(72 76 77 73) one side: 11 other side: 4
Finished face 11
local: 4(73 77 54 52) one side: 12 other side: 4
Finished face 12
local: 4(52 54 78 74) one side: 13 other side: 5
Finished face 13
local: 4(74 78 79 75) one side: 14 other side: 5
Finished face 14
local: 4(75 79 55 53) one side: 15 other side: 5
Finished face 15
local: 4(53 55 82 80) one side: 16 other side: 6
Finished face 16
local: 4(80 82 58 56) one side: 17 other side: 6
Finished face 17
local: 4(56 58 83 81) one side: 18 other side: 7
Finished face 18
local: 4(81 83 59 57) one side: 19 other side: 7
Finished face 19
Finished face 20
Finished face 21
Finished face 22
Finished face 23
Finished face 24
Finished face 25
Finished face 26
Finished face 27
Number of orphaned faces: master = 8 out of 8 slave = 20 out of 20
Retired 42 out of 42 points.
void slidingInterface::coupleInterface(polyTopoChange& ref) : Finished coupling sliding interface couple1
void slidingInterface::coupleInterface(polyTopoChange& ref) : Coupling sliding interface couple2
eeeeeeee
Processing slave edges
+ n-n-uun-n-n-
+ un-uun-uu
+ n-n-n-n-n-n-uun-n-
+ n-un-n-n-n-uuuu
+ n-n-n-n-n-uun-n-n-
+ n-n-un-n-uuun-u
+ n-n-n-n-un-u
+ n-n-uuun-n-
+ n-n-n-uun-n-
+ n-n-un-un-n-
+ un-n-un-n-u
+ un-un-n-un-
.+ un-n-n-un-n-un-n-
.+ un-n-un-n-un-n-n-n-n-n-
+ n-un-n-n-n-n-un-u
+ n-un-n-n-n-un-un-
.+ n-un-n-n-n-n-un-un-n-n-
.+ n-un-n-n-n-un-un-n-n-n-
+ n-n-un-n-n-un-n-u
+ n-n-un-n-un-un-n-
.+ n-n-un-n-n-n-n-n-un-n-u
.+ n-n-un-n-un-un-n-
+ n-n-n-n-n-uu
+ n-n-n-n-uuu
+ n-n-n-n-uun-
Enriched patch support OK. Slave faces: 8 Master faces: 4
local: 4(18 20 29 28) one side: 0 other side: 0
Finished face 0
local: 4(28 29 21 19) one side: 1 other side: 0
Finished face 1
local: 4(19 21 31 30) one side: 2 other side: 1
Finished face 2
local: 4(30 31 23 22) one side: 3 other side: 1
Finished face 3
local: 4(22 23 33 32) one side: 4 other side: 2
Finished face 4
local: 4(32 33 25 24) one side: 5 other side: 2
Finished face 5
local: 4(24 25 35 34) one side: 6 other side: 3
Finished face 6
local: 4(34 35 27 26) one side: 7 other side: 3
Finished face 7
Finished face 8
Finished face 9
Finished face 10
Finished face 11
Number of orphaned faces: master = 4 out of 4 slave = 8 out of 8
Retired 18 out of 18 points.
void slidingInterface::coupleInterface(polyTopoChange& ref) : Finished coupling sliding interface couple2
void slidingInterface::coupleInterface(polyTopoChange& ref) : Coupling sliding interface couple3
eeee
Processing slave edges
+ n-uun-n-n-n-
+ uuun-un-u
+ n-n-n-n-un-u
+ n-uuun-n-n-
+ n-n-uun-n-n-
+ n-un-un-n-n-
+ un-un-n-n-u
+ uun-n-un-n-
.+ un-un-n-n-u
.+ uun-n-un-n-
+ n-n-n-n-n-uu
+ n-n-n-n-uuu
+ n-n-n-n-uun-
Enriched patch support OK. Slave faces: 4 Master faces: 2
local: 4(10 12 17 16) one side: 0 other side: 0
Finished face 0
local: 4(16 17 13 11) one side: 1 other side: 0
Finished face 1
local: 4(11 13 19 18) one side: 2 other side: 1
Finished face 2
local: 4(18 19 15 14) one side: 3 other side: 1
Finished face 3
Finished face 4
Finished face 5
Number of orphaned faces: master = 2 out of 2 slave = 4 out of 4
Retired 10 out of 10 points.
void slidingInterface::coupleInterface(polyTopoChange& ref) : Finished coupling sliding interface couple3
void slidingInterface::coupleInterface(polyTopoChange& ref) : Coupling sliding interface couple4
ee
Processing slave edges
+ un-un-
+ un-uu
+ n-n-uu
+ un-n-u
+ n-uun-
+ uuun-
+ uun-n-
Enriched patch support OK. Slave faces: 2 Master faces: 1
local: 4(6 8 11 10) one side: 0 other side: 0
Finished face 0
local: 4(10 11 9 7) one side: 1 other side: 0
Finished face 1
Finished face 2
Number of orphaned faces: master = 1 out of 1 slave = 2 out of 2
Retired 6 out of 6 points.
void slidingInterface::coupleInterface(polyTopoChange& ref) : Finished coupling sliding interface couple4
void polyMesh::clearGeom() : clearing geometric data
void polyMesh::clearAddressing() : clearing topology
void polyMesh::setInstance(const fileName& inst) : Resetting file instance to "0"
void polyMesh::initMesh() : calculating faceCells
void polyMesh::initMesh() : unused faces detected. Number of used faces: 9807. Total number of faces: 9876
void polyMesh::initMesh() : unused points detected. Number of used points: 5118. Total number of points: 5220


--> FOAM FATAL ERROR : Error in point ordering: mixed used and unused points at the end of point list.
Number of used points: 5118 and point 5214 is used by a live face.

From function void polyMesh::initMesh(cellList& c) const
in file meshes/polyMesh/polyMeshCalcFaceCells.C at line 149.

FOAM aborting

Foam::error::printStack(Foam:stream&)
Foam::error::abort()
Foam::polyMesh::initMesh()
Foam::polyMesh::resetPrimitives(int, Foam::Field<foam::vector<double> > const&, Foam::List<foam::face> const&, Foam::List<int> const&, Foam::List<int> const&, Foam::List<int> const&, Foam::List<int> const&, bool)
Foam::polyTopoChanger::changeMesh(Foam::polyMesh&, Foam::polyTopoChange const&)
Foam::polyTopoChanger::changeMesh()
Foam::attachPolyTopoChanger::attach(bool)
blockMesh [0x80718c2]
__libc_start_main
__gxx_personality_v0
Aborted

[hjasak]

The mesh checker reports 12 open cells found. Max cell openness: 2.8744 and Failed 3 mesh geometry checks even before you try coupling the mesh interfaces. This is not good - you should make sure the blocked mesh is correct before you continue.

Hrv

[r2d2]

You´re right. Should fix whatever I messed up in the ordering of the nodes of the faces that I want to merge. To start with...
Will let you know whether I sort it out.
Cheers,
Radu

[r2d2]

Hi again,
I went back to the blockMeshDct file and re-defined correctly the order of the nodes in the patches to be merged. That meant that I got rid of the open cells but I am still left with the same sort of problem.
I will just quote the error messages from blockMesh generation:

Constructing polyMesh from cell and boundary shapes.
--> FOAM Warning :
From function polyMesh::polyMesh(... construct from shapes...)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 573
Found 121 undefined faces in mesh; adding to default patch.
void polyMesh::initMesh(cellList& c) : calculating faceCells
Number of non-orthogonality errors: 0. Number of severely non-orthogonal faces: 2.
--> FOAM Warning :
From function primitiveMesh::checkFaceSkewness(const bool report, labelHashSet* setPtr) const
in file meshes/primitiveMesh/primitiveMeshCheck.C at line 838
Large face skewness detected. Max skewness = 405.041 percent.
This may impair the quality of the result.
2 highly skew faces detected.
Failed 1 mesh geometry checks.
Failed some mesh checks.
Mesh OK

Check block mesh topology

Basic statistics
Number of internal faces : 64
Number of boundary faces : 154
Number of defined boundary faces : 154
Number of undefined boundary faces : 0

Checking patch -> block consistency

Creating block offsets

Creating merge list .

Creating points

Creating cells

Creating patches

Creating mesh from block mesh

Default patch type set to empty
Constructing polyMesh from cell and boundary shapes.
void polyMesh::initMesh(cellList& c) : calculating faceCells
Number of non-orthogonality errors: 0. Number of severely non-orthogonal faces: 66.
....etc, etc..merging part and ending with:
Processing slave edges
+ un-un-
+ un-uu
+ n-n-uu
+ un-n-u
+ n-uun-
+ uuun-
+ uun-n-
Enriched patch support OK. Slave faces: 2 Master faces: 1
local: 4(6 8 11 10) one side: 0 other side: 0
Finished face 0
local: 4(10 11 9 7) one side: 1 other side: 0
Finished face 1
Finished face 2
Number of orphaned faces: master = 1 out of 1 slave = 2 out of 2
Retired 6 out of 6 points.
void slidingInterface::coupleInterface(polyTopoChange& ref) : Finished coupling sliding interface couple4
void polyMesh::clearGeom() : clearing geometric data
void polyMesh::clearAddressing() : clearing topology
void polyMesh::setInstance(const fileName& inst) : Resetting file instance to "0"
void polyMesh::initMesh() : calculating faceCells
void polyMesh::initMesh() : unused faces detected. Number of used faces: 9807. Total number of faces: 9876
void polyMesh::initMesh() : unused points detected. Number of used points: 5118. Total number of points: 5220


--> FOAM FATAL ERROR : Error in point ordering: mixed used and unused points at the end of point list.
Number of used points: 5118 and point 5214 is used by a live face.

From function void polyMesh::initMesh(cellList& c) const
in file meshes/polyMesh/polyMeshCalcFaceCells.C at line 149.

FOAM aborting

Foam::error::printStack(Foam:stream&)
Foam::error::abort()
Foam::polyMesh::initMesh()
Foam::polyMesh::resetPrimitives(int, Foam::Field<foam::vector<double> > const&, Foam::List<foam::face> const&, Foam::List<int> const&, Foam::List<int> const&, Foam::List<int> const&, Foam::List<int> const&, bool)
Foam::polyTopoChanger::changeMesh(Foam::polyMesh&, Foam::polyTopoChange const&)
Foam::polyTopoChanger::changeMesh()
Foam::attachPolyTopoChanger::attach(bool)
blockMesh [0x80718c2]
__libc_start_main
__gxx_personality_v0
Aborted

Well..is that geometry check that gives me
trouble..and if so, is there a way to locate it?

[r2d2]

HI,
I think Hrvoje answered this, something like you cannot do two mesh modifiers at a time...as in the patches that you want to merge are intersecting..so you probably have to do one direction first and then with stitchMesh the remaining pair.
Am I right, Hrv?
If I find what I did sometime back I will post it here
Radu

[lord_kossity]

Hello Message Board,

I'm trying to use the function mergePatchPairs in order to solve the problem of hanging nodes.

The case can be described as follows:
Originally we have a closed channel with a square cylinder positioned vertically to the flow direction and is mounted on the front and back wall of the channel. In order to derive adequate results I'm defining several domains around the square cylinder becoming coarser from the inside to the outside domain.
___ ______________ ___
|___|_______1______|___|
|...|_2_|___3__|_4_|...|
|...|...|square|...|...|
|...|...|cylind|...|...|
|...|___|______|___|...|
|___|_6_|___7__|_8_|___|
|___|_______5______|___|

For this reason I'm defining a patch combined of 2,3,4 and merge it with 1 and also doing so for 6,7,8 and 5. That works fine!

BUT if I try to merge in the same principle in vertical direction, I receive the following error:

--> FOAM FATAL ERROR : Error in point ordering: mixed used and unused points at the end of point list.
Number of used points: 281841 and point 282051 is used by a live face.#0 Foam::error::printStack(Foam::-Ostream&) in "~/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"

What is going wrong here?
Is supposed, that there might be some problems on the corners that might be twice involved in merging (one time horizontal, other time vertical), therefore I changed the master to the slave patch and the other way round. but that didn't solve the problem either...

So please give me a hint!

Any help is appreciated!

Please excuse the bad picture

[nikhilmadduri]

hi everyone. i was trying to create a mesh around a solid cylinder.
can anyone suggest me which will be the best way? as im new to
openfoam, i started with quarter part of the whole mesh to be
generated containg the quarter of the cylinder. i did it in
blockmesh...while executing it, i was finding the following errors.
can anyone suggest me what might be the bug in here?

/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : blockMesh
/home/nikhil/OpenFOAM/nikhil-1.4.1/run/tutorials/icoFoam cavi
tyGrade
Date : Jun 24 2008
Time : 09:42:12
Host : localhost
PID : 22816
Root : /home/nikhil/OpenFOAM/nikhil-1.4.1/run/tutorials/icoFoam
Case : cavityGrade
Nprocs : 1
Create time

Reading block mesh description dictionary

Creating block mesh

Creating blockCorners

Creating curved edges

Creating blocks

Creating patches

Creating block mesh topology

Default patch type set to empty

Check block mesh topology

Basic statistics
Number of internal faces : 0
Number of boundary faces : 12
Number of defined boundary faces : 12
Number of undefined boundary faces : 0

Checking patch -> block consistency

Creating block offsets

Creating merge list . Creating points

Creating cells

Creating patches

Creating mesh from block mesh

Default patch type set to empty

Creating merge patch pairs

Adding point and face zones
Creating attachPolyTopoChanger
#0 Foam::error::printStack(Foam:stream&) in
"/home/nikhil/OpenFOAM/OpenFOAM-1
.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in
"/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/li
b/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe420]
#3 Foam::triangle, Foam::Vector const&>::ray(Foam:
:Vector const&, Foam::Vector const&, Foam::intersection::algorit
hm, Foam::intersection::direction) const in
"/home/nikhil/OpenFOAM/OpenFOAM-1.4.
1/lib/linuxGccDPOpt/libmeshTools.so"
#4 Foam::face::ray(Foam::Vector const&, Foam::Vector const&, Fo
am::Field > const&, Foam::intersection::algorithm, Foam::in
tersection::direction) const in
"/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/lib/linuxG
ccDPOpt/libOpenFOAM.so"
#5 Foam::List Foam::PrimitivePatch > const&>::projectPoints > const&>
>(Foam::Primitive Patch > const&> const&, Foam::Field > const&,
Foam::intersection::algorithm, Foam: :intersection::direction) const
in "/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/lib/lin
uxGccDPOpt/libdynamicMesh.so"
#6 Foam::slidingInterface::projectPoints() const in
"/home/nikhil/OpenFOAM/Open
FOAM-1.4.1/lib/linuxGccDPOpt/libdynamicMesh.so"
#7 Foam::slidingInterface::changeTopology() const in
"/home/nikhil/OpenFOAM/Ope
nFOAM-1.4.1/lib/linuxGccDPOpt/libdynamicMesh.so"
#8 Foam::polyTopoChanger::changeTopology() const in
"/home/nikhil/OpenFOAM/Open
FOAM-1.4.1/lib/linuxGccDPOpt/libdynamicMesh.so"
#9 Foam::polyTopoChanger::changeMesh(bool, bool) in
"/home/nikhil/OpenFOAM/Open
FOAM-1.4.1/lib/linuxGccDPOpt/libdynamicMesh.so"
#10 Foam::attachPolyTopoChanger::attach(bool) in
"/home/nikhil/OpenFOAM/OpenFOA
M-1.4.1/lib/linuxGccDPOpt/libdynamicMesh.so"
#11 main in
"/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOp
t/blockMesh"
#12 __libc_start_main in "/lib/i686/libc.so.6"
#13 Foam::regIOobject::readIfModified() in
"/home/nikhil/OpenFOAM/OpenFOAM-1.4.
1/applications/bin/linuxGccDPOpt/blockMesh"

regards,
nikhil.

[nikhilmadduri]

hi all.

the above problem i have posted was resolved but i have some other.

i have 4 blocks to be merged in a particular manner.
i have named the walls which are not to be merged with any other faces as 'patch0'(patch followed by a number representing the patch number). And i have named the patches which have to be merged as merge01-1 for example for the face which is involved in the merge between a face of block 0 and another face of block 1. and the last number '1' says that the face is belonging to block '1'. So this is the following blockMeshDict file obtained: (please look for problems in mergePatchPairs in particular. coz if i delete either second line or third line in mergePatchPairs, it was working fine at the cost of not being able to merge that deleted merger)

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

// FoamX Case Dictionary.

FoamFile
{
version 2.0;
format ascii;

root "/home/nikhil/OpenFOAM/nikhil-1.4.1/run/tutorials/icoFoam";
case "cavityGrade";
instance "constant/polyMesh";
local "";

class dictionary;
object blockMeshDict;
}

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

arguments "/home/nikhil/OpenFOAM/nikhil-1.4.1/run/tutorials/icoFoam/cavityGrade" off;

convertToMeters 0.001;

vertices
(
(0 0 0)
(10000 0 0)
(10000 10000 0)
(0 10000 0)
(0 0 9750)
(10000 0 9750)
(10000 10000 9750)
(0 10000 9750)
(0 0 9750)
(10000 0 9750)
(10000 10000 9750)
(0 10000 9750)
(0 0 10250)
(10000 0 10250)
(10000 10000 10250)
(0 10000 10250)
(10000 0 9750)
(20000 0 9750)
(20000 10000 9750)
(10000 10000 9750)
(10000 0 10250)
(20000 0 10250)
(20000 10000 10250)
(10000 10000 10250)
(10000 0 0)
(20000 0 0)
(20000 10000 0)
(10000 10000 0)
(10000 0 9750)
(20000 0 9750)
(20000 10000 9750)
(10000 10000 9750)
);

blocks
(
hex (0 1 2 3 4 5 6 7) (30 30 30) simpleGrading (1 1 1)
hex (8 9 10 11 12 13 14 15) (30 30 30) simpleGrading (1 1 1)
hex (16 17 18 19 20 21 22 23) (15 30 3) simpleGrading (1 1 1)
hex (24 25 26 27 28 29 30 31) (30 30 30) simpleGrading (1 1 1)
);

edges
(
line 0 1
line 3 2
line 7 6
line 4 5
line 0 3
line 1 2
line 5 6
line 4 7
line 0 4
line 1 5
line 2 6
line 3 7
line 8 9
line 11 10
line 15 14
line 12 13
line 8 11
line 9 10
line 13 14
line 12 15
line 8 12
arc 9 13 (9750 0 10000)
arc 10 14 (9750 10000 10000)
line 11 15
line 16 17
line 19 18
line 23 22
line 20 21
line 16 19
line 17 18
line 21 22
line 20 23
arc 16 20 (10250 0 10000)
line 17 21
line 18 22
arc 19 23 (10250 10000 10000)
line 24 25
line 27 26
line 31 30
line 28 29
line 24 27
line 25 26
line 29 30
line 28 31
line 24 28
line 25 29
line 26 30
line 27 31
);

patches
(
wall patch0
(
(0 1 5 4)
(0 4 7 3)
(3 7 6 2)
(1 0 3 2)
)
wall patch1
(
(9 13 12 8)
(8 12 15 11)
(11 15 14 10)
(10 14 13 9)
(12 13 14 15)
)
wall merge01_0
(
(4 5 6 7)
)
wall merge01_1
(
(9 8 11 10)
)
wall patch2
(
(16 20 23 19)
(19 23 22 18)
(18 22 21 17)
(17 21 20 16)
(20 21 22 23)
)
wall merge32_2
(
(16 19 18 17)
)
wall merge03_0
(
(2 6 5 1)
)
wall merge03_3
(
(24 28 31 27)
)
wall merge32_3
(
(30 31 28 29)
)
wall patch3
(
(25 29 28 24)
(26 30 29 25)
(27 31 30 26)
(25 24 27 26)
)
);

mergePatchPairs
(
(merge01_0 merge01_1)
(merge32_2 merge32_3)
(merge03_3 merge03_0)
);


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







and i was getting the following error in the terminal window. Cud anyone please suggest me what the problem might be:

/*---------------------------------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.4.1 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : blockMesh /home/nikhil/OpenFOAM/nikhil-1.4.1/run/tutorials/icoFoam cavityGrade
Date : Jun 30 2008
Time : 09:21:11
Host : localhost
PID : 23034
Root : /home/nikhil/OpenFOAM/nikhil-1.4.1/run/tutorials/icoFoam
Case : cavityGrade
Nprocs : 1
Create time


Reading block mesh description dictionary

Creating block mesh

Creating blockCorners

Creating curved edges

Creating blocks

Creating patches

Creating block mesh topology

Default patch type set to empty

Check block mesh topology

Basic statistics
Number of internal faces : 0
Number of boundary faces : 24
Number of defined boundary faces : 24
Number of undefined boundary faces : 0

Checking patch -> block consistency

Creating block offsets

Creating merge list .

Creating points

Creating cells

Creating patches

Creating mesh from block mesh

Default patch type set to empty

Creating merge patch pairs

Adding point and face zones
Creating attachPolyTopoChanger


--> FOAM FATAL ERROR : Error in point ordering: mixed used and unused points at the end of point list.
Number of used points: 88939 and point 88987 is used by a live face.#0 Foam::error::printStack(Foam::-Ostream&) in "/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::polyMesh::initMesh() in "/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#3 Foam::polyMesh::resetPrimitives(int, Foam::Field<foam::vector<double> > const&, Foam::List<foam::face> const&, Foam::List<int> const&, Foam::List<int> const&, Foam::List<int> const&, Foam::List<int> const&, bool) in "/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::polyTopoChange::changeMesh(Foam::polyMesh&, bool, bool) in "/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libdynamicMesh.so"
#5 Foam::polyTopoChanger::changeMesh(bool, bool) in "/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libdynamicMesh.so"
#6 Foam::attachPolyTopoChanger::attach(bool) in "/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libdynamicMesh.so"
#7 main in "/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/blockMesh"
#8 __libc_start_main in "/lib/i686/libc.so.6"
#9 Foam::regIOobject::readIfModified() in "/home/nikhil/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/blockMesh"


From function void polyMesh::initMesh()
in file meshes/polyMesh/polyMeshCalcFaceCells.C at line 153.

FOAM aborting




Thank you,
nikhil.

[mgz1985]

HI,

I am facing a problem while creating a mesh for a deflected airfoil. I have built the mesh in simple blocks. CheckMesh does not display any error but when I run icoFoam, I get the following error

This mesh contains patches of type empty but is not 1D or 2D
by virtue of the fact that the number of faces of this
empty patch is not divisible by the number of cells.

I had this error previously and it was due to defining an internal face at the boundaries. I have already checked for this error. It is not the case this time. Then individually all my blocks are fine and are 2-D. But when I run this file with blockMesh I get the message of 3-D mesh with refineMesh. For you view I am attaching my file with this mail.

Can somebody please help?

blockMeshDict