CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[Gmsh] boundaries with gmshToFoam‏

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By ericnutsch

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 9, 2009, 09:10
Default boundaries with gmshToFoam‏
  #1
New Member
 
ouafa
Join Date: Jul 2009
Posts: 15
Rep Power: 17
ouafa is on a distinguished road
hello,

I try to generate a mesh for use in OpenFoam by gmshToFoam conversion, but I alwas get the warning:

--> FOAM Warning :
From function polyMesh:olyMesh(... construct from shapes...)
in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 576
Found 28864 undefined faces in mesh; adding to default patch.
Finding faces of patch 0
Finding faces of patch 1
Finding faces of patch 2
--> FOAM Warning : Not using gmsh face 4(1446 14016 14017 1447) since zero vertex is not on boundary of polyMesh
.........

However i got polymesh directory with all files. but the faces or boundaries are not reading.
patches in boundary files are even taking nFace=0 even if startface#0 in some patches.
it seems to be arbitrary because i change many times the order and the way to define the physical surface and the result is different all the time.

here attached my script. and checkMesh outline.

Any idea what I am doing wrong? Thanks


Ouafa

lc = 5.;
N1 = 10;
N2 = N1;
N3 = 2*N1;
arctheta1 = 5;
arctheta2 = 10;
f = 1.1;
h1 = 10;
Nz1 = 10;
d1=7.0711;
d2=14.1421;
d3=16.2635;
d4=50;
//---------------------------------------------------------------------------
Point(1)={0, 0, 0, lc};
Point(2)={-d1, 0, d1, lc};
Point(3)={d1, 0, d1, lc};
Point(4)={d1, 0, -d1, lc};
Point(5)={-d1, 0, -d1, lc};
Point(6)={-d2, 0, d2, lc};
Point(7)={d2, 0, d2, lc};
Point(8)={d2, 0, -d2, lc};
Point(9)={-d2, 0, -d2, lc};
Point(10)={-d3, 0, d3, lc};
Point(11)={d3, 0, d3, lc};
Point(12)={d3, 0, -d3, lc};
Point(13)={-d3, 0, -d3, lc};
Point(14)={-d4, 0, d4, lc};
Point(15)={-d3, 0, d4, lc};
Point(16)={d3, 0, d4, lc};
Point(17)={d4, 0, d4, lc};
Point(18)={d4, 0, d3, lc};
Point(19)={d4, 0, -d3, lc};
Point(20)={d4, 0, -d4, lc};
Point(21)={d3, 0, -d4, lc};
Point(22)={-d3, 0, -d4, lc};
Point(23)={-d4, 0, -d4, lc};
Point(24)={-d4, 0, -d3, lc};
Point(25)={-d4, 0, d3, lc};
Point(26)={100, 0, d4, lc};
Point(27)={100, 0, d3, lc};
Point(28)={100, 0, -d3, lc};
Point(29)={100, 0, -d4, lc};
Point(30)={-100, 0, d4, lc};
Point(31)={-100, 0, d3, lc};
Point(32)={-100, 0, -d3, lc};
Point(33)={-100, 0, -d4, lc};
Point(34)={150, 0, d4, lc};
Point(35)={150, 0, d3, lc};
Point(36)={150, 0, -d3, lc};
Point(37)={150, 0, -d4, lc};
//---------------------------------------------------------------------------
Line(1) = {2,3};
Line(2) = {3,4};
Line(3) = {4,5};
Line(4) = {5,2};
Circle(5) = {6,1,7};
Circle(6) = {7,1,8};
Circle(7) = {8,1,9};
Circle(8) = {9,1,6};
Line(9 ) = {14,15};
Line(10) = {15,16};
Line(11) = {16,17};
Line(12) = {17,18};
Line(13) = {18,19};
Line(14) = {19,20};
Line(15) = {20,21};
Line(16) = {21,22};
Line(17) = {22,23};
Line(18) = {23,24};
Line(19) = {24,25};
Line(20) = {25,14};
Line(32) = {21,12};
Line(33) = {24,13};
Line(34) = {13,22};
Line(35) = {15,10};
Line(36) = {10,25};
Line(37) = {18,11};
Line(38) = {11,16};
Line(21) = {17,26};
Line(22) = {26,27};
Line(23) = {27,28};
Line(24) = {28,29};
Line(25) = {29,20};
Line(41) = {27,18};
Line(42) = {28,19};
Line(31) = {12,19};
Line(26) = {23,33};
Line(27) = {33,32};
Line(28) = {32,31};
Line(29) = {31,30};
Line(30) = {30,14};
Line(39) = {25,31};
Line(40) = {24,32};
Line(43) = {26,34};
Line(44) = {27,35};
Line(45) = {28,36};
Line(46) = {29,37};
Line(47) = {34,35};
Line(48) = {35,36};
Line(49) = {36,37};
Line(50) = {2,6};
Line(51) = {3,7};
Line(52) = {4,8};
Line(53) = {5,9};
Line(54) = {6,10};
Line(55) = {7,11};
Line(56) = {8,12};
Line(57) = {9,13};
Circle(58) = {10,1,11};
Circle(59) = {10,1,13};
Circle(60) = {13,1,12};
Circle(61) = {12,1,11};
//---------------------------------------------------------------------------
Line Loop(1) = {2,3,4,1};
Line Loop(2) = {-51,-1,50,5};
Line Loop(3) = {6,-52,-2,51};
Line Loop(4) = {52,7,-53,-3};
Line Loop(5) = {53,8,-50,-4};
Line Loop(6) = {-55,-5,54,58};
Line Loop(7) = {-61,-56,-6,55};
Line Loop(8) = {56,-60,-57,-7};
Line Loop(9) = {-8,57,-59,-54};
Line Loop(10) = {13,-31,61,-37};
Line Loop(11) = {-32,16,-34,60};
Line Loop(12) = {59,-33,19,-36};
Line Loop(13) = {-38,-58,-35,10};
Line Loop(14) = {12,37,38,11};
Line Loop(15) = {14,15,32,31};
Line Loop(16) = {34,17,18,33};
Line Loop(17) = {35,36,20,9};
Line Loop(18) = {22,41,-12,21};
Line Loop(19) = {23,42,-13,-41};
Line Loop(20) = {24,25,-14,-42};
Line Loop(21) = {-20,39,29,30};
Line Loop(22) = {-19,40,28,-39};
Line Loop(23) = {-18,26,27,-40};
Line Loop(24) = {47,-44,-22,43};
Line Loop(25) = {48,-45,-23,44};
Line Loop(26) = {49,-46,-24,45};
//---------------------------------------------------------------------------
s1 = newreg; Plane Surface(s1) = {1};
s2 = newreg; Plane Surface(s2) = {2};
s3 = newreg; Plane Surface(s3) = {3};
s4 = newreg; Plane Surface(s4) = {4};
s5 = newreg; Plane Surface(s5) = {5};
s6 = newreg; Plane Surface(s6) = {6};
s7 = newreg; Plane Surface(s7) = {7};
s8 = newreg; Plane Surface(s8) = {8};
s9 = newreg; Plane Surface(s9) = {9};
s10 = newreg; Plane Surface(s10) = {10};
s11 = newreg; Plane Surface(s11) = {11};
s12 = newreg; Plane Surface(s12) = {12};
s13 = newreg; Plane Surface(s13) = {13};
s14 = newreg; Plane Surface(s14) = {14};
s15 = newreg; Plane Surface(s15) = {15};
s16 = newreg; Plane Surface(s16) = {16};
s17 = newreg; Plane Surface(s17) = {17};
s18 = newreg; Plane Surface(s18) = {18};
s19 = newreg; Plane Surface(s19) = {19};
s20 = newreg; Plane Surface(s20) = {20};
s21 = newreg; Plane Surface(s21) = {21};
s22 = newreg; Plane Surface(s22) = {22};
s23 = newreg; Plane Surface(s23) = {23};
s24 = newreg; Plane Surface(s24) = {24};
s25 = newreg; Plane Surface(s25) = {25};
s26 = newreg; Plane Surface(s26) = {26};
//---------------------------------------------------------------------------
Transfinite Line{54,55,56,57}=arctheta1 Using Progression f;
Transfinite Line{50,51,52,53}=arctheta2 Using Progression 1/f;
Transfinite Line {11,37,31,15,9,36,33,17} = N1 Using Progression 1.0;
Transfinite Line {20,19,18,35,59,34,38,61,32,12,13,14,58,10,16,5,7, 60,61,6,8,1,2,3,4} = N1 Using Progression 1.0;
Transfinite Line {20,35,38,12} = N1 Using Progression 1.0;
Transfinite Line {19,13} = N1 Using Progression 1.0;
Transfinite Line {18,34,32,14} = N1 Using Progression 1.0;
Transfinite Line {21,41,42,25,30,39,40,26} = N3 Using Progression 1.0;
Transfinite Line {22,12,20,23,13,19,24,14,18} = N1 Using Progression 1.0;
Transfinite Line {43,44,45,46,30,39,40,26} = N3 Using Progression 1.0;
Transfinite Line {47,48,49,22,23,24,20,19,18,29,28,27} = N1 Using Progression 1.0;
Transfinite Surface {s1} = {2,3,4,5};
Transfinite Surface {s2} = {2,6,7,3};
Transfinite Surface {s3} = {3,7,8,4};
Transfinite Surface {s4} = {4,8,9,5};
Transfinite Surface {s5} = {5,9,6,2};
Transfinite Surface {s6} = {11,7,6,10};
Transfinite Surface {s7} = {11,7,8,12};
Transfinite Surface {s8} = {8,12,13,9};
Transfinite Surface {s9} = {9,13,10,6};
Transfinite Surface {s10} = {11,18,19,12};
Transfinite Surface {s11} = {12,21,22,13};
Transfinite Surface {s12} = {13,24,25,10};
Transfinite Surface {s13} = {11,10,15,16};
Transfinite Surface {s14} = {18,11,16,17};
Transfinite Surface {s15} = {12,19,20,21};
Transfinite Surface {s16} = {13,22,23,24};
Transfinite Surface {s17} = {10,25,14,15};
Transfinite Surface {s18} = {17,26,27,18};
Transfinite Surface {s19} = {18,27,28,19};
Transfinite Surface {s20} = {19,28,29,20};
Transfinite Surface {s21} = {14,25,31,30};
Transfinite Surface {s22} = {25,24,32,31};
Transfinite Surface {s23} = {23,24,32,33};
Transfinite Surface {s24} = {34,26,27,35};
Transfinite Surface {s25} = {36,35,27,28};
Transfinite Surface {s26} = {37,36,28,29};
Recombine Surface {s1:s26};
//---------------------------------------------------------------------------
i = 1;
For s In {s1:s26}
out1[] = Extrude {0,h1,0} {
Surface { s }; Layers{10}; Recombine;
};
ns[i] = out1[0];
fg[i] = out1[2];
fb[i] = out1[3];
fd[i] = out1[4];
fh[i] = out1[5];
nv[i] = newreg;
Physical Volume (nv[i]) = {out1[1]};
i += 1;
EndFor

