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

Reply via email to