[Gmsh] Error transfinite volume
David Colignon
David.Colignon at ulg.ac.be
Wed Nov 5 09:11:01 CET 2008
Hi Amit,
the ordering of the nodes in your "Transfinite Volume" commands is not correct. You should have:
Transfinite Volume{233}={1,2,3,4,101,102,103,104};
Transfinite Volume{235}={1,2,6,5,101,102,106,105};
Transfinite Volume{237}={2,3,7,6,102,103,107,106};
Transfinite Volume{239}={3,4,8,7,103,104,108,107};
Transfinite Volume{241}={4,1,5,8,104,101,105,108};
Cheers,
Dave
--
David Colignon, Ph.D.
Collaborateur Logistique du F.R.S.-FNRS
CÉCI - Consortium des Équipements de Calcul Intensif
ACE - Applied & Computational Electromagnetics
Sart-Tilman B28
Université de Liège
4000 Liège - BELGIQUE
Tél: +32 (0)4 366 37 32
Fax: +32 (0)4 366 29 10
WWW: http://hpc.montefiore.ulg.ac.be/
Agenda: http://www.google.com/calendar/embed?src=david.colignon%40gmail.com
Amit Itagi wrote:
> Hi,
>
> 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 ?
>
> Here is my "geo" file.
>
> Point(1) = {-1,-1,0,0.2};
> Point(2) = {1,-1,0,0.2};
> Point(3) = {1,1,0,0.2};
> Point(4) = {-1,1,0,0.2};
> Line(1) = {1,2};
> Line(2) = {2,3};
> Line(3) = {3,4};
> Line(4) = {4,1};
> Line Loop(5) = {4,1,2,3};
> Plane Surface(6) = {5};
> Point(5) = {-5,-5,0,1};
> Point(6) = {5,-5,0,1};
> Point(7) = {5,5,0,1};
> Point(8) = {-5,5,0,1};
> Line(5) = {5,6};
> Line(6) = {6,7};
> Line(7) = {7,8};
> Line(8) = {8,5};
> Line(9) = {1,5};
> Line(10) = {2,6};
> Line(11) = {3,7};
> Line(12) = {4,8};
> Line Loop(13) = {1,10,-5,-9};
> Plane Surface(14) = {13};
> Line Loop(15) = {2,11,-6,-10};
> Plane Surface(16) = {15};
> Line Loop(17) = {3,12,-7,-11};
> Plane Surface(18) = {17};
> Line Loop(19) = {4,9,-8,-12};
> Plane Surface(20) = {19};
> Point(101) = {-1,-1,10,0.2};
> Point(102) = {1,-1,10,0.2};
> Point(103) = {1,1,10,0.2};
> Point(104) = {-1,1,10,0.2};
> Line(101) = {101,102};
> Line(102) = {102,103};
> Line(103) = {103,104};
> Line(104) = {104,101};
> Line Loop(105) = {104,101,102,103};
> Plane Surface(106) = {105};
> Point(105) = {-5,-5,10,1};
> Point(106) = {5,-5,10,1};
> Point(107) = {5,5,10,1};
> Point(108) = {-5,5,10,1};
> Line(105) = {105,106};
> Line(106) = {106,107};
> Line(107) = {107,108};
> Line(108) = {108,105};
> Line(109) = {101,105};
> Line(110) = {102,106};
> Line(111) = {103,107};
> Line(112) = {104,108};
> Line Loop(113) = {101,110,-105,-109};
> Plane Surface(114) = {113};
> Line Loop(115) = {102,111,-106,-110};
> Plane Surface(116) = {115};
> Line Loop(117) = {103,112,-107,-111};
> Plane Surface(118) = {117};
> Line Loop(119) = {104,109,-108,-112};
> Plane Surface(120) = {119};
> Line(200) = {1,101};
> Line(201) = {2,102};
> Line(202) = {3,103};
> Line(203) = {4,104};
> Line(204) = {5,105};
> Line(205) = {6,106};
> Line(206) = {7,107};
> Line(207) = {8,108};
> Line Loop(208) = {1,201,-101,-200};
> Plane Surface(209) = {208};
> Line Loop(210) = {2,202,-102,-201};
> Plane Surface(211) = {210};
> Line Loop(212) = {3,203,-103,-202};
> Plane Surface(213) = {212};
> Line Loop(214) = {4,200,-104,-203};
> Plane Surface(215) = {214};
> Line Loop(216) = {9,204,-109,-200};
> Plane Surface(217) = {216};
> Line Loop(218) = {10,205,-110,-201};
> Plane Surface(219) = {218};
> Line Loop(220) = {11,206,-111,-202};
> Plane Surface(221) = {220};
> Line Loop(222) = {12,207,-112,-203};
> Plane Surface(223) = {222};
> Line Loop(224) = {5,205,-105,-204};
> Plane Surface(225) = {224};
> Line Loop(226) = {6,206,-106,-205};
> Plane Surface(227) = {226};
> Line Loop(228) = {7,207,-107,-206};
> Plane Surface(229) = {228};
> Line Loop(230) = {8,204,-108,-207};
> Plane Surface(231) = {230};
> Surface Loop(232) = {209,211,213,215,6,106};
> Volume(233) = {232};
> Surface Loop(234) = {14,225,114,217,219,209};
> Volume(235) = {234};
> Surface Loop(236) = {227,16,116,219,221,211};
> Volume(237) = {236};
> Surface Loop(238) = {18,229,118,223,221,213};
> Volume(239) = {238};
> Surface Loop(240) = {20,231,120,217,223,215};
> Volume(241) = {240};
>
> Transfinite Line{1,2,3,4}=5;
> Transfinite Line{5,6,7,8}=5;
> Transfinite Line{101,102,103,104}=5;
> Transfinite Line{9,10,11,12}=5 Using Progression 1.2;
> Transfinite Line{105,106,107,108}=5;
> Transfinite Line{109,110,111,112}=5 Using Progression 1.2;
> Transfinite Line{200:207}=5;
>
> Transfinite Surface{6}={1,2,3,4};
> Transfinite Surface{14}={1,2,6,5};
> Transfinite Surface{16}={2,3,7,6};
> Transfinite Surface{18}={3,4,8,7};
> Transfinite Surface{20}={4,1,5,8};
> Transfinite Surface{106}={101,102,103,104};
> Transfinite Surface{114}={101,102,106,105};
> Transfinite Surface{116}={102,103,107,106};
> Transfinite Surface{118}={103,104,108,107};
> Transfinite Surface{120}={104,101,105,108};
> Transfinite Surface{209}={1,2,102,101};
> Transfinite Surface{211}={2,3,103,102};
> Transfinite Surface{213}={3,4,104,103};
> Transfinite Surface{215}={4,1,101,104};
> Transfinite Surface{217}={1,5,105,101};
> Transfinite Surface{219}={2,6,106,102};
> Transfinite Surface{221}={3,7,107,103};
> Transfinite Surface{223}={4,8,108,104};
> Transfinite Surface{225}={5,6,106,105};
> Transfinite Surface{227}={6,7,107,106};
> Transfinite Surface{229}={7,8,108,107};
> Transfinite Surface{231}={8,5,105,108};
>
> Recombine Surface{6};
> Recombine Surface{14};
> Recombine Surface{16};
> Recombine Surface{18};
> Recombine Surface{20};
> Recombine Surface{106};
> Recombine Surface{114};
> Recombine Surface{116};
> Recombine Surface{118};
> Recombine Surface{120};
> Recombine Surface{209};
> Recombine Surface{211};
> Recombine Surface{213};
> Recombine Surface{215};
> Recombine Surface{217};
> Recombine Surface{219};
> Recombine Surface{221};
> Recombine Surface{223};
> Recombine Surface{225};
> Recombine Surface{227};
> Recombine Surface{229};Info :
> -------------------------------------------------------
> Info : Gmsh version : 2.0.7
> Info : Build OS : Linux
> Info : Build options : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB
> MATHEVAL ANN OCC
> Info : Build date : Mon Apr 2 11:09:24 CEST 2007
> Info : Build host : debian
> Info : Packager : geuzaine
> Info : Home directory : /home/amit/
> Info : Launch date : Tue Nov 4 15:50:00 2008
> Info : Command line : gmsh test.geo
> Info : -------------------------------------------------------
> Info : Reading 'test.geo'
> Info : Read 'test.geo'
> Info : Meshing 1D...
> Info : Meshing curve 1
> Info : Meshing curve 2
> Info : Meshing curve 3
> Info : Meshing curve 4
> Info : Meshing curve 5
> Info : Meshing curve 6
> Info : Meshing curve 7
> Info : Meshing curve 8
> Info : Meshing curve 9
> Info : Meshing curve 10
> Info : Meshing curve 11
> Info : Meshing curve 12
> Info : Meshing curve 101
> Info : Meshing curve 102
> Info : Meshing curve 103
> Info : Meshing curve 104
> Info : Meshing curve 105
> Info : Meshing curve 106
> Info : Meshing curve 107
> Info : Meshing curve 108
> Info : Meshing curve 109
> Info : Meshing curve 110
> Info : Meshing curve 111
> Info : Meshing curve 112
> Info : Meshing curve 200
> Info : Meshing curve 201
> Info : Meshing curve 202
> Info : Meshing curve 203
> Info : Meshing curve 204
> Info : Meshing curve 205
> Info : Meshing curve 206
> Info : Meshing curve 207
> Info : Mesh 1D complete (0.008001 s)
> Info : Mesh
> Info : Meshing 2D...
> Info : Meshing surface 6 (Plane)
> Info : Meshing surface 14 (Plane)
> Info : Meshing surface 16 (Plane)
> Info : Meshing surface 18 (Plane)
> Info : Meshing surface 20 (Plane)
> Info : Meshing surface 106 (Plane)
> Info : Meshing surface 114 (Plane)
> Info : Meshing surface 116 (Plane)
> Info : Meshing surface 118 (Plane)
> Info : Meshing surface 120 (Plane)
> Info : Meshing surface 209 (Plane)
> Info : Meshing surface 211 (Plane)
> Info : Meshing surface 213 (Plane)
> Info : Meshing surface 215 (Plane)
> Info : Meshing surface 217 (Plane)
> Info : Meshing surface 219 (Plane)
> Info : Meshing surface 221 (Plane)
> Info : Meshing surface 223 (Plane)
> Info : Meshing surface 225 (Plane)
> Info : Meshing surface 227 (Plane)
> Info : Meshing surface 229 (Plane)
> Info : Meshing surface 231 (Plane)
> Info : Mesh 2D complete (0 s)
> Info : Mesh
> Info : Meshing 3D...
> Info : Meshing volume 233
> Info : Meshing volume 235
> Error : Incompatible surface 14 in transfinite volume 235
> Info : Meshing volume 237
> Error : Incompatible surface 16 in transfinite volume 237
> Info : Meshing volume 239
> Error : Incompatible surface 18 in transfinite volume 239
> Info : Meshing volume 241
> Error : Incompatible surface 20 in transfinite volume 241
> Info : Mesh 3D complete (0 s)
> Info : Mesh
> Info : 337 vertices 544 elements
> Info : Writing '/home/amit/programs/FEM/myExamples/DTR/error.txt'
> Recombine Surface{231};
>
> Transfinite Volume{233}={1,2,3,4,101,102,103,104};
> Transfinite Volume{235}={1,2,5,6,101,102,105,106};
> Transfinite Volume{237}={2,3,6,7,102,103,106,107};
> Transfinite Volume{239}={3,4,7,8,103,104,107,108};
> Transfinite Volume{241}={4,1,8,5,104,101,108,105};
>
>
> I am using an older version of gmsh. Here is the error log.
> Info : -------------------------------------------------------
> Info : Gmsh version : 2.0.7
> Info : Build OS : Linux
> Info : Build options : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB
> MATHEVAL ANN OCC
> Info : Build date : Mon Apr 2 11:09:24 CEST 2007
> Info : Build host : debian
> Info : Packager : geuzaine
> Info : Home directory : /home/amit/
> Info : Launch date : Tue Nov 4 15:50:00 2008
> Info : Command line : gmsh test.geo
> Info : -------------------------------------------------------
> Info : Reading 'test.geo'
> Info : Read 'test.geo'
> Info : Meshing 1D...
> Info : Meshing curve 1
> Info : Meshing curve 2
> Info : Meshing curve 3
> Info : Meshing curve 4
> Info : Meshing curve 5
> Info : Meshing curve 6
> Info : Meshing curve 7
> Info : Meshing curve 8
> Info : Meshing curve 9
> Info : Meshing curve 10
> Info : Meshing curve 11
> Info : Meshing curve 12
> Info : Meshing curve 101
> Info : Meshing curve 102
> Info : Meshing curve 103
> Info : Meshing curve 104
> Info : Meshing curve 105
> Info : Meshing curve 106
> Info : Meshing curve 107
> Info : Meshing curve 108
> Info : Meshing curve 109
> Info : Meshing curve 110
> Info : Meshing curve 111
> Info : Meshing curve 112
> Info : Meshing curve 200
> Info : Meshing curve 201
> Info : Meshing curve 202
> Info : Meshing curve 203
> Info : Meshing curve 204
> Info : Meshing curve 205
> Info : Meshing curve 206Info :
> -------------------------------------------------------
> Info : Gmsh version : 2.0.7Info :
> -------------------------------------------------------
> Info : Gmsh version : 2.0.7
> Info : Build OS : Linux
> Info : Build options : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB
> MATHEVAL ANN OCC
> Info : Build date : Mon Apr 2 11:09:24 CEST 2007
> Info : Build host : debian
> Info : Packager : geuzaine
> Info : Home directory : /home/amit/
> Info : Launch date : Tue Nov 4 15:50:00 2008
> Info : Command line : gmsh test.geo
> Info : -------------------------------------------------------
> Info : Reading 'test.geo'
> Info : Read 'test.geo'
> Info : Meshing 1D...
> Info : Meshing curve 1
> Info : Meshing curve 2
> Info : Meshing curve 3
> Info : Meshing curve 4
> Info : Meshing curve 5
> Info : Meshing curve 6
> Info : Meshing curve 7
> Info : Meshing curve 8
> Info : Meshing curve 9
> Info : Meshing curve 10
> Info : Meshing curve 11
> Info : Meshing curve 12
> Info : Meshing curve 101
> Info : Meshing curve 102
> Info : Meshing curve 103
> Info : Meshing curve 104
> Info : Meshing curve 105
> Info : Meshing curve 106
> Info : Meshing curve 107
> Info : Meshing curve 108
> Info : Meshing curve 109
> Info : Meshing curve 110
> Info : Meshing curve 111
> Info : Meshing curve 112
> Info : Meshing curve 200
> Info : Meshing curve 201
> Info : Meshing curve 202
> Info : Meshing curve 203
> Info : Meshing curve 204
> Info : Meshing curve 205
> Info : Meshing curve 206
> Info : Meshing curve 207
> Info : Mesh 1D complete (0.008001 s)
> Info : Mesh
> Info : Meshing 2D...
> Info : Meshing surface 6 (Plane)
> Info : Meshing surface 14 (Plane)
> Info : Meshing surface 16 (Plane)
> Info : Meshing surface 18 (Plane)
> Info : Meshing surface 20 (Plane)
> Info : Meshing surface 106 (Plane)
> Info : Meshing surface 114 (Plane)
> Info : Meshing surface 116 (Plane)
> Info : Meshing surface 118 (Plane)
> Info : Meshing surface 120 (Plane)
> Info : Meshing surface 209 (Plane)
> Info : Meshing surface 211 (Plane)
> Info : Meshing surface 213 (Plane)
> Info : Meshing surface 215 (Plane)
> Info : Meshing surface 217 (Plane)
> Info : Meshing surface 219 (Plane)
> Info : Meshing surface 221 (Plane)
> Info : Meshing surface 223 (Plane)
> Info : Meshing surface 225 (Plane)
> Info : Meshing surface 227 (Plane)
> Info : Meshing surface 229 (Plane)
> Info : Meshing surface 231 (Plane)
> Info : Mesh 2D complete (0 s)
> Info : Mesh
> Info : Meshing 3D...
> Info : Meshing volume 233
> Info : Meshing volume 235
> Error : Incompatible surface 14 in transfinite volume 235
> Info : Meshing volume 237
> Error : Incompatible surface 16 in transfinite volume 237
> Info : Meshing volume 239
> Error : Incompatible surface 18 in transfinite volume 239
> Info : Meshing volume 241
> Error : Incompatible surface 20 in transfinite volume 241
> Info : Mesh 3D complete (0 s)
> Info : Mesh
> Info : 337 vertices 544 elements
> Info : Writing '/home/amit/programs/FEM/myExamples/DTR/error.txt'
> Info : Build OS : Linux
> Info : Build options : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB
> MATHEVAL ANN OCC
> Info : Build date : Mon Apr 2 11:09:24 CEST 2007
> Info : Build host : debian
> Info : Packager : geuzaine
> Info : Home directory : /home/amit/
> Info : Launch date : Tue Nov 4 15:50:00 2008
> Info : Command line : gmsh test.geo
> Info : -------------------------------------------------------
> Info : Reading 'test.geo'
> Info : Read 'test.geo'
> Info : Meshing 1D...
> Info : Meshing curve 1
> Info : Meshing curve 2
> Info : Meshing curve 3
> Info : Meshing curve 4
> Info : Meshing curve 5
> Info : Meshing curve 6
> Info : Meshing curve 7
> Info : Meshing curve 8
> Info : Meshing curve 9
> Info : Meshing curve 10
> Info : Meshing curve 11
> Info : Meshing curve 12
> Info : Meshing curve 101
> Info : Meshing curve 102
> Info : Meshing curve 103
> Info : Meshing curve 104
> Info : Meshing curve 105
> Info : Meshing curve 106
> Info : Meshing curve 107
> Info : Meshing curve 108
> Info : Meshing curve 109
> Info : Meshing curve 110
> Info : Meshing curve 111
> Info : Meshing curve 112
> Info : Meshing curve 200
> Info : Meshing curve 201
> Info : Meshing curve 202
> Info : Meshing curve 203
> Info : Meshing curve 204
> Info : Meshing curve 205
> Info : Meshing curve 206
> Info : Meshing curve 207
> Info : Mesh 1D complete (0.008001 s)
> Info : Mesh
> Info : Meshing 2D...
> Info : Meshing surface 6 (Plane)
> Info : Meshing surface 14 (Plane)
> Info : Meshing surface 16 (Plane)
> Info : Meshing surface 18 (Plane)
> Info : Meshing surface 20 (Plane)
> Info : Meshing surface 106 (Plane)
> Info : Meshing surface 114 (Plane)
> Info : Meshing surface 116 (Plane)
> Info : Meshing surface 118 (Plane)
> Info : Meshing surface 120 (Plane)
> Info : Meshing surface 209 (Plane)
> Info : Meshing surface 211 (Plane)
> Info : Meshing surface 213 (Plane)
> Info : Meshing surface 215 (Plane)
> Info : Meshing surface 217 (Plane)
> Info : Meshing surface 219 (Plane)
> Info : Meshing surface 221 (Plane)
> Info : Meshing surface 223 (Plane)
> Info : Meshing surface 225 (Plane)
> Info : Meshing surface 227 (Plane)
> Info : Meshing surface 229 (Plane)
> Info : Meshing surface 231 (Plane)
> Info : Mesh 2D complete (0 s)
> Info : Mesh
> Info : Meshing 3D...
> Info : Meshing volume 233
> Info : Meshing volume 235
> Error : Incompatible surface 14 in transfinite volume 235
> Info : Meshing volume 237
> Error : Incompatible surface 16 in transfinite volume 237
> Info : Meshing volume 239
> Error : Incompatible surface 18 in transfinite volume 239
> Info : Meshing volume 241
> Error : Incompatible surface 20 in transfinite volume 241
> Info : Mesh 3D complete (0 s)
> Info : Mesh
> Info : 337 vertices 544 elements
> Info : Writing '/home/amit/programs/FEM/myExamples/DTR/error.txt'
> Info : Meshing curve 207
> Info : Mesh 1D complete (0.008001 s)
> Info : Mesh
> Info : Meshing 2D...
> Info : Meshing surface 6 (Plane)
> Info : Meshing surface 14 (Plane)
> Info : Meshing surface 16 (Plane)
> Info : Meshing surface 18 (Plane)
> Info : Meshing surface 20 (Plane)
> Info : Meshing surface 106 (Plane)
> Info : Meshing surface 114 (Plane)
> Info : Meshing surface 116 (Plane)
> Info : Meshing surface 118 (Plane)
> Info : Meshing surface 120 (Plane)
> Info : Meshing surface 209 (Plane)
> Info : Meshing surface 211 (Plane)
> Info : Meshing surface 213 (Plane)
> Info : Meshing surface 215 (Plane)
> Info : Meshing surface 217 (Plane)
> Info : Meshing surface 219 (Plane)
> Info : Meshing surface 221 (Plane)
> Info : Meshing surface 223 (Plane)
> Info : Meshing surface 225 (Plane)
> Info : Meshing surface 227 (Plane)
> Info : Meshing surface 229 (Plane)
> Info : Meshing surface 231 (Plane)
> Info : Mesh 2D complete (0 s)
> Info : Mesh
> Info : Meshing 3D...
> Info : Meshing volume 233
> Info : Meshing volume 235
> Error : Incompatible surface 14 in transfinite volume 235
> Info : Meshing volume 237
> Error : Incompatible surface 16 in transfinite volume 237
> Info : Meshing volume 239
> Error : Incompatible surface 18 in transfinite volume 239
> Info : Meshing volume 241
> Error : Incompatible surface 20 in transfinite volume 241
> Info : Mesh 3D complete (0 s)
> Info : Mesh
> Info : 337 vertices 544 elements
>
>
>
> Thanks
>
> Rgds,
> Amit
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh