|
| 1 | +N=0.00002; |
| 2 | +dy=5e-5; |
| 3 | +//+ |
| 4 | +Point(1) = {0, 0, 0, N}; |
| 5 | +Point(2) = {0.01, 0, 0, 0.2*N}; |
| 6 | +Point(3) = {0.01, dy, 0, 0.2*N}; |
| 7 | +Point(4) = {0.0, dy, 0, N}; |
| 8 | +Point(9) = {0.0100+dy, 2*dy, 0, N}; |
| 9 | +Point(10) = {0.010, 2*dy, 0, N}; |
| 10 | +Point(11) = {0.0100+dy, 9*dy, 0, N}; |
| 11 | +Point(13) = {0.0100+8*dy, 0.00, 0, N}; |
| 12 | +Point(14) = {0.0100+8*dy, 2*dy, 0, N}; |
| 13 | +//+ |
| 14 | +Line(1) = {1, 4}; |
| 15 | +Line(2) = {4, 3}; |
| 16 | +Line(3) = {3, 2}; |
| 17 | +Line(4) = {2, 1}; |
| 18 | +Line(6) = {9, 11}; |
| 19 | +Line(10) = {13, 14}; |
| 20 | +Line(11) = {2, 13}; |
| 21 | +//+ |
| 22 | +Circle(5) = {3, 10, 9}; |
| 23 | +Circle(9) = {14, 9, 11}; |
| 24 | +//+ |
| 25 | +Transfinite Curve {3, 1} = 2 Using Progression 1; |
| 26 | +Transfinite Curve {4, 2} = 1000 Using Progression 1; |
| 27 | +//+ |
| 28 | +//+ |
| 29 | +Physical Curve("inlet", 12) = {1}; |
| 30 | +Physical Curve("symmetry_bottom", 13) = {4, 11}; |
| 31 | +Physical Curve("outlet", 14) = {9, 10}; |
| 32 | +Physical Curve("wall", 15) = {6, 5, 2}; |
| 33 | +//+ |
| 34 | +Curve Loop(1) = {2, 3, 4, 1}; |
| 35 | +//+ |
| 36 | +Plane Surface(1) = {1}; |
| 37 | +//+ |
| 38 | +Curve Loop(2) = {6, -9, -10, -11, -3, 5}; |
| 39 | +//+ |
| 40 | +Plane Surface(2) = {2}; |
| 41 | +//+ |
| 42 | +Transfinite Surface {1}; |
| 43 | +//+ |
| 44 | +Physical Surface("interior", 16) = {1, 2}; |
| 45 | +//+ |
| 46 | +Recombine Surface {1, 2}; |
0 commit comments