//---------------------------------------------------------------------------
Physical Surface("inlet") = {fg[1]};
Physical Surface("outlet") = {fd[1]};
Physical Surface("interface") = {s10:s20};
Physical Surface("jonction") = {s1:s9};
Physical Surface("wall") = {ns[1],ns[2],ns[3],ns[4],ns[5],ns[6],ns[7],ns[8],
ns[9],ns[10],ns[11],ns[12],ns[13],ns[14],ns[15],ns[16],
ns[17],ns[18],ns[19],ns[20],ns[21],ns[22],ns[23],ns[24],ns[25],ns[26]};







/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : checkMesh
Date : Jul 09 2009
Time : 15:36:04
Host : lambda31
PID : 13084
Case : /home/flurec/commun/ouafa/tmpgmsh
nProcs : 1
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create polyMesh for time = constant
Time = constant
Mesh stats
points: 135894
faces: 378092
internal faces: 349228
cells: 121220
boundary patches: 6
point zones: 0
face zones: 2
cell zones: 26
Number of cells of each type:
hexahedra: 121220
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0
Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Topological cell zip-up check OK.
Face vertices OK.
Face-face connectivity OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface
jonction 2565 2604 ok (not multiply connected)
interface 5111 5302 ok (not multiply connected)
inlet 0 0 ok (empty)
outlet 0 0 ok (empty)
wall 12122 12354 ok (not multiply connected)
defaultFaces 9066 9450 ok (not multiply connected)
Checking geometry...
Domain bounding box: (-100 0 -50) (150 50 50)
Boundary openness (-7.18902e-18 3.05372e-15 9.36245e-17) OK.
Max cell openness = 2.67949e-16 OK.
Max aspect ratio = 26.7476 OK.
Minumum face area = 0.312226. Maximum face area = 9.50492. Face area magnitudes OK.
Min volume = 1.56113. Max volume = 15.7656. Total volume = 1.25e+06. Cell volumes OK.
Mesh non-orthogonality Max: 41.0303 average: 5.01782
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.468149 OK.
All angles in faces OK.
Face flatness (1 = flat, 0 = butterfly) : average = 1 min = 1
All face flatness OK.
Mesh OK.
End
ouafa is offline   Reply With Quote

Old   January 6, 2010, 19:04
Default
  #2
Senior Member
 
Eric Nutsch
Join Date: Sep 2009
Location: Eugene, Oregon USA
Posts: 113
Rep Power: 17
ericnutsch is on a distinguished road
Send a message via Skype™ to ericnutsch
Did you set your boundary types in the boundary file?

"yourcasefoler" > constant > polyMesh > boundary
siefer92 likes this.
ericnutsch is offline   Reply With Quote

Old   January 22, 2010, 04:04
Default
  #3
New Member
 
ouafa
Join Date: Jul 2009
Posts: 15
Rep Power: 17
ouafa is on a distinguished road
thanks for UR replay. yes i made every patch in boundary file ok. but now i resolved the pb it was related to numerotation of faces before conversion to OF format.
ouafa is offline   Reply With Quote

Old   January 22, 2010, 12:48
Default
  #4
Senior Member
 
Eric Nutsch
Join Date: Sep 2009
Location: Eugene, Oregon USA
Posts: 113
Rep Power: 17
ericnutsch is on a distinguished road
Send a message via Skype™ to ericnutsch
Yeah, since your post was in July i figured you had it figured out by now. Its always nice to hear the ending for reference sake

Thanks
ericnutsch is offline   Reply With Quote

Old   May 20, 2010, 14:46
Default
  #5
