Hi,<br><br>I am getting an error while meshing a transfinite volume. When I comment out the transfinite volume lines, there is no error. Could you kindly take a look ?<br><br>Here is my "geo" file.<br><br>Point(1) = {-1,-1,0,0.2};<br>
Point(2) = {1,-1,0,0.2};<br>Point(3) = {1,1,0,0.2};<br>Point(4) = {-1,1,0,0.2};<br>Line(1) = {1,2};<br>Line(2) = {2,3};<br>Line(3) = {3,4};<br>Line(4) = {4,1};<br>Line Loop(5) = {4,1,2,3};<br>Plane Surface(6) = {5};<br>Point(5) = {-5,-5,0,1};<br>
Point(6) = {5,-5,0,1};<br>Point(7) = {5,5,0,1};<br>Point(8) = {-5,5,0,1};<br>Line(5) = {5,6};<br>Line(6) = {6,7};<br>Line(7) = {7,8};<br>Line(8) = {8,5};<br>Line(9) = {1,5};<br>Line(10) = {2,6};<br>Line(11) = {3,7};<br>Line(12) = {4,8};<br>
Line Loop(13) = {1,10,-5,-9};<br>Plane Surface(14) = {13};<br>Line Loop(15) = {2,11,-6,-10};<br>Plane Surface(16) = {15};<br>Line Loop(17) = {3,12,-7,-11};<br>Plane Surface(18) = {17};<br>Line Loop(19) = {4,9,-8,-12};<br>
Plane Surface(20) = {19};<br>Point(101) = {-1,-1,10,0.2};<br>Point(102) = {1,-1,10,0.2};<br>Point(103) = {1,1,10,0.2};<br>Point(104) = {-1,1,10,0.2};<br>Line(101) = {101,102};<br>Line(102) = {102,103};<br>Line(103) = {103,104};<br>
Line(104) = {104,101};<br>Line Loop(105) = {104,101,102,103};<br>Plane Surface(106) = {105};<br>Point(105) = {-5,-5,10,1};<br>Point(106) = {5,-5,10,1};<br>Point(107) = {5,5,10,1};<br>Point(108) = {-5,5,10,1};<br>Line(105) = {105,106};<br>
Line(106) = {106,107};<br>Line(107) = {107,108};<br>Line(108) = {108,105};<br>Line(109) = {101,105};<br>Line(110) = {102,106};<br>Line(111) = {103,107};<br>Line(112) = {104,108};<br>Line Loop(113) = {101,110,-105,-109};<br>
Plane Surface(114) = {113};<br>Line Loop(115) = {102,111,-106,-110};<br>Plane Surface(116) = {115};<br>Line Loop(117) = {103,112,-107,-111};<br>Plane Surface(118) = {117};<br>Line Loop(119) = {104,109,-108,-112};<br>Plane Surface(120) = {119};<br>
Line(200) = {1,101};<br>Line(201) = {2,102};<br>Line(202) = {3,103};<br>Line(203) = {4,104};<br>Line(204) = {5,105};<br>Line(205) = {6,106};<br>Line(206) = {7,107};<br>Line(207) = {8,108};<br>Line Loop(208) = {1,201,-101,-200};<br>
Plane Surface(209) = {208};<br>Line Loop(210) = {2,202,-102,-201};<br>Plane Surface(211) = {210};<br>Line Loop(212) = {3,203,-103,-202};<br>Plane Surface(213) = {212};<br>Line Loop(214) = {4,200,-104,-203};<br>Plane Surface(215) = {214};<br>
Line Loop(216) = {9,204,-109,-200};<br>Plane Surface(217) = {216};<br>Line Loop(218) = {10,205,-110,-201};<br>Plane Surface(219) = {218};<br>Line Loop(220) = {11,206,-111,-202};<br>Plane Surface(221) = {220};<br>Line Loop(222) = {12,207,-112,-203};<br>
Plane Surface(223) = {222};<br>Line Loop(224) = {5,205,-105,-204};<br>Plane Surface(225) = {224};<br>Line Loop(226) = {6,206,-106,-205};<br>Plane Surface(227) = {226};<br>Line Loop(228) = {7,207,-107,-206};<br>Plane Surface(229) = {228};<br>
Line Loop(230) = {8,204,-108,-207};<br>Plane Surface(231) = {230};<br>Surface Loop(232) = {209,211,213,215,6,106};<br>Volume(233) = {232};<br>Surface Loop(234) = {14,225,114,217,219,209};<br>Volume(235) = {234};<br>Surface Loop(236) = {227,16,116,219,221,211};<br>
Volume(237) = {236};<br>Surface Loop(238) = {18,229,118,223,221,213};<br>Volume(239) = {238};<br>Surface Loop(240) = {20,231,120,217,223,215};<br>Volume(241) = {240};<br><br>Transfinite Line{1,2,3,4}=5;<br>Transfinite Line{5,6,7,8}=5;<br>
Transfinite Line{101,102,103,104}=5;<br>Transfinite Line{9,10,11,12}=5 Using Progression 1.2;<br>Transfinite Line{105,106,107,108}=5;<br>Transfinite Line{109,110,111,112}=5 Using Progression 1.2;<br>Transfinite Line{200:207}=5;<br>
<br>Transfinite Surface{6}={1,2,3,4};<br>Transfinite Surface{14}={1,2,6,5};<br>Transfinite Surface{16}={2,3,7,6};<br>Transfinite Surface{18}={3,4,8,7};<br>Transfinite Surface{20}={4,1,5,8};<br>Transfinite Surface{106}={101,102,103,104};<br>
Transfinite Surface{114}={101,102,106,105};<br>Transfinite Surface{116}={102,103,107,106};<br>Transfinite Surface{118}={103,104,108,107};<br>Transfinite Surface{120}={104,101,105,108};<br>Transfinite Surface{209}={1,2,102,101};<br>
Transfinite Surface{211}={2,3,103,102};<br>Transfinite Surface{213}={3,4,104,103};<br>Transfinite Surface{215}={4,1,101,104};<br>Transfinite Surface{217}={1,5,105,101};<br>Transfinite Surface{219}={2,6,106,102};<br>Transfinite Surface{221}={3,7,107,103};<br>
Transfinite Surface{223}={4,8,108,104};<br>Transfinite Surface{225}={5,6,106,105};<br>Transfinite Surface{227}={6,7,107,106};<br>Transfinite Surface{229}={7,8,108,107};<br>Transfinite Surface{231}={8,5,105,108};<br><br>Recombine Surface{6};<br>
Recombine Surface{14};<br>Recombine Surface{16};<br>Recombine Surface{18};<br>Recombine Surface{20};<br>Recombine Surface{106};<br>Recombine Surface{114};<br>Recombine Surface{116};<br>Recombine Surface{118};<br>Recombine Surface{120};<br>
Recombine Surface{209};<br>Recombine Surface{211};<br>Recombine Surface{213};<br>Recombine Surface{215};<br>Recombine Surface{217};<br>Recombine Surface{219};<br>Recombine Surface{221};<br>Recombine Surface{223};<br>Recombine Surface{225};<br>
Recombine Surface{227};<br>Recombine Surface{229};Info    : -------------------------------------------------------<br>Info    : Gmsh version   : 2.0.7<br>Info    : Build OS       : Linux<br>Info    : Build options  : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB MATHEVAL ANN OCC <br>
Info    : Build date     : Mon Apr  2 11:09:24 CEST 2007<br>Info    : Build host     : debian<br>Info    : Packager       : geuzaine<br>Info    : Home directory : /home/amit/<br>Info    : Launch date    : Tue Nov  4 15:50:00 2008<br>
Info    : Command line   : gmsh test.geo <br>Info    : -------------------------------------------------------<br>Info    : Reading 'test.geo'<br>Info    : Read 'test.geo'<br>Info    : Meshing 1D...<br>Info    : Meshing curve 1<br>
Info    : Meshing curve 2<br>Info    : Meshing curve 3<br>Info    : Meshing curve 4<br>Info    : Meshing curve 5<br>Info    : Meshing curve 6<br>Info    : Meshing curve 7<br>Info    : Meshing curve 8<br>Info    : Meshing curve 9<br>
Info    : Meshing curve 10<br>Info    : Meshing curve 11<br>Info    : Meshing curve 12<br>Info    : Meshing curve 101<br>Info    : Meshing curve 102<br>Info    : Meshing curve 103<br>Info    : Meshing curve 104<br>Info    : Meshing curve 105<br>
Info    : Meshing curve 106<br>Info    : Meshing curve 107<br>Info    : Meshing curve 108<br>Info    : Meshing curve 109<br>Info    : Meshing curve 110<br>Info    : Meshing curve 111<br>Info    : Meshing curve 112<br>Info    : Meshing curve 200<br>
Info    : Meshing curve 201<br>Info    : Meshing curve 202<br>Info    : Meshing curve 203<br>Info    : Meshing curve 204<br>Info    : Meshing curve 205<br>Info    : Meshing curve 206<br>Info    : Meshing curve 207<br>Info    : Mesh 1D complete (0.008001 s)<br>
Info    : Mesh<br>Info    : Meshing 2D...<br>Info    : Meshing surface 6 (Plane)<br>Info    : Meshing surface 14 (Plane)<br>Info    : Meshing surface 16 (Plane)<br>Info    : Meshing surface 18 (Plane)<br>Info    : Meshing surface 20 (Plane)<br>
Info    : Meshing surface 106 (Plane)<br>Info    : Meshing surface 114 (Plane)<br>Info    : Meshing surface 116 (Plane)<br>Info    : Meshing surface 118 (Plane)<br>Info    : Meshing surface 120 (Plane)<br>Info    : Meshing surface 209 (Plane)<br>
Info    : Meshing surface 211 (Plane)<br>Info    : Meshing surface 213 (Plane)<br>Info    : Meshing surface 215 (Plane)<br>Info    : Meshing surface 217 (Plane)<br>Info    : Meshing surface 219 (Plane)<br>Info    : Meshing surface 221 (Plane)<br>
Info    : Meshing surface 223 (Plane)<br>Info    : Meshing surface 225 (Plane)<br>Info    : Meshing surface 227 (Plane)<br>Info    : Meshing surface 229 (Plane)<br>Info    : Meshing surface 231 (Plane)<br>Info    : Mesh 2D complete (0 s)<br>
Info    : Mesh<br>Info    : Meshing 3D...<br>Info    : Meshing volume 233<br>Info    : Meshing volume 235<br>Error   : Incompatible surface 14 in transfinite volume 235<br>Info    : Meshing volume 237<br>Error   : Incompatible surface 16 in transfinite volume 237<br>
Info    : Meshing volume 239<br>Error   : Incompatible surface 18 in transfinite volume 239<br>Info    : Meshing volume 241<br>Error   : Incompatible surface 20 in transfinite volume 241<br>Info    : Mesh 3D complete (0 s)<br>
Info    : Mesh<br>Info    : 337 vertices 544 elements<br>Info    : Writing '/home/amit/programs/FEM/myExamples/DTR/error.txt'<br>Recombine Surface{231};<br><br>Transfinite Volume{233}={1,2,3,4,101,102,103,104};<br>
Transfinite Volume{235}={1,2,5,6,101,102,105,106};<br>Transfinite Volume{237}={2,3,6,7,102,103,106,107};<br>Transfinite Volume{239}={3,4,7,8,103,104,107,108};<br>Transfinite Volume{241}={4,1,8,5,104,101,108,105};<br><br><br>
I am using an older version of gmsh. Here is the error log.<br>Info    : -------------------------------------------------------<br>Info    : Gmsh version   : 2.0.7<br>Info    : Build OS       : Linux<br>Info    : Build options  : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB MATHEVAL ANN OCC <br>
Info    : Build date     : Mon Apr  2 11:09:24 CEST 2007<br>Info    : Build host     : debian<br>Info    : Packager       : geuzaine<br>Info    : Home directory : /home/amit/<br>Info    : Launch date    : Tue Nov  4 15:50:00 2008<br>
Info    : Command line   : gmsh test.geo <br>Info    : -------------------------------------------------------<br>Info    : Reading 'test.geo'<br>Info    : Read 'test.geo'<br>Info    : Meshing 1D...<br>Info    : Meshing curve 1<br>
Info    : Meshing curve 2<br>Info    : Meshing curve 3<br>Info    : Meshing curve 4<br>Info    : Meshing curve 5<br>Info    : Meshing curve 6<br>Info    : Meshing curve 7<br>Info    : Meshing curve 8<br>Info    : Meshing curve 9<br>
Info    : Meshing curve 10<br>Info    : Meshing curve 11<br>Info    : Meshing curve 12<br>Info    : Meshing curve 101<br>Info    : Meshing curve 102<br>Info    : Meshing curve 103<br>Info    : Meshing curve 104<br>Info    : Meshing curve 105<br>
Info    : Meshing curve 106<br>Info    : Meshing curve 107<br>Info    : Meshing curve 108<br>Info    : Meshing curve 109<br>Info    : Meshing curve 110<br>Info    : Meshing curve 111<br>Info    : Meshing curve 112<br>Info    : Meshing curve 200<br>
Info    : Meshing curve 201<br>Info    : Meshing curve 202<br>Info    : Meshing curve 203<br>Info    : Meshing curve 204<br>Info    : Meshing curve 205<br>Info    : Meshing curve 206Info    : -------------------------------------------------------<br>
Info    : Gmsh version   : 2.0.7Info    : -------------------------------------------------------<br>Info    : Gmsh version   : 2.0.7<br>Info    : Build OS       : Linux<br>Info    : Build options  : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB MATHEVAL ANN OCC <br>
Info    : Build date     : Mon Apr  2 11:09:24 CEST 2007<br>Info    : Build host     : debian<br>Info    : Packager       : geuzaine<br>Info    : Home directory : /home/amit/<br>Info    : Launch date    : Tue Nov  4 15:50:00 2008<br>
Info    : Command line   : gmsh test.geo <br>Info    : -------------------------------------------------------<br>Info    : Reading 'test.geo'<br>Info    : Read 'test.geo'<br>Info    : Meshing 1D...<br>Info    : Meshing curve 1<br>
Info    : Meshing curve 2<br>Info    : Meshing curve 3<br>Info    : Meshing curve 4<br>Info    : Meshing curve 5<br>Info    : Meshing curve 6<br>Info    : Meshing curve 7<br>Info    : Meshing curve 8<br>Info    : Meshing curve 9<br>
Info    : Meshing curve 10<br>Info    : Meshing curve 11<br>Info    : Meshing curve 12<br>Info    : Meshing curve 101<br>Info    : Meshing curve 102<br>Info    : Meshing curve 103<br>Info    : Meshing curve 104<br>Info    : Meshing curve 105<br>
Info    : Meshing curve 106<br>Info    : Meshing curve 107<br>Info    : Meshing curve 108<br>Info    : Meshing curve 109<br>Info    : Meshing curve 110<br>Info    : Meshing curve 111<br>Info    : Meshing curve 112<br>Info    : Meshing curve 200<br>
Info    : Meshing curve 201<br>Info    : Meshing curve 202<br>Info    : Meshing curve 203<br>Info    : Meshing curve 204<br>Info    : Meshing curve 205<br>Info    : Meshing curve 206<br>Info    : Meshing curve 207<br>Info    : Mesh 1D complete (0.008001 s)<br>
Info    : Mesh<br>Info    : Meshing 2D...<br>Info    : Meshing surface 6 (Plane)<br>Info    : Meshing surface 14 (Plane)<br>Info    : Meshing surface 16 (Plane)<br>Info    : Meshing surface 18 (Plane)<br>Info    : Meshing surface 20 (Plane)<br>
Info    : Meshing surface 106 (Plane)<br>Info    : Meshing surface 114 (Plane)<br>Info    : Meshing surface 116 (Plane)<br>Info    : Meshing surface 118 (Plane)<br>Info    : Meshing surface 120 (Plane)<br>Info    : Meshing surface 209 (Plane)<br>
Info    : Meshing surface 211 (Plane)<br>Info    : Meshing surface 213 (Plane)<br>Info    : Meshing surface 215 (Plane)<br>Info    : Meshing surface 217 (Plane)<br>Info    : Meshing surface 219 (Plane)<br>Info    : Meshing surface 221 (Plane)<br>
Info    : Meshing surface 223 (Plane)<br>Info    : Meshing surface 225 (Plane)<br>Info    : Meshing surface 227 (Plane)<br>Info    : Meshing surface 229 (Plane)<br>Info    : Meshing surface 231 (Plane)<br>Info    : Mesh 2D complete (0 s)<br>
Info    : Mesh<br>Info    : Meshing 3D...<br>Info    : Meshing volume 233<br>Info    : Meshing volume 235<br>Error   : Incompatible surface 14 in transfinite volume 235<br>Info    : Meshing volume 237<br>Error   : Incompatible surface 16 in transfinite volume 237<br>
Info    : Meshing volume 239<br>Error   : Incompatible surface 18 in transfinite volume 239<br>Info    : Meshing volume 241<br>Error   : Incompatible surface 20 in transfinite volume 241<br>Info    : Mesh 3D complete (0 s)<br>
Info    : Mesh<br>Info    : 337 vertices 544 elements<br>Info    : Writing '/home/amit/programs/FEM/myExamples/DTR/error.txt'<br>Info    : Build OS       : Linux<br>Info    : Build options  : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB MATHEVAL ANN OCC <br>
Info    : Build date     : Mon Apr  2 11:09:24 CEST 2007<br>Info    : Build host     : debian<br>Info    : Packager       : geuzaine<br>Info    : Home directory : /home/amit/<br>Info    : Launch date    : Tue Nov  4 15:50:00 2008<br>
Info    : Command line   : gmsh test.geo <br>Info    : -------------------------------------------------------<br>Info    : Reading 'test.geo'<br>Info    : Read 'test.geo'<br>Info    : Meshing 1D...<br>Info    : Meshing curve 1<br>
Info    : Meshing curve 2<br>Info    : Meshing curve 3<br>Info    : Meshing curve 4<br>Info    : Meshing curve 5<br>Info    : Meshing curve 6<br>Info    : Meshing curve 7<br>Info    : Meshing curve 8<br>Info    : Meshing curve 9<br>
Info    : Meshing curve 10<br>Info    : Meshing curve 11<br>Info    : Meshing curve 12<br>Info    : Meshing curve 101<br>Info    : Meshing curve 102<br>Info    : Meshing curve 103<br>Info    : Meshing curve 104<br>Info    : Meshing curve 105<br>
Info    : Meshing curve 106<br>Info    : Meshing curve 107<br>Info    : Meshing curve 108<br>Info    : Meshing curve 109<br>Info    : Meshing curve 110<br>Info    : Meshing curve 111<br>Info    : Meshing curve 112<br>Info    : Meshing curve 200<br>
Info    : Meshing curve 201<br>Info    : Meshing curve 202<br>Info    : Meshing curve 203<br>Info    : Meshing curve 204<br>Info    : Meshing curve 205<br>Info    : Meshing curve 206<br>Info    : Meshing curve 207<br>Info    : Mesh 1D complete (0.008001 s)<br>
Info    : Mesh<br>Info    : Meshing 2D...<br>Info    : Meshing surface 6 (Plane)<br>Info    : Meshing surface 14 (Plane)<br>Info    : Meshing surface 16 (Plane)<br>Info    : Meshing surface 18 (Plane)<br>Info    : Meshing surface 20 (Plane)<br>
Info    : Meshing surface 106 (Plane)<br>Info    : Meshing surface 114 (Plane)<br>Info    : Meshing surface 116 (Plane)<br>Info    : Meshing surface 118 (Plane)<br>Info    : Meshing surface 120 (Plane)<br>Info    : Meshing surface 209 (Plane)<br>
Info    : Meshing surface 211 (Plane)<br>Info    : Meshing surface 213 (Plane)<br>Info    : Meshing surface 215 (Plane)<br>Info    : Meshing surface 217 (Plane)<br>Info    : Meshing surface 219 (Plane)<br>Info    : Meshing surface 221 (Plane)<br>
Info    : Meshing surface 223 (Plane)<br>Info    : Meshing surface 225 (Plane)<br>Info    : Meshing surface 227 (Plane)<br>Info    : Meshing surface 229 (Plane)<br>Info    : Meshing surface 231 (Plane)<br>Info    : Mesh 2D complete (0 s)<br>
Info    : Mesh<br>Info    : Meshing 3D...<br>Info    : Meshing volume 233<br>Info    : Meshing volume 235<br>Error   : Incompatible surface 14 in transfinite volume 235<br>Info    : Meshing volume 237<br>Error   : Incompatible surface 16 in transfinite volume 237<br>
Info    : Meshing volume 239<br>Error   : Incompatible surface 18 in transfinite volume 239<br>Info    : Meshing volume 241<br>Error   : Incompatible surface 20 in transfinite volume 241<br>Info    : Mesh 3D complete (0 s)<br>
Info    : Mesh<br>Info    : 337 vertices 544 elements<br>Info    : Writing '/home/amit/programs/FEM/myExamples/DTR/error.txt'<br>Info    : Meshing curve 207<br>Info    : Mesh 1D complete (0.008001 s)<br>Info    : Mesh<br>
Info    : Meshing 2D...<br>Info    : Meshing surface 6 (Plane)<br>Info    : Meshing surface 14 (Plane)<br>Info    : Meshing surface 16 (Plane)<br>Info    : Meshing surface 18 (Plane)<br>Info    : Meshing surface 20 (Plane)<br>
Info    : Meshing surface 106 (Plane)<br>Info    : Meshing surface 114 (Plane)<br>Info    : Meshing surface 116 (Plane)<br>Info    : Meshing surface 118 (Plane)<br>Info    : Meshing surface 120 (Plane)<br>Info    : Meshing surface 209 (Plane)<br>
Info    : Meshing surface 211 (Plane)<br>Info    : Meshing surface 213 (Plane)<br>Info    : Meshing surface 215 (Plane)<br>Info    : Meshing surface 217 (Plane)<br>Info    : Meshing surface 219 (Plane)<br>Info    : Meshing surface 221 (Plane)<br>
Info    : Meshing surface 223 (Plane)<br>Info    : Meshing surface 225 (Plane)<br>Info    : Meshing surface 227 (Plane)<br>Info    : Meshing surface 229 (Plane)<br>Info    : Meshing surface 231 (Plane)<br>Info    : Mesh 2D complete (0 s)<br>
Info    : Mesh<br>Info    : Meshing 3D...<br>Info    : Meshing volume 233<br>Info    : Meshing volume 235<br>Error   : Incompatible surface 14 in transfinite volume 235<br>Info    : Meshing volume 237<br>Error   : Incompatible surface 16 in transfinite volume 237<br>
Info    : Meshing volume 239<br>Error   : Incompatible surface 18 in transfinite volume 239<br>Info    : Meshing volume 241<br>Error   : Incompatible surface 20 in transfinite volume 241<br>Info    : Mesh 3D complete (0 s)<br>
Info    : Mesh<br>Info    : 337 vertices 544 elements<br><br><br><br>Thanks<br><br>Rgds,<br>Amit<br>