|
[Sponsors] |
October 23, 2015, 14:29 |
GmshTo Foam Warning
|
#1 |
New Member
vito
Join Date: Jun 2015
Posts: 22
Rep Power: 10 |
Hello Guys,
I have been using Gmsh for sometime to simulate flow past a cylinder and I keep getting this warning that I have posted below. I know that this has been posted separately in a thread and an answer was to check if all the external faces of the boundary have been defined or not. I have already checked that and ensured that all the boundary faces have been defined. I also tried visualizing this cell set called Internal and writing it to VTK and trying to view it in Parafoam. but unfortunately nothing turns up as in it is blank althiough the vtk file is quite huge..(by using the following command) foamToVTK -cellSet Internal My check mesh is ok and simulation runs fine but for the love of god I am not able to understand why this error is coming up. If some body can tell me what files they need me to upload I can do that and await suggestions. I am not ablr to upload the vtk file as it is too big. Also I am not able to find the default patch FOAM Warning : From function polyMesh:olyMesh(... construct from shapes...) in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 627 Found 220340 undefined faces in mesh; adding to default patch. |
|
October 23, 2015, 15:01 |
|
#2 |
Senior Member
|
Hi,
If you open gmshToFoam.C, go to line 869, you find commentary there: Code:
// Problem is that the orientation of the patchFaces does not have to // be consistent with the outwards orientation of the mesh faces. So // we have to construct the mesh in two stages: // 1. define mesh with all boundary faces in one patch // 2. use the read patchFaces to find the corresponding boundary face // and repatch it. Code:
faceListList boundaryFaces(patchFaces.size()); ... polyMesh mesh ( ... ); ... forAll(patchFaces, patchI) { const DynamicList<face>& pFaces = patchFaces[patchI]; Info<< "Finding faces of patch " << patchI << endl; forAll(pFaces, i) { const face& f = pFaces[i]; // Find face in pp using all vertices of f. label patchFaceI = findFace(pp, f); if (patchFaceI != -1) { label meshFaceI = pp.start() + patchFaceI; repatcher.changePatchID(meshFaceI, patchI); } ... } ... } |
|
October 24, 2015, 11:51 |
|
#3 |
New Member
vito
Join Date: Jun 2015
Posts: 22
Rep Power: 10 |
Thank you very Much alexeym.
So I guess all my faces found their owners .. Thank you much again. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] mesh airfoil NACA0012 | anand_30 | OpenFOAM Meshing & Mesh Conversion | 13 | March 7, 2022 17:22 |
[blockMesh] error message with modeling a cube with a hold at the center | hsingtzu | OpenFOAM Meshing & Mesh Conversion | 2 | March 14, 2012 09:56 |
[blockMesh] BlockMesh FOAM warning | gaottino | OpenFOAM Meshing & Mesh Conversion | 7 | July 19, 2010 14:11 |
Version 15 on Mac OS X | gschaider | OpenFOAM Installation | 113 | December 2, 2009 10:23 |
[blockMesh] Axisymmetrical mesh | Rasmus Gjesing (Gjesing) | OpenFOAM Meshing & Mesh Conversion | 10 | April 2, 2007 14:00 |