[algc]

Hi all,


I'm trying to open with a .geo file that I know works on the system on which it was created. However, I get a lot of errors when I try to work with it myself. WeDoes anyone know what the problem might be?


__________________________________________________ ______________


// Gmsh project created on Fri Oct 21 13:21:18 2022

//SetFactory("Built-in");
SetFactory("OpenCASCADE");

c= 1;
h = 3.36;
xLe = 0;
yLe = 0;
xStartDomain = -10*c;
xTEConnection = c*3.5+xLe;
xLastBoxMiddle = 4*c;
xLastBox = 6*c;
yMinA = -0.0778244;
yMaxA = 0.0619795;
t = yMaxA - yMinA;

xStartRegion2 = -4*c;
xStartRegion2Top = -2*c;
yGround = yMinA-h*c;

xStartGround = -10*c + xLe + (h*c -yMinA);
xEndGround = xLe + 11*c;

yOutletBottom = yGround;
yOutletTop = h*c - yMinA;
xOutlet = xEndGround;

nArc = 10;
nTop1x = 25;
nTop2x = 50;
nBot1x = nTop1x;
nBot2x = nTop2x;
ny = 100;
nWake1x = 30;
nWake2x = 15;
nProgress = 1.1;

Point(1) = {1, 0.0604, 0, 1};
Point(2) = {0.999605, 0.0605784, 0, 1};
Point(3) = {0.998419, 0.0610433, 0, 1};
Point(4) = {0.996446, 0.0616005, 0, 1};
Point(5) = {0.993688, 0.0619795, 0, 1};
Point(6) = {0.990149, 0.0619138, 0, 1};
Point(7) = {0.985835, 0.0612488, 0, 1};
Point(8) = {0.980752, 0.0600775, 0, 1};
Point(9) = {0.97491, 0.058799, 0, 1};
Point(10) = {0.968317, 0.0575758, 0, 1};
Point(11) = {0.960983, 0.0563594, 0, 1};
Point(12) = {0.95292, 0.0550487, 0, 1};
Point(13) = {0.944141, 0.0536418, 0, 1};
Point(14) = {0.93466, 0.0522031, 0, 1};
Point(15) = {0.924492, 0.0507849, 0, 1};
Point(16) = {0.913652, 0.0494166, 0, 1};
Point(17) = {0.902158, 0.0481248, 0, 1};
Point(18) = {0.890029, 0.0469341, 0, 1};
Point(19) = {0.877282, 0.0458759, 0, 1};
Point(20) = {0.863939, 0.0449858, 0, 1};
Point(21) = {0.850021, 0.0443008, 0, 1};
Point(22) = {0.835549, 0.043855, 0, 1};
Point(23) = {0.820547, 0.0436621, 0, 1};
Point(24) = {0.805038, 0.0437259, 0, 1};
Point(25) = {0.789046, 0.0440436, 0, 1};
Point(26) = {0.772597, 0.0445904, 0, 1};
Point(27) = {0.755718, 0.0453233, 0, 1};
Point(28) = {0.738434, 0.0461889, 0, 1};
Point(29) = {0.720773, 0.0471253, 0, 1};
Point(30) = {0.702763, 0.0480625, 0, 1};
Point(31) = {0.684432, 0.0489309, 0, 1};
Point(32) = {0.66581, 0.0497114, 0, 1};
Point(33) = {0.646925, 0.0504071, 0, 1};
Point(34) = {0.627808, 0.0510246, 0, 1};
Point(35) = {0.608489, 0.0515775, 0, 1};
Point(36) = {0.588999, 0.0520806, 0, 1};
Point(37) = {0.569368, 0.0525307, 0, 1};
Point(38) = {0.549627, 0.052908, 0, 1};
Point(39) = {0.529808, 0.053199, 0, 1};
Point(40) = {0.509941, 0.0534152, 0, 1};
Point(41) = {0.490059, 0.0535737, 0, 1};
Point(42) = {0.470192, 0.0536746, 0, 1};
Point(43) = {0.450373, 0.0537003, 0, 1};
Point(44) = {0.430632, 0.053639, 0, 1};
Point(45) = {0.411001, 0.0535024, 0, 1};
Point(46) = {0.391511, 0.0533079, 0, 1};
Point(47) = {0.372192, 0.0530609, 0, 1};
Point(48) = {0.353075, 0.052754, 0, 1};
Point(49) = {0.33419, 0.05238, 0, 1};
Point(50) = {0.315568, 0.0519341, 0, 1};
Point(51) = {0.297237, 0.051413, 0, 1};
Point(52) = {0.279227, 0.0508255, 0, 1};
Point(53) = {0.261566, 0.0502075, 0, 1};
Point(54) = {0.244282, 0.0495947, 0, 1};
Point(55) = {0.227403, 0.0489997, 0, 1};
Point(56) = {0.210954, 0.0484056, 0, 1};
Point(57) = {0.194962, 0.0477958, 0, 1};
Point(58) = {0.179453, 0.0471692, 0, 1};
Point(59) = {0.164451, 0.0465119, 0, 1};
Point(60) = {0.149979, 0.0457238, 0, 1};
Point(61) = {0.136061, 0.0448191, 0, 1};
Point(62) = {0.122718, 0.0438133, 0, 1};
Point(63) = {0.109971, 0.0426546, 0, 1};
Point(64) = {0.0978417, 0.0414712, 0, 1};
Point(65) = {0.086348, 0.0402678, 0, 1};
Point(66) = {0.0755083, 0.0388983, 0, 1};
Point(67) = {0.0653399, 0.0376177, 0, 1};
Point(68) = {0.0558587, 0.0362163, 0, 1};
Point(69) = {0.0470798, 0.0349255, 0, 1};
Point(70) = {0.0390171, 0.0336265, 0, 1};
Point(71) = {0.0316833, 0.0324608, 0, 1};
Point(72) = {0.0250899, 0.0312981, 0, 1};
Point(73) = {0.0192475, 0.0296216, 0, 1};
Point(74) = {0.0141653, 0.026429, 0, 1};
Point(75) = {0.00985132, 0.0229857, 0, 1};
Point(76) = {0.00631233, 0.0191972, 0, 1};
Point(77) = {0.00355396, 0.0142819, 0, 1};
Point(78) = {0.00158058, 0.00992171, 0, 1};
Point(79) = {0.000395301, 0.00402408, 0, 1};
Point(80) = {0, 0, 0, 1};
Point(81) = {0.000395301, -0.00382523, 0, 1};
Point(82) = {0.00158058, -0.00966363, 0, 1};
Point(83) = {0.00355396, -0.0141751, 0, 1};
Point(84) = {0.00631233, -0.0189984, 0, 1};
Point(85) = {0.00985132, -0.0227577, 0, 1};
Point(86) = {0.0141653, -0.0261142, 0, 1};
Point(87) = {0.0192475, -0.0290608, 0, 1};
Point(88) = {0.0250899, -0.0320986, 0, 1};
Point(89) = {0.0316833, -0.0354145, 0, 1};
Point(90) = {0.0390171, -0.0389275, 0, 1};
Point(91) = {0.0470798, -0.0426075, 0, 1};
Point(92) = {0.0558587, -0.0464556, 0, 1};
Point(93) = {0.0653399, -0.0502913, 0, 1};
Point(94) = {0.0755083, -0.0541613, 0, 1};
Point(95) = {0.086348, -0.0579993, 0, 1};
Point(96) = {0.0978417, -0.0616434, 0, 1};
Point(97) = {0.109971, -0.0651169, 0, 1};
Point(98) = {0.122718, -0.0683089, 0, 1};
Point(99) = {0.136061, -0.0711266, 0, 1};
Point(100) = {0.149979, -0.0736086, 0, 1};
Point(101) = {0.164451, -0.0755747, 0, 1};
Point(102) = {0.179453, -0.0768899, 0, 1};
Point(103) = {0.194962, -0.0776846, 0, 1};
Point(104) = {0.210954, -0.0778244, 0, 1};
Point(105) = {0.227403, -0.0773511, 0, 1};
Point(106) = {0.244282, -0.0764985, 0, 1};
Point(107) = {0.261566, -0.0755273, 0, 1};
Point(108) = {0.279227, -0.0745046, 0, 1};
Point(109) = {0.297237, -0.0733613, 0, 1};
Point(110) = {0.315568, -0.072032, 0, 1};
Point(111) = {0.33419, -0.0705279, 0, 1};
Point(112) = {0.353075, -0.0688931, 0, 1};
Point(113) = {0.372192, -0.0671575, 0, 1};
Point(114) = {0.391511, -0.0653079, 0, 1};
Point(115) = {0.411001, -0.0633244, 0, 1};
Point(116) = {0.430632, -0.0611936, 0, 1};
Point(117) = {0.450373, -0.0589077, 0, 1};
Point(118) = {0.470192, -0.0564646, 0, 1};
Point(119) = {0.490059, -0.0538821, 0, 1};
Point(120) = {0.509941, -0.0511821, 0, 1};
Point(121) = {0.529808, -0.048365, 0, 1};
Point(122) = {0.549627, -0.0454112, 0, 1};
Point(123) = {0.569368, -0.0423111, 0, 1};
Point(124) = {0.588999, -0.039096, 0, 1};
Point(125) = {0.608489, -0.035806, 0, 1};
Point(126) = {0.627808, -0.0324515, 0, 1};
Point(127) = {0.646925, -0.0290112, 0, 1};
Point(128) = {0.66581, -0.0254699, 0, 1};
Point(129) = {0.684432, -0.0218487, 0, 1};
Point(130) = {0.702763, -0.018181, 0, 1};
Point(131) = {0.720773, -0.0144874, 0, 1};
Point(132) = {0.738434, -0.0107618, 0, 1};
Point(133) = {0.755718, -0.00699652, 0, 1};
Point(134) = {0.772597, -0.00320282, 0, 1};
Point(135) = {0.789046, 0.000585036, 0, 1};
Point(136) = {0.805038, 0.00433399, 0, 1};
Point(137) = {0.820547, 0.00802701, 0, 1};
Point(138) = {0.835549, 0.011667, 0, 1};
Point(139) = {0.850021, 0.0152554, 0, 1};
Point(140) = {0.863939, 0.0187866, 0, 1};
Point(141) = {0.877282, 0.0222372, 0, 1};
Point(142) = {0.890029, 0.0255826, 0, 1};
Point(143) = {0.902158, 0.0288016, 0, 1};
Point(144) = {0.913652, 0.0318906, 0, 1};
Point(145) = {0.924492, 0.034874, 0, 1};
Point(146) = {0.93466, 0.0377273, 0, 1};
Point(147) = {0.944141, 0.0403881, 0, 1};
Point(148) = {0.95292, 0.0428671, 0, 1};
Point(149) = {0.960983, 0.0452231, 0, 1};
Point(150) = {0.968317, 0.0474328, 0, 1};
Point(151) = {0.97491, 0.0493574, 0, 1};
Point(152) = {0.980752, 0.0509158, 0, 1};
Point(153) = {0.985835, 0.0522553, 0, 1};
Point(154) = {0.990149, 0.0537617, 0, 1};
Point(155) = {0.993688, 0.0555086, 0, 1};
Point(156) = {0.996446, 0.0573218, 0, 1};
Point(157) = {0.998419, 0.0589173, 0, 1};
Point(158) = {0.999605, 0.060011, 0, 1};
//Line(1) = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158,1};



