I'm trying to mesh a simple airfoil case but gmsh "hangs" without producing a
result. I haven't used gmsh before but I cant seem to figure out what I'm doing
wrong. By the way, since my geometry is 2D I realize I can extrude the 2D to
make a 3D surface but other cases I want to run will be more 3D so I want to
do it this way.
// NACA0012 Aerofoil Shape and Spline fit
ps_lc = 0.005;
Point(1000) = { 1.00000000, 0.00000000, 0.00000000, ps_lc};
Point(1001) = { 0.98808000, -0.00221820, 0.00000000, ps_lc};
Point(1002) = { 0.97519000, -0.00441040, 0.00000000, ps_lc};
Point(1003) = { 0.96040000, -0.00668320, 0.00000000, ps_lc};
Point(1004) = { 0.94333000, -0.00894260, 0.00000000, ps_lc};
Point(1005) = { 0.92328000, -0.01153900, 0.00000000, ps_lc};
Point(1006) = { 0.90000000, -0.01447700, 0.00000000, ps_lc};
Point(1007) = { 0.87500000, -0.01754500, 0.00000000, ps_lc};
Point(1008) = { 0.85000000, -0.02052600, 0.00000000, ps_lc};
Point(1009) = { 0.82500000, -0.02342000, 0.00000000, ps_lc};
Point(1010) = { 0.80000000, -0.02623100, 0.00000000, ps_lc};
Point(1011) = { 0.77500000, -0.02895900, 0.00000000, ps_lc};
Point(1012) = { 0.75000000, -0.03160300, 0.00000000, ps_lc};
Point(1013) = { 0.72500000, -0.03416400, 0.00000000, ps_lc};
Point(1014) = { 0.70000000, -0.03663900, 0.00000000, ps_lc};
Point(1015) = { 0.67500000, -0.03902700, 0.00000000, ps_lc};
Point(1016) = { 0.65000000, -0.04132500, 0.00000000, ps_lc};
Point(1017) = { 0.62500000, -0.04352900, 0.00000000, ps_lc};
Point(1018) = { 0.60000000, -0.04563400, 0.00000000, ps_lc};
Point(1019) = { 0.57500000, -0.04763400, 0.00000000, ps_lc};
Point(1020) = { 0.55000000, -0.04952400, 0.00000000, ps_lc};
Point(1021) = { 0.52500000, -0.05129600, 0.00000000, ps_lc};
Point(1022) = { 0.50000000, -0.05294000, 0.00000000, ps_lc};
Point(1023) = { 0.47500000, -0.05444800, 0.00000000, ps_lc};
Point(1024) = { 0.45000000, -0.05580700, 0.00000000, ps_lc};
Point(1025) = { 0.42500000, -0.05700600, 0.00000000, ps_lc};
Point(1026) = { 0.40000000, -0.05803000, 0.00000000, ps_lc};
Point(1027) = { 0.37500000, -0.05886300, 0.00000000, ps_lc};
Point(1028) = { 0.35000000, -0.05948600, 0.00000000, ps_lc};
Point(1029) = { 0.32500000, -0.05987900, 0.00000000, ps_lc};
Point(1030) = { 0.30000000, -0.06001700, 0.00000000, ps_lc};
Point(1031) = { 0.27500000, -0.05987300, 0.00000000, ps_lc};
Point(1032) = { 0.25000000, -0.05941200, 0.00000000, ps_lc};
Point(1033) = { 0.22500000, -0.05859600, 0.00000000, ps_lc};
Point(1034) = { 0.20000000, -0.05737500, 0.00000000, ps_lc};
Point(1035) = { 0.17500000, -0.05568800, 0.00000000, ps_lc};
Point(1036) = { 0.15138000, -0.05359100, 0.00000000, ps_lc};
Point(1037) = { 0.13070000, -0.05128000, 0.00000000, ps_lc};
Point(1038) = { 0.11284000, -0.04885900, 0.00000000, ps_lc};
Point(1039) = { 0.09751400, -0.04640300, 0.00000000, ps_lc};
Point(1040) = { 0.08343700, -0.04377600, 0.00000000, ps_lc};
Point(1041) = { 0.07151400, -0.04121400, 0.00000000, ps_lc};
Point(1042) = { 0.06141100, -0.03874600, 0.00000000, ps_lc};
Point(1043) = { 0.05291100, -0.03641100, 0.00000000, ps_lc};
Point(1044) = { 0.04588300, -0.03426100, 0.00000000, ps_lc};
Point(1045) = { 0.04000000, -0.03227700, 0.00000000, ps_lc};
Point(1046) = { 0.03500000, -0.03043000, 0.00000000, ps_lc};
Point(1047) = { 0.03000000, -0.02840100, 0.00000000, ps_lc};
Point(1048) = { 0.02500000, -0.02614700, 0.00000000, ps_lc};
Point(1049) = { 0.02031000, -0.02376600, 0.00000000, ps_lc};
Point(1050) = { 0.01635200, -0.02148700, 0.00000000, ps_lc};
Point(1051) = { 0.01302200, -0.01930800, 0.00000000, ps_lc};
Point(1052) = { 0.01025700, -0.01724400, 0.00000000, ps_lc};
Point(1053) = { 0.00803040, -0.01534300, 0.00000000, ps_lc};
Point(1054) = { 0.00627150, -0.01362500, 0.00000000, ps_lc};
Point(1055) = { 0.00490630, -0.01210200, 0.00000000, ps_lc};
Point(1056) = { 0.00386450, -0.01077900, 0.00000000, ps_lc};
Point(1057) = { 0.00307410, -0.00964250, 0.00000000, ps_lc};
Point(1058) = { 0.00245990, -0.00864800, 0.00000000, ps_lc};
Point(1059) = { 0.00193200, -0.00768320, 0.00000000, ps_lc};
Point(1060) = { 0.00149900, -0.00678320, 0.00000000, ps_lc};
Point(1061) = { 0.00111400, -0.00586120, 0.00000000, ps_lc};
Point(1062) = { 0.00078050, -0.00491760, 0.00000000, ps_lc};
Point(1063) = { 0.00050290, -0.00395680, 0.00000000, ps_lc};
Point(1064) = { 0.00028350, -0.00297820, 0.00000000, ps_lc};
Point(1065) = { 0.00012610, -0.00199060, 0.00000000, ps_lc};
Point(1066) = { 0.00003160, -0.00099950, 0.00000000, ps_lc};
ss_lc = 0.005;
// This is the leading edge
Point(1067) = { 0.00000000, 0.00000000, 0.00000000, ss_lc};
Point(1068) = { 0.00003160, 0.00099950, 0.00000000, ss_lc};
Point(1069) = { 0.00012610, 0.00199060, 0.00000000, ss_lc};
Point(1070) = { 0.00028350, 0.00297820, 0.00000000, ss_lc};
Point(1071) = { 0.00050290, 0.00395680, 0.00000000, ss_lc};
Point(1072) = { 0.00078050, 0.00491760, 0.00000000, ss_lc};
Point(1073) = { 0.00111400, 0.00586120, 0.00000000, ss_lc};
Point(1074) = { 0.00149900, 0.00678320, 0.00000000, ss_lc};
Point(1075) = { 0.00193200, 0.00768320, 0.00000000, ss_lc};
Point(1076) = { 0.00245990, 0.00864800, 0.00000000, ss_lc};
Point(1077) = { 0.00307410, 0.00964250, 0.00000000, ss_lc};
Point(1078) = { 0.00386450, 0.01077900, 0.00000000, ss_lc};
Point(1079) = { 0.00490630, 0.01210200, 0.00000000, ss_lc};
Point(1080) = { 0.00627150, 0.01362500, 0.00000000, ss_lc};
Point(1081) = { 0.00803040, 0.01534300, 0.00000000, ss_lc};
Point(1082) = { 0.01025700, 0.01724400, 0.00000000, ss_lc};
Point(1083) = { 0.01302200, 0.01930800, 0.00000000, ss_lc};
Point(1084) = { 0.01635200, 0.02148700, 0.00000000, ss_lc};
Point(1085) = { 0.02031000, 0.02376600, 0.00000000, ss_lc};
Point(1086) = { 0.02500000, 0.02614700, 0.00000000, ss_lc};
Point(1087) = { 0.03000000, 0.02840100, 0.00000000, ss_lc};
Point(1088) = { 0.03500000, 0.03043000, 0.00000000, ss_lc};
Point(1089) = { 0.04000000, 0.03227700, 0.00000000, ss_lc};
Point(1090) = { 0.04588300, 0.03426100, 0.00000000, ss_lc};
Point(1091) = { 0.05291100, 0.03641100, 0.00000000, ss_lc};
Point(1092) = { 0.06141100, 0.03874600, 0.00000000, ss_lc};
Point(1093) = { 0.07151400, 0.04121400, 0.00000000, ss_lc};
Point(1094) = { 0.08343700, 0.04377600, 0.00000000, ss_lc};
Point(1095) = { 0.09751400, 0.04640300, 0.00000000, ss_lc};
Point(1096) = { 0.11284000, 0.04885900, 0.00000000, ss_lc};
Point(1097) = { 0.13070000, 0.05128000, 0.00000000, ss_lc};
Point(1098) = { 0.15138000, 0.05359100, 0.00000000, ss_lc};
Point(1099) = { 0.17500000, 0.05568800, 0.00000000, ss_lc};
Point(1100) = { 0.20000000, 0.05737500, 0.00000000, ss_lc};
Point(1101) = { 0.22500000, 0.05859600, 0.00000000, ss_lc};
Point(1102) = { 0.25000000, 0.05941200, 0.00000000, ss_lc};
Point(1103) = { 0.27500000, 0.05987300, 0.00000000, ss_lc};
Point(1104) = { 0.30000000, 0.06001700, 0.00000000, ss_lc};
Point(1105) = { 0.32500000, 0.05987900, 0.00000000, ss_lc};
Point(1106) = { 0.35000000, 0.05948600, 0.00000000, ss_lc};
Point(1107) = { 0.37500000, 0.05886300, 0.00000000, ss_lc};
Point(1108) = { 0.40000000, 0.05803000, 0.00000000, ss_lc};
Point(1109) = { 0.42500000, 0.05700600, 0.00000000, ss_lc};
Point(1110) = { 0.45000000, 0.05580700, 0.00000000, ss_lc};
Point(1111) = { 0.47500000, 0.05444800, 0.00000000, ss_lc};
Point(1112) = { 0.50000000, 0.05294000, 0.00000000, ss_lc};
Point(1113) = { 0.52500000, 0.05129600, 0.00000000, ss_lc};
Point(1114) = { 0.55000000, 0.04952400, 0.00000000, ss_lc};
Point(1115) = { 0.57500000, 0.04763400, 0.00000000, ss_lc};
Point(1116) = { 0.60000000, 0.04563400, 0.00000000, ss_lc};
Point(1117) = { 0.62500000, 0.04352900, 0.00000000, ss_lc};
Point(1118) = { 0.65000000, 0.04132500, 0.00000000, ss_lc};
Point(1119) = { 0.67500000, 0.03902700, 0.00000000, ss_lc};
Point(1120) = { 0.70000000, 0.03663900, 0.00000000, ss_lc};
Point(1121) = { 0.72500000, 0.03416400, 0.00000000, ss_lc};
Point(1122) = { 0.75000000, 0.03160300, 0.00000000, ss_lc};
Point(1123) = { 0.77500000, 0.02895900, 0.00000000, ss_lc};
Point(1124) = { 0.80000000, 0.02623100, 0.00000000, ss_lc};
Point(1125) = { 0.82500000, 0.02342000, 0.00000000, ss_lc};
Point(1126) = { 0.85000000, 0.02052600, 0.00000000, ss_lc};
Point(1127) = { 0.87500000, 0.01754500, 0.00000000, ss_lc};
Point(1128) = { 0.90000000, 0.01447700, 0.00000000, ss_lc};
Point(1129) = { 0.92328000, 0.01153900, 0.00000000, ss_lc};
Point(1130) = { 0.94333000, 0.00894260, 0.00000000, ss_lc};
Point(1131) = { 0.96040000, 0.00668320, 0.00000000, ss_lc};
Point(1132) = { 0.97519000, 0.00441040, 0.00000000, ss_lc};
Point(1133) = { 0.98808000, 0.00221820, 0.00000000, ss_lc};
// Put a Spline fit through the appropriate points - note that we could
// also use BSpline but you end up with a different shape
Spline(100) = {1000:1133,1000};
Point(1300) = { 1.00000000, 0.00000000, 1.00000000, ps_lc};
Point(1301) = { 0.98808000, -0.00221820, 1.00000000, ps_lc};
Point(1302) = { 0.97519000, -0.00441040, 1.00000000, ps_lc};
Point(1303) = { 0.96040000, -0.00668320, 1.00000000, ps_lc};
Point(1304) = { 0.94333000, -0.00894260, 1.00000000, ps_lc};
Point(1305) = { 0.92328000, -0.01153900, 1.00000000, ps_lc};
Point(1306) = { 0.90000000, -0.01447700, 1.00000000, ps_lc};
Point(1307) = { 0.87500000, -0.01754500, 1.00000000, ps_lc};
Point(1308) = { 0.85000000, -0.02052600, 1.00000000, ps_lc};
Point(1309) = { 0.82500000, -0.02342000, 1.00000000, ps_lc};
Point(1310) = { 0.80000000, -0.02623100, 1.00000000, ps_lc};
Point(1311) = { 0.77500000, -0.02895900, 1.00000000, ps_lc};
Point(1312) = { 0.75000000, -0.03160300, 1.00000000, ps_lc};
Point(1313) = { 0.72500000, -0.03416400, 1.00000000, ps_lc};
Point(1314) = { 0.70000000, -0.03663900, 1.00000000, ps_lc};
Point(1315) = { 0.67500000, -0.03902700, 1.00000000, ps_lc};
Point(1316) = { 0.65000000, -0.04132500, 1.00000000, ps_lc};
Point(1317) = { 0.62500000, -0.04352900, 1.00000000, ps_lc};
Point(1318) = { 0.60000000, -0.04563400, 1.00000000, ps_lc};
Point(1319) = { 0.57500000, -0.04763400, 1.00000000, ps_lc};
Point(1320) = { 0.55000000, -0.04952400, 1.00000000, ps_lc};
Point(1321) = { 0.52500000, -0.05129600, 1.00000000, ps_lc};
Point(1322) = { 0.50000000, -0.05294000, 1.00000000, ps_lc};
Point(1323) = { 0.47500000, -0.05444800, 1.00000000, ps_lc};
Point(1324) = { 0.45000000, -0.05580700, 1.00000000, ps_lc};
Point(1325) = { 0.42500000, -0.05700600, 1.00000000, ps_lc};
Point(1326) = { 0.40000000, -0.05803000, 1.00000000, ps_lc};
Point(1327) = { 0.37500000, -0.05886300, 1.00000000, ps_lc};
Point(1328) = { 0.35000000, -0.05948600, 1.00000000, ps_lc};
Point(1329) = { 0.32500000, -0.05987900, 1.00000000, ps_lc};
Point(1330) = { 0.30000000, -0.06001700, 1.00000000, ps_lc};
Point(1331) = { 0.27500000, -0.05987300, 1.00000000, ps_lc};
Point(1332) = { 0.25000000, -0.05941200, 1.00000000, ps_lc};
Point(1333) = { 0.22500000, -0.05859600, 1.00000000, ps_lc};
Point(1334) = { 0.20000000, -0.05737500, 1.00000000, ps_lc};
Point(1335) = { 0.17500000, -0.05568800, 1.00000000, ps_lc};
Point(1336) = { 0.15138000, -0.05359100, 1.00000000, ps_lc};
Point(1337) = { 0.13070000, -0.05128000, 1.00000000, ps_lc};
Point(1338) = { 0.11284000, -0.04885900, 1.00000000, ps_lc};
Point(1339) = { 0.09751400, -0.04640300, 1.00000000, ps_lc};
Point(1340) = { 0.08343700, -0.04377600, 1.00000000, ps_lc};
Point(1341) = { 0.07151400, -0.04121400, 1.00000000, ps_lc};
Point(1342) = { 0.06141100, -0.03874600, 1.00000000, ps_lc};
Point(1343) = { 0.05291100, -0.03641100, 1.00000000, ps_lc};
Point(1344) = { 0.04588300, -0.03426100, 1.00000000, ps_lc};
Point(1345) = { 0.04000000, -0.03227700, 1.00000000, ps_lc};
Point(1346) = { 0.03500000, -0.03043000, 1.00000000, ps_lc};
Point(1347) = { 0.03000000, -0.02840100, 1.00000000, ps_lc};
Point(1348) = { 0.02500000, -0.02614700, 1.00000000, ps_lc};
Point(1349) = { 0.02031000, -0.02376600, 1.00000000, ps_lc};
Point(1350) = { 0.01635200, -0.02148700, 1.00000000, ps_lc};
Point(1351) = { 0.01302200, -0.01930800, 1.00000000, ps_lc};
Point(1352) = { 0.01025700, -0.01724400, 1.00000000, ps_lc};
Point(1353) = { 0.00803040, -0.01534300, 1.00000000, ps_lc};
Point(1354) = { 0.00627150, -0.01362500, 1.00000000, ps_lc};
Point(1355) = { 0.00490630, -0.01210200, 1.00000000, ps_lc};
Point(1356) = { 0.00386450, -0.01077900, 1.00000000, ps_lc};
Point(1357) = { 0.00307410, -0.00964250, 1.00000000, ps_lc};
Point(1358) = { 0.00245990, -0.00864800, 1.00000000, ps_lc};
Point(1359) = { 0.00193200, -0.00768320, 1.00000000, ps_lc};
Point(1360) = { 0.00149900, -0.00678320, 1.00000000, ps_lc};
Point(1361) = { 0.00111400, -0.00586120, 1.00000000, ps_lc};
Point(1362) = { 0.00078050, -0.00491760, 1.00000000, ps_lc};
Point(1363) = { 0.00050290, -0.00395680, 1.00000000, ps_lc};
Point(1364) = { 0.00028350, -0.00297820, 1.00000000, ps_lc};
Point(1365) = { 0.00012610, -0.00199060, 1.00000000, ps_lc};
Point(1366) = { 0.00003160, -0.00099950, 1.00000000, ps_lc};
ss_lc = 0.005;
// This is the leading edge
Point(1367) = { 0.00000000, 0.00000000, 1.00000000, ss_lc};
Point(1368) = { 0.00003160, 0.00099950, 1.00000000, ss_lc};
Point(1369) = { 0.00012610, 0.00199060, 1.00000000, ss_lc};
Point(1370) = { 0.00028350, 0.00297820, 1.00000000, ss_lc};
Point(1371) = { 0.00050290, 0.00395680, 1.00000000, ss_lc};
Point(1372) = { 0.00078050, 0.00491760, 1.00000000, ss_lc};
Point(1373) = { 0.00111400, 0.00586120, 1.00000000, ss_lc};
Point(1374) = { 0.00149900, 0.00678320, 1.00000000, ss_lc};
Point(1375) = { 0.00193200, 0.00768320, 1.00000000, ss_lc};
Point(1376) = { 0.00245990, 0.00864800, 1.00000000, ss_lc};
Point(1377) = { 0.00307410, 0.00964250, 1.00000000, ss_lc};
Point(1378) = { 0.00386450, 0.01077900, 1.00000000, ss_lc};
Point(1379) = { 0.00490630, 0.01210200, 1.00000000, ss_lc};
Point(1380) = { 0.00627150, 0.01362500, 1.00000000, ss_lc};
Point(1381) = { 0.00803040, 0.01534300, 1.00000000, ss_lc};
Point(1382) = { 0.01025700, 0.01724400, 1.00000000, ss_lc};
Point(1383) = { 0.01302200, 0.01930800, 1.00000000, ss_lc};
Point(1384) = { 0.01635200, 0.02148700, 1.00000000, ss_lc};
Point(1385) = { 0.02031000, 0.02376600, 1.00000000, ss_lc};
Point(1386) = { 0.02500000, 0.02614700, 1.00000000, ss_lc};
Point(1387) = { 0.03000000, 0.02840100, 1.00000000, ss_lc};
Point(1388) = { 0.03500000, 0.03043000, 1.00000000, ss_lc};
Point(1389) = { 0.04000000, 0.03227700, 1.00000000, ss_lc};
Point(1390) = { 0.04588300, 0.03426100, 1.00000000, ss_lc};
Point(1391) = { 0.05291100, 0.03641100, 1.00000000, ss_lc};
Point(1392) = { 0.06141100, 0.03874600, 1.00000000, ss_lc};
Point(1393) = { 0.07151400, 0.04121400, 1.00000000, ss_lc};
Point(1394) = { 0.08343700, 0.04377600, 1.00000000, ss_lc};
Point(1395) = { 0.09751400, 0.04640300, 1.00000000, ss_lc};
Point(1396) = { 0.11284000, 0.04885900, 1.00000000, ss_lc};
Point(1397) = { 0.13070000, 0.05128000, 1.00000000, ss_lc};
Point(1398) = { 0.15138000, 0.05359100, 1.00000000, ss_lc};
Point(1399) = { 0.17500000, 0.05568800, 1.00000000, ss_lc};
Point(1400) = { 0.20000000, 0.05737500, 1.00000000, ss_lc};
Point(1401) = { 0.22500000, 0.05859600, 1.00000000, ss_lc};
Point(1402) = { 0.25000000, 0.05941200, 1.00000000, ss_lc};
Point(1403) = { 0.27500000, 0.05987300, 1.00000000, ss_lc};
Point(1404) = { 0.30000000, 0.06001700, 1.00000000, ss_lc};
Point(1405) = { 0.32500000, 0.05987900, 1.00000000, ss_lc};
Point(1406) = { 0.35000000, 0.05948600, 1.00000000, ss_lc};
Point(1407) = { 0.37500000, 0.05886300, 1.00000000, ss_lc};
Point(1408) = { 0.40000000, 0.05803000, 1.00000000, ss_lc};
Point(1409) = { 0.42500000, 0.05700600, 1.00000000, ss_lc};
Point(1410) = { 0.45000000, 0.05580700, 1.00000000, ss_lc};
Point(1411) = { 0.47500000, 0.05444800, 1.00000000, ss_lc};
Point(1412) = { 0.50000000, 0.05294000, 1.00000000, ss_lc};
Point(1413) = { 0.52500000, 0.05129600, 1.00000000, ss_lc};
Point(1414) = { 0.55000000, 0.04952400, 1.00000000, ss_lc};
Point(1415) = { 0.57500000, 0.04763400, 1.00000000, ss_lc};
Point(1416) = { 0.60000000, 0.04563400, 1.00000000, ss_lc};
Point(1417) = { 0.62500000, 0.04352900, 1.00000000, ss_lc};
Point(1418) = { 0.65000000, 0.04132500, 1.00000000, ss_lc};
Point(1419) = { 0.67500000, 0.03902700, 1.00000000, ss_lc};
Point(1420) = { 0.70000000, 0.03663900, 1.00000000, ss_lc};
Point(1421) = { 0.72500000, 0.03416400, 1.00000000, ss_lc};
Point(1422) = { 0.75000000, 0.03160300, 1.00000000, ss_lc};
Point(1423) = { 0.77500000, 0.02895900, 1.00000000, ss_lc};
Point(1424) = { 0.80000000, 0.02623100, 1.00000000, ss_lc};
Point(1425) = { 0.82500000, 0.02342000, 1.00000000, ss_lc};
Point(1426) = { 0.85000000, 0.02052600, 1.00000000, ss_lc};
Point(1427) = { 0.87500000, 0.01754500, 1.00000000, ss_lc};
Point(1428) = { 0.90000000, 0.01447700, 1.00000000, ss_lc};
Point(1429) = { 0.92328000, 0.01153900, 1.00000000, ss_lc};
Point(1430) = { 0.94333000, 0.00894260, 1.00000000, ss_lc};
Point(1431) = { 0.96040000, 0.00668320, 1.00000000, ss_lc};
Point(1432) = { 0.97519000, 0.00441040, 1.00000000, ss_lc};
Point(1433) = { 0.98808000, 0.00221820, 1.00000000, ss_lc};
Spline(200) = {1300:1433,1300};
// Now add in points to max a box.
edge_lc = 0.2;
Point(1200) = { 10, 5, 0, edge_lc};
Point(1201) = { 10, -5, 0, edge_lc};
Point(1202) = { -10, -5, 0, edge_lc};
Point(1203) = { -10, 5, 0, edge_lc};
Point(1204) = { 10, 5, 1, edge_lc};
Point(1205) = { 10, -5, 1, edge_lc};
Point(1206) = { -10, -5, 1, edge_lc};
Point(1207) = { -10, 5, 1, edge_lc};
//Connect the points to form lines
Line(201) = {1204, 1207};
Line(202) = {1207, 1206};
Line(203) = {1206, 1205};
Line(204) = {1205, 1204};
Line(205) = {1200, 1203};
Line(206) = {1202, 1202};
Line(207) = {1203, 1202};
Line(208) = {1202, 1201};
Line(209) = {1201, 1200};
Line(210) = {1204, 1200};
Line(211) = {1205, 1201};
Line(212) = {1206, 1202};
Line(213) = {1207, 1203};
//Connect the lines to form surfaces
Line Loop(214) = {201, 202, 203, 204};
Line Loop(216) = {205, 207, 208, 209};
Line Loop(219) = {201, 213, -205, -210};
Line Loop(221) = {213, 207, -212, -202};
Line Loop(223) = {203, 211, -208, -212};
Line Loop(225) = {211, 209, -210, -204};
Line Loop(227) = {200};
Line Loop(229) = {100};
Plane Surface(220) = {219};
Plane Surface(222) = {221};
Plane Surface(224) = {223};
Plane Surface(226) = {225};
Plane Surface(231) = {214, 227};
Plane Surface(232) = {216, 229};
//Create Physical surfaces
Physical Surface("top") = {220};
Physical Surface("front") = {231};
Physical Surface("back") = {232};
Physical Surface("inlet") = {222};
Physical Surface("bottom") = {224};
Physical Surface("outlet") = {226};
//Create volume
Surface Loop(233) = {231, 220, 222, 232, 224, 226};
Volume(234) = {233};
Physical Volume("air") = {234};
_______________________________________________
gmsh mailing list
[email protected]
http://www.geuz.org/mailman/listinfo/gmsh