Member
 
Robin Gilbert
Join Date: Jan 2010
Posts: 66
Rep Power: 16
robingilbert is on a distinguished road
hey i have the same problem. how did u guys solve it??
robingilbert is offline   Reply With Quote

Old   May 20, 2010, 15:01
Default
  #6
Senior Member
 
Eric Nutsch
Join Date: Sep 2009
Location: Eugene, Oregon USA
Posts: 113
Rep Power: 17
ericnutsch is on a distinguished road
Send a message via Skype™ to ericnutsch
Hello robingilbert,

Quote:
Originally Posted by ouafa View Post
i resolved the problem, it was related to numerotation of faces before conversion to OF format.
I had the same problem in the past, and fixed mine by double checking the boundary conditions:

-----
Did you set your boundary types in the boundary file?
"yourcasefoler" > constant > polyMesh > boundary
-----

If that doesn't help post back with your specific problem.
ericnutsch is offline   Reply With Quote

Old   May 21, 2010, 05:00
Default
  #7
Member
 
Robin Gilbert
Join Date: Jan 2010
Posts: 66
Rep Power: 16
robingilbert is on a distinguished road
Hello Eric,

Thank you for your reply. This is what i did:
i made a geometry in gmsh and converted to openfoam using gmsh2tofoam.
i wanted to convert some of the internal faces to patch so as to apply boundary conditions on them. so i used the command splitMesh. ( i did it the proper way i suppose : ie. edited the boundary file, added 2 new patches with nFaces 0 ) when i ran a checkMesh it showed the error: boundary openness .
I am kind of stuck there. please advice. is there any other better (robust) way to apply boundary conditions on internal walls??
robingilbert is offline   Reply With Quote

Old   May 21, 2010, 12:43
Default
  #8
Senior Member
 
Eric Nutsch
Join Date: Sep 2009
Location: Eugene, Oregon USA
Posts: 113
Rep Power: 17
ericnutsch is on a distinguished road
Send a message via Skype™ to ericnutsch
I have run OF1.6 and OF1.5dev and I have use thed "gmshToFoam" command exclusively. gmsh2tofoam must be an old old version...

I have also never used splitMesh. I define all of my boundary names in gmsh with "Physical Surface", "Physical Volume", etc.

Once the mesh is created and all the boundaries are named you save the mesh, run gmshToFoam, then go in and change the boundary and "0" directory files to correspond with your mesh.

Hope that helps!
ericnutsch is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[mesh manipulation] Problem with using createPatchDict to define cyclic boundaries kaszt OpenFOAM Meshing & Mesh Conversion 0 April 1, 2016 21:18
2D Rotating Detonation Engine Periodic Boundaries whitet86 FLUENT 2 June 25, 2015 11:04
Setting Flow/Pressure Boundaries in Floworks Eran FloEFD, FloWorks & FloTHERM 3 August 11, 2009 04:23
periodic boundaries - flow through a net PK FLUENT 0 July 12, 2007 11:58
problems replacing old boundaries Jared Siemens 4 August 5, 2005 19:36


All times are GMT -4. The time now is 07:21.