Point(160) = {xStartGround,yGround,0, 1}; //startGround
Point(161) = {xEndGround,yGround,0,1}; //endGround

Point(162) = {xOutlet, yOutletBottom, 0, 1}; //outletBottom
Point(163) = {xOutlet, yOutletTop, 0, 1}; //outletTop

Point(164) = {xLe,yGround,0,1}; //LE Projected on the ground
Point(165) = {xStartGround,yOutletTop,0,1}; //LE Projected to top of domain

Point(166) = {xLe,yLe,0,1}; //LE
Point(167) = {xStartDomain,yLe,0,1}; //Left side of the domain
Point(168) = {xStartGround,yLe,0,1}; //Middle of front circle

Point(169) = {xStartDomain,yOutletTop,0,1};
Point(170) = {xStartDomain,yOutletBottom,0,1};
Point(171) = {xStartRegion2Top,yOutletTop,0,1};
Point(172) = {xStartRegion2,yOutletBottom,0,1};
//lines between domain points
Line(2) = {160,161};
Line(3) = {162,163};
//+
//+
//+
Circle(4) = {160,168,165};
//+
Line(5) = {165, 163};
//+

//+
Line(6) = {168, 170};
//+
Line(7) = {168, 169};


//+

//+
BooleanUnion{ Line{7}; Delete; }{ Line{4}; Delete; }
//+
Recursive Delete {
Line{7}; Line{8};
}
//+
BooleanUnion{ Line{6}; Delete; }{ Line{9}; Delete; }
//+
Recursive Delete {
Line{11};
}
//+
Recursive Delete {
Line{12};
}

Point(187) = {xStartDomain,yOutletBottom,0,1};
Line(8) = {75, 173};
//+
Line(9) = {85, 187};
Line(38) = {173,187};
//+

Line(11) = {51, 171};
//+
Line(12) = {109, 172};

Point(176) = {xTEConnection,yOutletTop,0,1};
Point(177) = {xTEConnection,yOutletBottom,0,1};
//+
Line(15) = {1, 176};
//+
Line(16) = {1, 177};

Point(178) = {xLastBox,yOutletBottom,0,1};
Point(179) = {xLastBox,yOutletTop,0,1};
Point(180) = {xLastBox,yLe,0,1};
Point(181) = {xOutlet,yLe,0,1};
Point(186) = {xLastBoxMiddle,yLe,0,1};
//+
Line(17) = {1, 186};
//+
Line(18) = {186, 181};
//+


Line(19) = {186, 179};
//+
Line(20) = {186, 178};
//+ airfoil splines
Spline(21) = {75, 74, 73, 72, 71, 70, 69, 68, 67, 66, 65, 64, 63, 62, 61, 60, 59, 58, 57, 56, 55, 54, 53, 52, 51}; //top part 1
//+
Spline(22) = {75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85}; //arc
//+
Spline(23) = {85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109}; //bottom part 1
//+
Spline(24) = {51, 50, 49, 48, 47, 46, 45, 44, 43, 42, 41, 40, 39, 38, 37, 36, 35, 34, 33, 32, 31, 30, 29, 28, 27, 26, 25, 24, 23, 22, 21, 20, 19, 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1}; //top part 2
//+
Spline(25) = {109, 110, 111, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 1}; //bottom part 2

