Hi Christophe,
thanks a lot!
Actually, it seems that I wasted a lot of time to create a hex mesh for
a room airflow setup with a jet just in the middle of one wall.
I attached my approach. It looks like, that I could have done it a lot
easier and with much less transfinite surfaces/volumnes using 'Extrude'.
Though, do you think that it could be a problem to expand the inlet
section just like in my attached approach using extrude? The middle
section should expand, whereas the other part of the room would become
smaller in the x-direction...
Regards!
Fabian
Christophe Geuzaine wrote:
Fabian wrote:
Hi Christophe,
Christophe Geuzaine wrote:
Fabian wrote:
Hi,
I just saw the geometry and wonder, if there is a chance to extrude
the
sure: cf. attached file.
Oh, thanks! Am I able to declare any progression or bump to the extrude
volume as well? Would be great.
You can specify the number of elements and the relative height of each
extruded layer. For example:
Extrude {0, 0, 1} {
Surface{9, 1};
Layers{{1,1,1,1,1,1}, {0.05,0.15,0.5,0.85,0.95,1}};
Recombine;
}
You could thus generate a progression with something like this:
n = 6;
r = 2;
a = (r - 1) / (r^n - 1);
one[0] = 1;
layer[0] = a;
For i In {1:n-1}
one[i] = 1;
layer[i] = layer[i-1] + a * r^i;
EndFor
Extrude {0, 0, 1} { Surface{9, 1}; Layers{one[], layer[]}; Recombine; }
Fabian
surfaces and create a transfinite volume?
It's probably not working due to more than six surfaces building the
volumen!? Would be nice, if this would work, though...
Fabian
Christophe Geuzaine wrote:
Durufle Marc wrote:
Hi,
I have two questions :
1- How do you create hybrid meshes in 2-D (quad-dominant with a few
triangles) ?
Until now, I am using the Recombine option, but often I get
non-convex quadrilaterals.
In the picture attached to this message, you can see that two
elements (in black) are non-convex
quads and that's a problem ... So, often I change the characterisc
length in the hope that I will have only
fine quadrilaterals. Is it possible to require that all
quadrilaterals are convex ?
Yes. In your file you use
Recombine Surface {1,s0[0]} = 100;
which means that you allow quads with angles between 90-100 and
90+100 degrees (which does not make sense).
Use
Recombine Surface {1,s0[0]} = 45;
instead, which will allow angles between 90-45=45 and 90+45=135
degrees. (45 is the default if you don't specify anything, i.e.,
with just "Recombine Surface {1,s0[0]}".)
We are working on an automatic full quad algorithm, but it's not
ready yet. In the meantime for simple geometries you could also use
Transfinite meshes to generate a full quad mesh (see attached file).
2- In 3-D, when I try "Recombine Volume", it doesn't work. So, how
is it possible with gmsh
to create 3-D hexahedral-dominant meshes (with some tets, pyramids
and wedges) ?
Non. Your only choices for generating hex/prism/pyramid meshes are
extruded (see "Extrude Surface") or transfinite (see "Transfinite
Volume") grids.
Thank you for your help.
PS : I am using gmsh 2.24
------------------------------------------------------------------------
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
------------------------------------------------------------------------
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
//========== PARAMETER MESH ===============================//
// Lauflaenge des Strahls (x-Richtung) //
nodesLaenge=10;
progLaenge=1.01;
// Strahl etc. //
nodesStrahl=10;
progStrahl=1;
// Hoehe Konus innen //
nodesKonInnen=5;
progKonInnen=1.05;
// Hoehe Konus aussen //
nodesKonAussen=5;
progKonAussen=1.04;
// Hoehe Raum //
nodesRaum=5;
progRaum=1.01;
// Einlass //
nodesEinlassLaenge=5;
progEinlassLaenge=1;
// Auslass //
nodesAuslassHoehe=5;
progAuslassHoehe=1;
//========== PARAMETER GEOMETRIE ==========================//
// Box //
gesLaenge=3.0;
gesBreite=1.8;
gesHoehe=1.1;
mitteBreite=gesBreite/2;
mitteHoehe=gesHoehe/2;
nullLaenge=0.0;
nullBreite=0.0;
nullBreite_v2=0.2;
nullHoehe=0.0;
// Strahl //
rStr1Klein=0.0185;
rStr1Gross=0.1;
rStr2Klein=0.1;
rStr2Gross=0.3;
bStrBoxKlein=rStr1Klein/3;
hStrBoxKlein=rStr1Klein/3;
bStrBoxGross=rStr1Gross/3;
hStrBoxGross=rStr1Gross/3;
angl = 45. * (Pi/180.) ;
// Einlass //
einlassLaenge=0.04;
// Auslass //
auslassHoehe=0.01;
//========== EXAMPLES =====================================//
//Point(1) = {0, 0, 0, lc0} ;
//Point(4) = {r_far, 0, 0, lc2} ;
//Point(23) = {r_far*Cos(angl), 0, r_far*Sin(angl), lc2} ;
//Point(8) = {0, 0, r_far, lc2} ;
//Point(14) = {r_int/2, 0, 0, lc2} ;
//Point(16) = {r_int/2, 0, r_int/2, lc2} ;
//Point(18) = {0, 0, r_int/2, lc2} ;
//========== DATA GEOMETRIE ==============================//
Point(1) = { nullLaenge , mitteBreite , mitteHoehe ,1e+22};
Point(2) = { nullLaenge , mitteBreite-bStrBoxKlein , mitteHoehe ,1e+22};
Point(3) = { nullLaenge , mitteBreite-rStr1Klein , mitteHoehe ,1e+22};
Point(4) = { nullLaenge , mitteBreite-rStr2Klein , mitteHoehe ,1e+22};
Point(5) = { nullLaenge , nullBreite , mitteHoehe ,1e+22};
Point(6) = { nullLaenge , mitteBreite , mitteHoehe+hStrBoxKlein ,1e+22};
Point(7) = { nullLaenge , mitteBreite , mitteHoehe+rStr1Klein ,1e+22};
Point(8) = { nullLaenge , mitteBreite , mitteHoehe+rStr2Klein ,1e+22};
Point(9) = { nullLaenge , mitteBreite , gesHoehe ,1e+22};
Point(10) = { nullLaenge , mitteBreite-bStrBoxKlein , mitteHoehe+hStrBoxKlein
,1e+22};
Point(11) = { nullLaenge , mitteBreite-rStr1Klein*Cos(angl),
mitteHoehe+rStr1Klein*Sin(angl), 1e+22} ;
Point(12) = { nullLaenge , mitteBreite-rStr2Klein*Cos(angl),
mitteHoehe+rStr2Klein*Sin(angl), 1e+22} ;
Point(13) = { nullLaenge , nullBreite , gesHoehe ,1e+22};
Point(14) = { nullLaenge , mitteBreite , mitteHoehe-hStrBoxKlein ,1e+22};
Point(15) = { nullLaenge , mitteBreite , mitteHoehe-rStr1Klein ,1e+22};
Point(16) = { nullLaenge , mitteBreite , mitteHoehe-rStr2Klein ,1e+22};
Point(17) = { nullLaenge , mitteBreite , nullHoehe+auslassHoehe ,1e+22};
Point(18) = { nullLaenge , mitteBreite-bStrBoxKlein , mitteHoehe-hStrBoxKlein
,1e+22};
Point(19) = { nullLaenge , mitteBreite-rStr1Klein*Cos(angl),
mitteHoehe-rStr1Klein*Sin(angl), 1e+22} ;
Point(20) = { nullLaenge , mitteBreite-rStr2Klein*Cos(angl),
mitteHoehe-rStr2Klein*Sin(angl), 1e+22} ;
Point(21) = { nullLaenge , nullBreite , nullHoehe+auslassHoehe ,1e+22};
Point(22) = { gesLaenge , mitteBreite , mitteHoehe ,1e+22};
Point(23) = { gesLaenge , mitteBreite-bStrBoxGross , mitteHoehe ,1e+22};
Point(24) = { gesLaenge , mitteBreite-rStr1Gross , mitteHoehe ,1e+22};
Point(25) = { gesLaenge , mitteBreite-rStr2Gross , mitteHoehe ,1e+22};
Point(26) = { gesLaenge , nullBreite , mitteHoehe ,1e+22};
Point(27) = { gesLaenge , mitteBreite , mitteHoehe+hStrBoxGross ,1e+22};
Point(28) = { gesLaenge , mitteBreite , mitteHoehe+rStr1Gross ,1e+22};
Point(29) = { gesLaenge , mitteBreite , mitteHoehe+rStr2Gross ,1e+22};
Point(30) = { gesLaenge , mitteBreite , gesHoehe ,1e+22};
Point(31) = { gesLaenge , mitteBreite-bStrBoxGross , mitteHoehe+hStrBoxGross
,1e+22};
Point(32) = { gesLaenge , mitteBreite-rStr1Gross*Cos(angl),
mitteHoehe+rStr1Gross*Sin(angl), 1e+22} ;
Point(33) = { gesLaenge , mitteBreite-rStr2Gross*Cos(angl),
mitteHoehe+rStr2Gross*Sin(angl), 1e+22} ;
Point(34) = { gesLaenge , nullBreite , gesHoehe ,1e+22};
Point(35) = { gesLaenge , mitteBreite , mitteHoehe-hStrBoxGross ,1e+22};
Point(36) = { gesLaenge , mitteBreite , mitteHoehe-rStr1Gross ,1e+22};
Point(37) = { gesLaenge , mitteBreite , mitteHoehe-rStr2Gross ,1e+22};
Point(38) = { gesLaenge , mitteBreite , nullHoehe+auslassHoehe ,1e+22};
Point(39) = { gesLaenge , mitteBreite-bStrBoxGross , mitteHoehe-hStrBoxGross
,1e+22};
Point(40) = { gesLaenge , mitteBreite-rStr1Gross*Cos(angl),
mitteHoehe-rStr1Gross*Sin(angl), 1e+22} ;
Point(41) = { gesLaenge , mitteBreite-rStr2Gross*Cos(angl),
mitteHoehe-rStr2Gross*Sin(angl), 1e+22} ;
Point(42) = { gesLaenge , nullBreite , nullHoehe+auslassHoehe ,1e+22};
Point(43) = { nullLaenge-einlassLaenge , mitteBreite , mitteHoehe ,1e+22};
Point(44) = { nullLaenge-einlassLaenge , mitteBreite-bStrBoxKlein , mitteHoehe
,1e+22};
Point(45) = { nullLaenge-einlassLaenge , mitteBreite-rStr1Klein , mitteHoehe
,1e+22};
Point(46) = { nullLaenge-einlassLaenge , mitteBreite , mitteHoehe+hStrBoxKlein
,1e+22};
Point(47) = { nullLaenge-einlassLaenge , mitteBreite , mitteHoehe+rStr1Klein
,1e+22};
Point(48) = { nullLaenge-einlassLaenge , mitteBreite-bStrBoxKlein ,
mitteHoehe+hStrBoxKlein ,1e+22};
Point(49) = { nullLaenge-einlassLaenge , mitteBreite-rStr1Klein*Cos(angl),
mitteHoehe+rStr1Klein*Sin(angl), 1e+22} ;
Point(50) = { nullLaenge-einlassLaenge , mitteBreite , mitteHoehe-hStrBoxKlein
,1e+22};
Point(51) = { nullLaenge-einlassLaenge , mitteBreite , mitteHoehe-rStr1Klein
,1e+22};
Point(52) = { nullLaenge-einlassLaenge , mitteBreite-bStrBoxKlein ,
mitteHoehe-hStrBoxKlein ,1e+22};
Point(53) = { nullLaenge-einlassLaenge , mitteBreite-rStr1Klein*Cos(angl),
mitteHoehe-rStr1Klein*Sin(angl), 1e+22} ;
Point(54) = { nullLaenge , mitteBreite , nullHoehe ,1e+22};
Point(55) = { nullLaenge , nullBreite , nullHoehe ,1e+22};
Point(56) = { gesLaenge , mitteBreite , nullHoehe ,1e+22};
Point(57) = { gesLaenge , nullBreite , nullHoehe ,1e+22};
//--------------------------------------------------//
Line(1) = {10,6};
Line(2) = {10,2};
Line(3) = {2,18};
Line(4) = {18,14};
Circle(5) = {7,1,11};
Circle(6) = {11,1,3};
Circle(7) = {3,1,19};
Circle(8) = {19,1,15};
Circle(9) = {8,1,12};
Circle(10) = {12,1,4};
Circle(11) = {4,1,20};
Circle(12) = {20,1,16};
Line(13) = {9,13};
Line(14) = {13,5};
Line(15) = {5,21};
Line(16) = {21,17};
Line(17) = {1,6};
Line(18) = {6,7};
Line(19) = {7,8};
Line(20) = {8,9};
Line(21) = {1,14};
Line(22) = {14,15};
Line(23) = {15,16};
Line(24) = {16,17};
Line(25) = {1,2};
Line(26) = {2,3};
Line(27) = {3,4};
Line(28) = {4,5};
Line(29) = {10,11};
Line(30) = {11,12};
Line(31) = {12,13};
Line(32) = {18,19};
Line(33) = {19,20};
Line(34) = {20,21};
Line(35) = {27,31};
Line(36) = {31,23};
Line(37) = {23,39};
Line(38) = {39,35};
Circle(39) = {28,22,32};
Circle(40) = {32,22,24};
Circle(41) = {24,22,40};
Circle(42) = {40,22,36};
Circle(43) = {29,22,33};
Circle(44) = {33,22,25};
Circle(45) = {25,22,41};
Circle(46) = {41,22,37};
Line(47) = {30,34};
Line(48) = {34,26};
Line(49) = {26,42};
Line(50) = {42,38};
Line(51) = {22,27};
Line(52) = {27,28};
Line(53) = {28,29};
Line(54) = {29,30};
Line(55) = {22,35};
Line(56) = {35,36};
Line(57) = {36,37};
Line(58) = {37,38};
Line(59) = {22,23};
Line(60) = {23,24};
Line(61) = {24,25};
Line(62) = {25,26};
Line(63) = {31,32};
Line(64) = {32,33};
Line(65) = {33,34};
Line(66) = {39,40};
Line(67) = {40,41};
Line(68) = {41,42};
Line(69) = {1,22};
Line(70) = {6,27};
Line(71) = {7,28};
Line(72) = {8,29};
Line(73) = {9,30};
Line(74) = {14,35};
Line(75) = {15,36};
Line(76) = {16,37};
Line(77) = {17,38};
Line(78) = {2,23};
Line(79) = {10,31};
Line(80) = {18,39};
Line(81) = {3,24};
Line(82) = {11,32};
Line(83) = {19,40};
Line(84) = {4,25};
Line(85) = {12,33};
Line(86) = {20,41};
Line(87) = {5,26};
Line(88) = {13,34};
Line(89) = {21,42};
Line(90) = {46,48};
Line(91) = {48,44};
Line(92) = {44,52};
Line(93) = {52,50};
Circle(94) = {47,43,49};
Circle(95) = {49,43,45};
Circle(96) = {45,43,53};
Circle(97) = {53,43,51};
Line(98) = {43,44};
Line(99) = {44,45};
Line(100) = {43,46};
Line(101) = {46,47};
Line(102) = {43,50};
Line(103) = {50,51};
Line(104) = {48,49};
Line(105) = {52,53};
Line(106) = {43,1};
Line(107) = {46,6};
Line(108) = {47,7};
Line(109) = {50,14};
Line(110) = {51,15};
Line(111) = {48,10};
Line(112) = {44,2};
Line(113) = {52,18};
Line(114) = {45,3};
Line(115) = {49,11};
Line(116) = {53,19};
Line(117) = {21,55};
Line(118) = {42,57};
Line(119) = {38,56};
Line(120) = {17,54};
Line(121) = {54,55};
Line(122) = {55,57};
Line(123) = {57,56};
Line(124) = {56,54};
//--------------------------------------------------//
Line Loop(125) = {17,-1,2,-25};
Plane Surface(126) = {125};
Line Loop(127) = {21,-4,-3,-25};
Plane Surface(128) = {127};
Line Loop(129) = {18,5,-29,1};
Plane Surface(130) = {129};
Line Loop(131) = {29,6,-26,-2};
Plane Surface(132) = {131};
Line Loop(133) = {26,7,-32,-3};
Plane Surface(134) = {133};
Line Loop(135) = {32,8,-22,-4};
Plane Surface(136) = {135};
Line Loop(137) = {19,9,-30,-5};
Plane Surface(138) = {137};
Line Loop(139) = {6,27,-10,-30};
Plane Surface(140) = {139};
Line Loop(141) = {27,11,-33,-7};
Plane Surface(142) = {141};
Line Loop(143) = {33,12,-23,-8};
Plane Surface(144) = {143};
Line Loop(145) = {20,13,-31,-9};
Plane Surface(146) = {145};
Line Loop(147) = {31,14,-28,-10};
Plane Surface(148) = {147};
Line Loop(149) = {28,15,-34,-11};
Plane Surface(150) = {149};
Line Loop(151) = {34,16,-24,-12};
Plane Surface(152) = {151};
Line Loop(153) = {16,120,121,-117};
Plane Surface(154) = {153};
Line Loop(155) = {59,-36,-35,-51};
Plane Surface(156) = {155};
Line Loop(157) = {55,-38,-37,-59};
Plane Surface(158) = {157};
Line Loop(159) = {35,63,-39,-52};
Plane Surface(160) = {159};
Line Loop(161) = {36,60,-40,-63};
Plane Surface(162) = {161};
Line Loop(163) = {37,66,-41,-60};
Plane Surface(164) = {163};
Line Loop(165) = {38,56,-42,-66};
Plane Surface(166) = {165};
Line Loop(167) = {39,64,-43,-53};
Plane Surface(168) = {167};
Line Loop(169) = {40,61,-44,-64};
Plane Surface(170) = {169};
Line Loop(171) = {41,67,-45,-61};
Plane Surface(172) = {171};
Line Loop(173) = {42,57,-46,-67};
Plane Surface(174) = {173};
Line Loop(175) = {43,65,-47,-54};
Plane Surface(176) = {175};
Line Loop(177) = {44,62,-48,-65};
Plane Surface(178) = {177};
Line Loop(179) = {45,68,-49,-62};
Plane Surface(180) = {179};
Line Loop(181) = {46,58,-50,-68};
Plane Surface(182) = {181};
Line Loop(183) = {119,-123,-118,50};
Plane Surface(184) = {183};
Line Loop(185) = {73,47,-88,-13};
Plane Surface(186) = {185};
Line Loop(187) = {14,87,-48,-88};
Plane Surface(188) = {187};
Line Loop(189) = {49,-89,-15,87};
Plane Surface(190) = {189};
Line Loop(191) = {89,118,-122,-117};
Plane Surface(192) = {191};
Line Loop(193) = {20,73,-54,-72};
Plane Surface(194) = {193};
Line Loop(195) = {72,-53,-71,19};
Plane Surface(196) = {195};
Line Loop(197) = {71,-52,-70,18};
Plane Surface(198) = {197};
Line Loop(199) = {17,70,-51,-69};
Plane Surface(200) = {199};
Line Loop(201) = {74,-55,-69,21};
Plane Surface(202) = {201};
Line Loop(203) = {74,56,-75,-22};
Plane Surface(204) = {203};
Line Loop(205) = {76,-57,-75,23};
Plane Surface(206) = {205};
Line Loop(207) = {76,58,-77,-24};
Plane Surface(208) = {207};
Line Loop(209) = {124,-120,77,119};
Plane Surface(210) = {209};
Line Loop(211) = {72,43,-85,-9};
Ruled Surface(212) = {211};
Line Loop(213) = {84,-44,-85,10};
Ruled Surface(214) = {213};
Line Loop(215) = {84,45,-86,-11};
Ruled Surface(216) = {215};
Line Loop(217) = {86,46,-76,-12};
Ruled Surface(218) = {217};
Line Loop(219) = {82,-39,-71,5};
Ruled Surface(220) = {219};
Line Loop(221) = {81,-40,-82,6};
Ruled Surface(222) = {221};
Line Loop(223) = {83,-41,-81,7};
Ruled Surface(224) = {223};
Line Loop(225) = {83,42,-75,-8};
Ruled Surface(226) = {225};
Line Loop(227) = {79,-35,-70,-1};
Plane Surface(228) = {227};
Line Loop(229) = {69,59,-78,-25};
Plane Surface(230) = {229};
Line Loop(231) = {80,38,-74,-4};
Plane Surface(232) = {231};
Line Loop(233) = {80,-37,-78,3};
Plane Surface(234) = {233};
Line Loop(235) = {79,36,-78,-2};
Plane Surface(236) = {235};
Line Loop(237) = {77,-50,-89,16};
Plane Surface(238) = {237};
Line Loop(239) = {122,123,124,121};
Plane Surface(240) = {239};
Line Loop(241) = {110,-8,-116,97};
Ruled Surface(242) = {241};
Line Loop(243) = {116,-7,-114,96};
Ruled Surface(244) = {243};
Line Loop(245) = {114,-6,-115,95};
Ruled Surface(246) = {245};
Line Loop(247) = {115,-5,-108,94};
Ruled Surface(248) = {247};
Line Loop(249) = {113,4,-109,-93};
Plane Surface(250) = {249};
Line Loop(251) = {98,112,-25,-106};
Plane Surface(252) = {251};
Line Loop(253) = {111,1,-107,90};
Plane Surface(254) = {253};
Line Loop(255) = {113,-3,-112,92};
Plane Surface(256) = {255};
Line Loop(257) = {2,-112,-91,111};
Plane Surface(258) = {257};
Line Loop(259) = {22,-110,-103,109};
Plane Surface(260) = {259};
Line Loop(261) = {109,-21,-106,102};
Plane Surface(262) = {261};
Line Loop(263) = {17,-107,-100,106};
Plane Surface(264) = {263};
Line Loop(265) = {108,-18,-107,101};
Plane Surface(266) = {265};
Line Loop(267) = {102,-93,-92,-98};
Plane Surface(268) = {267};
Line Loop(269) = {100,90,91,-98};
Plane Surface(270) = {269};
Line Loop(271) = {93,103,-97,-105};
Plane Surface(272) = {271};
Line Loop(273) = {105,-96,-99,92};
Plane Surface(274) = {273};
Line Loop(275) = {91,99,-95,-104};
Plane Surface(276) = {275};
Line Loop(277) = {90,104,-94,-101};
Plane Surface(278) = {277};
Line Loop(279) = {89,-68,-86,34};
Ruled Surface(280) = {279};
Line Loop(281) = {31,88,-65,-85};
Ruled Surface(282) = {281};
Line Loop(283) = {87,-62,-84,28};
Plane Surface(284) = {283};
Line Loop(285) = {27,84,-61,-81};
Plane Surface(286) = {285};
Line Loop(287) = {78,60,-81,-26};
Plane Surface(288) = {287};
Line Loop(289) = {86,-67,-83,33};
Ruled Surface(290) = {289};
Line Loop(291) = {85,-64,-82,30};
Ruled Surface(292) = {291};
Line Loop(293) = {82,-63,-79,29};
Ruled Surface(294) = {293};
Line Loop(295) = {80,66,-83,-32};
Ruled Surface(296) = {295};
Line Loop(297) = {113,32,-116,-105};
Ruled Surface(298) = {297};
Line Loop(299) = {111,29,-115,-104};
Ruled Surface(300) = {299};
Line Loop(313) = {112,26,-114,-99};
Plane Surface(314) = {313};
//--------------------------------------------------//
Surface Loop(301) = {186,194,146,176,282,212};
Volume(302) = {301};
Surface Loop(303) = {188,148,178,214,282,284};
Volume(304) = {303};
Surface Loop(305) = {190,180,150,284,280,216};
Volume(306) = {305};
Surface Loop(307) = {208,182,152,238,280,218};
Volume(308) = {307};
Surface Loop(309) = {154,210,240,192,184,238};
Volume(310) = {309};
Surface Loop(311) = {272,260,242,136,298,250};
Volume(312) = {311};
Surface Loop(315) = {274,244,134,298,314,256};
Volume(316) = {315};
Surface Loop(317) = {314,246,276,132,300,258};
Volume(318) = {317};
Surface Loop(319) = {266,248,278,300,130,254};
Volume(320) = {319};
Surface Loop(321) = {270,264,258,126,254,252};
Volume(322) = {321};
Surface Loop(323) = {268,262,256,252,250,128};
Volume(324) = {323};
Surface Loop(325) = {174,206,144,226,218,290};
Volume(326) = {325};
Surface Loop(327) = {168,196,138,212,292,220};
Volume(328) = {327};
Surface Loop(329) = {214,170,140,286,292,222};
Volume(330) = {329};
Surface Loop(331) = {216,172,142,290,286,224};
Volume(332) = {331};
Surface Loop(333) = {130,160,198,228,294,220};
Volume(334) = {333};
Surface Loop(335) = {136,226,296,166,204,232};
Volume(336) = {335};
Surface Loop(337) = {162,294,222,236,288,132};
Volume(338) = {337};
Surface Loop(339) = {134,164,296,288,224,234};
Volume(340) = {339};
Surface Loop(341) = {126,156,200,228,236,230};
Volume(342) = {341};
Surface Loop(343) = {158,202,230,234,232,128};
Volume(344) = {343};
//========== DATA MESH ===================================//
// Lauflaenge des Strahls //
Transfinite Line
{88,73,87,72,85,84,71,82,79,81,70,78,69,80,74,83,75,86,89,122,76,124,77} =
nodesLaenge Using Progression progLaenge;
// Strahl //
Transfinite Line
{1,17,25,2,3,4,21,5,6,7,8,100,90,91,92,93,102,98,94,95,96,97,9,10,11,12,13,14,15,16,121,50,123,49,48,47,43,44,45,46,39,40,41,42,35,36,37,38,51,59,55}
= nodesStrahl Using Progression progStrahl;
// Hoehe Konus innen //
Transfinite Line {52,63,60,66,56,18,29,26,32,22,101,104,99,105,103} =
nodesKonInnen Using Progression progKonInnen;
// Hoehe Konus aussen //
Transfinite Line {19,30,27,33,23,53,64,61,67,57} = nodesKonAussen Using
Progression progKonAussen;
// Hoehe Raum //
Transfinite Line {20,31,28,34,24,54,65,62,68,58} = nodesRaum Using Progression
progRaum;
// Einlass //
Transfinite Line {106,107,108,109,110,113,112,111,115,116,114} =
nodesEinlassLaenge Using Progression progEinlassLaenge;
// Auslass //
Transfinite Line {117,120,118,119} = nodesAuslassHoehe Using Progression
progAuslassHoehe;
//--------------------------------------------------//
Transfinite Surface {126} = {1,2,10,6};
Transfinite Surface {128} = {14,18,1,2};
Transfinite Surface {130} = {6,7,11,10};
Transfinite Surface {132} = {10,11,3,2};
Transfinite Surface {134} = {19,18,2,3};
Transfinite Surface {136} = {18,19,15,14};
Transfinite Surface {138} = {7,8,12,11};
Transfinite Surface {140} = {11,3,4,12};
Transfinite Surface {142} = {3,4,20,19};
Transfinite Surface {144} = {15,19,20,16};
Transfinite Surface {146} = {12,8,9,13};
Transfinite Surface {148} = {12,4,5,13};
Transfinite Surface {150} = {20,4,5,21};
Transfinite Surface {152} = {16,20,21,17};
Transfinite Surface {154} = {21,55,54,17};
Transfinite Surface {156} = {27,31,23,22};
Transfinite Surface {158} = {22,23,39,35};
Transfinite Surface {160} = {31,27,28,32};
Transfinite Surface {162} = {31,23,24,32};
Transfinite Surface {164} = {40,24,23,39};
Transfinite Surface {166} = {35,36,40,39};
Transfinite Surface {168} = {32,33,29,28};
Transfinite Surface {170} = {32,24,25,33};
Transfinite Surface {172} = {40,41,25,24};
Transfinite Surface {174} = {41,40,36,37};
Transfinite Surface {176} = {33,34,30,29};
Transfinite Surface {178} = {34,33,25,26};
Transfinite Surface {180} = {42,26,25,41};
Transfinite Surface {182} = {41,37,38,42};
Transfinite Surface {184} = {57,42,38,56};
Transfinite Surface {186} = {9,13,34,30};
Transfinite Surface {188} = {5,13,34,26};
Transfinite Surface {190} = {21,5,26,42};
Transfinite Surface {192} = {55,21,42,57};
Transfinite Surface {194} = {8,9,30,29};
Transfinite Surface {196} = {7,8,29,28};
Transfinite Surface {198} = {6,7,28,27};
Transfinite Surface {200} = {27,22,1,6};
Transfinite Surface {202} = {22,35,14,1};
Transfinite Surface {204} = {15,14,35,36};
Transfinite Surface {206} = {36,37,16,15};
Transfinite Surface {208} = {16,37,38,17};
Transfinite Surface {210} = {54,17,38,56};
Transfinite Surface {212} = {33,29,8,12};
Transfinite Surface {214} = {12,33,25,4};
Transfinite Surface {216} = {4,20,41,25};
Transfinite Surface {218} = {41,37,16,20};
Transfinite Surface {220} = {32,28,7,11};
Transfinite Surface {222} = {3,11,32,24};
Transfinite Surface {224} = {24,40,19,3};
Transfinite Surface {226} = {15,19,40,36};
Transfinite Surface {228} = {31,27,6,10};
Transfinite Surface {230} = {23,22,1,2};
Transfinite Surface {232} = {39,35,14,18};
Transfinite Surface {234} = {23,39,18,2};
Transfinite Surface {236} = {2,10,31,23};
Transfinite Surface {238} = {42,38,17,21};
Transfinite Surface {240} = {54,55,57,56};
Transfinite Surface {242} = {51,53,19,15};
Transfinite Surface {244} = {53,45,3,19};
Transfinite Surface {246} = {11,3,45,49};
Transfinite Surface {248} = {7,47,49,11};
Transfinite Surface {250} = {50,52,18,14};
Transfinite Surface {252} = {1,2,44,43};
Transfinite Surface {254} = {6,10,48,46};
Transfinite Surface {256} = {2,18,52,44};
Transfinite Surface {258} = {10,2,44,48};
Transfinite Surface {260} = {51,50,14,15};
Transfinite Surface {262} = {50,43,1,14};
Transfinite Surface {264} = {1,6,46,43};
Transfinite Surface {266} = {46,47,7,6};
Transfinite Surface {268} = {43,44,52,50};
Transfinite Surface {270} = {46,48,44,43};
Transfinite Surface {272} = {50,52,53,51};
Transfinite Surface {274} = {52,44,45,53};
Transfinite Surface {276} = {48,49,45,44};
Transfinite Surface {278} = {46,48,49,47};
Transfinite Surface {280} = {20,21,42,41};
Transfinite Surface {282} = {34,33,12,13};
Transfinite Surface {284} = {26,25,4,5};
Transfinite Surface {286} = {3,4,25,24};
Transfinite Surface {288} = {2,3,24,23};
Transfinite Surface {290} = {41,40,19,20};
Transfinite Surface {292} = {11,12,33,32};
Transfinite Surface {294} = {10,11,32,31};
Transfinite Surface {296} = {40,39,18,19};
Transfinite Surface {298} = {18,52,53,19};
Transfinite Surface {300} = {48,49,11,10};
Transfinite Surface {314} = {3,2,44,45};
//--------------------------------------------------//
Transfinite Volume{344} = {23, 22, 35, 39, 18, 2, 1, 14};
Transfinite Volume{342} = {31, 27, 22, 23, 2, 10, 6, 1};
Transfinite Volume{340} = {23, 39, 40, 24, 3, 2, 18, 19};
Transfinite Volume{338} = {32, 31, 23, 24, 3, 11, 10, 2};
Transfinite Volume{336} = {40, 39, 35, 36, 15, 19, 18, 14};
Transfinite Volume{334} = {32, 28, 27, 31, 10, 11, 7, 6};
Transfinite Volume{332} = {3, 4, 20, 19, 40, 24, 25, 41};
Transfinite Volume{330} = {12, 11, 3, 4, 25, 33, 32, 24};
Transfinite Volume{328} = {12, 8, 7, 11, 32, 33, 29, 28};
Transfinite Volume{326} = {20, 19, 15, 16, 37, 41, 40, 36};
Transfinite Volume{324} = {2, 1, 14, 18, 52, 44, 43, 50};
Transfinite Volume{322} = {10, 6, 1, 2, 44, 48, 46, 43};
Transfinite Volume{320} = {11, 49, 47, 7, 6, 10, 48, 46};
Transfinite Volume{318} = {11, 10, 2, 3, 45, 49, 48, 44};
Transfinite Volume{316} = {18, 2, 44, 52, 53, 19, 3, 45};
Transfinite Volume{312} = {18, 52, 50, 14, 15, 19, 53, 51};
Transfinite Volume{310} = {54, 56, 57, 55, 21, 17, 38, 42};
Transfinite Volume{308} = {21, 17, 38, 42, 41, 20, 16, 37};
Transfinite Volume{306} = {41, 20, 4, 25, 26, 42, 21, 5};
Transfinite Volume{304} = {12, 33, 25, 4, 5, 13, 34, 26};
Transfinite Volume{302} = {33, 12, 8, 29, 30, 34, 13, 9};
//--------------------------------------------------//
Recombine Surface {168, 160, 156, 176, 158, 162, 164, 166, 170, 172, 174, 178,
182, 180, 184, 198, 288, 204, 200, 228, 236, 230, 202, 234, 232, 196, 286, 206,
194, 186, 212, 220, 292, 138, 294, 130, 126, 222, 128, 132, 296, 224, 134, 136,
226, 214, 140, 146, 282, 290, 142, 216, 144, 218, 266, 248, 254, 264, 300, 278,
258, 252, 270, 262, 246, 256, 250, 314, 268, 276, 260, 298, 244, 274, 272, 242,
208, 284, 188, 148, 210, 280, 238, 240, 152, 190, 150, 192, 154};
//Physical Surface(345) = {188, 190, 148, 146, 150, 152, 138, 140, 142, 144,
154, 248, 246, 278, 276, 270, 268, 272, 274, 244, 242, 192, 240, 266, 264, 262,
2
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh