1 // Parametres geometriques
2 r1 = 200*Cos(Pi/4)/1000;
3 r2 = 200*Cos(Pi/4)/1000;
10 // Parametres de maillage
12 rCells = 10/2; rRatio = 0.85;
14 S1Cells = 30/2; S1ratio = 1;
16 S2Cells = 35/2; S2ratio = 0.95;
18 S3Cells = 20/2; S3ratio = 1;
20 Point(1) = {r1, r1, h4};
21 Point(2) = {r1, r1, h3};
22 Point(3) = {r2, r2, h2};
23 Point(4) = {r2, r2, h1};
24 Point(5) = {0, 0, h1};
25 Point(6) = {0, 0, h2};
26 Point(7) = {0, 0, h3};
27 Point(8) = {0, 0, h4};
40 Line Loop(11) = {1, 2, -3, -4};
41 Ruled Surface(12) = {11};
42 Line Loop(13) = {5, -6, -7, 3};
43 Ruled Surface(14) = {13};
44 Line Loop(15) = {8, -9, -10, 6};
45 Ruled Surface(16) = {15};
46 Transfinite Line {1, 3, 6, 9} = rCells Using Progression rRatio;
47 Transfinite Line {4, 2} = S1Cells Using Progression S1ratio;
48 Transfinite Line {7, 5} = S2Cells Using Progression S2ratio;
49 Transfinite Line {10, 8} = S3Cells Using Progression S3ratio;
50 Transfinite Surface {12} = {8, 1, 2, 7};
51 Transfinite Surface {14} = {7, 2, 3, 6};
52 Transfinite Surface {16} = {6, 3, 4, 5};
53 Recombine Surface {12, 14, 16};
54 Extrude {{0, 0, 1}, {0, 0, 0}, Pi/2} {
59 Extrude {{0, 0, 1}, {0, 0, 0}, Pi/2} {
64 Extrude {{0, 0, 1}, {0, 0, 0}, Pi/2} {
65 Surface{84, 101, 118};
69 Extrude {{0, 0, 1}, {0, 0, 0}, Pi/2} {
70 Surface{135, 152, 169};
74 Physical Surface("entree") = {126, 75, 24, 177};
75 Physical Surface("S1") = {28, 181, 130, 79};
76 Physical Surface("S2") = {93, 42, 193, 144};
77 Physical Surface("S3") = {110, 59, 205, 161};
78 Physical Surface("fond") = {113, 62, 208, 164};
79 Physical Volume("fluide") = {4, 7, 10, 1, 5, 8, 11, 2, 9, 12, 3, 6};