//+
//+
Recursive Delete {
Line{5}; Line{2};
}
//+
Line(26) = {165, 171};
//+
Line(27) = {171, 176};
//+
Line(28) = {176, 179};
//+
Line(29) = {179, 163};
//+
Line(30) = {160, 172};
//+
Line(31) = {172, 177};
//+
Line(32) = {177, 178};
//+
Line(33) = {178, 162};
//+
//+
Point(182) = {xStartGround-1.5*c,yOutletTop,0,1};
Point(183) = {xStartGround+2*c,yOutletTop,0,1};
Point(184) = {xStartGround-1.5*c,yOutletBottom,0,1};
Point(185) = {xStartGround-2*c,yOutletBottom,0,1};

//

//+
Bezier(34) = {173, 182, 183, 171};
//+
Bezier(35) = {175, 184, 185, 172};
//+
Recursive Delete {
Line{10}; Line{26}; Line{13}; Line{30};
}
//+
Recursive Delete {
Line{3};
}
//+
Line(36) = {181, 163};
//+
Line(37) = {181, 162};
//+

Recursive Delete {
Line{14}; Line{35};
}
//+
Line(39) = {187, 172};

//+
Point(188) = {-8, 1.6, 0, 1.0};
//+
Spline(40) = {187, 188, 171};
//+
BooleanUnion{ Line{8}; Delete; }{ Line{40}; Delete; }
//+
//+
Recursive Delete {
Line{41}; Line{34}; Line{38};
}
//Line(44) = {75, 189};

