[Gmsh] Error transfinite volume
Amit Itagi
aitagi at gmail.com
Wed Nov 5 15:26:11 CET 2008
David,
Thanks for the corrections. This works.
What is the rule for node ordering (in the transfinite volume/surface
statements) ? Is this documented ?
Thanks
Rgds,
Amit
On Wed, Nov 5, 2008 at 3:11 AM, David Colignon <David.Colignon at ulg.ac.be>wrote:
> 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
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20081105/a1cc9eb3/attachment.html>