|
[Sponsors] |
November 21, 2006, 08:22 |
BlockMeshmergePatchPairs
|
#1 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,906
Rep Power: 33 |
For some reason, your patch merging has failed. Try switching on debug for:
polyTopoChanger polyMesh enrichedPatch slidingInterface and see what it says. Hrv
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk |
|
November 21, 2006, 08:29 |
Sorry about the possibly stupi
|
#2 |
Member
Radu Mustata
Join Date: Mar 2009
Location: Zaragoza, Spain
Posts: 99
Rep Power: 17 |
Sorry about the possibly stupid question...but how do I do that?
Radu |
|
November 21, 2006, 08:36 |
Look at ~/.OpenFOAM-1.3/contro
|
#3 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,906
Rep Power: 33 |
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
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk |
|
November 21, 2006, 09:08 |
Hrv,
Did what you said and th
|
#4 |
Member
Radu Mustata
Join Date: Mar 2009
Location: Zaragoza, Spain
Posts: 99
Rep Power: 17 |
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 |
|
November 21, 2006, 09:23 |
The mesh checker reports 12 op
|
#5 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,906
Rep Power: 33 |
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
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk |
|
November 21, 2006, 10:38 |
You´re right. Should fix whate
|
#6 |
Member
Radu Mustata
Join Date: Mar 2009
Location: Zaragoza, Spain
Posts: 99
Rep Power: 17 |
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 |
|
November 22, 2006, 07:59 |
Hi again,
I went back to the
|
#7 |
Member
Radu Mustata
Join Date: Mar 2009
Location: Zaragoza, Spain
Posts: 99
Rep Power: 17 |
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? |
|
May 16, 2007, 03:30 |
HI,
I think Hrvoje answered
|
#8 |
Member
Radu Mustata
Join Date: Mar 2009
Location: Zaragoza, Spain
Posts: 99
Rep Power: 17 |
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 |
|
April 14, 2008, 03:04 |
Hello Message Board,
I'm tr
|
#9 |
Member
Andreas Dietz
Join Date: Mar 2009
Location: Munich
Posts: 79
Rep Power: 17 |
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 |
|
June 24, 2008, 01:30 |
hi everyone. i was trying to c
|
#10 |
New Member
nikhil babu madduri
Join Date: Mar 2009
Posts: 17
Rep Power: 17 |
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. |
|
June 30, 2008, 01:56 |
hi all.
the above problem i
|
#11 |
New Member
nikhil babu madduri
Join Date: Mar 2009
Posts: 17
Rep Power: 17 |
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. |
|
August 15, 2008, 07:36 |
HI,
I am facing a problem w
|
#12 |
Senior Member
mayank gupta
Join Date: Mar 2009
Posts: 110
Rep Power: 17 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|