Transfinite Line {42} = nArc Using Bump 0.3; //arcDomain
Transfinite Line {43} = nTop1x Using Progression 1; //BezierTop
Transfinite Line {27} = nTop2x Using Progression 1; //TopAirfoilDomain
Transfinite Line {28} = nWake1x Using Progression 1; //TopWake1Domain
Transfinite Line {29} = nWake2x Using Progression 1; //TopWake2Domain
Transfinite Line {36} = ny Using Progression nProgress; //OutletUpperDomain
Transfinite Line {37} = ny Using Progression nProgress; //OutletLowerDomain
Transfinite Line {33} = nWake2x Using Progression 1; //BottomWake2Domain
Transfinite Line {32} = nWake1x Using Progression 1; //BottomWake1Domain
Transfinite Line {31} = nBot2x Using Progression 1; //BottomAirfoilDomain
Transfinite Line {39} = nBot1x Using Progression 1; //BezierBottom
Transfinite Line {9} = ny Using Progression nProgress; //ArcLower
Transfinite Line {40} = ny Using Progression nProgress; //ArcUpper
Transfinite Line {11} = ny Using Progression nProgress; //AirfoilUp
Transfinite Line {15} = ny Using Progression nProgress; //TeUp
Transfinite Line {19} = ny Using Progression nProgress; //Wake1Up
Transfinite Line {20} = ny Using Progression nProgress; //Wake1Down
Transfinite Line {16} = ny Using Progression nProgress; //TeDown
Transfinite Line {12} = ny Using Progression nProgress; //airfoilDown
Transfinite Line {22} = nArc Using Progression 1; //AirofilArc
Transfinite Line {21} = nTop1x Using Progression 1; //AirfoilTop1
Transfinite Line {24} = nTop2x Using Progression 1; //AirfoilTop2
Transfinite Line {23} = nBot1x Using Progression 1; //AirfoilBottom1
Transfinite Line {25} = nBot2x Using Progression 1; //AirfoilBottom2
Transfinite Line {17} = nWake1x Using Progression nProgress; //wake1X
Transfinite Line {18} = nWake2x Using Progression nProgress; //wake2X
//+
Line Loop(1) = {22, 9, -42, -40};
//+
Plane Surface(1) = {1};
//+
Line Loop(2) = {21, 11, -40, -43};
//+
Plane Surface(2) = {2};
//+
Line Loop(10) = {24, 15, -27, -11};
//+
Plane Surface(3) = {10};
//+
Line Loop(4) = {17, 19, -28, -15};
//+
Plane Surface(4) = {4};
//+
Line Loop(5) = {18, 36, -29, -19};
//+
Plane Surface(5) = {5};
//+
Line Loop(6) = {18, 37, -33, -20};
//+
Plane Surface(6) = {6};
//+
Line Loop(7) = {17, 20, -32, -16};
//+
Plane Surface(7) = {7};
//+
Line Loop(8) = {25, 16, -31, -12};
//+
Plane Surface(8) = {8};
//+
Line Loop(9) = {23, 12, -39, -9};
//+
Plane Surface(9) = {9};
//+
Transfinite Surface {1};
//+
Transfinite Surface {2};
//+
Transfinite Surface {3};
//+
Transfinite Surface {4};
//+
Transfinite Surface {5};
//+
Transfinite Surface {6};
//+
Transfinite Surface {7};
//+
Transfinite Surface {8};
//+
Transfinite Surface {9};
//+
Recombine Surface {1, 2, 3, 4, 5, 6, 7, 8, 9};

//+


//+

//+
//+
Extrude {0, 0, 1} {
Surface{2}; Surface{1}; Surface{9}; Line{42}; Line{43}; Surface{3}; Line{27}; Line{11}; Line{40}; Line{15}; Line{28}; Surface{5}; Line{19}; Line{29}; Line{36}; Line{37}; Surface{6}; Line{18}; Line{33}; Line{20}; Surface{7}; Line{16}; Line{32}; Line{31}; Line{12}; Line{9}; Line{39}; Surface{4}; Surface{8}; Layers{1}; Recombine;
}
//+
Physical Volume("Fluid") = {3, 2, 1, 4, 9, 7, 8, 5, 6};
//+
Physical Surface("Inlet") = {21, 43, 37, 33, 17, 12, 25, 40, 29};
//+
Physical Surface("Outlet") = {28, 32};
//+
Physical Surface("Side") = {14, 26, 41, 31, 35, 39, 44, 22, 18, 1, 9, 2, 3, 8, 4, 7, 6, 5};
//+
Physical Surface("Airfoil") = {42, 19, 23, 10, 15};
//+
__________________________________________________ _________________